Effect of extraction treatment on upper airway dimensions in patients with bimaxillary skeletal protrusion relative to their vertical skeletal pattern.

OBJECTIVE: To investigate dimensional changes in regional pharyngeal airway spaces after premolar extraction in bimaxillary skeletal protrusion (BSP) patients according to vertical skeletal pattern, and to further identify dentoskeletal risk factors to predict posttreatment pharyngeal changes.
METHODS: Fifty-five adults showing BSP treated with microimplant anchorage after four premolar extractions were included in this retrospective study. The subjects were divided into two groups according to the mandibular plane steepness: hyperdivergent (Frankfort horizontal plane to mandibular plane [FH-MP] ≥ 30) and nonhyperdivergent groups (FH-MP < 30). The control group consisted of 20 untreated adults with skeletal Class I normodivergent pattern and favorable profile. Treatment changes in cephalometric variables were evaluated and compared. The association between posttreatment changes in the dentoskeletal and upper airway variables were analyzed using linear regression analysis.
RESULTS: The BSP patients showed no significant decrease in the pharyngeal dimensions to the lower level in comparison with controls, except for middle airway space (MAS, p < 0.01). The upper airway variable representing greater decrease in the hyperdivergent group than in the nonhyperdivergent group was the MAS (p < 0.01). Posttreatment changes in FH-MP had negative correlation with changes in MAS (β = -0.42, p < 0.01) and inferior airway space (β = -0.52, p < 0.01) as a result of multivariable regression analysis adjusted for sagittal skeletal relationship.
CONCLUSIONS: Decreased pharyngeal dimensions after treatment in BSP patients showed no significant difference from the normal range of pharyngeal dimensions. However, the glossopharyngeal airway space may be susceptible to treatment when vertical dimension increased in hyperdivergent BSP patients.

INTRODUCTION

Bimaxillary protrusion (BMP) is defined as a condition wherein the upper and lower incisors are protrusive or proclined, leading to lip protrusion, which mostly corresponds to bimaxillary dentoalveolar protrusion.[1] In a specific sense, bimaxillary skeletal protrusion (BSP) is characterized by protrusive upper and lower jaws with or without incisor proclination.[2] Most patients with skeletal and/or dentoalveolar protrusion are treated with microimplant anchorage for maximum anterior retraction after extraction of four premolars with the aim of esthetic improvement.[3-5] In particular, BSP patients require greater retrusion of both A-point and B-point areas concomitant with incisor retraction than BMP patients. However, there may be concerns regarding the effects of a decreased oral cavity with a retruded anterior alveolus and subsequent parapharyngeal tissue displacement on the airway dimensions and respiratory function.[6,7] On the other hand, the pharyngeal airway dimensions of BSP patients are presumed to be initially large enough to maintain normal airway patency even after extraction treatment, albeit with no scientific evidence.[8]

Previous studies have investigated the pharyngeal dimensional changes after premolar extraction treatment in BMP patients. A few studies found significant constriction of the velopharyngeal space behind the soft palate, glossopharyngeal space behind the tongue base, and/or hypopharyngeal space below the epiglottis base after extraction treatment,[9,10] whereas one study reported no significant posttreatment changes at any level of the pharyngeal spaces.[11] These contradictory findings might be due to the differences in sample characteristics between a skeletal Class I BMP with a well-developed chin and skeletal Class II BMP with a retruded chin. These two patterns tend to show different initial pharyngeal morphologies and dimensions, which may respond differently to the dentoskeletal treatment changes.[12] Moreover, no previous study has differentiated the BSP from BMP, which requires a different amount and pattern of anterior retraction of teeth and apical bone, to demonstrate the posttreatment pharyngeal airway changes.

Notably, the vertical skeletal pattern may affect the treatment outcome of the pharyngeal dimension in BSP patients, because a hyperdivergent vertical pattern with steep mandibular and occlusal planes is known to be a predisposing factor for pharyngeal narrowing and collapsibility.[13-15] One unique study reported no significant difference in the pharyngeal changes between hyperdivergent and nonhyperdivergent groups in BMP patients treated with premolar extraction.[10] More importantly, an increase in the vertical dimension by clockwise mandibular rotation after treatment may contribute to pharyngeal narrowing; however, the effects of vertical skeletal changes on the pharyngeal dimensions in BSP patients has not been discussed so far.

Existing evidence does not support the notion that extraction treatment may induce a detrimental effect on respiratory function or obstructive sleep apnea (OSA).[16] Although pharyngeal narrowing may be noted after treatment in two-dimensional or three-dimensional radiographic images, most patients show the potential for functional adaptation to normal function.[17,18] However, some orthodontic patients may have risks or symptoms of OSA in relation to the hidden neuromuscular pathophysiology, which might increase the susceptibility to the pharyngeal change after orthodontic treatment.[19] Therefore, orthodontists should reflect the possible pharyngeal responses to various types of orthodontic treatments in daily practice. It is worth questioning whether the posttreatment pharyngeal airway dimensions would be stable in BSP patients despite the greater retraction of the anterior teeth and alveolar bone due to their initially larger patency than that in the normal skeletal pattern, or it might be unfavorably changed by some dentoskeletal risk factors. The objective of this retrospective study was to investigate the dimensional changes in regional pharyngeal airway spaces after premolar extraction treatment in BSP patients according to their vertical skeletal pattern. Further, any dentoskeletal risk factors that might be used to predict posttreatment pharyngeal changes in BSP patients were also evaluated.

