Anatomic development of the upper airway during the first five years of life: A three-dimensional imaging study.

PURPOSE: Normative data on the growth and development of the upper airway across the sexes is needed for the diagnosis and treatment of congenital and acquired respiratory anomalies and to gain insight on developmental changes in speech acoustics and disorders with craniofacial anomalies.
METHODS: The growth of the upper airway in children ages birth to 5 years, as compared to adults, was quantified using an imaging database with computed tomography studies from typically developing individuals. Methodological criteria for scan inclusion and airway measurements included: head position, histogram-based airway segmentation, anatomic landmark placement, and development of a semi-automatic centerline for data extraction. A comprehensive set of 2D and 3D supra- and sub-glottal measurements from the choanae to tracheal opening were obtained including: naso-oro-laryngo-pharynx subregion volume and length, each subregion's superior and inferior cross-sectional-area, and antero-posterior and transverse/width distances.
RESULTS: Growth of the upper airway during the first 5 years of life was more pronounced in the vertical and transverse/lateral dimensions than in the antero-posterior dimension. By age 5 years, females have larger pharyngeal measurement than males. Prepubertal sex-differences were identified in the subglottal region.
CONCLUSIONS: Our findings demonstrate the importance of studying the growth of the upper airway in 3D. As the lumen length increases, its shape changes, becoming increasingly elliptical during the first 5 years of life. This study also emphasizes the importance of methodological considerations for both image acquisition and data extraction, as well as the use of consistent anatomic structures in defining pharyngeal regions.

Introduction

The upper airway, a virtual conduit as characterized by Marcus et al. [1], has an anatomic boundary defined by other tissues (bony, cartilaginous and soft) while serving the functions of respiration, food ingestion (mastication and deglutition), as well as vocalization/speech, hence the function-based terms ‘aerodigestive tract’, ‘vocal tract’, or more comprehensively the ‘aerodigestive and vocal tract’. During the course of development, especially from infancy to early childhood, the upper airway undergoes drastic changes in size, shape and mechanical properties due to the restructuring of its anatomical sub-components, such as the descent of the larynx and the hyoid bone [2-4]. The anatomic growth process persists while adapting to the various functional needs and demands during maturation. As posited by current theory on craniofacial growth, the development of the upper airway is shaped by both genetic as well as intrinsic and extrinsic epigenetic factors, such as function, mechanical forces, and trauma [5-13].

The lack of knowledge regarding the growth and development of the upper airway, defined as the air conduit from the level of the nose to the carina, was addressed in a workshop by the National Heart, Lung, and Blood Institute (NHLBI) in 2009 with a large team of clinicians and scientists from diverse fields in healthcare and the biological sciences [1]. The outcome was a comprehensive set of research guidelines on various aspects of the upper airway, each with a set of priorities relevant to clinical disorders of upper airway functions. Among the priorities was the need to study the developmental changes of the upper airway anatomy and function during childhood (neonatal to puberty) across sexes and ethnicities and to provide normative values of the upper airway. Normative data are needed to better understand common respiratory disorders such as obstructive sleep apnea syndrome (OSAS), as well as a number of other congenital and acquired respiratory anomalies [1]. Furthermore, normative data can provide additional insight on developmental speech acoustics [14, 15], as well as speech disorders, particularly where craniofacial anomalies are present [16, 17]. As listed in Table 1, a large number of studies have examined the upper airway anatomy using different modalities, methodologies, airway regions, and age ranges. Table 1 summarizes the studies to date that have examined the typical development of the aerodigestive and vocal tract from the choanae or the soft palate superiorly to the epiglottis or the trachea inferiorly. A subset of studies listed have factored in growth and/or sex in their data analysis. Most studies have employed imaging to obtain quantitative measurements, including linear, angular and/or area measurements, based on the midsagittal or axial slices, as well as volumetric measurements. However, only a very limited number of studies have assessed multidimensional volumetric measurements during early childhood. Among the 34 studies listed in Table 1, only 17 studies included linear, area and volumetric measurements, and fewer than half of those studies controlled for head position during or after data acquisition. Of the 12 studies summarized in Table 1 that examined the pre-pubertal age group, the majority obtained measurements in 2D that were collected primarily from radiographic images using mid-sagittal, axial, or coronal visualization planes; an approach frequently used to assess the upper airway, as it is cost effective and less time-consuming to process. However, this approach does not provide accurate representation of the complex airway morphology, as it overlooks information of lateral dimensions [18, 19]. Two of those 12 studies [20, 21] quantified the prepubertal airway in 3D but only Abramson et al. [20] covered the entire prepubertal period from birth to 5 years and assessed sexual dimorphism. Neither of those retrospective studies reported controlling for head position or using it as an inclusion criterion.

Table 1

Summary of studies on typical upper airway development.

Methodology		Study			Age Cohort (yrs)				Assessment/Control		Measurement Type					Airway Region
Type	Modality	Author (year) a	n = Total (M/F)	Age	Pre-Pub (0–5)	Peri-Pub(6–10)	Pubertal (11–17)	Post-Pub (18+)	Sexual Dimorph-ism	Head Position	Dist-ance	Region Length	CSA (*)	Volume	Center-line
2D	CT	Ronen et al. (2007)	n = 69 (38/31)	4–10, 14–19	+	+	+	+	+		YES					PNS—hyoid superior
	CT, MRI	Vorperian et al. (2009)	n = 605 (327/278)	0–19	+	+	+	+	+		YES	YES				Choanae—glottis
	MRI	Fitch & Giedd (1999)	n = 129 (76/53)	2–25	+	+	+	+	+		YES	YES				PNS—C4ai
	MRI	Litman et al. (2002)	n = 16 (9/7)	10mos-7yrs	+	+				+	YES		YES (a)			Oral-pharyngeal
	MRI	Machata et al. (2010)	n = 138 (60/78)	0–6	+	+				+	YES		YES (a)			Soft palate—tongue
	MRI	Yi et al. (2017)	n = 521 (296/225)	0–6	+	+			+		YES					Soft palate—tip of epiglottis
	MRI	Fregosi et al. (2003)	n = 18 (13/5)	7–12		+	+				YES		YES (a)	YES		Choanae—epiglottis base
	X-ray	Barbier et al. (2015)	n = 966 (494/470)	0–25	+	+	+	+	+		YES	-	-	-	-	PNS—tip of epiglottis
	X-ray	Goncalves et al. (2011)	n = 390 (195/195)	6–18		+	+	+	+		YES					Oral-pharyngeal
	X-ray	Sheng et al. (2009)	n = 239 (107/132)	7–27		+	+	+	+		YES					PNS—mandible/pharynx intersection
	X-ray	Mislik et al. (2014)	n = 880 (458/422)	6–17		+	+				YES					PNS—epiglottis base
	X-ray	Daraze et al. (2017)	n = 117 (48/69)	21–25				+	+	+	YES					PNS—tip of epiglottis
	Fluoroscopy	Rommel et al. (2003)	n = 23 (14/9)	0–4	+				+		YES	YES				Choanae—epiglottis base
3D	CBCT	Chiang et al. (2012)	n = 387 (173/214)	8–18		+	+	+	+	+	YES		YES (o)	YES	YES	PNS—C4ai
	CBCT	Jiang et al. (2014)	n = 254 (119/135)	6–18		+	+	+	+	+	YES	YES	YES (x)	YES		Choanae—hyoid superior
	CBCT	Lenza et al. (2010)	n = 34 (14/20)	11–56		+	+	+			YES		YES (x)	YES		Choanae—tip of epiglottis
	CBCT	Schendel et al. (2012)	n = 1300 (571/729)	6–60		+	+	+	+		YES		YES (a)	YES		PNS—C4as
	CBCT	Anandarajah et al. 2017	n = 105 (44/61)	7–13		+	+						YES (a)	YES		PNS—epiglottis base
	CBCT	Masoud et al. (2020)	n = 81 (32/49)	7–17		+	+			+	YES	YES	YES (a)	YES		Choanae—epiglottis base
	CBCT	Yanagita et al. (2017)	n = 61 (0/12)	10–16		+	+			+	YES	YES	YES (a)	YES		Choanae—epiglottis base
	CBCT	Alves et al. (2012)	n = 50 (27/23)	8–10		+				+	YES		YES (a)	YES		PNS—tip of epiglottis
	CBCT	Claudino et al. (2013)	n = 54 (21/33)	13–20			+	+	+	+	YES	YES	YES (a)	YES		Choanae—hyoid inferior
	CT	Abramson et al. (2009)	n = 46 (31/15)	4mos-64yrs	+	+	+	+	+		YES		YES (a)	YES		PNS—epiglottis base
	CT	Li et al. (2011)	n = 281 (141/140)	10–18		+	+	+		+	YES	YES	YES (a)	YES		PNS—epiglottis base
	CT	Kim et al. (2011)	n = 73 (64/9)	14–72			+	+	+	+	YES		YES (a)	YES		PNS—hyoid superior
	CT	Gibelli et al. (2020)	n = 80 (40/40)	21–86				+	+		YES			YES		Choanae—glottis
	CT	Inamoto et al. (2015)	n = 107 (30/77)	23–77				+	+	+	YES			YES		PNS—Glottis
	CT	Shigeta et al. (2008)	n = 38 (19/19)	24–84				+	+		YES	YES		YES		PNS—tip of epiglottis
	MRI	Smitthimedhin et al. (2018)	n = 96 (53/43)	0.5-46wks	+									YES		Choanae—trachea
	MRI	Welch et al. (2002)	n = 12 (0/12)	30–64				+		+	YES		YES (o)	YES	YES	PNS—tip of epiglottis
Other	APh	Leboulanger et al. (2011)	n = 59 (28/31)	2–18	+	+	+	+	+	+			YES			Oral-pharyngeal
	APh	Martin et al. (1997)	n = 114 (60/54)	16–74			+	+	+				YES	YES		Oral-pharyngeal
	APh, Spiro	Brooks et al. (1992)	n = 175 (77/98)	19–64				+	+				YES	YES		Airway
	APh, Spiro	Brown et al. (1986)	n = 24 (14/10)	26–42				+	+				YES	YES		Airway

