Trajectories of Prenatal Alcohol Consumption in An Australian Population-Based Cohort of Pregnant Women - A Data Driven Approach

Accurate information on dose, frequency and timing of maternal alcohol consumption is critically important when investigating fetal risks from prenatal alcohol exposure. Identication of distinct alcohol use behaviours can also assist in developing directed public health interventions to prevent adverse child outcomes, including Fetal Alcohol Spectrum Disorder. We aimed to determine group-based trajectories of time-specic, unit-level, alcohol consumption using data from 1458 pregnant women in the Asking Questions about Alcohol in Pregnancy (AQUA) longitudinal study in Melbourne, Australia. Six alcohol consumption trajectories were identied incorporating four timepoints across gestation. Labels were assigned based on consumption in trimester one and whether alcohol use was continued throughout pregnancy: abstained (33.8%); low discontinued (trimester one) (14.4%); moderate discontinued (11.7%); low sustained (13.0%); moderate sustained (23.5%); and high sustained (3.6%). Median weekly consumption in trimester one ranged from 3 grams (low discontinued) to 184 grams of absolute alcohol (high sustained). Alcohol use after pregnancy recognition decreased dramatically for all sustained drinking trajectories, indicating some awareness of fetal risk. Improved understanding of factors that contribute to alcohol consumption in pregnancy in specic sub-populations is critical when developing prevention strategies. Maternal characteristics associated with different trajectories can inform the strategic development of such initiatives.


Introduction
Our understanding of the dose-response relationship between prenatal alcohol exposure (PAE) and adverse child outcomes, especially at low levels of exposure, is limited despite several decades of research. Over the past 10 years, several systematic reviews have attempted to consolidate the available evidence, but meta-analytic inference has been hampered by methodological limitations, and heterogeneity in measurement and reporting of alcohol consumption during pregnancy [1][2][3].
Although objective laboratory measures of PAE exist, including testing newborn meconium for fatty acid ethyl esters and maternal hair for ethyl glucuronide, these have high positivity thresholds, low sensitivity and speci city, and do not capture the timing or duration of PAE [4]. Consequently, researchers continue to rely on maternal self-report, either retrospective, or more reliably, collected during the prenatal period. Highly variable methods are then used to quantify PAE, mostly based on risk groupings with pre-de ned cut off scores [5]. An important consideration is that any potential non-linear dose-response relationships between prenatal alcohol consumption and fetal effects will be obscured by arbitrary ordinal or nominal measures of exposure. If the relationship is non-linear, there may be a threshold of effect. Understanding the dose-relationship is most important when delivering harm reduction messages, especially to women for whom abstinence is di cult to achieve. Alternatively, if the risk to the fetus is known to increase rapidly with only small increases in the level of PAE, it would be important to emphasise the need for abstinence.
Further, increasing use of meta-analysis to consolidate research evidence has highlighted the lack of detailed, comparable measures of PAE which can be aggregated across studies [5,6]. Ideally, these measures should provide detailed information on metric units (e.g. grams of absolute alcohol) that characterise consumption patterns at speci c pregnancy timepoints.
The 'Asking Questions about Alcohol in Pregnancy' longitudinal cohort study (AQUA) was speci cally designed to capture common prenatal drinking patterns in order to investigate potential effects on long-term developmental outcomes of children in the general population [7]. The PAE classi cations originally used in the study were categorised to re ect real-life maternal drinking, taking into account the dose, pattern and timing of alcohol consumption during pregnancy and based on the 'composite' method described by O'Leary et al. in 2010 [8]. Brie y, classi cation was based on a composite of continuous measures, i.e. the total number of grams of absolute alcohol consumed per week, and the maximum amount of absolute alcohol consumed per occasion of drinking. Units of alcohol (grams) were calculated from maternal self-report of the type and amount of alcohol consumed, using a detailed pictorial drinks guide at four timepoints during pregnancy. This enabled different patterns of PAE to be discerned, particularly low-level consumption, episodic binge drinking, special occasion drinking, and cessation of alcohol consumption upon pregnancy recognition, usually between ve and seven weeks of gestation [9].
Developments in classifying temporal continuous data allow for maternal alcohol consumption patterns to be identi ed more objectively. Methods such as group-based trajectory modelling (GBTM) are increasingly employed in clinical research to describe the course of an outcome or behaviour over time [10]. In alcohol research, GBTM can be used as a tool to measure consumption trajectories arising directly from the source data, without the need for pre-determined classi cation. This allows analysis of temporal drinking patterns and provides more nuanced results than aggregate or presence/absence drinking information [5,11,12]. This is especially important when investigating the potential relationship between low-level alcohol use and adverse child outcomes [6]. Moreover, trajectory modelling using unit-level and temporal alcohol consumption data has the potential to provide more detailed information on actual drinking patterns than predetermined alcohol consumption categories (i.e. low/moderate/high levels).
The main aim of this paper was to determine the different alcohol consumption trajectories in a general pregnant population using a data driven approach. Another aim was to describe salient maternal characteristics predictive of these patterns, which may both assist in targeting prevention approaches in different populations and identify potential confounding factors to be considered when investigating the casual relationship between PAE and child outcomes.