MATERIALS AND METHODS

Subjects

Seventy-two adult patients who were diagnosed with BSP and treated at the Department of Orthodontics, Kyung Hee University Dental Hospital were identified. Fifty-five patients (10 males and 45 females) with a mean age of 23.4 ± 5.21 years (range, 18.2 to 34.9 years) who met the following inclusion criteria were included: (1) maxillary and mandibular protrusion as indicated by A-point and B-point protrusion, respectively (Nper-A ≥ 2 mm, Nper-B ≥ –3 mm); (2) arch length discrepancy ≤ 4 mm; (3) extraction of four first premolars; (4) microimplant anchorage for maximum anterior retraction; (5) pre- and posttreatment lateral cephalograms taken in a natural head position, and in a breath-holding position after the end of expiration to obtain reproducible pharyngeal caliber[20]; and (6) no evidence of adenoid and/or tonsillar hypertrophy in the cephalogram. The exclusion criteria were as follows: (1) chronic medication; (2) history of pharyngeal soft tissue surgery; (3) obesity with body mass index over 25 kg/m2; and (4) bidentoalveolar protrusion with no skeletal protrusion. This study was approved by the Institutional Review Board of Kyung Hee University Dental Hospital (IRB No. KH-DT19015).

The subjects with BSP were divided into two groups according to the initial mandibular plane angle (FH-MP): the hyperdivergent group (n = 27) with an FH-MP of 30° or more; and the nonhyperdivergent group (n = 28) with FH-MP less than 30°. A control group was set up to basically evaluate whether the pretreatment and posttreatment pharyngeal dimensions in the BSP group showed significant differences, respectively. The control group consisted of 20 untreated adults (8 males and 12 females; mean age, 25.1 ± 4.8 years) with skeletal Class I normovergent patterns and a favorable facial profile and incisor inclination. All subjects were treated using 0.022 × 0.028-inch (in) self-ligating brackets (Clippy-C; Tomy, Tokyo, Japan). Maximum anterior retraction was performed using sliding mechanics with 0.019 × 0.025-in stainless steel wires and microimplant anchorage on all quadrants.

Pretreatment (T0) and posttreatment (T1) lateral cephalograms were evaluated by one expert orthodontist at two-week intervals. Cephalometric measurements were consisted of 17 dentoskeletofacial and 14 upper airway variables (Figure 1). Dimensions of the pharyngeal airway were evaluated as follows (Figure 2): the nasopharyngeal space above the palatal plane was represented by PNS-Ad1 and PNS-Ad2; the velopharyngeal space behind the soft palate was represented by superior posterior airway space (SPAS); and the glossopharyngeal space behind the tongue base was described as middle airway space (MAS) and inferior airway space (IAS).

Figure 1

Dentoskeletal (A) and pharyngeal (B) cephalometric landmarks and measurements.

S, sella; Na, nasion; Po, porion; Or, orbitale; ANS, anterior nasal spine; PNS, posterior nasal spine; A, point A; B, point B; Pog, pogonion; Me, menton; Go, gonion; Sn, subnasale; UL, upper lip; LL, lower lip; Pog’; soft-tissue pogonion; Pt, pterygomaxillary point; Ad1, adenoid point 1; Ad2, adenoid point 2; Ba, basion; TT, tongue tip; H, hyoid bone; RGN, retrognathion; C3, third vertebra; P, tip of soft palate; R, midpoint of S-Ba line; Sb, belly of soft palate; Eb, base of epiglottis.

See footnote of Table 1 for the definition of measurements.

Figure 2

Cephalometric pharyngeal measurements to describe each pharyngeal airway section: PNS-Ad2 and PNS-Ad1 for the nasopharynx; superior posterior airway space (SPAS) for the velopharynx behind the soft palate; middle airway space (MAS) and inferior airway space (IAS) for the glossopharynx behind the tongue base; and vertical airway length (VAL) for the pharyngeal airway length.[11]

See Figure 1 and footnote of Table 1 for the definition of other landmarks or measurements.

Statistical analysis

To assess the error in cephalometric analysis, measurement sets were obtained twice with a two-week interval and analyzed by a paired t-test. The mean intra-observer differences were less than 0.45 mm and 0.64° in linear and angular measurements, respectively. No significant differences were found in any of the variables. Pearson correlation coefficients for all variables represented a significant range of 0.93 to 0.97, showing the reliability of the measurements.

A Shapiro–Wilk test confirmed the normal distribution of the data. Overall treatment change (ΔT) was analyzed using a paired t-test, and intergroup comparison was performed using an independent t-test. Univariable linear regression analyses were conducted to detect any correlation between posttreatment dentoskeletal changes and upper airway changes. Then, multivariable linear regression analysis was performed to determine whether the posttreatment changes in vertical skeletal divergency (FH-MP) affect the changes in upper airway variables when adjusted for the sagittal skeletal relationship (ANB angle). SPSS software, version 22.0 (IBM Corp., Armonk, NY, USA) was used, and all data were expressed as mean and standard deviation (p < 0.05).