Summary of studies on typical upper airway development listed per methodology in the first column, using measurement type/dimension: 2D (two-dimensional), indicates that measurements were collected from radiographic images representing either the mid-sagittal, axial or coronal visualization planes; 3D (three-dimensional), indicates that measurements were collected utilizing multiplanar visualizations (axial, coronal and/or sagittal planes) and 3D representation of the upper airway; Other indicates that measurements were made using non-imaging techniques. The first column also lists study modality including imaging (CT, CBCT, MRI, X-ray) or non-imaging (Acoustic Pharyngometry (APh), Spirometer (Spiro)) techniques. The second column lists alphabetically per modality the study author(s) with year of publication in parentheses, sample size (n =) with total male/female (M/F) numbers specified in parentheses, and age range examined. Age of study participants is also classified using pubertal age cohorts in the third column, followed by assessment or control of sex-differences and head position in the fourth column. The fifth column lists measurement reported including overall pharyngeal distance, and pharyngeal region/subregion measurements including: length, cross-sectional area* (CSA), volume and centerline length. The final column lists the defined superior to inferior anatomical boundaries in each study.

* CSA specifications include the automatic or manual calculation of the area of a plane of the pharyngeal airway in: the axial plane (a), a plane that is orthogonal to the pharyngeal centerline (o), or a plane connecting two or more specific anatomical landmarks (x).

a Cited literature: [20–53].

Since the upper airway is a lumen, attention must be paid to a number of methodological considerations, given their potential effect on various pharyngeal measurements. Methodological procedures known to affect pharyngeal measurements include head/neck position (flexion/extension), body position (upright/supine), and sedation [37, 54, 55]. While most studies to date have accounted for one or more of these confounders, it is difficult to compare findings across studies unless all confounders have been addressed. For example, Inamoto et al. [47], reported significant sex differences between the adult male and female laryngopharynx, but Gibelli et al. [46], who also used CT but did not control for head position, reported no sex differences. Additionally, variations in the anatomical boundary and subregion borders of the upper airway morphology (nasopharynx, oropharynx, and laryngopharynx/hypopharynx) as defined by different studies, summarized in Table 1 (final column) and Table 2, further complicates the ease and feasibility of comparing findings across studies. Thus, to ensure an accurate and reliable assessment of the developmental changes of this cavity, it is critical to use standardized imaging procedures, well-defined anatomical regions, and established airway data extraction protocols.

Table 2

Summary of boundaries of the pharyngeal regions.

Main Pharyngeal Subdivisions	Anatomical boundaries	Definitions b
Nasopharynx	Superior Border	Nares (Ayappa & Rapoport, 2002)
		Nasal cavity (Adewale, 2009)
		End of nasal septum (Netter, 2019)
		Choanae (Arens et al., 2004; Logan et al., 2017; Schuenke et al., 2010)
	Inferior Border	Hard palate (Ayappa & Rapoport, 2002)
		Soft palate (Laird et al., 2019; Moore et al., 2006; Standring et al., 2017)
		Level of soft palate (Arens et al., 2004)
		Above soft palate (Adewale, 2009; Gu et al., 2016)
		Inferior/Lower border of soft palate (Netter, 2019; Logan et al., 2017)
Oropharynx	Superior Border	Soft palate (Ayappa & Rapoport, 2002)
		Uvula (Schuenke et al., 2010)
	Inferior Border	Epiglottis (Ayappa & Rapoport, 2002; Schuenke et al., 2010)
		Superior/Upper border of epiglottis (Gu et al., 2016; Laird et al., 2019; Moore et al., 2006; Netter, 2019; Standring et al., 2017)
		Larynx (Arens et al., 2004)
Laryngopharynx/Hypopharynx	Superior Border	Posteriorlateral to the larynx (Arens et al., 2004)
		Base of tongue (Ayappa & Rapoport, 2002)
		Tip of epiglottis (Adewale, 2009)
		Superior/Upper border of epiglottis (Logan et al., 2017)
	Inferior Border	Larynx (Ayappa & Rapoport, 2002)
		Opening of esophagus (Netter, 2019)
		Cricoid cartilage (Schuenke et al., 2010)
		Inferior/Lower border of cricoid cartilage (Adewale, 2009; Gu et al., 2016; Laird et al., 2019; Logan, 2017; Moore et al., 2006; Standring et al., 2017)

Summary of pharyngeal regions’ boundaries as defined in relevant anatomy textbooks and published papers. Note the lack of consistency and/or specificity in the anatomical boundaries for each region.

b Cited literature: [56–65].

This study aims to systematically study the anatomic development of the upper airway, specifically the structural changes from the choanae to the tracheal opening (inferior border of the cricoid cartilage) from birth to 5 years, as compared to adults. To acquire normative data of the anatomy that subserve the aerodigestive and speech functions, we used an imaging database with computed tomography (CT) studies from typically developing individuals to obtain a comprehensive set of two-dimensional (2D) and three-dimensional (3D) measurements quantifying the growth of the upper airway. Our comprehensive set of methodological criteria included control of head position, histogram-based upper airway segmentation, placement of anatomic landmarks, and development of a semi-automatic method to determine lumen centerline. In addition, to better understand the resonance/acoustic characteristics of the vocal tract, this study aimed to examine the nature of the developmental changes of the upper airway dimensions and to determine if there are sex differences in the upper airway dimensions during the pre-pubertal period. We hypothesize all pediatric airway dimensions to be substantially smaller than adult dimensions. We also hypothesize sex differences in both children and adults.

Materials and methods

a. Image acquisition/dataset

Using imaging studies performed at the University of Wisconsin Hospital and Clinics (UWHC), our Vocal Tract Development Lab (VTLab) has curated a lifespan retrospective database of more than 2000 head and neck CT scans to study the anatomic growth and development of the oral and pharyngeal structures. This database was established following approval of the University of Wisconsin-Madison Institutional Review Board (IRB) and anonymized accordingly. All CT imaging studies, performed in the supine body position, were acquired using CT scanners manufactured by General Electric Medical Systems or Siemens and stored in Digital Imaging and Communications in Medicine (DICOM) format. Additional details on this imaging database and image acquisition are provided in Kelly et al. [66], Miller et al. [67, 68] and Vorperian et al. [23].

To ensure the adequacy of imaging studies selected for this study, the VTLab imaging database was reviewed for typically developing cases between the ages 0–5 years (pediatric) and 20–30 years (adults) who were imaged for conditions that do not affect typical growth. A total of 410 (208 Males (M), and 202 Females (F)) CT imaging studies that included 264 (161M, 103F) pediatric scans and 146 (47M, 99F) adult scans, from 276 (143M, 133F) individuals (195 [115M, 80F] children; and 81 [28M, 53F] adults), were inspected for cases that met the following inclusion criteria: (1) slice thickness ≤ 2.5mm, (2) 14-22cm field-of-view (FOV), (3) 512x512 matrix size, (4) no movements or dental artifacts affecting the view of pharynx structure, and (5) neutral or flexed head position as confirmed using Miller et al.’s [68] head position classification protocol. While all extreme flexion/extension cases were excluded, including all sedation cases, neutral-flexed head position cases were not excluded given that the larger infant head is prone to being flexed in the supine position. The total yield of cases that met the inclusion criteria for this study’s dataset included 61 (32M, 29F) pediatric cases from 78 imaging studies (41M, 37F), and 17 (9M, 8F) adult cases from 72 (39M, 33F) imaging studies. The individuals whose images were used included 56 (31M, 25F) children, and 16 (8M, 8F) adults. Age specific demographics are presented in Table 3.

Table 3

Distribution of male and female cases per age group.

Group (Age range (yr;mos))	M	F	Total
<1 (00;00–00;11)	4	5	9
1 (01;00–01;11)	4	5	9
2 (02;00–02;11)	6	7	13
3 (03;00–03;11)	9	8	17
4 (04;00–04;11)	6	7	13
5 (20;00–30;00)	8	9	17

