ABSTRACT
Up to 40% of genetic and rare disorders (RD) present facial dysmorphologies. Visual assessment of facial gestalt is commonly used for clinical diagnosis, health management and treatment monitoring. Quantitative approaches to facial phenotypes are more objective and provide first diagnoses of RD with relatively high accuracy, but are mainly based on populations of European descent, disregarding the influence of population ancestry. Here we assessed the facial phenotypes associated to four genetic disorders in a Latino-American population from Colombia. We recorded the coordinates of 18 facial landmarks in 2D images from 79 controls 51 pediatric individuals diagnosed with Down (DS), Morquio (MS), Noonan (NS) and Neurofibromatosis type 1 (NF1) syndromes. We quantified facial differences using Euclidean Distance Matrix Analysis (EDMA) and assessed the diagnostic accuracy of Face2gene, an automatic deep learning algorithm with widespread use in the clinical practice.
Quantitative comparisons indicated that individuals diagnosed with DS and MS were associated with the most severe phenotypes, with 58.2% and 65.4% of facial traits significantly different as compared to controls. The percentage decreased to 47.7% in NS and to 11.4% in NF1. Each syndrome presented a characteristic pattern of facial dysmorphology, supporting the potential of facial biomarkers for disorder diagnosis. However, our results detected population-specific traits in the Colombian population as compared to the facial gestalt described in literature for DS, NS and NF1. When clinical diagnosis based on genetic testing was used to verify the diagnosis based on 2D facial pictures, our results showed that Face2Gene accuracy was very high in DS, moderate in NS and NF1, and very low in MS, with low gestalt similarity scores in highly admixed individuals. Our study underscores the added value of precise quantitative comparison of facial dysmorphologies in genetic and rare disorders and the need to incorporate populations with diverse contributions of Amerindian, African and European ancestry components to further improve automatic diagnostic methods.
INTRODUCTION
According to the Online Mendelian Inheritance in Man (OMIM) databank, there are more than 10,000 genetic and rare diseases (RD) affecting 7% of the world’s population (Nguengang Wakap et al., 2020; Viteri et al., 2020), which corresponds to approximately 500 million people. Although as a whole genetic and RD are a significant cause of morbidity and mortality in the pediatric population (Suárez-Obando, 2018), by separate each disorder affects a very reduced number of people. Depending on the country, the prevalence to consider a disease as rare ranges from 1 affected individual in 50,000 people to 1 in 200,000. This low prevalence has limited the research on rare disorders.
Currently, there is limited knowledge on the etiology of these disorders, a reduced percentage of diseases (20%) presents a known molecular basis associated to a detailed phenotype description, and treatment is only available for 0.04% of RD (Suárez-Obando, 2018). As orphan diseases, many RD are chronic and incurable, representing severe and debilitating conditions (Cortés, 2015). The diagnosis and management of genetic RD is currently a clinical challenge (Schieppati et al., 2008). Precise and early diagnosis is crucial for individuals and their families to get effective care and to reduce disease progression. However, due to the limited knowledge and complexity of these pathologies, diagnosis may take several years (Bannister et al., 2020) and people may experience a diagnostic odyssey with delays in their correct treatment and management (González-Lamuño & García Fuentes, 2008). For most rare diseases there are no reliable biomarkers for an early diagnosis (Gülbakan et al., 2016).
Among the wide constellation of clinical symptoms associated to genetic and rare disorders, craniofacial dysmorphologies emerge as potential biomarkers (Gurovich et al., 2019; Hallgrímsson et al., 2020). These phenotypes are highly prevalent (Bannister et al., 2020; Viteri et al., 2020) and are commonly used for diagnosis, management and treatment monitoring of genetic and RD (Bannister et al., 2020). Up to 40% of these disorders present characteristic craniofacial phenotypes, as in Down, Morquio, Noonan, Apert, Rett, Fragile X, Williams-Beuren and Treacher-Collins and velocardiofacial syndromes, as well as other conditions such as microcephaly, holoprosencephaly, palate/lip cleft and other 2,000 rare genetic disorders (Farrera et al., 2019; Hallgrímsson et al., 2020).
The genetic and environmental factors associated with these disorders alter the complex process that orchestrates facial morphogenesis during pre- and postnatal development, inducing facial dysmorphologies. Facial development is highly regulated by multiple signaling pathways (Martínez-Abadías et al., 2011; Richtsmeier & Flaherty, 2013; Hallgrímsson et al., 2015; Marcucio et al., 2015), including Fibroblast Growth Factor (FGF), Hedgehog (HH), Wingless (WNT) and Transforming Growth Factor Beta (TGF-β) and Bone Morphogenetic Proteins (BPMs). Disruptions in the regulation of any of these signaling pathways can lead to facial dysmorphogenesis (Kouskoura et al., 2011).
The facial patterns associated with each disorder are unique but vary within and among diagnostics, from subtle facial anomalies to severe malformations (Jones et al., 2021). In the clinical practice, craniofacial dysmorphology is commonly assessed through qualitative visual assessment and basic anthropometric measurements. However, this approach may not capture with maximum precision the anatomical complexity of the facial traits associated with these disorders. Qualitative descriptions of facial phenotypes are sometimes based on general terms such as coarse face, large and bulging head; saddle-like, flat bridged nose with broad, fleshy tip; or malformed teeth (Aase, 1990; Johannes et al., 2003; Reardon et al., 2007). Accurate identification of dysmorphic features for diagnosis depends on the clinician’s expertise and highly trained dysmorphologists are required to recognize the facial “gestalt” characteristic of many dysmorphic syndromes, especially for rarest disorders (Reardon et al., 2007).
Recent research seeks to incorporate into the clinical diagnosis of RD the use of objective and quantitative tools to assess facial phenotypes (Hammond et al., 2004; Hammond, 2007; Hammond and Suttie, 2012; Hurst, 2018; Köhler et al., 2019; Agbolade et al., 2020). Automated systems have been developed to improve and to accelerate the diagnostic process (Gurovich et al., 2019, Hallgrímsson et al., 2020, Hsieh et al., 2021). Within the clinical practice, Face2Gene is the most commonly used system (FDNA Inc., https://www.face2gene.com/), a community-driven phenotyping platform trained over 17,000 people representing more than 200 syndromes (Gurovich et al., 2019). Face recognition is performed on 2D images that can be collected with any type of digital camera or phone, without previous training. Syndrome classification is achieved using DeepGestalt, a cascade Deep Convolutional Neural Network (DCNN)-based method that achieved 91% top-10 accuracy in identifying the correct syndrome (Gurovich et al., 2019).
Diagnostic approaches based on 3D photogrammetry have been developed more recently (Hammond et al., 2004; Hammond, 2007; Hallgrímsson et al., 2020). The advantage of 3D facial models is that they are more efficient than 2D images in capturing the complexity of facial phenotypes, but their widespread use is limited because the photographic equipment required for generating 3D models is usually not commonly available in the clinical practice. Hallgrímsson et al. (2020) analyzed 3D facial models from 7,057 subjects including subjects with 396 different syndromes, relatives and unrelated unaffected subjects (https://www.facebase.org/). Deep phenotyping based on quantitative 3D facial imaging and machine learning presented a balanced accuracy of 73% for syndrome diagnosis (Hallgrímsson et al., 2020).
Automated methods have thus demonstrated high potential to facilitate the diagnosis of facial dysmorphic syndromes (Gurovich et al., 2019; Bannister et al., 2020; Hallgrímsson et al., 2020; Hsieh et al., 2021). These tools present high accuracy diagnosis in European and North American populations in which the software predictions are mostly based. However, the tools have not been thoroughly tested in populations with different ancestry, and it is not well understood the how facial phenotypes associated with genetic and RD might be influenced by the complex patterns of population ancestry characterizing human populations.
Population ancestry in facial dysmorphologies: a long-disregarded factor
Facial shape shows wide variation across world-wide human populations (Xiong et al., 2019). Facial differences between populations detected in the shape of the forehead, brow ridges, eyes, nose, cheeks, mouth and jaw (Qiao et al., 2018) result from divergent evolutionary and adaptive histories of human populations occurred during the evolution of Homo sapiens over the last 200,000 years. Continuous migration and admixture keep shaping the facial phenotypes of human populations. Admixed populations display a variety of craniofacial morphologies, resembling the parental groups or presenting novel phenotypes (Martínez-Abadías et al., 2006), depending on dominance and epistatic interactions between alleles fixed or predominant in each parental group (Quinto-Sánchez et al., 2015). Therefore, the evolutionary and population dynamics of human populations results in genetic and phenotypic patterns that surrogate population ancestry (Ruiz-Linares et al., 2014; Sheehan and Nachman, 2014; Quinto-Sánchez et al., 2015) and can modulate the facial phenotypes associated to disease.
Few studies to date have analyzed the craniofacial phenotypes associated with genetic and RD in populations of non-European descent (Kruszka, et al., 2017; Kruszka, Addissie, et al., 2017; Kruszka, Sobering, et al., 2017; Dowsett et al., 2019), leaving African, Asian and Latin-American populations often disregarded and underrepresented. Thus, there is no reliable representation of facial phenotypes in genetic and rare diseases in populations of non-European descent. To develop quantitative approaches that efficiently diagnose genetic and RD in populations from all over the world, it is crucial to account for the influence of population ancestry on facial variation.
To contribute covering this gap, here we assessed the facial dysmorphologies associated to prevalent RD in a Latin-American population from the Southwest of Colombia. Latin-Americans are fascinating cases of hybrid/admixed populations that evolved on relatively short periods of time (Quinto-Sánchez et al., 2015; Mendoza-Revilla et al., 2021). Peopling of the Americas likely started 33,000 years ago (Ardelean et al., 2020; Becerra-Valdivia & Higham, 2020) by migration waves coming from North and South East Asia, and probably Australia (Castro E Silva et al., 2021). The first settlers arrived at North America and migrated towards southern regions following coastal and continental routes (González-José et al., 2003). Amerindian populations established all over the continent and adapted to a variety of environments over thousands of years. During the last 600 years, admixture with European and African populations has further shaped the genetic ancestry of Latin-American populations (Salzano & Bortolini, 2002; Salzano & Sans, 2014). In particular, the population from the region of Cali is the result of diverse migratory processes throughout the history of the Colombian settlement (Urrea-Giraldo & Álvarez, 2017). Admixture with the indigenous Amerindian population began in the sixteenth century with the arrival of the Spanish colonizers. In the eighteenth century, large colonial settlements of slaves brought from Africa were established in Cali for the exploitation of sugar cane and changed the population structure of Valle del Cauca. Nowadays, the population of Cali is characterized by indigenous and mestizo communities, with Amerindian and African ancestry components predominating over the European ancestry contribution (Urrea-Giraldo & Álvarez, 2017).
In this study, we compared the facial phenotypes associated to four genetic and RD, including Down syndrome, Mucopolysaccharidosis type IVA metabolic disorder known as Morquio syndrome, and two types of RASopathies, Noonan syndrome (NS) and neurofibromatosis type 1 (NF). The facial phenotype of these syndromes has not been previously characterized in Latin-American or Colombian populations and differences between populations with different ancestry backgrounds have not been assessed (Kruszka et al., 2017; Dowsett et al., 2019). Here, we quantitatively assessed the facial phenotypes associated to these syndromes and tested the hypothesis whether the diagnostic accuracy of automated methods commonly used in the clinical practice is affected by the admixed ancestry of this Colombian population of non-European descent.
MATERIALS AND METHODS
Participant recruitment for photographic sessions
The sample comprised 130 individuals from 0 to 18 years old from Valle del Cauca, a Southwest region in Colombia. The cohort included 79 age matched controls and 51 pediatric individuals diagnosed with Down, Morquio, Noonan and Neurofibromatosis type 1 syndromes that were recruited from the clinical genetics consultation at Hospital-Fundación Valle del Lili in Cali (Colombia), a tertiary health reference center for these genetic and rare disorders. See Table 1 for details on sample composition. In all participants, clinical diagnoses were confirmed by molecular genetic testing.
Human cohort by diagnostic group. Total number of participants as well as separated by sex (Male (M) / Female (F)) are provided.
Down syndrome (DS, OMIM 190685) was first selected because it is one of the most common genetic disorders, and previous studies have shown that the clinical manifestations associated with DS vary across ethnicities (Kruszka et al., 2017). DS presents a worldwide prevalence of 14 per 10,000 live births, with life expectancy increasing from 25 to 60 years in developing countries (Glasson et al., 2002; Roper & Reeves, 2006; Patterson, 2009; Aivazidis et al., 2017). In most Latino-American regions, the real incidence of patients with DS remains unknown and is usually underreported. A cross-sectional study in Brazil reported a DS birth rate of 4 for every 10,000 live births (Laignier et al., 2021); whereas in Colombia several studies have reported a prevalence rate between 1 for every 1,000 to 5 per 10,000 live births (Hernández Ramírez & Manrique Hernández, 2006; Valencia Arana et al., 2008). DS is an aneuploidy caused by trisomy of chromosome 21 and is the leading genetic cause of intellectual disability (Roper & Reeves, 2006). DS is associated with craniofacial dysmorphologies that impair vital functions such as breathing, eating, and speaking. DS craniofacial phenotype is described as brachycephalic, with maxillary hypoplasia leading to facial flatness, depressed nasal bridge and reduced airway passages (Starbuck et al., 2013); dysplastic ears with lobe absence; and eyes with oblique palpebral fissures, epicanthal folds, strabismus and nystagmus (Jones et al., 2021; Korayem & Bakhadher, 2019). Oral alterations include open mouth, cleft lip, lingual furrows and protrusion, macroglossia, micrognathia and narrow palate (Hennequin et al., 1999; Oliveira et al., 2008).
Within RD, we included Morquio syndrome type A (MS, OMIM 253000) because Colombia presents one of the highest prevalence of MS in the world, probably as a result of founder effects. The worldwide prevalence, including other Latin-American countries, ranges from 1 case per 75,000 to 200,000 live births, whereas in Colombia the prevalence rises up to 0.68 per 100,000 live births (Pachajoa et al., 2021). Morquio syndrome is a subtype of Mucopolysaccharidosis disorder caused by more than 180 autosomal recessive mutations in the GALNS gene (Herrera et al., 2017) that alter the metabolism of the extracellular matrix glycosaminoglycans (Sawamoto et al., 2020). Keratan and chondroitin sulfate alterations cause irreparable damage to leukocytes and fibroblasts and the typical alterations of MS thus involve the supporting tissue and the osteoarticular system (Suárez-Guerrero et al., 2016). Individuals with MS display skeletal abnormalities such as skeletal dysplasia, short stature and trunk, kyphoscoliosis, pectus carinatum, genu valgum, and joint hyperlaxity (Ortiz-Quiroga et al., 2018). Oral diseases often include periodontal disease, malocclusions, caries and premature tooth loss (Herrera et al., 2017). Individuals with MS show coarse facies with an excessively rapid growth of the head (Suárez-Guerrero et al., 2013). Craniofacial features include a prominent forehead, hypertelorism, prognathism, wide mouth and nose, depressed nasal bridge, plump cheeks and lips with an oversize tongue (Suárez-Guerrero et al., 2013). Dysmorphologies associated with MS vary among individuals, with typical severe phenotypes, although less severe forms have been described as mild or attenuated phenotypes (Suárez-Guerrero et al., 2016).
Finally, we also included in the analyses two RASopathies, Noonan syndrome and Neurofibromatosis type 1, which are prevalent in Valle del Cauca and present altered craniofacial development by genetic mutations that cause Ras/MAPK pathway dysregulation. Noonan syndrome (NS, OMIM 163950) is the most common type of RASopathy and is a rare genetically heterogeneous autosomal dominant disorder caused by mutations in either the PTPN11, SOS 1, KRAS, BRAF or RAF1 genes. The worldwide prevalence of NS is 1 per 1,000 to 1 per 2,500 live births (Hernández-Martín & Torrelo, 2011). Individuals with NS display facial features such as hypertelorism, epicanthic folds, strabismus, downward slanting palpebral fissures, ptosis, high arched palate, deeply grooved philtrum with high peaks of upper lip vermillion border, midfacial hypoplasia and micrognathia, broad flat nose, low-set posteriorly rotated ears, curly/sparse/coarse hair and short webbed neck (Athota et al., 2020).
Neurofibromatosis type 1 (NF1, OMIM 162200) is an autosomal dominantly inherited neurocutaneous disorder caused by a mutation in the neurofibromin gene. NF1 has a worldwide incidence of 1 in 2,500 to 1 in 3,000 individuals. The clinical manifestations in NF1 are variable, even within the same family group, and the timing of the onset has a major influence (Hernández-Martín & Torrelo, 2011). Regarding craniofacial traits, individuals with NF1 present macrocephaly, facial asymmetry caused by dysplasia of the sphenoid wings (Khosrotehrani et al., 2003), as well as bone deformities caused by plexiform neurofibromas, enlarged mandibular canal, retrognathic mandible and maxilla, and short cranial base (Visnapuu et al., 2018).
To assess the facial phenotypes associated with these disorders, individuals diagnosed with DS, MS, NS and NF1 and age matched controls were recruited for photographic sessions at educational and research centers in Cali (Colombia) in 2021. The photographic material was taken under the protocol approved by ethics committee “Human Research Ethics Committee of the Icesi University” with Approval Act No. 309”. To photograph the participants and to record relevant clinical information, we obtained informed consent from the participants or from their parents or legal guardians in the case of minor children, in accordance with national guidelines and regulations.
Facial image acquisition and anatomical landmark collection
Facial shape was captured from 2D images taken using a professional digital camera (SONY Alpha 58 +18-55) that was attached to a tripod and placed at one-meter distance in front of the participants. To capture a natural facial gesture, the images were acquired in an upright position with facial neutral expression.
To measure facial shape of each individual and to detect the traits associated with each disorder, we recorded the 2D coordinates of a set of 18 anatomical facial landmarks (Figure 1 and Supplementary Table 1). Landmarks were acquired using an automatic facial landmark detection procedure adapted from the open-source software library Dlib (King, 2009). From the set of landmarks registered by Dlib, 15 landmarks directly matched our configuration of 18 facial landmarks (Figure 1). Three additional landmarks were approximated through direct computations between the landmarks coordinates automatically returned by Dlib: the glabella was computed as the midpoint point between the innermost points located in the eyebrows, and the palpebrale inferius landmarks of the right and left eyes were computed as the midpoint between the two central lower eyelid landmarks. The validity of the data was assessed by comparing the coordinates of landmarks automatically detected by Dlib to the coordinates of landmarks manually collected by an expert facial morphologist in a set of 20 facial images. The average deviation measured as the root mean square error (RMSE) was 2 mm, which falls within acceptable measurement error standards (Stull et al., 2014).
Anatomical position of facial landmarks used in morphometric and statistical analyses to quantify dysmorphologies associated to genetic and rare disorders Down, Morquio, Noonan syndromes and Neurofibromatosis type 1 in a Colombian population.
After inspection of landmark coordinates for gross error, anatomical landmark coordinates were used to conduct morphometric analyses.
Quantification of facial phenotypes
We used Euclidean distance matrix analysis (EDMA) to describe the facial phenotype associated to each syndrome and to detect specific facial traits that were significantly different between controls and individuals diagnosed with DS, MS, NS and NF. EDMA is a robust morphometric method for assessing local differences between samples (Lele & Richtsmeier, 1991). EDMA detects which linear distances significantly differ between pairwise sample contrasts and compares patterns of significant differences across samples.
First, to account for size differences between subjects, the 2D coordinates of the facial landmarks of each subject were scaled by their centroid size, estimated as the square root of the sum of squared distances of all the landmarks from their centroid (Rohlf & Slice, 1990). After scaling, as EDMA represents shape as a matrix of linear distances between all possible pairs of landmarks, a total of 153 unique facial measurements were calculated for each individual. Linear distances were compared for each group of DS, MS, NS and NF syndromes with control individuals by performing a two-tailed two-sample shape contrasts on all unique inter-landmark linear distances from each sample. Statistical significance was assessed using a non-parametric bootstrap test with 10,000 resamples. EDMA statistically evaluated the number of significant local linear distances in each two-sample comparison based on confidence interval testing (α = 0.10). Relative differences between patients and controls were computed as (mean distance in controls -mean distance in patients) / mean distance in controls. To pinpoint specific local shape differences and to reveal the unique morphological pattern of variation associated with each disorder, the ten longest and shortest significant relative differences were plotted on facial figures.
Facial dysmorphology score
To confirm that the results were not random due to the small sample sizes available in rare diseases, we combined the results from EDMA with an iterative bootstrapping method that further assessed whether the facial dysmorphologies associated to each syndrome were statistically significant (Starbuck et al, 2021). First, we estimated from the EDMA results a facial dysmorphology score (FDS) as the percentage of significantly different distances between patient and control groups. Then, we ran simulations with random samples of controls and patients generated by iterative bootstrapping to assess the statistical significance of the patterns revealed by EDMA. For each disorder, we first created subsamples of N randomly chosen controls (where N is the total number of patients available in the sample). Then, using a subsampling approach, we automatically generated random pseudo-subsamples containing a known number of patients (namely M). This procedure was repeated with increasing numbers of patients and resulted in a series of staggered pseudo sub-samples that contained from M = 0 to M = N patients. In each of these simulations, we computed an EDMA analysis and an FDS score. The results from each round of random groups were separately represented in histograms. Finally, we compared the distribution of FDS random values with the FDS observed in the whole sample. The P-value assessing the statistical significance of the comparison (α = 0.05) was computed as the ratio between the number of simulations containing no patients that provided a higher FDS than the observed FDS and the total number of simulations.
Face2Gene diagnostic assessment
To assess the accuracy of automated diagnostic methods in the Colombian sample, we compared the clinical diagnosis based on genetic testing with the diagnosis estimated by the Face2Gene technology (FDNA Inc., Boston, MA, USA) using the frontal facial 2D images. Following Gurovich et al. (2019), we assessed the top-one and top-five accuracies for each disorder, estimated as the percentage of cases where the Face2Gene model predicted the correct syndrome as the first result or within the five first results from the sorted list of probable diagnosis. We also calculated these accuracies expanding the diagnostic range to the disorder family.
Moreover, we evaluated the similarity between the facial Gestalt models used by Face2Gene for syndrome classification and the Colombian individuals with clinical diagnosis. For each individual, we selected the first diagnostic prediction that matched their clinical and genetic diagnosis and recorded the Gestalt similarity. We classified the level of similarity between the individual and the corresponding Gestalt model into seven categories, including “very low”, “low”, “low-medium”, “medium”, “medium-high”, “high”, and “very high” Gestalt similarity, using the “gestalt level” barplot provided by Face2Gene.
RESULTS
EDMA analyses showed that each syndrome presented a characteristic facial phenotype. In individuals with Down syndrome, all facial structures including the eyes, nose and mouth presented significant differences as compared to controls. Overall, DS was associated with wider but shorter facial traits (Figure 2A).
A) EDMA results. Dotted lines represent facial measurements significantly different in control and patient groups. Blue lines indicate measurements that are shorter in patients as compared to controls, whereas orange lines represent measurements that are longer in patients. B) Iterative bootstrapping tests based on facial dysmorphology scores (FDS). Histograms represent the simulation results for each random group separately, which contain an increasing number of patients, from no patients (M = 0) to all patients (M = N). The blue line shows the FDS score obtained with the complete sample of control and patients. * Statistically significant P-value.
Results showed 6.5% increase of relative distance between the midpoint between the eyebrows (glabella) and the most inferior medial point of the lower right eyelid (palpelabre inferius), and 7.5% increase between the right palpelabre inferius and the outer commissure of the right eyes (exocanthion), indicating hypertelorism. Additionally, in this Colombian sample, people with DS exhibited longer measurements in the buccal portion, with a 6-8% increase of mouth width as measured from the crista philtri to the chelions (Figure 2A). However, the midfacial and nasal regions were reduced (Figure 2A). People with DS presented a 6-8% reduction in measurements of midfacial height, with the largest difference detected as a 9.7% reduction of the distance between the tip and the root of the nose (Figure 2A). The facial dysmorphology score (FDS) indicated that up to 58.2% of facial traits were significantly different in people with DS (Figure 2B).
The facial pattern associated with Morquio syndrome was also characterized by wider and shorter midfacial traits, as observed in Down syndrome. However, facial dysmorphologies were more abundant and severe in MS than in DS, with 65.4% of facial traits significantly different in diagnosed individuals and higher percentages of relative change (Figure 3 A, B).
For details see legend in Figure 2.
The most affected regions were the midface and the nose, whereas the mouth was the least affected. Individuals with MS presented hypertelorism, with 14% increase in the distance between the midpoint between the eyebrows (glabella) and the inner commissures of the left and right eyes (endocanthions), even larger compared to the differences observed in DS. Individuals with MS also showed larger distances in the base of the nose, with 14-19% increase in the distance from the tip of the nose to the insertion of the right and left alar bases (subalare) as compared to controls. Mouth width was also increased in MS. On the contrary, midfacial heights measuring the distance between the eyes and the nose were significantly reduced from 10 to 16% in individuals with MS (Figure 3A).
In Noonan syndrome, facial dysmorphologies were abundant and concentrated in the orbital and nasal regions. EDMA detected significantly increased distances in the upper face, but decreased distances in the midface (Figure 4A). Patients presented a lower position of the eyes, with 9 to 13% increased distances between the glabella or sellion and the landmarks located in the eyes. The mouth also showed a more inferior position, with 8-10% increased relative distance between the tip of the nose and the superior lip, but the shape of the mouth did not show large differences between patients and controls. The reduction of midfacial heights in individuals with NS ranged from 5 to 11%, with a similar magnitude as in DS (Figure 4A). FDS indicated that 47.7% of facial traits were significantly different in NS (Figure 4B).
For details see legend in Figure 2.
Neurofibromatosis type 1 was associated with minor facial dysmorphologies, which were less abundant and less severe than the previous syndromes (Figure 5A). Individuals with NS only presented 11.4% of significantly different facial traits as compared to controls and the percentages of relative change were low, mostly ranging from 1 to 5% (Figure 5A, B). The largest difference was a 10% increase in facial distance between the glabella and the labiale superius (Figure 5A). Along with larger distances in the midline of the face, EDMA detected reduced distances on the right and left sides of the face, with shorter distances from the right and left chelion to the eye landmarks, the endocanthion and the palpebrale inferius. Hypertelorism was not present in individuals with NF1 (Figure 5A). In NF1, the FDS score was not significant (Figure 5B), indicating that the facial dysmorphology pattern associated with NF1 is so subtle that overall is not larger than facial differences that can be randomly detected within control individuals.
For details see legend Figure 2.
Nevertheless, the simulation tests demonstrated that the facial dysmorphologies measured as FDS scores on individuals diagnosed with Down, Morquio and Noonan syndromes were significant and different from random comparisons in control individuals. Few simulations resulted in a higher FDS than the FDS obtained with the complete sample (Figures 2B, 3B, 4B, first row and blue line). Moreover, in DS, MS and NS, facial dysmorphology scores increased as larger numbers of diagnosed individuals were included in the simulations (Figures 2B, 3B, 4B, middle rows), confirming the severity of the facial dysmorphologies associated to these syndromes. Finally, the simulations comparing all recruited diagnosed individuals with random subsamples of control individuals (Figure 2B, 3B, 4B, last row) indicated that FDS scores can range widely from 10 to 80%, underscoring the biasing effects of control samples.
After quantifying the facial dysmorphologies associated to DS, MS, NS and NF1 in the Colombian sample, we tested the accuracy of the diagnosis provided by the automatic diagnostic algorithms of Face2Gene. We assessed the correspondence between the estimated Face2Gene diagnosis based on facial frontal 2D images with the diagnosis made by clinicians based on genetic testing.
Face2Gene estimated Down syndrome diagnosis with top-1 accuracy of 100%, as DS diagnosis was listed as the first diagnosis in all individuals with an average Gestalt similarity of 6.2 (Table 2). In 79% of individuals, DS diagnosis was associated with very high Gestalt similarity values, but in 21% of individuals the Gestalt similarity was lower, ranging from medium-high to very low values (Table 2, Figure 6).
% Of cases matching the genetic diagnosis are provided for each syndrome, as well as Gestalt similarity values.
Histograms based on top-5 accuracy Face2Gene predictions within family disorder.
On the contrary, the top-1 accuracy of Morquio syndrome was 0%, as Face2Gene never listed the specific diagnostic of mucopolysaccharidosis type IVA (MPSIVA) as a first prediction (Table 2). Although Face2Gene could not identify the specific type of MS, the automatic diagnostic algorithms associated the facial dysmorphologies with a diagnosis related with mucopolysaccharidosis disorders in 36.4% of cases, with a medium-high average Gestalt similarity of 5.6 (Table 2). When the first 5 diagnostic predictions were considered, the top-5 accuracy raised to 45.4% for exact MPSIVA diagnosis and to 100% for mucopolysaccharidosis disorders, but with a low-medium Gestalt similarity (Table 2, Figure 6).
The top-1 accuracy of Face2Gene for Noonan syndrome was 66.7%, with a medium-high average Gestalt similarity of 5.2 (Table 2). Top-5 accuracy increased to 77.8% for exact NS diagnosis, and to 88.9% when considering Noonan Syndrome-Like Disorder diagnoses, with a medium Gestalt similarity of 4.4 and wide variation among individuals (Table 2, Figure 6).
Finally, in Neurofibromatosis type 1, Face2Gene presented a top-1 accuracy of 8.3% associated with a very low Gestalt similarity of 1 (Table 2). When diagnoses within the RASopathies disorder family were considered, the top-1 accuracy raised to 50%, as 5 out of 12 individuals were diagnosed as Noonan syndrome (Table 2). The top-5 accuracy was 66.6% and was associated with low Gestalt similarity values of between 1 and 2 in 87.5% of individuals (Table 2, Figure 6).
DISCUSSION
Our analyses provided an accurate quantitative comparison of facial dysmorphologies in Down, Morquio and Noonan syndromes, as well as in Neurofibromatosis type 1, in a Latin-American population from Colombia. An objective and highly detailed description of the facial phenotype is a major improvement over qualitative descriptions of the complex facial dysmorphologies associated with these genetic disorders. We quantified local facial trait differences presented in people diagnosed with these disorders as compared to age matched controls of the same population, localizing the largest statistically significant facial dysmorphologies.
Our results indicated differential facial patterns associated with each disorder, with major significant dysmorphologies in Down, Morquio and Noonan syndromes, and minor facial dysmorphologies associated with NF1. Different types of genetic alterations, which ranged from aneuploidy and overall genetic imbalance as in DS; to point genetic mutations affecting different processes, such as the metabolism of mucopolysaccharides in MS; or signaling pathways, such as the RAS/MAPK pathway in NS and NF1, significantly affected the facial phenotypes. In genetic and rare diseases, genetic alterations deviate the signaling pathways regulating normal facial development (Jones et al., 2021; Terrazas et al., 2017), altering normal morphogenesis and growth during pre- and postnatal development (Kouskoura et al., 2011).
Population-specific facial traits in Colombian individuals with rare disorders
Overall, the facial patterns associated to each syndrome in the Colombian Latin-American population coincide with the descriptions reported in the literature (Khosrotehrani et al., 2003; Suárez-Guerrero et al., 2013; Starbuck et al., 2013; Athota et al., 2020). However, there are specific local traits that differ, suggesting that facial traits associated to disease might be modulated by population ancestry and influenced by different evolutionary and adaptive histories of human populations (Kruszka et al., 2017).
In Down syndrome, we found a facial dysmorphologies that are consistent with the patterns reported in the literature for populations of European descent. Our analyses detected differences in linear facial measurements that correspond to typical DS traits such as hypertelorism, maxillary hypoplasia and shorter and wider faces associated to a brachycephalic head (Jones et al., 2021; Korayem & Bakhadher, 2019). Results also suggested other characteristic traits of DS, such as midfacial retrusion and depressed nasal bridge (Starbuck et al., 2013). Open mouth and macroglossia (Hennequin et al., 1999; Oliveira et al., 2008) were also observed during the photographic sessions in the participants of our study. However, in contrast to European and North American populations (Starbuck et al., 2021), in the Colombian population we detected that the mouth was wider in individuals diagnosed with DS as compared to euploid controls. This difference may be caused by unnatural facial gestures of the participants when asked to close the mouth during the photo shoot, or by facial differences associated to ancestry. Indeed, Kruszka et al., (2017) analyzed individuals diagnosed with DS in diverse populations and showed craniofacial differences between individuals from different ethnicities (Africans, Asians, and Latin Americans), demonstrating that ancestry is a relevant factor when assessing craniofacial variation associated to rare diseases.
In Morquio syndrome, the facial dysmorphologies observed in Colombian individuals were consistent with traits reported in the literature, which included hypertelorism, prognathism, wide nose and wide mouth (Hung et al., 2016 Suárez-Guerrero et al., 2013; Herrera et al., 2017). In our pediatric cohort from Colombia, Morquio syndrome was associated with the most severe facial dysmorphologies. Considering that keratan and chondroitin sulfate alterations associated with MS accumulate over life and cause irreparable damage to the osteoarticular system, facial dysmorphologies associated with MS are expected to increase with age due to an excessively rapid growth of the head, becoming even more severe in adult individuals (Suárez-Guerrero et al., 2013). Further research is required to assess this hypothesis and to test whether pharmacological treatments can slow down the progression of the disease and reduce the facial dysmorphologies associated with MS. This is especially relevant in Colombia, a country with one of the highest prevalence of MS in the world (Pachajoa et al., 2021).
In Noonan syndrome, we detected hypertelorism, downward slanting palpebral fissures, and midfacial hypoplasia in the Colombian population, as reported in populations of European descent (Athota et al., 2020). Moreover, our results quantified changes in the position of the mouth in Colombian individuals diagnosed with NS not reported before (Lores et al., 2020).
The facial pattern associated with NF1 in individuals from Colombia was also compatible with some typical traits of NF1, such as midface hypoplasia (Hernández-Martín & Torrelo, 2011). However, our results did not underscore facial asymmetry or hypertelorism as prominent facial differences between diagnosed individuals and controls in the Colombian population (Hernández-Martín & Torrelo, 2011).
Comparative quantitative studies from different world regions are not usually available for most genetic and rare disorders, and reference data for diagnosis is mainly based on phenotypes defined on populations of European descent. In fact, almost no images of individuals of Latin American origin are included in reference medical texts such as Smith’s Recognizable Patterns of Human Malformation (Jones et al., 2021). However, previous studies have reported that in Noonan and other genetic syndromes, such as Turner, 22q11.2 and Cornelia de Lange syndromes presented distinctive facial traits that were population specific, with clinical features that were significantly different in Africans, Asians, and Latin Americans (Kruszka et al., 2017). Our results in a Colombian population further support this evidence, highlighting the need to extend the analyses to populations from all over the world to achieve a complete and more accurate phenotypic representation of genetic and RD to optimize the diagnostic potential of facial biomarkers in the clinical practice.
Reduced accuracy diagnosis in a Colombian population with diverse ancestry
Deep learning algorithms such as Face2gene have shown potential as a reliable and precise tool for genetic diagnosis by image recognition (Gurovich et al., 2019; Hsieh et al., 2021; Park et al., 2022; Pascolini et al., 2022). In the Colombian cohort analyzed here, Face2gene diagnosed with 100% accuracy Down syndrome, which is one of the most common and more accurately represented genetic disorder. In Noonan syndrome, the percentage of top1-accuracy was lower (66.7%), but still correctly identified the disorder in the majority of individuals when considering the top5-accuracy within Noonan syndrome-like disorders (88.9%). However, in Morquio syndrome, despite being associated with the most severe facial dysmorphologies, the top1-accuracy for exact diagnosis of Mucopolysaccharidosis type IVA was 0% and a low percentage of cases (36.4%) were diagnosed with a Mucopolysaccharidosis-like syndrome as the first prediction.
Interestingly, there was a wide range of variation in gestalt similarity for most disorders, even for Down syndrome. Visual inspection of cases suggested that most individuals associated with the lowest gestalt similarity scores were highly admixed, with diverse contributions of Amerindian, African and European ancestry components. Future analyses need to investigate the role of population ancestry and further assess the reliability and validity of automatic diagnostic tools in admixed populations from non-European descent. As evidenced in this study, where syndromes such as MS, NS, NF1 were associated with low accuracy scores versus more frequent syndromes such as DS, this is critical in rare syndromes with heterogenous clinical presentation and phenotype, where clinical diagnosis is a challenging process (Schieppati et al., 2008; Bannister et al., 2020) that may take several years, leading to the so-called diagnostic odyssey (González-Lamuño & García Fuentes, 2008).
Precise and early diagnosis of genetic and rare disorders is crucial for adequate health care and clinical management. Without a diagnosis, individuals and their families must proceed without basic information regarding their health and future developmental outcomes (Bannister et al., 2020). Even though gene-based technologies have greatly improved diagnostic procedures (Agbolade et al., 2020), the mutations causing many rare diseases are still unknown and access to genetic testing is limited (Suárez-Obando, 2018). Genetic consultations may become a long process, and broad molecular testing such as exome and genome sequencing represent a high expense that is not affordable for all families and health care systems, especially in low-medium income countries (González-Lamuño & García Fuentes, 2008).
In this context, faster, non-invasive and low-cost diagnostic methods based on facial phenotypes emerge as complementary tools for providing earlier first reliable diagnoses (Gurovich et al., 2019; Agbolade et al., 2020; Hallgrímsson et al., 2020; Hsieh et al., 2021), especially in low-middle income countries with low budget and resources for molecular testing. Qualitative visual assessment of craniofacial dysmorphology is frequently employed for diagnosis, clinical management and treatment monitoring (Jones et al., 2021). Experts in dysmorphologies can identify the facial “gestalt” distinctive of many dysmorphic syndromes (Jones et al., 2021). However, this facial assessment relies on the expertise of the clinician and is very challenging, as there is no clear one-to-one correspondence between disorders and facial dysmorphologies. Different genetic mutations can cause the same syndrome or similar phenotypes, whereas the same mutation can induce different phenotypes (Aldridge et al., 2010; Martínez-Abadías et al., 2011). In addition, within the same rare disease there may be several subtypes, and symptoms may vary even within individuals of the same genetic disorder and the same family (Suárez-Obando, 2018). This complex biology generates confusion at the time of diagnosis and warrants the development of efficient, objective and reliable diagnostic methods.
Computer-assisted phenotyping can overcome these pitfalls and provide widely accessible technologies for quick syndrome screening (Bannister et al., 2020). In this automated approach, methods can be based on 2D or 3D images (Gurovich et al., 2019; Hallgrímsson et al., 2020; Hsieh et al., 2021). The advantage of 2D methods is that data collection is easy and can be readily translated into the clinical practice, as physicians can take facial images even with simple digital cameras or smartphones. The collection of 3D models is more sophisticated and requires specialized equipment but provides more accurate phenotype descriptions by incorporating the depth dimension.
To further improve the methods of craniofacial assessment to diagnose individuals with genetic syndromes and RD that exhibit facial dysmorphologies, it is crucial to assess the large morphological variation displayed by human populations in facial phenotypes. Factors such as age, sex and ancestry should be accounted for in diagnostic methods. Clinical manifestations in some genetic disorders usually begin at an early age, with two thirds of patients expressing symptoms before the second year of birth (Suárez-Obando, 2018), but in other disorders facial dysmorphologies develop later, during postnatal development. Male and female faces present sexual dimorphism at adulthood (Enlow & Hans, 1996) and diseases can differently affect the facial phenotype depending on sex differences (Martínez-Abadías et al., 2021). As suggested here, differences in population ancestry can also significantly influence the facial phenotype and the potential for clinical diagnosis.
Therefore, recruitment of participants must be expanded to include as many individuals with RD as possible, together with large comparative samples of age-matched controls, from both sexes and from all the populations in the world representing varied ancestries. For instance, the population in Southwestern Colombia is characterized by high levels of admixture from people with Native American, African, and European ancestry (Urrea-Giraldo & Álvarez, 2017; Adhikari et al., 2017). Including the morphological variation of faces from such different ancestry backgrounds is key to pinpoint the facial dysmorphologies associated with diseases in worldwide diverse populations (Conley et al., 2017). Our simulation analyses further highlighted the importance of maximizing the recruitment of diagnosed and control individuals, as results considerably changed depending on the cohort and sample sizes
CONCLUSION
Our results underscore that facial phenotypes associated with genetic and rare disorders are influenced by population ancestry (Kruszka et al., 2017; Dowsett et al., 2019) and that diverse genetic background variation can modulate the phenotypic response to disease, affecting the accuracy of current tools of clinical diagnosis. In the future, deep learning algorithms including a high variety of populations with different ancestry backgrounds will optimize the precision and accuracy of diagnosis in an unbiased approach. Such predictive models will support clinicians in decision-making across the world.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
ACKNOWLEDGEMENTS
We are grateful for the voluntary collaboration of all participants, including children and their families. We are thankful to Colegio Ecológico Scout and Universidad Icesi for granting us permission to organize the photographic sessions. We thank Max Rubert for technical photographic assistance. We acknowledge support from Proyecto COL0012168-1097 Interfacultades-ICESI, and 2017SGR1630 Research group (AGAUR).
LITERATURE CITED