RESULTS

Overall treatment changes in the BSP patients in comparison with control samples (Table 1)

Before treatment, BSP patients showed no significant differences in pharyngeal dimensions from the normal control group due to the large interindividual variations, except for nasopharyngeal space (PNS-Ad1, p < 0.05). More protruded tongue posture was noted in the BSP patients (TT-PTV, p < 0.001). In assessments of the dentoskeletal changes after treatment, significant reductions in SNA, SNB, Nper-A, and Nper-B (p < 0.001) were observed along with significant reductions in U1-FH and L1-MP (p < 0.001). IIA and OB (p < 0.001) significantly increased while OJ (p < 0.001) decreased.

Regarding the pharyngeal and parapharyngeal tissue changes after treatment, SPAS (p < 0.001) behind the belly of the soft palate, MAS (p < 0.001) behind the soft palate tip, and IAS (p < 0.001) at the level of the mandibular plane significantly decreased. In accordance with these findings, the soft palate inclination to the palatal plane (SPI, p < 0.001) increased and tongue position was retruded (TT-PTV, p < 0.001) and lowered (TGH, p < 0.01). Nasopharyngeal spaces and vertical airway length showed no significant changes. Notably, however, the pharyngeal airway variable that showed reduction compared to the control samples was only MAS (p < 0.01), in relation to the retroclined soft palate (p < 0.05) and the lowered tongue (p < 0.05). Nonetheless, the mean difference in MAS between the posttreatment dimensions in BSP patients (11.11 ± 2.85 mm) and the normal dimensions in the controls (12.64 ± 1.05 mm) was 1.53 mm, showing little clinical relevance.

Comparisons between hyperdivergent and nonhyperdivergent groups (Tables 2 and 3)

Before treatment (Table 2), the hyperdivergent group showed a smaller SPAS (p < 0.01) and IAS (p < 0.05) than the nonhyperdivergent group. However, the hyperdivergent group showed no significant difference in the pharyngeal dimensions at all levels from the normal controls, while nonhyperdivergent BSP patients exhibited rather larger pharyngeal dimensions than the controls: PNS-Ad1 (p < 0.05); SPAS (p < 0.01); and IAS (p < 0.05).

Table 2

Comparison of pretreatment upper airway variables between hyperdivergent and nonhyperdivergent groups

Upper airwayvariable (T0)		Normal samples(n = 20)			Hyperdivergent group(n = 27)				Nonhyperdivergent group(n = 28)				Intergroup comparison
		Mean	Standard deviation		Mean	Standard deviation	p-value†		Mean	Standard deviation	p-value†		p-value‡
Hyoid	H-RGN (mm)	33.35	2.64		35.47	4.04	0.131		36.71	5.37	0.056		0.397
	MPH (mm)	8.88	3.13		9.00	4.43	0.937		6.75	5.50	0.236		0.144
	H-PTV (mm)	−0.35	2.81		−3.47	5.28	0.171		−0.33	4.33	0.989		0.082
	H-C3Me (mm)	1.34	4.69		−2.82	4.54	0.024*		−0.02	4.88	0.434		0.053
Tongue	TGH (mm)	1.91	0.72		4.63	2.36	0.025*		2.16	2.47	0.315		0.020*
	TT-PTV (mm)	50.25	2.53		50.34	4.22	0.007		52.46	4.00	0.100		0.108
Soft palate	SPL (mm)	31.46	2.03		32.55	3.19	0.067		32.45	3.29	0.108		0.913
	SPI (°)	127.10	2.83		129.71	6.74	0.151		128.25	3.98	0.389		0.403
Pharyngeal airway	PNS-Ad1 (mm)	25.10	2.54		27.63	4.48	0.097		28.18	3.53	0.012*		0.642
	PNS-Ad2 (mm)	26.89	1.50		26.26	4.93	0.610		27.20	4.36	0.939		0.744
	SPAS (mm)	13.30	1.57		12.92	2.52	0.659		15.54	2.48	0.009**		0.001**
	MAS (mm)	12.64	1.05		12.47	2.98	0.826		13.45	2.65	0.182		0.246
	IAS (mm)	11.16	1.30		11.68	3.35	0.548		13.11	3.24	0.011*		0.023*
	VAL (mm)	63.65	8.50		64.87	3.39	0.670		63.14	3.18	0.855		0.303

*p < 0.05, **p < 0.01.

†Independent t-test was performed for the intergroup comparisons with the control group.

‡Independent t-test was performed for the intergroup comparisons between the hyperdivergent group and nonhyperdivergent group.

After treatment (Table 3), changes in the dentoskeletal and facial soft tissue variables showed no significant intergroup differences. For the upper airway variables, both groups showed reductions in SPAS, MAS, IAS in relation to the decrease in TT-PTV and increase in SPI. However, a significant difference in the pharyngeal decrease between the two groups was found only in the MAS (p < 0.01); the hyperdivergent group showed a greater reduction in MAS at the level of the soft palatal tip than the nonhyperdivergent group.

Table 3