Distribution of male (M) and female (F) cases per age group. Age groups specified in years; months (group <1 includes cases birth (00;00) to 11 months (00;11); group 1 includes cases 1 year (01;00) to 1 year 11 months (01;11) etc., and group 5 adults ages 20-to-30 years.

b. Image reconstruction

The standard reconstruction kernel was the preferred CT reconstruction algorithm, and was available for the majority of the imaging studies. For cases/imaging studies processed without the standard kernel, imaging features of the standard kernel were simulated by processing the soft kernel with an unsharp enhance filter using a kernel size of 5x5, or by processing the bone kernel with a low pass filter using a kernel size of 3x3. Next, the software Analyze 12.0 [69] was used to reconstruct CT images from DICOM format into 3D volume.

A histogram-based threshold method was applied to the reconstructed CT volume in order to identify the intensity in Hounsfield Unit (HU) that allows an optimal representation of the airway. Guided by the technique of Nakano et al. [70], per image, we used the midpoint between the air threshold peak (-1000 HU) and soft tissue peak (+100 HU to +300 HU), as the applied upper threshold intensity to segment the airway. The range of upper thresholds used in this study was between -556 HU and -445 HU. The Volume Render and Volume Edit modules were then used to visualize and segment the 3D pharynx model from the reconstructed CT volume. Using the identified threshold value, the airway region studied was restricted inferiorly at the first tracheal ring (lower limit of the cricoid cartilage), and superiorly at the choanae. The resulting 3D pharynx model was saved in Analyze Object Map format [.obj].

c. Anatomic landmarks and variables

As depicted in Figs 1 and 2, and listed with descriptions in Table 4, a set of 26 anatomic landmarks that included 20 pharyngeal, 4 maxillary, and 2 reference landmarks were manually placed on each of the 78 3D pharynx models to quantify upper airway growth. The set of landmarks selected were carefully determined following a thorough review of landmarks and airway variables examined in studies to date [20, 24, 36, 47, 54, 71]. Landmark placement entailed using the Volume Render module in Analyze 12.0 [69], to manually place each of the 26 landmarks by overlaying them on their respective CT images while using the axial, coronal and sagittal planes to guide accuracy of landmark placement. The landmarks were similarly saved in Analyze Object Map format. To ensure reliability in landmark placement, prior to data collection, two researchers modeled and placed landmarks on six cases. The differences in resulting measurements, calculated from the raters’ landmarks, had an average relative error (ARE) that was less than or equal to 5% between researchers. The landmarks were then used to establish a data extraction protocol, described in the following section, that generates pharyngeal cross sections perpendicular to the centerline and calculates landmark-based measurements. The comprehensive set of 30 pharyngeal variables measured, as listed and defined in Table 5 below, are described in the section on variable measurements.

Fig 1

Illustration of airway regions and measurements.

The airway was examined using landmark-derived planes orthogonal to the centerline, as described in text. The four airway regions bounded by five cross sectional areas, a-to-d as depicted in the right panel, were quantified developmentally using the following measurements: volume, region length, cross-sectional area (CSA), anterior-posterior distance, and lateral width—as defined in Table 5. The airway regions above the glottis (d; Table 5, definition 19), included the following pharyngeal regions: a. Nasopharynx (blue; definition 1); b. Oropharynx (red; definition 6); c. Laryngopharynx (cyan; definition 11); and the airway below the glottis consisted of the subglottal region (magenta; definition 22). The pharynx (Table 5, definition 16) consisted of all three supraglottal regions.

Fig 2

3D airway model (choanae to trachea) of a 4-year 8-month old typically developing male as visualized in MATLAB.

Blue dots represent the 26 anatomic landmarks listed in Table 4. The centerline is shown in green, and CSAs closest to each of the anatomic landmarks are shown in red.

Table 4

Description/definition of the 26 anatomic landmarks.

#	Description of landmarks	Landmark Name	Abbreviation
	Pharynx Landmarks
1	The point of attachment of the vocal folds with the thyroid cartilage at the anterior commissure, 2–3 mm below the thyroid notch of the larynx.	Glottis Anterior	(ga)
2	In the axial plane at the level of the glottis as determined by the tear-shaped glottal area (GA), the most posterior point of the laryngo-pharynx between the lateral-most sides of the vocal folds attached to the arytenoid cartilages.	Glottis Posterior	(gp)
3	Most superior and posterior point of the epiglottis in the midsagittal plane.	Epiglottis Superior	(epS)
4	Attachment of the epiglottis with the hyoepiglottic ligament, visualized as the most anteroinferior point of contact in the midsagittal plane.	Epiglottis Base	(epBse)
5	The most posterior point on the pharyngeal wall at the level of the epiglottis base in the midsagittal plane.	Epiglottis Base Posterior	(epBsePo)
6	Most inferior point of the left piriform sinus.	Piriform Sinus Inferior Left	(PSInL)
7	Most inferior point of the right piriform sinus.	Piriform Sinus Inferior Right	(PSInR)
8	The midpoint (most ‘curved’ point) of the left aryepiglottic fold. Approximately halfway between the base and tip of the epiglottis.	Piriform Sinus Superior Left	(PSSuL)
9	The midpoint (most ‘curved’ point) of the right aryepiglottic fold. Approximately halfway between the base and tip of the epiglottis.	Piriform Sinus Superior Right	(PSSuR)
10	Most inferior point on the left vallecula.	Vallecula Inferior Left	(VaInL)
11	Most inferior point on the right vallecula.	Vallecula Inferior Right	(VaInR)
12	Most anterior point of the anterior pharyngeal wall at the level of the velum tip as visualized in the midsagittal plane.	Velum Anterior	(VeAn)
13	Most posterior point of the posterior pharyngeal wall at the level of the velum tip, as visualized in the midsagittal plane.	Velum Posterior	(VePo)
14	Midpoint between VeAn and VePo (landmarks 12 and 13 as defined above). Midpoint was calculated based on VePo and VeAn landmark coordinates.	Midpoint between VeAn and VePo	(MidVe)
15	The most supero-posterior point of the velum.	Velum Back	(VeBa)
16	The inferior tip of the velum.	Velum End	(VeEnd)
17	The most anterior point of the pharynx at the level of the PNS in the axial plane as guided by the midsagittal plane of the pharynx.	Nasopharynx Anterior	(NpxAn)
18	Most posterior point of the pharyngeal wall at the level of the PNS in the axial plane as guided by the midsagittal plane of the pharynx.	Nasopharynx Posterior	(NpxPo)
19	Midpoint between NpxAn and NpxPo (landmarks 17 and 18). Midpoint was calculated based on NpxAn and NpxPo landmark coordinates).	Midpoint between NpxAn and NpxPo	(NpxMid)
20	Landmark placed on the soft tissue between the posterior vomer bone and the nasal crest of the palatine bone.	Posterior Nasal Septum	(NasalS)
Reference Landmarks
21	The most anterior point of the Anterior Nasal Spine.	Anterior Nasal Spine	(ANS)
22	The most posterior point of the Posterior Nasal Spine.	Posterior Nasal Spine	(PNS)
Maxilla Landmarks
23	The most postero-inferior point of the maxillary alveolar bone in the midsagittal plane (landmark location is between the first incisors)	Alveolar bone of incisor	(ABI)
24	The most posterior point of the incisive canal in the midsagittal plane of the maxilla.	Posterior edge of incisive canal	(PIC)
25	The intersection between the transverse palatine suture and the median palatine suture.	Palatine Sutures intersection	(PALS)
26	Midpoint between PIC and PALS (landmarks 24 and 25) along the median palatine suture on the maxilla. Coordinates calculated based on PIC and PALS landmark coordinates x, y, z.	Maxilla midpoint	(MMax)

Description/definition of the 26 anatomic landmarks (pharynx, reference, and maxilla), listed from the inferior to the superior regions of the airway as displayed in Fig 1. Landmark placement entailed use of multiplanar views (at least two of the sagittal, axial, and coronal planes, or all three) for accuracy. These landmarks were used in defining study variables and extracting the quantitative measurements of the upper airway as specified in Table 5.

Table 5

All upper airway variables examined.

Measure-ment #	Variable Description	Variable Name	(Abbreviation)
Nasopharynx
1	Orthogonal volume of the region bound by the intersections of the centerline with the palatal plane (ANS-PNS) superiorly and inferiorly with the tip of the velum.	Nasopharynx Volume	Nasopharynx
2	The curvilinear segment length along the centerline of the Nasopharynx Volume.	Nasopharynx Length	NasopharynxL
3	Cross sectional area (CSA) of the superior border of Nasopharynx Volume.	Nasopharynx Area	NasopharynxArea
4	The distance between the most anterior and posterior points along the midline of the superior border of the Nasopharynx Volume.	Nasopharynx Anterior-Posterior Distance	NasopharynxAPDist
5	The distance between the most lateral left and right points along the midline of the superior border of the Nasopharynx Volume.	Nasopharynx Width	NasopharynxWidth
Oropharynx
6	Orthogonal volume of the region bounded superiorly by the tip of the velum, and inferiorly by the midpoint of the aryepiglottic folds.	Oropharynx Volume	Oropharynx
7	The curvilinear segment length along the centerline of the Oropharynx Volume.	Oropharynx Length	OropharynxL
8	CSA of the superior border of the Oropharynx Volume.	Oropharynx Area	OropharynxArea
9	The distance between the most anterior and posterior points along the midline of the superior border of the Oropharynx Volume.	Oropharynx Anterior-Posterior Distance	OropharynxAPDist
10	The distance between the most lateral left and right points along the midline of the superior border of the Oropharynx Volume.	Oropharynx Width	OropharynxWidth
Laryngopharynx
11	Orthogonal volume of the region bounded by the orthogonal planes where the centerline intersects superiorly with the tip of the midpoint of the aryepiglottic folds, and inferiorly with the glottis.	Laryngopharynx Volume	Laryngopharynx
12	The curvilinear segment length along the centerline of the Laryngopharynx Volume.	Laryngopharynx Length	LaryngopharynxL
13	CSA of the superior orthogonal border/plane of Laryngopharynx Volume.	Laryngopharynx Area	LaryngopharynxArea
14	The distance between the most anterior and posterior points along the midline of the most superior border of the Laryngopharynx Volume.	Laryngopharynx Antero-posterior Distance	LaryngopharynxAPDist
15	The distance between the most lateral left and right points along the midline of the most superior border of the Laryngopharynx Volume.	Laryngopharynx Width	LaryngopharynxWidth
Pharynx
16	Orthogonal volume of the supraglottal region bounded superiorly by the intersections of centerline with the palatal plane, and inferiorly by the glottis.	Total Pharynx Volume	PharynxVolume
17	The curvilinear length along the centerline of the Pharynx Volume.	Pharynx Length	PharynxLength
18	CSA of the superior border/plane of the Subglottal Volume.	Glottis Area	GlottisArea
19	The distance between the most anterior and posterior points along the midline of the most superior border of the Subglottal Volume.	Glottis Anterior-Posterior Distance	GlottisAPDist
20	The distance between the most lateral left and right points along the midline of the most superior border of the Subglottal Volume.	Glottis Width	GlottisWidth
26	The curvilinear distance extending from the posterior aspect of the maxillary incisor teeth (ABI, landmark 23) through the oral and pharyngeal cavities to the level of the glottis. (i.e. traditional VTL measures starting at incisors).	Vocal Tract Lengthi	VTLengthi
27	Curvilinear length of the velum from the PNS to the tip of the velum (VeEnd).	Velum Length	Velum Length
28	Length of the left piriform sinus measured as the 3D distance from the midpoint of the left aryepiglottic fold and the most inferior aspect of the left piriform sinus.	Piriform Sinus Length Left	PSLengthLeft
29	Length of the right piriform sinus measured as the 3D distance from the midpoint of the right aryepiglottic fold and the most inferior aspect of the right piriform sinus.	Piriform Sinus Length Right	PSLengthRight
30	Average length of the piriform sinus (average of PSLengthLeft and PSLengthRight if both exist, or the length of either if one is missing).	Average Piriform Sinus Length	AveragePSLength
Subglottal
21	Orthogonal volume of the region bounded superiorly by the glottis, and inferiorly by the opening of the trachea.	Subglottal Volume	Subglottal
22	The curvilinear segment length along the centerline of the Subglottal Volume	Subglottal Length	SubglottalL
23	CSA of the most inferior border of the Subglottal Volume, below the cricoid bone at the level of the first tracheal ring	Trachea Area	TracheaArea
24	The distance between the most anterior and posterior points along the midline of the most inferior border of the Subglottal Volume.	Trachea Anterior-Posterior Distance	TracheaAPDist
25	The distance between the most lateral left and right points along the midline of the most inferior border of the Subglottal Volume.	Trachea Width	TracheaWidth

