Abstract
BACKGROUND: Characterizing maternal glucose sampling over the course of the entire pregnancy is an important step toward improvement in prediction of adverse birth outcome, such as preterm birth, for women with type 1 diabetes mellitus (T1DM). OBJECTIVES: To characterize the relationship between the gestational glycemic profile and risk of preterm birth using a joint modeling approach. METHODS: A joint model was developed to simultaneously characterize the relationship between a longitudinal outcome (daily blood glucose sampling) and an event process (preterm birth). A linear mixed effects model using natural cubic splines was fitted to predict the longitudinal submodel. Covariates included mother's age at last menstrual period, age at diabetes onset, body mass index, hypertension, retinopathy, and nephropathy. Various association structures (value, value plus slope, and area under the curve) were examined before selecting the final joint model. We compared the joint modeling approach to the time-dependent Cox model (TDCM). RESULTS: A total of 16,480 glucose readings over gestation (range: 50-260 days) with 32 women (28%) having preterm birth was included in the study. Mother's age at last menstrual period and age at diabetes onset were statistically significant (beta = 1.29, 95% CI 1.10, 1.72; beta = 0.84, 95% CI 0.62, 0.98) for the longitudinal submodel, reflecting that older women tended to have higher mean blood glucose and those with later diabetes onset tended to have a lower mean blood glucose level. The presence of nephropathy was statistically significant in the event submodel (beta = 2.29, 95% CI 1.05, 4.48). Cumulative association parameterization provided the best joint model fit. The joint model provided better fit compared to the time-dependent Cox model (DIC (JM) = 19,895; DIC (TDCM) = 19,932). CONCLUSION: The joint model approach was able to simultaneously characterize the glycemic profile and assess the risk of preterm birth and provided additional insights and a better model fit compared to the time-dependent Cox model.