Comparison of posttreatment changes between hyperdivergent and nonhyperdivergent groups

Treatment changes invariable (ΔT)		Hyperdivergent group (n = 27)				Nonhyperdivergent group (n = 28)				Intergroup comparison
		Mean	Standard deviation	p-value†		Mean	Standard deviation	p-value†		p-value‡
Skeletal	SNA (°)	−1.63	1.32	0.000 ***		−1.39	0.98	0.000***		0.726
	SNB (°)	−1.24	1.55	0.003 **		−0.89	0.99	0.000***		0.358
	ANB (°)	−0.39	1.85	0.511		−0.50	0.96	0.011*		0.436
	Nper-A (mm)	−1.68	1.23	0.000 ***		−1.80	0.98	0.000***		0.272
	Nper-B (mm)	−1.45	2.07	0.007 **		−1.59	1.78	0.000***		0.803
	Nper-Pog (mm)	0.26	1.39	0.419		−0.82	1.81	0.024*.		0.033*
	FH-MP (°)	−0.39	1.28	0.259		0.45	1.31	0.084		0.793
	FH-OP (°)	0.03	2.71	0.967		1.18	2.26	0.010*		0.121
	LFH (mm)	−0.11	1.59	0.777		0.39	1.34	0.105		0.346
Dental	U1-FH (°)	−11.00	4.95	0.000 ***		−12.68	3.69	0.000***		0.190
	Nper-U1 (mm)	−9.52	3.28	0.000 ***		−10.32	3.11	0.000 ***		0.165
	L1-MP (°)	−9.42	2.38	0.000 ***		−9.16	2.23	0.000***		0.704
	Nper-L1 (mm)	−7.88	2.22	0.000 ***		−7.35	2.05	0.000 ***		0.256
	IIA (°)	23.95	10.83	0.000 ***		22.98	8.20	0.000***		0.730
	OJ (mm)	−1.58	1.20	0.057		−1.18	1.61	0.001**		0.178
	OB (mm)	1.08	1.55	0.007 **		0.71	1.50	0.018*		0.424
Soft tissue	SnV-UL (mm)	−2.66	1.61	0.000 ***		−2.38	1.43	0.000***		0.530
	SnV-LL (mm)	−4.18	1.96	0.000 ***		−4.38	2.70	0.000***		0.853
	SnV-Pog’ (mm)	0.83	1.74	0.051		−0.20	2.17	0.636		0.092
Hyoid	H-RGN (mm)	−0.87	4.64	0.425		−0.80	4.03	0.301		0.960
	MPH (mm)	1.39	2.39	0.020 *		0.50	2.49	0.209		0.282
	H-PTV (mm)	0.79	4.02	0.403		−0.05	3.09	0.932		0.423
	H-C3Me (mm)	0.87	2.49	0.146		0.13	2.16	0.762		0.283
Tongue	TGH (mm)	1.50	2.55	0.019 *		0.71	2.14	0.066		0.285
	TT-PTV (mm)	−5.00	3.26	0.000 ***		−5.02	3.71	0.000***		0.987
Soft palate	SPL (mm)	1.18	2.86	0.088		1.30	2.83	0.022*		0.888
	SPI (°)	1.89	2.83	0.009 **		2.48	1.97	0.000***		0.295
Pharyngeal airway	PNS-Ad1 (mm)	−1.21	2.52	0.051		0.21	2.53	0.657		0.064
	PNS-Ad2 (mm)	−0.42	1.49	0.292		0.11	1.52	0.712		0.135
	SPAS (mm)	−1.95	1.43	0.000 ***		−1.96	0.82	0.000***		0.383
	MAS (mm)	−2.50	1.07	0.000 ***		−1.57	1.13	0.000***		0.007**
	IAS (mm)	−2.16	1.65	0.000 ***		−1.46	1.29	0.000***		0.113
	VAL (mm)	1.58	2.33	0.508		0.75	2.89	0.180		0.123

*p < 0.05, **p < 0.01, ***p < 0.001.

†Paired t-test was performed for intragroup comparisons of treatment changes (ΔT).

‡Independent t-test was performed for intergroup comparisons.

See footnote of Table 1 for the definition of measurements.

Relationship between posttreatment changes in dentoskeletal and upper airway variables (Tables 4 and 5)

In the univariable linear regression analysis (Table 4), posttreatment retrusion of the tongue tip (TT-PTV) showed a positive correlation with the retrusion amount of A-point (β = 1.25; p < 0.01), and posttreatment retroclination of the soft palate (SPI) was correlated with the retrusion amount of B-point (β = –0.49; p < 0.01) and pogonion point (β = –0.58; p < 0.01). On the other hand, posttreatment changes in the pharyngeal airway spaces showed no significant correlation with the retraction amounts of A-point and B-point as well as with the distance of upper and lower incisor retraction. Instead, reductions in SPAS, MAS, IAS were associated with the increases in ANB and FM-MP, respectively. Accordingly, multivariable linear regression analysis was performed using the mandibular plane change as an independent variable adjusted for the changes in the sagittal skeletal relationship (ANB). As a consequence, the posttreatment change in FH-MP showed a negative correlation with the change in MAS (β = –0.42; p < 0.01) and IAS (β = –0.52; p < 0.01).