The 30 upper airway variables examined. Measurements extracted for each region include: the orthogonal volume, curvilinear/centerline volume-length, the orthogonal superior or inferior cross-sectional area (CSAs) of each of the subregions, as well as the anterior-posterior distance (APDist) and lateral width (Width) of each CSAs (measurements 1-to-25). All planar measurements are orthogonal to the centerline. The sum of the nasopharynx, oropharynx and laryngopharynx subregions’ length and volume was used to calculate total pharynx volume and pharynx length (measurements 16 and 17). See text for additional vocal tract (VT) measurements (measurements 26–30).

d. Pharynx centerline and data extraction protocol

A semi-automatic, centerline-based data extraction pipeline was developed to extract quantitative measurements from the 3D pharynx in MATLAB (The MathWorks, Natick, MA). First, the built-in marching-cube algorithm in MATLAB was used to generate 3D meshes of the pharynx model to serve as input to the pipeline [72]. This pipeline adapted the implicit fairing diffusion method to smooth the 3D pharynx meshes iteratively while preserving the intrinsic geometric properties [73]. Next, a level-contour-based centroid-extraction method was applied on the smoothed pharynx, obtaining a set of coordinates along the tubular center of the pharynx [74-76]. These coordinates were further interpolated and smoothed with the B-spline de Boor algorithm, generating a centerline representative of the center of the airway lumen [77, 78]. This centerline was then used as input to an in-house written script that calculated planes orthogonal (i.e., perpendicular) to the line segment formed by each centerline coordinate and its subsequent centerline coordinate. Finally, the intersections between the orthogonal planes and the 3D meshes were extracted as boundary vertices. With the boundary vertices, cross sectional areas (CSAs) as well as additional variable measurements were calculated along the centerline. See Fig 2 for an illustration of the 3D pharynx model and the cross sections.

e. Variable measurements

A total of 30 airway variables, as listed and defined in Table 5, were measured by the above-described protocol using planes orthogonal to the centerline. The variables extracted are described below and include overall pharyngeal length and volume, modified vocal tract length (VTLength), velum length, and piriform sinuses length measurements, as well as measurements from the following four subregions: (i) Nasopharynx, (ii) Oropharynx, (iii) Laryngopharynx, and (iv) Subglottal. See Fig 1. Each subregion was isolated using its respective ‘landmark-derived planes’ orthogonal to the centerline using the following boundary definitions: The nasopharynx region was defined as an orthogonal volume bound by the intersection of the centerline with the palatal plane–formed by the anterior nasal spine (ANS) and posterior nasal spine (PNS) landmarks–superiorly, and with the tip of the velum inferiorly. The oropharynx region was defined as an orthogonal volume bound by the orthogonal planes at the tip of the velum superiorly, and by the aryepiglottic fold inferiorly. The laryngopharynx region, that includes the piriform sinuses, was defined as the orthogonal volume bound by the orthogonal planes at the midpoint of the aryepiglottic folds superiorly (the most curved point at approximate halfway between the base and tip of the aryepiglottic folds), and by the glottis inferiorly. The subglottal region was bound superiorly by the inferior boundary of the laryngopharynx region, and inferiorly by the first axial slice displaying the first tracheal ring. The first tracheal ring was used as a guide to the inferior border of the cricoid cartilage since the unossified cricoid cartilage in pediatric cases was difficult to delineate on the CT images [79].

Measurements extracted, as defined in Table 5, included for each region: the orthogonal volume, curvilinear/centerline volume-length, the orthogonal superior or inferior cross-sectional areas (CSAs, five total) of each of the subregions, as well as the anterior-posterior distance (APDist) and lateral width (Width) of each CSAs. The sum of the nasopharynx, oropharynx and laryngopharynx supraglottal subregions was used to calculate pharynx volume and pharynx length (Table 5, measurements 16 and 17). Additional vocal tract (VT) measurements (Table 5, measurements 26-to-30) included: Vocal Tract Lengthincisor (VTL); calculated as the curvilinear distance extending from the posterior border of the maxillary incisor (seen as the most anterior landmark in Fig 1) to the glottis, representing the vocal tract portion of the upper airway starting at the incisor i.e. excluding the lip and teeth region. Velum length; calculated as the curvilinear distance extending from the PNS to the tip of the velum (VeEnd). Piriform sinus length (PSLength); left, right, and average PSLength, measured using the defined anatomical landmarks (PSSuL/R and PSInL/R).

f. Statistical analysis

All statistical analyses were performed in R. A linear mixed-effect model was used to capture sex-specific growth in young children and allow for developmental comparison with adult pharyngeal morphology. This model, using the lmerTest package for mixed-effects in R, accounted for the repeat scans included in the dataset from individuals with multiple visits. The model was specified as follows: with “Adult” a dummy variable for adult subjects, “PediatricAge” the age for non-adult subjects (0 for adults), and α a random per-subject effect.

Outliers were first excluded using the model, by removing data points with residuals exceeding 2.576 of standard deviation, as described in [23, 71]. The model was then refitted on log-scale for each of the variables to assess for growth trends and sex differences.

Likelihood ratio test (LRT) was conducted to assess overall age effect in the first 5 years of life. To assess sex-differences, Wald test was performed at three time points: age groups <1 year, 5 years and adults. Tests were conducted at a nominal significance level of α = 0.05; in Table 6, significance at the stringent Bonferroni corrected level (< .0004) was also indicated. Finally, using point estimate of modeled means, percent growth at age 5 years was calculated using data at age group <1 year, and adults for the purpose of gaining insight on upper airway growth type as described by Scammon [80]. Scammon determined two primary postnatal growth types, neural and general growth types, or their combination, to characterize growth of head and neck structures. Furthermore, he noted that while all primary growth types are characterized by a period of rapid growth during infancy, by early childhood neural growth type achieves greater than two-thirds of the adult size, while somatic growth type barely achieves a quarter of the adult size [80].

Table 6

Likelihood ratio test and Wald test results.

Measure-ments	Variable	Age Effect LRT		Sex Effect—Wald Test											Percent Growth at Age 5
Types		p value				p value				p value			p value
		Pediatric		M	F	Age <1		M	F	Age 5	M	F	Adult		M	F
Volume (mm 3 )	Nasopharynx	0.0006	***	807.28	927.40	0.6125		2453.77	2105.76	0.5155	5221.76	4255.25	0.3591		51.96	47.51
	Oropharynx	0.0000	****	304.42	185.39	0.1546		879.79	1338.51	0.1570	7460.16	4927.30	0.1485		22.28	47.01
	Laryngopharynx	0.0000	****	265.04	296.85	0.6296		1064.72	1188.46	0.5845	6168.93	4215.44	0.0464	*	31.79	41.35
	PharynxVolume	0.0000	****	1450.67	1480.53	0.9229		4597.95	5257.92	0.4566	19742.95	13573.47	0.0288	*	33.80	48.10
	Subglottal	0.0000	****	132.25	154.83	0.3167		670.82	722.98	0.5761	5998.06	3463.51	0.0000	****	27.99	36.94
Volume Segment Length (mm)	NasopharynxL	0.0010	**	20.68	24.28	0.0800		30.15	29.74	0.8622	35.92	32.71	0.2105		62.14	64.75
	OropharynxL	0.0004	****	6.26	3.80	0.0635		10.35	15.14	0.0972	37.05	30.24	0.3612		13.30	42.89
	LaryngopharynxL	0.0017	**	13.11	13.32	0.8246		16.75	16.97	0.8364	27.60	22.01	0.0002	****	25.12	41.93
	PharynxLength	0.0000	****	40.54	42.09	0.4439		58.20	61.78	0.1530	101.65	84.68	0.0000	****	28.90	46.23
	SubglottalL	0.0000	****	7.47	10.07	0.0106	*	14.66	15.12	0.7570	28.00	26.44	0.5538		35.04	30.86
Area (mm 2 )	NasopharynxArea	0.0868		54.65	39.11	0.2446		74.51	83.92	0.6260	248.52	218.61	0.5852		14.80	34.08
	OropharynxArea	0.0434	*	56.37	53.05	0.8100		100.62	94.33	0.7651	171.08	140.44	0.3344		45.28	53.18
	LaryngopharynxArea	0.0014	**	48.66	56.11	0.5462		111.26	123.17	0.6120	344.29	276.37	0.2512		30.85	39.50
	GlottisArea	0.0000	****	11.28	5.65	0.0598		42.56	49.90	0.6103	121.90	74.74	0.1079		41.20	74.77
	TracheaArea	0.0000	****	23.20	24.41	0.6389		53.31	51.56	0.7178	238.62	150.92	0.0000	****	23.37	30.49