Methods
Data from 1458 (92.9%) women who took part in the AQUA study were used in this analysis. [9] We excluded 112 women (7.1%) who were lifetime abstainers because our target population was women who would normally drink some alcohol but who may abstain during pregnancy.
Prenatal alcohol consumption data Detailed information on the quantity and frequency of alcohol consumption was collected via questionnaires delivered at:1) recruitment (< 18 weeks' gestation); 2) 25 weeks' gestation; and 3) 35 weeks' gestation. Together these provided data on alcohol consumption in the three months pre-pregnancy, post conception but pre-pregnancy recognition, and for each trimester of pregnancy. The mean (SD) gestational age at pregnancy recognition was 4.9 (1.5) weeks [9].
Women were provided with a pictorial drinks guide showing common types and volumes of alcoholic drinks including red and white wine, champagne, beer, cider, spirits, alcoholic sodas, pre-mixed spirits, port, sherry, and cocktails. This drinks guide was developed with input from a focus group study [13]. Women were asked to use the drinks guide to identify their 'usual' pattern of drinking, with provision for up to ve types of drinks. For each beverage identi ed, they were asked how often they usually drank this type of alcohol (less than once per month, 1-2 days per month, 1-2 days per week, 3-4 days per week, 5 or more days per week) and how many drinks they usually consumed on each occasion (less than one drink, 1-2 drinks, 3-4 drinks, 5-6 drinks, 7 or more drinks). Women were then asked if there were any 'special occasions' (or di cult times) when they consumed more alcohol than usual, the frequency of these occasions, the drink types, and the number of drinks per occasion. If a woman reported consuming seven or more drinks on any occasion, she was asked to provide the maximum number. Estimates from 'special occasions' were combined with information from 'usual' alcohol consumption to calculate a maximum weekly intake [9].
The amount of alcohol consumed per week was derived from the number and types of drink reported by women, which were converted to standard drinks to calculate the amount of absolute alcohol in grams (gAA) consumed. One standard drink in Australia is equal to 10 gAA. A binge episode was de ned as consumption of at least ve standard drinks (50gAA) per drinking occasion.

Maternal characteristics
Other alcohol behaviour Pre-pregnancy binge drinking for the three months prior to pregnancy was dichotomised as "yes, at least one binge episode" or "no binge drinking". Women were also asked about their drinking history, including how old they were when they rst started drinking regularly and the age when they rst became intoxicated from drinking alcohol (de ned as slurred speech, unsteady on their feet, or blurred vision). Responses were dichotomised into whether or not women were at least 18 years-old at the time to re ect the legal drinking age in Australia. To gauge possible individual variation in alcohol metabolism, women were asked if, prior to their pregnancy, they felt the effects of alcohol very quickly, quickly, normally, slowly, or very slowly.