Table 4

Association between posttreatment dentoskeletal changes and upper airway changes tested by univariable linear regression analysis

Variable		ΔH-RGN(mm)	ΔMPH(mm)	ΔH-PTV(mm)	ΔH-C3Me(mm)	ΔTGH(mm)	ΔTT-PTV(mm)	ΔSPL(mm)	ΔSPI(°)	ΔPNS-Ad1 (mm)	ΔPNS-Ad2 (mm)	ΔSPAS(mm)	ΔMAS(mm)	ΔIAS(mm)	ΔVAL(mm)
ΔSNA (°)	β	0.03	−0.35	−0.36	−0.05	0.26	0.81	0.20	0.46	0.32	0.47	−0.20	−0.16	−0.23	−0.10
	p	0.960	0.282	0.434	0.879	0.395	0.078	0.594	0.138	0.358	0.055	0.170	0.297	0.237	0.790
ΔSNB (°)	β	0.35	−0.60	0.27	−0.24	0.24	0.58	−0.18	−0.30	0.08	0.26	0.25	0.17	0.25	−0.35
	p	0.486	0.059	0.524	0.392	0.382	0.167	0.593	0.291	0.789	0.148	0.052	0.231	0.157	0.276
ΔANB (°)	β	−0.27	0.26	−0.46	0.16	−0.02	0.07	0.28	0.55	0.14	0.10	−0.34	−0.25	−0.36	0.22
	p	0.556	0.338	0.221	0.516	0.925	0.856	0.356	0.027*	0.612	0.531	0.003**	0.049*	0.022*	0.445
ΔNper-A (mm)	β	−0.24	−0.33	−0.02	−0.07	−0.11	1.25	−0.18	0.22	0.27	0.52	−0.21	−0.14	−0.10	−0.31
	p	0.680	0.321	0.964	0.835	0.738	0.007**	0.643	0.503	0.448	0.109	0.162	0.375	0.633	0.405
ΔNper-B (mm)	β	0.02	−0.27	0.32	0.11	0.12	0.19	−0.11	−0.49	0.02	0.17	0.09	0.07	0.11	−0.15
	p	0.956	0.167	0.238	0.532	0.531	0.490	0.614	0.006**	0.906	0.148	0.310	0.490	0.353	0.491
ΔNper-Pog (mm)	β	0.22	−0.14	0.06	0.12	−0.04	0.16	0.04	−0.58	−0.03	−0.12	0.13	0.14	0.17	−0.26
	p	0.553	0.505	0.837	0.548	0.851	0.599	0.873	0.003**	0.897	0.367	0.167	0.168	0.168	0.255
ΔFH-MP (°)	β	−0.69	0.19	−0.01	−0.20	−0.01	0.18	0.02	0.48	−0.15	0.20	−0.24	−0.45	−0.59	0.45
	p	0.156	0.505	0.986	0.454	0.975	0.658	0.963	0.077	0.609	0.256	0.053	0.001**	0.000***	0.141
ΔFH-OP (°)	β	0.03	0.24	−0.22	0.14	0.11	0.06	0.09	0.27	−0.01	0.10	−0.07	0.02	−0.07	0.13
	p	0.920	0.101	0.298	0.310	0.432	0.760	0.605	0.053	0.948	0.269	0.285	0.760	0.443	0.427
ΔLFH (mm)	β	−0.02	−0.14	−0.60	−0.45	0.12	0.09	−0.03	0.45	0.28	0.18	−0.05	−0.20	−0.24	0.16
	p	0.963	0.575	0.095	0.055	0.630	0.809	0.927	0.061	0.297	0.244	0.641	0.093	0.113	0.556
ΔNper-U1 (mm)	β	−0.02	0.05	0.10	−0.04	−0.01	0.10	−0.07	−0.06	0.11	−0.04	0.02	0.00	0.00	0.09
	p	0.898	0.551	0.424	0.588	0.918	0.420	0.489	0.491	0.231	0.462	0.647	0.926	0.949	0.351
ΔNper-L1 (mm)	β	0.14	0.01	0.00	−0.17	−0.01	0.46	−0.25	0.01	0.20	−0.04	−0.06	0.07	0.06	0.06
	p	0.612	0.943	0.989	0.251	0.927	0.043*	0.170	0.955	0.239	0.710	0.380	0.339	0.566	0.720
ΔIIA (°)	β	−0.11	−0.07	0.06	0.01	−0.03	−0.05	0.07	−0.03	−0.05	0.04	−0.02	−0.03	−0.04	−0.04
	p	0.090	0.079	0.271	0.764	0.490	0.360	0.098	0.468	0.233	0.103	0.268	0.100	0.084	0.342

Univariable linear regression analysis was performed. Δ indicates the amount of treatment-induced change.

*p < 0.05, **p < 0.01, ***p < 0.001.

See footnote of Table 1 for the definition of measurements.

DISCUSSION