Anterior-Posterior Distance (mm)	NasopharynxAPDist	0.2087		5.41	3.61	0.0793		5.89	6.07	0.8801	14.47	12.40	0.4132		5.37	27.93
	OropharynxAPDist	0.0094	**	6.47	7.49	0.4237		10.44	12.38	0.2730	12.82	11.55	0.4807		62.44	120.23
	LaryngopharynxAPDist	0.3702		7.22	8.23	0.3171		8.91	9.03	0.8998	14.48	12.69	0.2206		23.26	18.01
	GlottisAPDist	0.0000	****	5.71	5.23	0.2835		9.78	10.48	0.3147	19.93	14.35	0.0000	****	28.66	57.59
	TracheaAPDist	0.0001	****	5.47	5.27	0.5954		7.52	6.85	0.1295	15.97	12.30	0.0000	****	19.48	22.50
Width (mm)	NasopharynxWidth	0.0086	**	11.37	10.35	0.5211		16.12	16.78	0.7464	22.31	23.03	0.7880		43.42	50.69
	OropharynxWidth	0.0163	*	8.31	8.81	0.7408		14.65	12.61	0.3252	19.60	15.60	0.1167		56.20	55.91
	LaryngopharynxWidth	0.0000	****	6.39	5.75	0.3324		10.35	12.00	0.1166	20.99	18.35	0.1341		27.11	49.58
	GlottisWidth	0.0000	****	2.58	1.76	0.1098		6.23	6.73	0.7073	8.92	7.16	0.2674		57.58	92.06
	TracheaWidth	0.0000	****	5.20	5.38	0.6007		8.52	8.75	0.6373	17.64	14.67	0.0005	***	26.72	36.27
VT Additional Length (mm)	VTLength	0.0000	****	78.88	76.92	0.4097		108.35	111.13	0.3312	159.81	141.01	0.0000	****	36.42	53.38
	VelumLength	0.0000	****	25.19	25.13	0.9533		33.61	34.40	0.5550	42.41	38.92	0.0206	*	48.85	67.22
	PSLengthLeft	0.3470		6.13	5.97	0.9279		9.01	8.26	0.6893	17.68	14.13	0.2472		24.99	28.02
	PSLengthRight	0.1376		6.01	5.79	0.8925		10.18	9.41	0.7050	20.74	15.64	0.1262		28.35	36.69
	AveragePSLength	0.0768		5.70	5.95	0.8571		9.80	8.76	0.5465	19.33	14.94	0.1376		30.09	31.31

Likelihood ratio test (LRT) results for age effect, and sex-effect using the Wald test at age <1 year, 5 years, and adults. Significant differences are denoted with asterisk (* < .05; ** < .01; *** < .001; **** < .0004 Bonferroni corrected value). Also, percent of adult size at age 5 years (final column) for each of the 30 variables as listed by measurement type (column 1), and airway sub-regions (column 2). Percentages are based on the point estimate of modelled means, see Fig 3 for 95% confidence intervals. The regions are listed superior to inferior with supra-laryngeal (above glottis) measurements listed first. Refer to Fig 1 to visualize variables and subregions, and Table 5 for variable definitions.

Fig 3

Measurements extracted for pharyngeal region (3A. Nasopharynx; 3B. Oropharynx; 3C. Laryngopharynx; 3D. Pharynx; and 3E. Subglottal). Measurements extracted for each pharyngeal region (3A to 3E) is depicted in top panel image on left, with measurements as defined in Table 5 and consisting of: Top panel; the orthogonal volume, and the curvilinear/centerline volume-length. Bottom panel; the orthogonal cross-sectional area (CSA) (superior 3A, or inferior 3B–3E), its anterior-posterior distance (APDist), and lateral width (Width). Plots include measurements for male in blue filled square symbols, and for female in red shaded circle symbols. Pediatric data include linear fits for males (blue solid line) and females (red dashed line). Point estimate of modeled means and confidence intervals are plotted for adult data, and at ages 0 and 5 years respectively for males (purple) and female (magenta). The second Y-axis reflects the percent growth for males (blue, inwards tick orientation) and females (black, outwards tick orientation). Significance for sex differences at birth, age 5 years and/or adults are denoted with single, double, and triple asterisk(s) at the interval plots using the nominal α < 0.05, < 0.01, and < 0.001 levels respectively; the numeric p values are displayed in Table 6.

Results

Measurements extracted for males and females are displayed in Fig 3A–3C for each pharyngeal subregion (naso- oro- laryngo- pharynx), Fig 3D for the entire pharynx (supraglottal region), and Fig 3E for the subglottal region with sex-specific linear fits and confidence intervals at age <1 year, 5 years, and adults. The plots also include a second y-axis depicting the percent growth of adult size. Statistical analysis results are also summarized numerically in Table 6. Significance at the .05, < .01, < .001, and Bonferroni corrected < .0004 levels are marked with one, two, three and four asterisks respectively in Table 6.

a. Age effect

As expected, likelihood ratio test results confirmed that all airway measurements for total and subregion volume, length (including VTL), and width (lateral) exhibited statistically significant growth in size during the first 5 years of life (<1 year to 5 years) (Table 6). Also, four of the five CSAs examined displayed significant growth in size, except for the CSA at the level of the nasopharynx (p = 0.0868). Similarly, the linear anterior-posterior distance (APDist) measurements displayed significant growth except for the two measurements at the level of the nasopharynx and the laryngopharynx (p = 0.2087 and 0.3702 respectively). Limited growth was noted for average piriform sinus length (AveragePSLength; p = 0.0768), but growth in Velum Length was highly significant (p < .0001).

Compared to the mature adult airway, both male and female pediatric upper airway dimensions by age 5 years were significantly smaller (with higher percent growth in females than males as discussed below). However, one exception was the oropharynx anterior-posterior distance (APDist) and to some extent oropharynx width, where by age 5 years, children had essentially attained their adult size (see Fig 3B).

b. Sex effect

The Wald test performed on pediatric data indicated that only subglottal length (SubglottalL) showed significance (p = .0106) at age <1 year with females’ mean length being longer than males (see Table 7 mean (s.d.); M = 8.06 (3.02) mm; F = 10.56 (3.17) mm). No statistical significance was detected for any other variable at age group <1 year or age group 5 years. However, sex differences–though not significant–were noted at age group <1 year in the volume-length of the nasopharynx (p = 0.08) and the oropharynx (p = 0.0635) subregions with females having longer measurement than males; also, differences in nasopharynx APDist (p = 0.0793) and glottis area (p = 0.0598) with males having larger measurements than females. By age 5 years, insignificant differences in the volume-length of the oropharynx subregion was noted (p = 0.0972) with females having slightly longer measurements than males (see Table 7; M = 10.43 (3.57) mm; F = 14.28 (7.21) mm).

Table 7

Age-specific mean (standard deviation) of the different measurement types for each pharyngeal region.

Types	Variables	Sex	Age Groups (Age range (yr;mos))
			<1(00;00–00;11)	1(01;00–01;11)	2(02;00–02;11)	3(03;00–03;11)	4(04;00–04;11)	5(20;00–30;00)
Volume (mm 3 )	Nasopharynx	M	1027.75 (681.71)	1082.45 (521.43)	2184.46 (879.28)	1709.06 (629.70)	2097.12 (699.82)	5372.79 (1543.29)
		F	1022.04 (853.53)	1488.75 (335.17)	2100.80 (1073.60)	1456.65 (640.11)	2232.27 (588.07)	4483.83 (1516.60)
	Oropharynx	M	380.64 (140.79)	353.09 (162.78)	727.05 (335.33)	1021.23 (751.86)	754.96 (614.13)	8025.12 (2397.84)
		F	299.93 (250.23)	482.34 (304.17)	613.47 (280.78)	796.34 (266.50)	1359.76 (1123.66)	5179.36 (1811.89)
	Laryngopharynx	M	412.27 (154.59)	333.90 (96.72)	607.52 (405.70)	866.48 (380.82)	868.67 (317.73)	6588.71 (2894.12)
		F	328.35 (179.26)	571.61 (174.32)	666.87 (321.22)	873.74 (489.16)	1009.18 (285.87)	4262.30 (681.45)
	PharynxVolume	M	1820.66 (823.81)	1769.44 (621.23)	3519.03 (1227.26)	3596.77 (1582.71)	3720.76 (1282.02)	19986.62 (3401.72)
		F	1650.32 (963.36)	2542.70 (573.50)	3381.15 (1451.68)	3552.66 (1430.01)	4601.21 (881.81)	13925.49 (3407.88)
	Subglottal	M	181.69 (88.19)	231.53 (67.30)	316.81 (93.86)	396.33 (133.49)	617.19 (203.74)	6023.72 (1161.30)
		F	179.72 (47.69)	264.70 (34.89)	321.03 (80.81)	526.63 (107.58)	533.77 (93.47)	3572.99 (979.31)
Volume Segment Length (mm)	NasopharynxL	M	22.19 (4.35)	23.12 (5.24)	26.32 (5.36)	27.61 (5.13)	27.88 (3.48)	36.13 (4.16)
		F	22.39 (4.64)	30.18 (2.29)	29.50 (4.07)	26.26 (3.37)	29.01 (2.32)	33.03 (4.65)
	OropharynxL	M	6.71 (1.87)	7.95 (2.02)	9.34 (4.08)	9.58 (5.01)	10.43 (3.57)	38.57 (8.49)
		F	7.57 (6.03)	5.72 (3.70)	6.51 (1.57)	11.32 (3.32)	14.28 (7.21)	30.79 (6.18)
	LaryngopharynxL	M	13.52 (2.14)	14.81 (1.67)	14.14 (2.25)	15.65 (1.40)	16.56 (1.75)	27.52 (2.79)
		F	13.79 (2.92)	14.33 (2.28)	15.31 (1.62)	15.94 (2.14)	16.53 (2.69)	22.08 (1.81)
	PharynxLength	M	42.42 (5.27)	45.75 (5.42)	49.09 (4.24)	52.85 (6.31)	54.80 (4.26)	101.83 (6.48)
		F	43.59 (6.95)	50.22 (3.20)	51.32 (3.15)	53.53 (4.60)	59.81 (4.55)	84.87 (6.16)
	SubglottalL	M	8.06 (3.02)	9.74 (1.85)	12.41 (3.18)	10.96 (2.15)	14.20 (2.62)	28.17 (2.69)
		F	10.56 (3.17)	11.97 (1.97)	11.68 (2.19)	14.40 (1.91)	13.68 (2.01)	26.65 (3.68)