Statistical analysis
Group-based Trajectory modelling (GBTM), a specialised form of nite mixture modelling that does not require complete data across all time points [10] was used to investigate prenatal alcohol consumption trajectories. For a hypothesised number of underlying latent groups, it uses maximum-likelihood estimation to: identify distinctive clusters of individuals who follow similar trajectories for an outcome; outline the shape of each trajectory and size of each group; and pro le the characteristics of individuals within trajectory groups. GBTM allows analysis of factors in uencing group membership through the inclusion of time-invariant predictors. Analyses were conducted in Stata/ICv15.1 (StataCorp LLC) using the traj plugin.
Model selection involved 2 stages: (1) identi cation of the optimal number of trajectory groups and (2) determination of preferred polynomial orders specifying the shape of the identi ed trajectories. Best tting models were determined for two to six groups (models of seven groups and above failed to converge) and then compared using the Bayesian Information Criterion (an increase of BIC [ΔBIC] > 2), model parsimony, entropy closer to 1, and t with prior theory [14,15]. Participants were then classi ed into trajectory groups according to their maximum posterior probability of group membership. Model t was further assessed by calculating group mean posterior probability and odds of correct classi cation (OCC) [10].
Total alcohol intake per week in grams was modelled over four time points, roughly classi ed in terms of gestational week: prior to pregnancy recognition = 5 weeks, trimester one = 13 weeks, trimester two = 25 weeks, and trimester three = 38 weeks. To accommodate the commonly skewed distribution of grams of alcohol outcomes [12,16], a square-root transformation was applied to adjust for non-normality.
Binge episodes were also recorded at each time point. Almost 1 in 5 women had a binge episode during pregnancy, but 99.4% of these occurred at the 'prior to pregnancy recognition' time point [9]. Therefore, presence of one or more binge episodes was included as a dichotomised, time-static predictor of group membership in GBTM.
Following identi cation of a best model t, we investigated the association of several non-pregnancy related maternal alcohol use characteristics with group membership using chi-squared tests, making planned comparisons between different alcohol consumption groups and abstainers.
Multivariate logistic regression was used to examine associations between maternal characteristics and group membership as compared to abstinent women (control). Unadjusted and adjusted odds ratios (controlling for all characteristics signi cantly related to any of the drinking patterns) were calculated. For predictor variables with more than two categories (maternal age, educational attainment, household income and pre-pregnancy body mass index), p-values from likelihood ratio tests were used to evaluate the predictors. Alpha was set to 0.05 for all analyses.