The present study elucidated the influence of premolar extraction treatment with absolute anchorage on the dimensional changes in the upper airway in BSP patients, reflecting the influence of the vertical skeletal pattern. The BSP patients showed significant reductions in velopharyngeal (SPAS) and glossopharyngeal airway spaces (MAS and IAS) after treatment, with no significant changes in nasopharyngeal space (PNS-Ad1 and PNS-Ad2) and vertical airway length (Figure 2). However, the important finding is that the posttreatment pharyngeal narrowing does not indicate an increased risk for functional problems such as snoring or OSA, because the reduced airway dimensions after treatment in the BSP patients mostly corresponded to the normal dimensions in the control samples (Table 1). Uniquely, the posttreatment pharyngeal space at the level of the soft palatal tip (MAS) was smaller than the normal MAS value by 1.53 mm, which may not be clinically relevant. When the pharyngeal dimensional changes were compared between the two groups, the hyperdivergent BSP group showed a greater reduction in the MAS than the nonhyperdivergent BSP group by 0.93 mm (Table 3). In addition, the posttreatment increase in the mandibular plane angle might affect the reductions of MAS and IAS in BSP patients (Table 5). This study demonstrates that the pharyngeal dimensions of BSP patients tended to be resistant to the anteroposterior changes with maximum anterior retraction, while the glossopharyngeal space might become susceptible to the vertical increase after extraction treatment.

Table 1

Overall treatment changes in dentoskeletal and upper airway variables in BSP patients

Variable		Normal samples(n = 20)			BSP patients(n = 55)
		Mean	Standard deviation		Pretreatment (T0)				Posttreatment (T1)				Treatment change (ΔT)
					Mean	Standard deviation	p-value†		Mean	Standard deviation	p-value†		Mean	Standard deviation	p-value‡(ΔT)
Skeletal	SNA (°)	80.68	2.60		82.85	2.60	0.016***		81.36	2.57	0.585		−1.49	1.13	0.000***
	SNB (°)	78.57	2.68		78.69	2.74	0.892		77.66	2.80	0.334		−1.03	1.24	0.000***
	ANB (°)	2.11	0.75		4.16	1.93	0.000***		3.70	2.05	0.000***		−0.46	1.37	0.029*
	Nper-A (mm)	0.41	1.37		3.73	1.57	0.000***		2.28	1.54	0.001**		−1.46	1.09	0.000***
	Nper-B (mm)	−4.49	1.36		−1.52	2.95	0.000***		−3.05	2.93	0.021*		−1.53	1.88	0.000***
	Nper-Pog (mm)	−3.27	1.38		−2.18	3.65	0.113		−2.56	2.99	0.248		−0.38	1.72	0.178
	FH-MP (°)	25.85	2.48		28.11	4.74	0.036*		28.53	4.70	0.013*		0.43	1.29	0.058
	FH-OP (°)	8.93	1.99		9.35	3.10	0.666		10.06	2.84	0.215		0.71	2.49	0.056
	LFH (mm)	68.92	2.45		70.09	4.85	0.100		70.36	4.60	0.140		0.28	1.44	0.203
Dental	U1-FH (°)	115.87	4.68		121.81	5.23	0.001**		109.81	5.18	0.001**		−12.00	4.28	0.000***
	L1-MP (°)	91.63	2.35		99.04	4.50	0.000***		89.78	3.85	0.132		−9.27	2.27	0.000***
	IIA (°)	126.46	4.57		111.50	6.97	0.000***		134.87	6.38	0.000***		23.37	9.25	0.000***
	OJ (mm)	2.95	0.67		3.94	1.52	0.002**		3.00	0.77	0.591		−0.94	1.49	0.000***
	OB (mm)	2.13	0.92		1.93	1.53	0.757		2.79	0.94	0.078		0.86	1.51	0.000***
Facial profile	SnV-UL (mm)	5.13	1.62		7.14	2.06	0.004**		4.65	1.78	0.414		−2.49	1.49	0.000***
	SnV-LL (mm)	2.85	2.35		5.55	2.77	0.004**		1.26	2.90	0.096		−4.30	2.40	0.000***
	SnV-Pog’ (mm)	−3.72	3.19		−5.42	4.03	0.198		−5.20	3.24	0.176		0.22	2.05	0.467
Hyoid	H-RGN (mm)	33.35	2.64		36.21	4.87	0.065		35.38	4.40	0.200		−0.83	4.24	0.186
	MPH (mm)	8.88	3.13		7.66	5.16	0.457		8.52	4.73	0.813		0.86	2.46	0.051
	H-PTV (mm)	−0.35	2.81		−1.60	4.93	0.422		−1.31	4.65	0.514		0.29	3.48	0.571
	H-C3Me (mm)	1.34	4.69		−1.15	4.89	0.132		−0.72	4.54	0.222		0.43	2.31	0.212
Tongue	TGH (mm)	1.91	0.72		2.95	2.41	0.335		3.98	2.96	0.018*		1.03	2.32	0.005**
	TT-PTV (mm)	50.25	2.53		53.61	4.18	0.000***		48.60	3.77	0.775		−5.01	3.50	0.000***
Soft palate	SPL (mm)	31.46	2.03		34.49	3.21	0.094		35.74	3.66	0.100		1.26	2.81	0.104
	SPI (°)	127.10	2.83		128.84	5.26	0.295		131.09	5.68	0.028*		2.24	2.34	0.000***