Area (mm 2 )	NasopharynxArea	M	59.01 (22.48)	58.16 (24.12)	93.06 (39.75)	60.69 (31.72)	89.03 (41.34)	256.15 (92.03)
		F	47.54 (20.78)	52.03 (18.20)	71.78 (31.60)	67.30 (34.31)	99.45 (40.13)	227.33 (69.79)
	OropharynxArea	M	67.66 (22.94)	51.51 (23.51)	118.30 (54.47)	94.23 (51.27)	87.63 (21.59)	176.34 (42.19)
		F	51.52 (32.64)	96.06 (13.02)	79.93 (40.44)	74.60 (19.79)	91.80 (14.02)	145.07 (41.00)
	LaryngopharynxArea	M	62.05 (22.29)	53.06 (19.45)	93.10 (44.07)	117.89 (49.90)	81.66 (33.04)	355.86 (78.68)
		F	55.90 (24.64)	92.19 (12.97)	103.54 (31.96)	93.93 (28.32)	123.42 (53.04)	278.90 (39.83)
	GlottisArea	M	19.29 (5.79)	14.66 (4.33)	18.23 (9.49)	30.01 (8.85)	42.59 (12.77)	128.42 (53.84)
		F	11.58 (4.61)	11.88 (10.45)	21.27 (5.55)	30.63 (7.69)	36.82 (5.11)	80.05 (27.61)
	TracheaArea	M	27.24 (10.18)	31.82 (5.90)	34.70 (7.73)	42.59 (6.76)	47.83 (6.40)	238.74 (31.36)
		F	28.22 (10.35)	30.88 (3.93)	31.71 (4.79)	45.02 (7.12)	43.96 (5.21)	152.23 (22.33)
Anterior-Posterior Distance (mm)	NasopharynxAPDist	M	5.40 (2.21)	6.25 (2.25)	7.94 (3.11)	4.83 (1.69)	6.68 (2.00)	14.62 (4.34)
		F	4.35 (1.22)	4.11 (0.87)	4.43 (1.73)	5.41 (2.37)	6.94 (3.49)	12.77 (3.50)
	OropharynxAPDist	M	7.02 (2.85)	7.59 (3.33)	9.92 (2.74)	9.81 (4.82)	9.63 (2.78)	13.28 (3.44)
		F	7.11 (1.53)	9.39 (1.60)	10.11 (1.86)	11.77 (2.36)	10.99 (3.71)	11.95 (3.33)
	LaryngopharynxAPDist	M	7.65 (1.47)	7.08 (0.74)	9.20 (2.96)	9.26 (2.20)	7.81 (1.23)	14.63 (1.80)
		F	7.52 (1.37)	9.63 (1.95)	10.32 (2.58)	8.92 (1.23)	9.04 (2.85)	12.83 (2.03)
	GlottisAPDist	M	7.09 (1.27)	6.08 (0.61)	7.03 (1.68)	8.47 (0.90)	9.60 (1.26)	19.98 (1.60)
		F	5.41 (1.20)	6.60 (1.26)	7.72 (0.69)	8.82 (0.60)	9.08 (0.55)	14.44 (1.68)
	TracheaAPDist	M	5.79 (1.07)	6.16 (0.85)	6.36 (1.23)	7.07 (0.69)	7.06 (0.72)	16.01 (1.21)
		F	5.75 (1.19)	5.45 (0.46)	5.74 (0.70)	6.81 (0.98)	6.40 (0.52)	12.33 (0.94)
Width (mm)	NasopharynxWidth	M	12.54 (4.58)	12.16 (3.97)	15.64 (2.86)	14.66 (4.48)	15.73 (3.65)	22.50 (3.02)
		F	10.90 (3.06)	12.84 (3.02)	14.23 (2.90)	13.82 (3.71)	16.99 (3.18)	23.27 (3.40)
	OropharynxWidth	M	10.06 (3.98)	8.09 (3.56)	13.82 (2.63)	12.90 (4.34)	13.51 (3.42)	19.77 (2.85)
		F	7.60 (2.29)	12.26 (1.51)	12.79 (4.53)	10.73 (4.15)	12.76 (1.20)	15.72 (2.15)
	LaryngopharynxWidth	M	8.43 (1.55)	6.29 (1.49)	8.26 (1.43)	8.41 (1.48)	10.67 (2.04)	21.08 (2.04)
		F	5.96 (1.56)	7.44 (1.65)	8.73 (1.50)	10.01 (1.24)	10.58 (2.04)	18.52 (2.68)
	GlottisWidth	M	3.60 (0.79)	3.17 (0.94)	3.58 (1.28)	4.95 (1.54)	6.10 (1.21)	9.13 (3.02)
		F	2.78 (0.47)	2.45 (1.81)	3.65 (0.67)	4.90 (1.01)	5.66 (0.76)	7.57 (2.50)
	TracheaWidth	M	5.62 (1.37)	6.23 (0.49)	6.74 (0.96)	7.31 (0.83)	8.10 (0.49)	17.66 (1.18)
		F	5.57 (1.00)	6.63 (0.52)	6.63 (0.85)	7.76 (0.41)	7.89 (0.63)	14.72 (1.27)
VT Additional Length (mm)	VTLength	M	80.97 (5.11)	87.12 (3.85)	94.53 (6.37)	98.09 (6.00)	103.03 (5.02)	161.06 (7.03)
		F	78.63 (7.79)	89.66 (4.34)	93.46 (2.80)	100.11 (5.99)	103.86 (6.34)	141.14 (6.53)
	VelumLength	M	26.49 (1.46)	26.96 (2.40)	30.32 (3.54)	29.86 (3.32)	32.97 (1.15)	42.60 (4.21)
		F	24.46 (1.72)	29.62 (1.73)	30.92 (2.29)	30.65 (2.01)	32.51 (0.83)	39.00 (2.67)
	PSLengthLeft	M	6.87 (0.85)	6.63 (2.93)	8.41 (3.31)	8.35 (2.33)	8.90 (2.99)	17.96 (5.27)
		F	5.35 (3.99)	10.73 (3.25)	8.89 (4.41)	8.46 (3.01)	8.26 (3.61)	14.43 (3.16)
	PSLengthRight	M	6.84 (1.50)	5.69 (1.55)	8.54 (3.38)	10.49 (1.85)	8.76 (3.06)	20.71 (3.23)
		F	6.47 (4.12)	10.93 (3.63)	7.65 (4.70)	9.30 (2.87)	9.34 (2.71)	15.81 (2.55)
	AveragePSLength	M	6.53 (1.21)	5.75 (1.69)	8.47 (3.31)	9.26 (1.70)	8.82 (2.82)	19.33 (3.86)
		F	5.91 (4.04)	10.83 (3.44)	8.27 (4.53)	8.58 (3.04)	8.87 (2.76)	15.12 (2.61)

Age-specific mean (standard deviation) of the different measurement types for each of the variables examined with M/F denoting the average Male/Female measurements. Age groups as specified in Table 3.