Results
Group-based trajectory modelling (GBTM) GBTM was conducted to determine the best-tting models for two to six groups. In the best-tting models for each number of groups, we found no higherorder polynomial effects of time. All groups followed a linear or intercept-only trajectory. Each best-tting model per group number is detailed in Table 1. Comparatively, the AIC and BIC showed that the ve-and six-group models were a better t than models with four groups or less. Optimal BIC was found in the ve-group model. This is illustrated by the difference in BIC: ΔBIC was 114.72 for the ve-group model whereas the six-group model showed poorer t with a ΔBIC of -3.03. Entropy was best in two-and four-group models, but acceptable (> 0.8) in models with ve or less groups. With all t information considered, the ve-group model was determined the best t (Supp Fig. 1).
However, the GBTM found no mathematical difference between pregnancy abstainers and women who had consumed some alcohol at the 'prior to pregnancy recognition' time point. To accommodate this theoretical implication, whole of pregnancy abstainers were 'forced' into their own group, resulting in a nal six group solution. Trajectories of the ve-group model with abstainer/control separation (group six) are illustrated in Fig. 1. Group means with 95% con dence intervals (CI -shaded area) are presented on a logarithmic scale to illustrate group differences at the lower end of alcohol consumption. The same six alcohol trajectories are also presented on a normal y-scale in Supp Fig. 2.
Examination of the trajectories, which are summarised in Table 2, resulted in six groups which have been named as follows:  Pre-pregnancy recognition binge episodes were con ned to the three moderate and high drinking trajectories with the occasional post pregnancy-recognition binge episode in the moderate and high sustained trajectories only. ( Association between group membership and non-pregnancy related alcohol use behaviour Pregnant women with a moderate to high alcohol consumption trajectory were less likely to report that they felt the effects of alcohol quickly than controls or women with a low consumption trajectory. Women with moderate to high consumption were also more likely to have experienced their rst alcohol intoxication before the Australian legal drinking age of 18 years and to have had at least one binge drinking episode in the three months before pregnancy. (Table 3).

Association between group membership and demographic and pregnancy-related characteristics
Multivariate analysis revealed no discernible difference between controls and the low discontinued trajectory in any of the characteristics investigated (Table  4). Compared to controls, women in all other alcohol consumption trajectory groups were two to seven times more likely to be Caucasian (e.g. low sustained (AOR 2.32 (95%CI 1.40-3.85)) and moderate sustained (AOR 7.13 (95%CI 3.99-12.73)). Cigarette smoking in pregnancy was associated with all moderate to high drinking trajectories, e.g. moderate sustained (AOR 4.05 (95%CI 2.60-6.31)) and high sustained (AOR 4.28 (95%CI 1.92-9.54)). Women in the moderate discontinued trajectory were more likely to have an unplanned pregnancy (AOR 2.99 (95%CI 1.92-4.67)) and be pregnant with their rstborn (AOR 2.09 (95%CI 1.38-3.16)). Women with low sustained group membership were less likely to be primipara (AOR 0.61 (95%CI 0.42-0.91)). A sustained alcohol consumption pattern was more likely in women in their early to mid-thirties and a high sustained level was more likely in women aged 35 years or more. Increasing household income was associated with moderate to high sustained group membership.

Discussion
Group-based trajectory modelling of continuous data on grams of absolute alcohol consumed per week during pregnancy identi ed ve distinct trajectories of prenatal alcohol consumption. In this population-based cohort of pregnant women, the group that consumed one to two standard drinks, once or twice per month until they became aware that they were pregnant (14.4%), was mathematically indistinguishable from the group that abstained from alcohol (33.8%). A second group of women who discontinued their alcohol consumption at some point during the rst trimester averaged around three standard drinks per week until then. Women in this moderate-discontinued group were more likely to have an unplanned pregnancy, be primiparous and smoke cigarettes. Cigarette smoking was also associated with moderate and high sustained alcohol use, as was higher maternal age and household income. Pre-pregnancy and early pregnancy binge episodes were common in the moderate and high-level groups, regardless of whether alcohol use was discontinued or sustained. Other characteristics of all moderate and high-level groups were a history of underage intoxication and higher self-reported tolerance for the effects of alcohol.
Importantly, all sustained drinking trajectories showed a dramatic decrease in levels of alcohol use after pregnancy recognition, potentially indicating a degree of awareness of the potential harms to the unborn child.
Improved understanding of the factors that contribute to alcohol consumption in pregnancy in speci c sub-populations is critical when developing prevention initiatives. The maternal characteristics associated with different trajectories may inform the strategic development of such initiatives and interventions. Here, the most important trajectory we identi ed for targeted prevention interventions is the moderate-sustained group, which comprised almost a quarter of all women in our study. Women following this consumption pattern are clearly not responding to existing public health messaging advising abstinence.
Qualitative research exploring the reasons for alcohol use in pregnancy has shown that while most women are aware that abstinence is recommended, there is a general perception that the risk of harm from occasional alcohol use is low. This usually results from con icting advice by maternity clinicians, the women's own observations of the behaviour of family and friends, plus the lack of convincing research evidence on harm from low-level consumption patterns [17]. Consequently, some women make individual decisions about the perceived quantity of alcohol that is without risk of harm, even if they received best practice health messages advising abstinence. The women in our study who followed a moderate-sustained alcohol use pattern had often been drinking alcohol regularly from an early age, and although they reduced their intake following pregnancy recognition, alcohol consumption may be well-established and normal aspect of their social environment. Further, a perception of not being easily affected by alcohol may contribute to feeling that some level of alcohol consumption is unlikely to affect the unborn baby. This sizeable group of pregnant women will require sophisticated health messages that acknowledge the uncertainties around risk of harm from low-level or occasional alcohol use, but also emphasise the importance of maximising health outcomes for their baby through abstinence. Brief psychosocial interventions have established bene ts in women with heavier alcohol consumption and although the evidence of effectiveness is not as strong for pregnant women identi ed as consuming low-levels of alcohol, behaviour change techniques such as tailored information about consequences and fostering positive social support, or goal setting appear to increase abstinence rates [18]. It may be that abstinence rates among pregnant women will only improve if maternity service systems consider the different social and cultural contexts which in uence women's drinking choices.
Consideration could be given to encouraging positive involvement from partners, family and friends, and attention to the provision of clear and consistent messages about the bene ts of abstaining from alcohol use in pregnancy.
We previously classi ed drinking patterns during pregnancy according to pre-determined cut-off levels based on the 2001 Australian National Health & Medical Research Council Alcohol Guidelines: Health Risks and Bene ts [19]. These guidelines, which were revoked in 2009, stated that pregnant women should consider not drinking alcohol, but if they chose to drink, they should have less than seven standard drinks (< 70g absolute alcohol) over the course of a week, and no more than two standard drinks (≤ 20g absolute alcohol) per day. To date, we have classi ed pregnant women in the AQUA study who followed this drinking pattern as "low level" drinkers.
However, GBTM showed that most women in the lowest drinking trajectory (low discontinued) consumed a minimal amount of alcohol, less than one standard drink (< 10g absolute alcohol) per week, compared with our original classi cation which included women who consumed up to almost seven standard drinks (< 70g absolute alcohol) per week. This distinction may prove invaluable when investigating the potential harmful effects on the unborn child of various prenatal alcohol exposure patterns.
This study is not the rst to use GBTM to describe maternal alcohol consumption patterns. In an earlier Australian study Tran et al. identi ed three trajectories from six pre-speci ed frequency and quantity questions asked at four different time points, ranging from the period before pregnancy to six months postpartum in a longitudinal pre-birth cohort of 6,597 Australian women that commenced in 1981 [20]. The three trajectories comprised women who abstained or drank minimally (53%), those who uctuated at an average of about 0.37 glasses per day (39%) and those who drank at a higher level of about 2.5 glasses per day before pregnancy but dropped their alcohol intake to about 0.6 glasses per day during pregnancy. A major difference between these data and the present study is that even with a smaller sample size, we identi ed an additional two trajectories of women who discontinued alcohol use at, or soon after, pregnancy recognition. This difference most likely re ects changes in community awareness and advice on alcohol abstinence given by maternity providers in the 30 years between studies.
A more recent analysis by Dukes et al. of about 11,692 women taking part in the Safe Passage Study in 2007 identi ed ve trajectories more akin to those we found in the AQUA study [21]. Although the levels of consumption and the timing of cessation differed in their population, the ve groups included one abstinent/minimal use group, two that discontinued and two that continued some alcohol consumption throughout pregnancy, like the current study. In 2019, Bandoli et al. published an analysis of ve GBTM-derived alcohol consumption trajectories from 471 pregnant women and their potential association with infant growth and early development [11]. The authors reported an association between the highest consumption trajectory and de cits in infant birth weight and length and psychomotor development at six to 12 months of age. However, the evidence generated from the study is limited given its small sample size and inclusion of only 24 participants in the highest consumption trajectory.
Most importantly, both Dukes and Bandoli et al. reported speci c metrics that characterised each trajectory. These were presented as a daily average; either as the number of standard drinks de ned as 14gAA (Dukes [21]) or in ounces of absolute alcohol (Bandoli [11]). Use of similar reporting methods across studies will improve our ability to integrate our results with those from other studies going forward. This is critical to accumulate robust research evidence to better predict which prenatal alcohol exposure patterns are most strongly associated with particular adverse child outcomes.

Strengths and limitations
A strength of this study is the detail of the alcohol measures and the focus on the most frequent prenatal patterns of alcohol consumption (low, moderate and discontinued) rather than heavy and sustained alcohol use. This focus has played a key role in the study's high participation and low attrition rates over the course of the women's pregnancy, but most importantly, in providing alcohol consumption data of the highest quality possible [7,13,22]. Although we measured these data prospectively, and thus optimised our ability to measure frequency, dose and timing of exposure, the use of self-reported questionnaires runs the risk of reporting bias. However, our focus group research showed that if questions on alcohol in pregnancy are appropriately contextualised and include an option to report unusual drinking episodes, this encourages more accurate reporting [13], a nding which appears con rmed by the high number of binge episodes reported in response to the special occasion question [9].
In our nal GBTM model we found one small group (high sustained, n = 52), which contained fewer than the suggested minimum in any one group, being < 5% of the sample [10]. We acknowledge this limitation but believe that this nding is a true re ection of the small number of consistently high drinkers in this population-based cohort of pregnant women and that it is imperative to describe such an important clinical group.
Another strength of this study lies in the ability of GBTM to directly classify the continuous source data without the need for arbitrary cut-offs. However, GBTM is an application of nite mixture modelling, which assumes that the study population is composed of distinct groups de ned by their trajectory membership.
This theoretical assumption may be compromised when using non-research data and there may be women who could be assigned to more than one trajectory. We considered entropy as a measure of classi cation accuracy in our nal model selection and found this to indicate a high degree of precision in the assignment of individuals to their most likely group.

Conclusion
GBTM-derived trajectories of prenatal alcohol consumption can re ect real-life maternal drinking patterns because they preserve the timing, quantity, and frequency of consumption derived directly from unit-level source data. Understanding these distinct consumption trajectories and their associated maternal characteristics can assist in identifying antenatal populations for targeted alcohol cessation approaches. The trajectories also provide a discerning classi cation method for investigating causal relationships between prenatal alcohol exposure and child outcomes. Further, an inherent ability to mathematically de ne the underlying unit-level consumption patterns of each trajectory may reduce heterogeneity in exposure classi cation across studies, thereby improving the ability to aggregate data in future meta-analyses.

Declarations
Data availability The AQUA datasets analysed during the current study are not publicly available due to institutional ethics requirements restricting access to study investigators but are available from the corresponding author on reasonable request.

Figure 1
Trajectory groups for prenatal alcohol consumption based on the ve-group model with abstainer/control separation

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download. Mugglietalsupplementarymaterial.docx