Soft palate	SPL (mm)	31.46	2.03		34.49	3.21	0.094		35.74	3.66	0.100		1.26	2.81	0.104
	SPI (°)	127.10	2.83		128.84	5.26	0.295		131.09	5.68	0.028*		2.24	2.34	0.000***
Pharyngeal airway	PNS-Ad1 (mm)	25.10	2.54		27.96	3.90	0.025*		27.60	3.77	0.082		−0.36	2.60	0.344
	PNS-Ad2 (mm)	26.89	1.50		26.82	4.57	0.928		26.71	4.52	0.824		−0.11	1.51	0.632
	SPAS (mm)	13.30	1.57		14.48	2.79	0.183		12.52	2.92	0.399		−1.96	1.09	0.000***
	MAS (mm)	12.64	1.05		13.05	2.80	0.429		11.11	2.85	0.005**		−1.95	1.19	0.000***
	IAS (mm)	11.16	1.30		12.53	3.33	0.053		10.79	3.12	0.539		−1.74	1.47	0.000***
	VAL (mm)	63.65	8.50		63.84	5.06	0.946		64.93	4.95	0.642		1.09	2.68	0.058

BSP, bimaxillary skeletal protrusion, SNA, sella-nasion-A point; SNB, sella-nasion-B point; ANB, A point-nasion-B point; Nper, nasion perpendicular plane; FH, Frankfort horizontal plane, MP, mandibular plane; OP, occlusal plane; LFH, lower facial height; U1, upper incisor; L1, lower incisor; IIA, interincisal angle; OJ, overjet; OB, overbite; SnV, subnasale vertical line; H-RGN, distance between H (the most superior and anterior point of hyoid bone) and RGN (the most protrusive point of retrognathion); MPH, the shortest distance from the mandibular plane (Me-Go) to H; H-PTV, distance between H and the perpendicular line from Pt to FH plane; H-C3Me, distance between H and C3 (anteroinferior limit of the third vertebra)-Me line; TGH, the shortest distance from the tongue dorsum to hard palate; TT-PTV, distance between the tongue tip and the perpendicular line from Pt to FH plane; SPL, soft palatal length, PNS-P (tip of soft palate); SPI, soft palate inclination, angle between the ANS-PNS and PNS-P lines; PNS-Ad1, lower nasopharyngeal airway space, width of the airway along the PNS-Ba line; PNS-Ad2, upper nasopharyngeal airway space, width of the airway along the PNS-R (midpoint of S-Ba line) line; SPAS, superior posterior airway space corresponding to the velopharynx, width of the airway along the parallel line to the palatal plane through Sb (belly of the soft palate); MAS, middle airway space corresponding to the upper glossopharynx, width of the airway along the parallel line to the palatal plane through P; IAS, inferior airway space, corresponding to the lower glossopharynx, width of the airway along the parallel line to the palatal plane through the intersection of the mandible and airway; VAL, vertical airway length, distance from PNS to Eb (base of epiglottis).

*p < 0.05, **p < 0.01, ***p < 0.001.

†Independent t-test was performed for intergroup comparisons with the control group at T0 and T1, respectively.

‡Paired t-test was performed to assess treatment changes (ΔT) in BSP patients.

See Figure 1 for the definition of other landmarks.

Table 5

Posttreatment changes in upper airway variables according to the change in mandibular plane steepness, tested by multivariable linear regression analysis adjusted for sagittal skeletal classification

Dependent variables	Independent variable
	ΔFH-MP (°)
	β-estimate	CI		p-value
ΔH-RGN (mm)	−0.70	−1.80	0.39	0.202
ΔMPH (mm)	0.09	−0.56	0.73	0.791
ΔH-PTV (mm)	0.27	−0.63	1.17	0.553
ΔH-C3Me (mm)	−0.35	−0.95	0.25	0.246
ΔTGH (mm)	0.00	−0.61	0.62	0.990
ΔTT-PTV (mm)	0.18	−0.74	1.11	0.689
ΔSPL (mm)	−0.15	−0.88	0.58	0.683
ΔSPI (°)	0.27	−0.32	0.85	0.363
ΔPNS-Ad1 (mm)	−0.28	−0.96	0.40	0.411
ΔPNS-Ad2 (mm)	0.19	−0.21	0.58	0.344
ΔSPAS (mm)	−0.10	−0.36	0.16	0.443
ΔMAS (mm)	−0.42	−0.69	−0.14	0.004**
ΔIAS (mm)	−0.52	−0.85	−0.19	0.003**
ΔVAL (mm)	0.44	−0.26	1.13	0.211

Multivariable linear regression analysis was performed using the mandibular plane angle (FH-MP) change as an independent variable adjusted for the ANB angle. Δ indicates the amount of treatment change.

CI, confidence interval.

**p < 0.01.

See footnote of Table 1 for the definition of measurements.