As for adults, sexual dimorphism was highly significant for overall pharynx length and VTL (p < .0001; M = 161.06 (7.03) mm; F = 141.14 (6.53) mm), and significant for pharynx volume (p < .05; M = 53.73 (16.43) cm3; F = (44.84 (15.17) cm3), with males having larger measurements than females (see Tables 6 and 7). Sexual dimorphism was also present in laryngopharynx length (p < .0004) and volume (p < .05), as well as subglottal volume (p < .0004), with highly significant differences in Glottis APDist (p < .0004) superiorly and CSA at tracheal ring inferiorly (p < .0004; including differences in tracheal APDist (p < .0004) and width p = .0005) with males having larger measurements than females.

c. Percent growth

As displayed in second y-axes of Fig 3 plots, and listed in Table 6 (final column), percent growth assessment based on modeled point estimates, revealed that by age 5 years, female upper airway measurements were closer than male measurements to the adult mature size in 26 out of the 30 upper airway variables examined. (See the tabulated average data per age group in Table 7.) Female data revealed 9 out of 30 variables to have reached over 50% of adult size nasopharynx volume-length (64.75%), nasopharynx width (50.69%), oropharynx area (53.18%), oropharynx APDist (120.23%), oropharynx width (55.91%), glottis area (74.77%), glottis APDist (57.59%), glottis width (92.05%), and VTL (53.38%). In contrast, males had only 5 of the 30 variables reach over 50% of the adult mature size: nasopharynx volume (51.96%), nasopharynx volume-length (62.14%), oropharynx APDist (62.44%), oropharynx width (56.20%), and glottis width (57.58%). The only four measurements where both male and female growth reached over 50% of their respective adult size were: nasopharynx length, oropharynx APDist, oropharynx width, and glottis width.

Discussion

This study addresses a void in normative data on the upper airway during the first 5 years of life. After developing a protocol that controls for variables that can affect measurement accuracy (e.g., head position, sedation, threshold for airway segmentation), CT studies from 61 typically developing pediatric and 17 adults were used to quantify the multidimensional growth of the airway systematically with respect to age and sex. Our findings are novel in that, to our knowledge, this is the first study that examines the birth to 5 years age range, as compared to adults, using a comprehensive set of 2D and 3D measurements from the choanae to below the cricoid ring (opening to the trachea), including: supra- and sub-glottal volume and length, naso-oro-laryngo-pharynx subregion volume and length, each subregion’s superior and inferior CSA, and their antero-posterior and transverse/width distances. Additionally, the data were collected using a protocol that included a well-defined and established threshold for airway segmentation, and a semi-automatic centerline that we developed for the extraction of accurate measurements to quantify the upper airway using the natural anatomic orientation of airflow for respiration and speech production.

The use of a centerline, as an added methodological consideration, is critical for obtaining accurate measurements of the airway. As summarized in Table 1, two studies [35, 49] have used a centerline to quantify the airway, but Chiang et al. [35] is the only study to date that performed a centerline-based technique to obtain quantitative data on the growth and development of the airway. However, their measurements stopped at the level of the epiglottis, and they did not include the pre-pubertal age group.

The upper airway subregion dimensions are sensitive to altered head and tongue posture, particularly for 3D assessment but also for 2D measurements as Gurani et al. [54] point out. Given the need for a valid method to classify head position of imaging studies, our laboratory first developed a reliable protocol that uses 14 landmarks to account for both head and neck positions [68]. We therefore first employed this protocol for the selection of cases with a neutral head position for inclusion in this study, then applied the centerline protocol.

Given all the methodological considerations we accounted for, the attrition rate of cases included in this study from the imaging studies available in our database was high. We retained only 19% of the cases reviewed. Given this rigorous approach to control for positioning and other potential confounders, we expect our findings to reliably reflect typical airway growth. Furthermore, we anticipate that the inclusion of additional cases in future studies, using the above-described airway data extraction protocol, will further strengthen present findings and observations.

a. Age differences

Our findings reflect persistent positive increase in size for all variables examined during the first 5 years of life for all measurements in all subregions as displayed in Fig 3, with the means per age group summarized in Table 7. The age effect of the likelihood ratio test confirmed the significant growth in size for 24 out of the 30 variables studied with some variables displaying more rapid and extensive growth than others. Such findings provide insight on the proportional and relational growth of upper airway dimensions with age during anatomic restructuring (e.g., hyo-laryngeal descent).

As expected, all pediatric airway dimensions were substantially smaller than adult dimensions except oropharynx APDist and to some extent oropharynx width (see Fig 3). Abramson et al. [20] also found that volume, CSA and transverse measurements–but not AP dimensions–of the pediatric naso-oropharynx airway were significantly smaller than adult measurements. This rapid and early maturation in oropharynx APDist dimension corresponds to our previous research findings where growth of oral structures in the horizontal plane, in line with neural growth, achieved most of their growth towards the adults size by age 5 years [23, 66, 71]. Alternatively, though unlikely, it is possible that hypertrophy of lingual tonsils interfered with lumen APDist measurements. As for growth in width, while not as fully developed as AP dimensions by age 5 years, oropharynx lateral dimensions had reached over 55% of the adult, suggesting that growth in this subregion undergoes a combination of neural and general growth types, which Scammon had noted is present in the growth of structures in the neck region.

Volume, volume-length and width measurements increased with age for all pharyngeal subregions, consistent with prior studies on upper airway development in infant and pre-pubertal children [20, 24]. The CSA measurements in this study were extracted from anatomical landmarks representing the superior and inferior borders of the pharyngeal subregions. This is in contrast to the typical approach of measuring minimum or maximum CSA to examine sites of constriction for assessment of patients’ risk for OSAS [39, 40, 42, 43, 45, 49, 81]. Furthermore, we used oblique planes–orthogonal to the centerline–which cannot be compared to studies that used the axial plane, as in most of the above listed studies.

A factor that further complicates comparisons, including within-study cross-sectional comparisons, is the hypertrophy of tonsils in young children that follow a lymphoid growth type. In particular, nasopharyngeal tonsils referred to as adenoids, where hypertrophy is the highest in 4-6-year-old children [82]. Keeping these issues in mind, among the five CSA measurements in this study, the nasopharynx region was the only site that did not have a significant age-effect. Similarly, the APDist measurements in the nasopharynx and laryngopharynx were the only sites that did not have significant age effect between the ages <1 year and 5 years. Such findings could be attributed to adenoid hypertrophy, typically occurring between the ages 2–6 years, that diminish airway dimensions [83-85]. The decrease of mean nasopharynx CSA and APDist measurements per age group can be noted in Table 7, with changes most evident between the ages 2-to-4 years in this study.

In all age groups, the APDist dimensions at the nasopharynx, oropharynx and trachea (Fig 1A, 1B and 1E) were smaller than the width/transverse measurements, but larger than width/transverse measurements at the glottis (Fig 1D, and Table 7). Such findings are in line with Abramson et al. [20] who reported significant upper airway growth along the transverse dimension with age, where the airway becomes more elliptical in shape. Similarly, Machata et al. [26] using MRI studies of children ages 0–6 years, reported smaller anteroposterior dimensions than transverse dimensions for all ages at the level of the soft palate, the base of the tongue, and the tip of the epiglottis. In contrast, the laryngopharynx APDist dimensions (Fig 1C), were larger than width/transverse measurements from birth to age 3 years, but became smaller than width/transverse measurements at age 3 and beyond, which likely contributed to the absence of age effect for APDist.

Changes in APDist versus width dimensions could be attributed to the cartilaginous composition of the larynx. The laryngopharyngeal cross section in this study was designed to capture its surrounding structures—the aryepiglottic folds on each side, the laryngeal vestibule anteriorly, and posteriorly by arytenoid cartilages, corniculate cartilages and the interarytenoid fold. Cartilage ossification is usually not observed until past age 20 years, and the pediatric laryngopharynx region is often described to be “featureless” and difficult to assess using commonly acquired medical images [79]. We employed an established method for airway segmentation that uses image-specific airway thresholds, and therefore are confident that our data is reflective of airway development. Since we used landmarks on the aryepiglottic folds that connect to the piriform sinuses on each side of the cavity, aditus of larynx, the laryngopharynx width/transverse measurement in this study excluded the piriform sinuses (see Fig 1C). The significant age effect along this transverse dimension is therefore truly reflective of the strong lateral growth in the laryngopharyngeal region.

The CSA and transverse width measurements at the level of the glottis are smaller than the area and width at the first tracheal ring immediately inferior to the cricoid (i.e., subglottal region); however, the average APDist measurement of the glottis is larger than the APDist dimension at the level of the first tracheal ring at ages <1 year, 5 years, and in adults. This finding is consistent with Luscan et al. [86], who concluded that “the cricoid has a round shape regardless of the child’s age.” Indeed, the mean APDist and width measurements were very close or similar for all ages at the level of the first tracheal ring proximal to the inferior border of the cricoid, and comparable to the cricoid outlet’s (to trachea) anteroposterior and transverse interior diameters of Liu et al. [87]. While growth trends were comparable, our measurements were closer to those of Liu et al [87] than to those of Luscan et al. [86], and indicate the importance of methodological considerations, including having well-defined data extraction protocols such as the determination of an appropriate threshold level (HU) to segment the airway.

Growth in VTLi was significant during the first 5 years of life, confirming an increase of about 3 cm, which is consistent with VTL findings to date [23] and reflects that this modified measure captured growth in both the oral and pharyngeal portions of the VT. The measurements in this study were smaller than what has been reported previously, which is to be expected given the modified length measure had an onset at the posterior margin of the incisors in lieu of the typical anterior margin of the lips. Findings of a significant age effect on velum length were comparable to values reported by Perry et al. [88] and Yi et al. [27]. Closer examination of the developmental data on pharyngeal length and pharyngeal volume revealed a close relationship particularly after about age 2. Before age 2, the growth rate was slightly more pronounced in length than in volume, likely due to the drastic anatomic restructuring of the skeletal framework in the region of the pharyngeal cavity, including hyo-laryngeal descent, and rapid neural growth in length of the second cervical spine (C2) [68].

As for the piriform sinuses, our findings revealed a borderline average PSLength age effect (p = .077) with average measurements per age group summarized in Table 7. The pediatric average PSLength measurements ranged from 6 to 9 mm at ages <1 year to 5 years, whereas the adult average PSLength measurements ranged from 1.5 to 1.9 cm. While the development of piriform sinuses has not been examined to date, and therefore comparative measurements were not available, adult PSLength measurements were comparable to the 1.6 to 2 cm piriform sinus depth measurements of Dang and Honda [89-91]. This similarity was despite the fact that our PSLength measurement extended from the most inferior aspect of the piriform sinuses to the midpoint of the aryepiglottic folds, which is beyond the arytenoid apex plane used by Dang and Honda [89]. This could be in part due to differences in imaging modality used (CT vs MRI) and/or segmentation thresholding levels used to obtain reliable airway measures, particularly given the small size of this region of interest. Additional factors include methodological differences (oblique vs, axial plane) in obtaining measurements, as well as the height of participants, which has been shown to be related to vocal tract length [24] and pharyngeal dimensions [47]. The piriform sinuses play an important role during swallowing by diverting liquids around the aditus of the larynx and into the esophagus. They also affect speech acoustics and attenuate the vocal tract resonant frequencies in adults [89, 92–94] by an estimated range of 5% of formant frequencies [89]. Thus, detailed developmental data on the piriform sinuses would provide needed normative data and could help provide insight on pediatric dysphagia. Furthermore, such data can be used to implement modeling [95] to systematically examine the effect of the piriform sinuses on the resonances of the developing vocal tract, particularly given the intriguing findings that formant frequencies reportedly remain stable during the first 24-to-36 months of life [14, 15, 96], despite documented increases in vocal tract length [2, 23, 24].

In summary, the upper airway dimensions reveal persistent growth during the first 5 years of life, with some dimensions growing at a faster pace than others. Growth in the vertical and transverse/lateral dimensions are more pronounced than growth in the AP dimension.

b. Sex differences

As depicted in Table 6, sexual dimorphism was present in a number of supra- and sub-glottal variables in adults. However, while there was evidence towards sexual dimorphism in all three supra-glottal regions for a number of variables at age <1 year (specifically, larger nasopharynxL and oropharynxL in females; also, larger nasopharynx APDist and glottis area in males), with the larger oropharynxL persisting in females at age 5 years, none of these supra-glottal or pharyngeal variables were significant in children.

As for the subglottal region, only subglottal volume-length displayed significant sexual dimorphism at age <1 year (with males shorter than females), but not at age 5 years or in adults (see Table 6). To our knowledge, this specific subglottal volume-length measurement has not been examined, despite its importance in procedures like tracheotomy [97], and laryngotracheal infections/diseases including SIDS [98] where the incidence is higher in males [99]. However, two studies have performed distance measurement in this region, specifically anterior commissure to first tracheal ring [100], and vocal folds to the cricoid cartilage [101]. Contrary to present findings, Khadivi and colleagues, who used laryngoscopy to collect subglottal length data from 82 adults (57 males and 25 females), documented significant sexual dimorphism. While our measurements for pediatric subglottal length are comparable to the normative values reported by Sirisopana and colleagues, from the CT scans of 56 children (29 males, 27 females), they unfortunately neither assessed for sex differences, nor reported sex-specific measurements given their primary focus on tracheal tube design.

Despite methodological differences, the absence of sexual dimorphism in pediatric upper airway data for most of the measurements analyzed in this study was consistent with past studies [22, 24, 29]. Barbier et al. [29] did not find sex difference in pre-pubertal data but suggested that sexual dimorphism in VTL emerged during puberty. Among studies reporting regional upper airway normative data for the age range between 0-to-5 years, Abramson et al. [20] found no difference in naso-pharyngeal airway size or shape between the sexes in children, but reported longer airway length in post-pubertal males. Jeans et al. [83], using lateral cephalometric radiographs to study the nasopharyngeal airway area–comparable to our nasopharyngeal region–, however, reported mild decreases in nasopharyngeal area in both 3-to-5-year-old males and 3-to-6-year-old females. Sex-differences in the same region using an anteroposterior distance measure in the midsagittal plane of medical imaging studies (MRI & CT), referred to as oropharyngeal-width, have similarly been noted to display evidence, albeit not significant, towards sexual dimorphism in 3-to-4-year-old children with males having larger width measurements [71]. Linder-Aronson et al. [85] noted that nasopharynx airway depth/AP dimension in males were consistently larger than females throughout ages 3-to-16 years. Sexual dimorphism in the pharyngeal portion of the VT in ages 8-to-19 years, has also been reported by Vorperian et al. [71], with the vertical nasopharyngeal length being longer in females than males and the vertical posterior cavity length being longer in males than females. In contrast, Yi et al. [27] reported no sex differences in any of their linear dimensions using MRI in infants and children up to 72 months. Rommel et al. [34], who used curvilinear length drawn on 2D X-ray images to assess naso-oropharynx segments, found no sex difference in children as young as 0-to-4 years. Griscom [102] similarly found no significant sex differences in trachea dimensions until late in adolescence. Definitive prepubertal sexual dimorphism of the pharynx thus cannot be confirmed with studies available to date.

Detecting sex differences is a difficult task given the critical methodological considerations outlined in our methods section and the importance of having a large number of participants per age group. Statistical analysis methods can overcome differences in growth rate between males and females, such as implementing continuous-window comparisons across age [e.g., 71]. This latter approach was particularly effective in unveiling sexual dimorphism that does not persist during the course of development, since growth in females outpaces males during early development, but growth in males begins to outpace females during the peripubertal period, with sexual dimorphism emerging during puberty.

Sexual dimorphism in adults, however, was mostly present and aligned with research findings to date in pharynx volume, pharynx length [47], VTL [71], velum length [103], glottis APDist [47], subglottal volume [102], and tracheal dimensions [86]. The finding that subglottal volume-length, the only measurement that displayed significant sexual dimorphism at age <1 year old was not sexually dimorphic at age 5 years is not surprising, given growth rate differences in males versus females as noted above. However, the absence of differences in adults is likely due to both methodological differences and the limited number of adults in this study (n = 17), given our stringent inclusion criteria.

Although present findings revealed significant prepubertal sex-differences only in the subglottal region, findings in this study and others as noted above, in the naso-oro-pharyngeal region, provide sufficient justification to further examine this issue using a larger number of cases particularly given the above noted critical methodological considerations, and the nuance of growth rate differences between the two sexes. Such a conclusion is further supported by auditory-perceptual and acoustic findings where Bloom et al. [104] reported that adults accurately identified 3-month old infants’ vocalizations as boy vs girl. The only acoustic difference was the feature of nasality with girls’ vocalizations being more nasal than boys. Furthermore, several studies have reported sex differences in vowel formants (i.e., vocal tract resonant frequencies) in children as young as 3 or 4 years of age [105-107].

The growth pattern of anatomical structures in the craniofacial and upper airway region are known to be non-uniform with the primary growth types in this region being neural, general/somatic, and lymphoid. Since we did not include adenoid or tonsil measurements, we will limit this discussion to the first two types. Both neural and general growth types display rapid growth during the first few years of life. Scammon [80] summarized schematically, however, that by age 5 years, the percent of the adult mature size reached was drastically different for the neural (~80%) versus the general (< 40%) growth types. He also noted that growth in the neck region could be a combination of both neural and general primary growth types [80]. The final column of Table 6 presents the model-based point-estimates of percent growth of the adult mature size by age 5 years for all 30 variables examined in this study.

The general finding that by age 5 years, female upper airway dimensions were larger than males is not surprising since typically females have a faster growth rate during childhood and reach the adult mature size sooner than males. Based on findings to date, structures in the upper airway were expected to follow a mostly somatic growth type or a composite growth of somatic and neural growth types [23, 108, 109]. Despite differences in VTL versus VTLi measurements (where the onset of the former at the anterior margin of the lips and the latter is at the posterior margin of the incisors), the general growth findings in this study are in line with the expected growth trend indicating that VTLi growth type is predominantly hybrid somatic/neural in females (53.38%) and somatic in males (36.42%) [23]. Similarly, pharynx volume, pharynx volume-length and all other pharyngeal subdivision volume and volume-length results except for the nasopharynx subregion confirmed the predominant somatic growth types at age 5 years, in line with the reported growth patterns for pharyngeal cavity length and VT vertical data [23]. As for the nasopharyngeal region, the expected hybrid somatic/neural or combination growth type, was indeed reflected in present findings where both male and female nasopharynx volume and volume-length measurements ranged between 47% and 65% of adult size by age 5 years. The oropharynx AP dimension exceeded 60% for males and 85% for females, suggestive of a more pronounced hybrid somatic/neural growth type for the males and neural growth type for the females. These findings are in line with Vorperian et al. [23] where structures in oral region in the horizontal plane reached maturation earlier than structures in pharyngeal region in the vertical plane. Percent growth analysis (in Table 6, last column) also reflect the presence of sex-specific differences in the combination of growth types for upper airway structures. Since multiple factors contribute to growth, it is possible that such sex-specific differences further contributed to the difficulty in detecting sexual dimorphism during early childhood.

Conclusions, study limitations, and future direction for research

This study, using CT studies, provides data quantifying the 3D growth of the upper airway with minimal methodological concerns for an age group with scant normative data. Findings confirm persistent growth of the upper airway during the first 5 years of life with growth in the vertical and transverse/lateral dimensions having a faster pace and greater prominence than growth of anteroposterior dimension. Findings also reveal that at age 5 years, females have larger airway dimensions than males. Such findings confirm the importance of studying sex-specific developmental changes of the upper airway in 3D. A better understanding of pharyngeal functions and disorders will require further, more detailed examination of the developmental changes in pharyngeal length versus volume and the piriform sinuses, as well as re-examination of prepubertal sex-differences.

Our painstaking efforts to optimize accurate and reliable upper airway measurements limited the sample size. Specifically, the attrition rate of only retaining retrospective imaging studies with a neutral head position was high given the head position protocol we applied [68]. Also, despite the imaging protocol at the University Hospital to maintain the head in midline, it is likely that some rotation was present. However, given that the airway is functional during imaging (breathing and swallowing), it is difficult to determine source of asymmetry as noted with the piriform sinuses, where we resorted to averaging them in this study. Furthermore, airway anatomic measurements can be affected by breathing phase where significant effect of breathing phase in the oro-laryngo-pharyngeal region has been reported [34]. Such concerns, including the above suggested assessments, could be addressed by replicating this study using a larger sample size, and ideally increasing the age range to cover the entire developmental period particularly ages 5-to-20 years.

Aside from using a larger sample, having relevant demographic information such as height, weight, and race would be valuable to include. Given the retrospective nature of this study, relevant demographic information was not available for all cases. Our imaging database however, was representative of regional Dane County demographics, with growth between the 10th and 95th percentiles. We believe our present findings are representative of typical growth since the natural variability in craniofacial dimensions within individual races is related to the natural variability or variations within the different racial/ethnic groups [110].

Aside from the use of 3D anatomic landmarks and establishing standardized procedures to minimize, if not eliminate, methodological limitations on measurement accuracy and reliability, it is also imperative to establish standardized definitions of the pharyngeal subregions using well-defined anatomical boundaries. This will facilitate accurate representation and comparison of the developmental morphology of the aerodigestive and vocal tract, within and across disciplines, including comparison between different imaging modalities. This will undoubtedly enrich our understanding of the growth of a region that serves multiple life-functions since each modality has its strengths and limitations. For example, while MRI provides more accurate information on the growth of lymphoid tissue in the pharyngeal region, CBCT could address concerns on the effect of body position and gravity on soft-tissue structures for obtaining reliable airway dimensions. The feasibility of comparing developmental findings across disciplines, including the individual and relational growth of structures that provide the skeletal framework of the aerodigestive and vocal tract, will facilitate the understanding of upper airway pathophysiology, improve surgical planning such as estimation of laryngeal mask airway size or endotracheal tube diameter, evaluation of pharyngeal collapsibility in early childhood in the assessment of OSAS, and other upper airway anomalies including swallowing difficulties and speech disorders. Such information would also facilitate the advancement of developmental models to assess various typical and atypical functions related to airflow, swallowing, and speech production.

Data used in this study.

Variable measurements as described in Table 5.

(CSV)

Click here for additional data file.