Previous studies investigated the pharyngeal airway changes after extraction treatment mostly in BMP patients, with no special consideration of double jaw protrusion. Although it has been difficult to draw a qualified conclusion due to the limited number of studies and heterogeneity across studies, a systematic review[21] summarized that premolar extraction followed by a large retraction of the anterior teeth might lead to upper airway constriction in BMP patients. However, previous findings have never been supported by comparisons with normal pharyngeal dimensions obtained from control samples with normal skeletal and dental occlusion. It is clinically important to understand what the posttreatment upper airway dimension implies in relation to the dentoskeletal changes considering the normal range of airway patency, rather than limiting the interpretation just on the absolute amount of airway change. The present study could show that the retruded tongue position and reduced velopharyngeal spaces after treatment in BSP patients fell within normal ranges as in the control samples, despite the maximum retraction of the upper and lower incisors and anterior alveolus.[15,17,22]

Previous studies have attempted to identify dentoskeletal predictive factors for the pharyngeal airway changes after extraction treatment, mostly in the horizontal direction. Wang et al.[10] suggested a correlation between the retraction distance of the lower incisors and reductions in velopharyngeal and glossopharyngeal spaces in BMP patients in a cephalometric study, but their measurements focused on dental parameters with no definite description of the maxillary and mandibular jaw positions. Chen et al.[9] reported a correlation between the retraction distance of the upper incisors and hypopharyngeal constriction in a cone-beam computed tomographic (CBCT) study, but their samples were irrelevant to the BSP population. Recently, Zheng et al.[23] found correlations among pressure, airflow, and pharyngeal cross-sectional reductions in BMP patients by performing a computational fluid dynamics study, but they did not address dentoskeletal risk factors. Interestingly, the present study suggested that the vertical skeletal change with clockwise mandibular rotation might be a risk factor for glossopharyngeal narrowing in BSP patients, as supported by multivariable linear regression analysis using the FH-MP as an independent variable adjusted for the sagittal skeletal relationship. The clinical significance of this finding is that it would be more important to avoid the increasement of the vertical dimension than to compromise the amount of incisor retraction in the treatment of BSP patients with respiratory concerns, in order to achieve both esthetic and functional respiratory treatment goals.

Regarding the mechanism underlying the pharyngeal narrowing after maximum retraction, the backward displacement of the tongue base, soft palate, and hyoid bone have been suspected.[15,21] Backward displacement of the tongue along with a decrease in the oral cavity pushes the soft palate backward via the communicating action of palatoglossus muscles, leading to velopharyngeal and glossopharyngeal airway reductions in the horizontal direction.[19,22] The impact of hyoid positional changes on the upper airway has been a topic of debate among the studies since it depended on the presence or absence of compensatory positional adaptation to prevent encroachment of the tongue into the pharynx.[7,24] In the present study, posttreatment retrusion of the tongue tip (TT-PTV) showed a positive correlation with the retrusion amount of A-point, and posttreatment retroclination of the soft palate (SPI) was correlated with the retrusion amount of B-point and pogonion point. On the other hand, posttreatment changes in the pharyngeal airway spaces and hyoid bone showed no significant correlation with the amount of dentoalveolar retraction. Notably, the positional changes in parapharyngeal soft tissues showed too little reliability or reproducibility in lateral cephalography to have clinical significance for the posttreatment pharyngeal changes. Our understanding of the physiologic mechanism of functional adaptation and restoration of the tissues needs to be strengthened.

It is still unclear whether the pharyngeal narrowing with parapharyngeal tissue displacement after extraction treatment might bring about stable functional adaptation in the long term. Nevertheless, in consideration of various phenotypic causes of OSA, including the anatomical and nonanatomical factors,[25] it should be noted that even BSP patients might have nonanatomical OSA risks irrelevant to the craniofacial condition. In this context, upper airway dimensions and respiratory or sleep functions need to be maintained under control by orthodontists, especially when maximum anterior retraction treatment is planned for the hyperdivergent BSP patients with pharyngeal constriction and respiratory symptoms.

The present study included the following fundamental limitations: (1) morphological image analysis of pharyngeal airway taken in the awake state and upright position; (2) difficulty in standardization of cephalometry to represent consistent parapharyngeal soft tissue positions; (3) no inclusion of three-dimensional analysis of the skeletal pattern and pharyngeal volume; (4) the absence of long-term treatment changes matched with follow-up functional investigations. Although the pharyngeal changes in the transverse dimension could not be reflected in the present analysis, our study had the significance of attempting pharyngeal airway prediction in relation to the sagittal and vertical dentoskeletal treatment-induced changes by using the routinely employed screening tool of lateral cephalometry. It is still clinically possible to see the dentoskeletal pattern, pharyngeal morphology and dimension, and parapharyngeal soft tissues at a glance in planning orthodontic treatment. Further studies are ongoing to evaluate the relationship between altered upper airway dimension after orthodontic treatment and the resultant effects on respiratory or sleep function, based on a three-dimensional CBCT analysis related to polysomnographic studies.

The BSP patients showed significant reductions in the velopharyngeal and glossopharyngeal airway spaces after extraction treatment; however the reduced posttreatment dimensions corresponded to the normal dimensions in the control samples.

The glossopharyngeal space might be susceptible after treatment in the hyperdivergent BSP group, especially when the mandibular plane angle is increased during treatment.

Control of the vertical dimensions during maximum anterior retraction would be important to achieve both esthetic and respiratory functional improvements in BSP patients.