A latent class analysis to identify socio-economic and health risk profiles among mothers of young children predicting longitudinal risk of food insecurity

Background The purpose of the current study was to use a social determinants of health (SDOH) framework and latent class analysis (LCA) to identify risk classes among mothers with young children. The risk classes were then used to predict food insecurity severity and stability/change of food insecurity over time. Method The secondary data from the Fragile Families and Child Wellbeing Study (n = 2,368; oversampled for non-marital births) was used in this study. Household food insecurity was assessed using the 18-items USDA Food Security Survey. A seventeen-item inventory of educational, economic stability, incarceration (i.e. social context), neighborhood safety (i.e. neighborhood and built environment), health and health care, and substance use behaviors at baseline/Year-1 were included to identify SDOH risk indicators in the LCA. Covariate-adjusted multinomial logistic regression models were used to examine the relation between risk classes at Year-1 and the severity of food insecurity at Year-3 and stability/change of food insecurity between Year-3 and Year -5. Results LCA identified five risk classes: High utility and medical hardship (Class 1), high housing and employment hardship, high substance use, and incarceration (Class 2), high housing and medical hardship, poor health, and health care (Class 3), high employment hardship and low-income (Class 4) and low-risk (Class 5). The Class 1, Class 2 and Class 3 had greater odds of low food security and very low food security at Year-3 compared to Class 4. In addition, compared to Class 4, Class 1, Class 2 and Class 3 had greater odds unstable food insecurity and persistent food insecurity over time. Conclusions LCA could be used to identify distinctive family system risk profiles predictive of food insecurity. The generated risk profiles could be used by health care providers as an additional tool to identify families in need for resources to ensure household food security.


Introduction
Households with limited access to adequate healthy and nutritious food due to lack of money and other resources are considered to be food insecure [1]. According to the US Department of Agriculture [1] in 2019, 13.6% of US households with children experienced food insecurity, and the numbers are expected to increase with the current inflation concerns. Household food insecurity is a major public health concern as it is linked to a range of negative health consequences among adults. Food insecurity experiences are associated with the risk of depression, mood or anxiety disorders, overweight/obesity, hypertension, diabetes, metabolic syndrome and poorer overall health status [2][3][4][5][6][7]. Literature also indicates that there is a bidirectional relationship between food insecurity and depression among mothers residing in urban and rural settings [8,9]. While children are usually buffered from directly experiencing food insecurity [10], toddlers score lower on cognitive and health assessments when they reside with adults who experience food insecurity [11]. Children in food-insecure households are more likely to report poor general health, experience chronic health conditions, including asthma, eczema or other skin allergies, depressive symptoms and poorer mental health, and acute health conditions such as cold and stomach problems compared with peers in food-secure households [12][13][14][15][16][17][18][19][20][21]. Therefore, it is important to identify households at risk of food insecurity as earlier as possible and provide adequate resources to prevent/control these potential adverse health conditions.
However, these risk factors do not exist in isolation. The social vulnerability perspective suggest that individuals are unable to escape negative outcomes associated with experiencing food insecurity due to the multiple socio-economic disadvantaged characteristics that they are exposed to [44]. Specifically, low-income families are more likely to experience multiple overlapping risk factors at the same time [40,[45][46][47]. Studies using cumulative risk approach [48] (accumulation of risk factors rather than any one particular risk factor that influences a particular outcome) have provided evidence of accumulation of multiple family risk factors to be associated with greater risk of food insecurity. Hernandez [32] reported that exposure to cumulative financial strain was related to experiencing marginal food security over food security among non-poor households. In addition, poor health and risky health behaviors index was associate with food insecurity experiences over marginal food security experiences among domiciled populations [32] and associated with food insecurity experiences over food secure experiences among adult experiencing homelessness [49]. O'Reilly, Hager, Harrington, et al. [50] created a cumulative risk index based on caregivers' endorsement of five demographic, behavioral, and psychosocial risk factors. Each additional cumulative risk index factor was related to a 54% increase in the odds of experiencing food insecurity. One limitation of using cumulative approach is that based on the variables used to create the cumulative risk index, it may not capture the impacts of risk factors across domains and the relative strength of association of factors within a domain.
Latent class analysis (LCA) is an approach that could be used to identify the risk profiles of households experiencing food insecurity, that could be used as a screener to early identify families at risk of food insecurity. LCA is a statistical method that uses observed categorical/continuous variables to identify subpopulations that have similar characteristics among a greater population that has varying characteristics [51]. Therefore, LCA can be beneficial in distinguishing subpopulations that are at greater risk of food insecurity and understand the complex challenges experienced by these families. In clinical settings, the characteristic profiles could be used to identify families that may need additional recourses to ensure household food security/prevent insecurity. Identifying risk profiles can be useful for both treatment and prevention of food insecurity, particularly with the examination of multiple levels of food security (low, marginal, high) and those with fluctuating/persistent/intermittent FI. Further, the identified profiles will be beneficial in designing programs or revising policies to meet the specific needs of the at-risk groups to overcome or prevent future food hardships.
Within the food security literature there is currently only one US-based study that used the LCA approach to identify subgroups (latent classes) at risk of food insecurity [50]. Using five risk factors-employment status, education, having a smoker in the household, stress, and depression-in the LCA, they have identified three latent classes among low-income African American urban households with adolescent daughters: 1) high stress/depression, 2) low education/low stress and depression, and 3) low risk overall. Compared to the low risk class, other two groups had greater risk of food insecurity [50].

Current study
According to current literature there is one study that uses the LCA approach to identify latent risk classes based on 5 risk factors among mothers with children and assess whether the risk classes are associated with the risk of food insecurity [50]. The literature would benefit from research that uses a broader framework that considers individuals' physical, economic, and social environment to create latent risk indicators among a sample of racially and ethnically diverse mothers. A Social Determinants of Health (SDOH) framework [52] considers five domains (1. educational, 2. economic stability, 3. social context, 4. neighborhood and built environment, and 5. health, health care, and substance use) that influences food insecurity. Consequently, there are no studies that then use a longitudinal study design to examine the SDOH latent risk indicators as predictors of the severity of food insecurity or stability/change in food insecurity. To address this research gap, the current study expands upon established characteristics from a prior research study that described the social vulnerability of lowincome families [50] to identify risk classes using LCA among urban mothers with young children using a 17-item inventory that encompass the five SDOH domains. Using a longitudinal study design that covers four years, the risk classes are then used to predict the severity of food insecurity 2 years later and to predict the stability/change of food insecurity over a 2-4-year period.

Participants
The current study utilizes secondary data from the Fragile Families and Child Wellbeing study (FFCWS). The sample and the study design of the FFCWS is reported elsewhere [53]. Briefly, the study follows a cohort of 4898 children born in the US between 1998 and 2000 and their parents/primary caregivers. The study sample included children from 20 large cities (cities with populations of 200,000 or more) and was oversampled for non-marital births (76% children born to unmarried mothers). The core study included maternal and paternal interviews at the child's birth (baseline) and follow-up interviews via telephone when the child was approximately age one (Year-1), three (Year-3), five (Year-5) and nine (Year-9). The core study at age  included an interview with the child and the primary caregiver. In addition, an in-home study was conducted at Year-3 and -5 that included measures of household food insecurity. All participants in the Year-3/5 core study were invited for the Year-3/5 in-home study. The FFCWS was conducted according to the guidelines described in the Declaration of Helsinki and was approved by Princeton University and Columbia University Institutional Review Boards. Written informed consent was obtained from all participants by the FFCWS research team. The current study was approved by the Institutional Review Board at University of Texas Health Science Center.
The current study utilized data from the baseline, Year-1, Year-3 and Year-5 surveys of the FFCWS. The analytic sample selection is illustrated in Fig 1. The current study sample included only the mothers who completed the core survey at Year-1 and the in-home surveys at Year-3 and Year-5. In addition, mothers with missing data for the household food security at either Years -3 and -5 and those with missing covariates were excluded from the current study sample. Attrition analysis was conducted to compare the current study sample to excluded participants using either chi-square tests (categorical variables) or t-tests (continuous variables). The final study sample included 2,348 mothers. Compared to the mothers excluded from the study sample, the mothers included were significantly more likely to be non-Hispanic white (23% vs. 20%, χ 2 = 7.51, p = 0.006) or black (50% vs. 45%, χ 2 = 11.53, p = 0.001), US born (88% vs. 78%, χ 2 = 76.06, p<0.001), have public health insurance (60% vs. 54%, χ 2 = 16.89, p<0.001) and report utility hardships (42% vs. 37%, χ 2 = 8.87, p = 0.003) and depression (17% vs. 14%, χ 2 = 4.44, p = 0.035). Mothers in the study sample reported significantly greater rates of stable food security during Years-3-5 compared to the excluded participants (76% vs. 66%, χ 2 = 5.35, p = 0.021).

Measures
Food insecurity. Household food insecurity at Year-3 and Year-5 was assessed using the 18-items United States Department of Agriculture's (USDA) Food Security Survey. The survey items inquire about household food availability and the expenditure on food during the past year (e.g., "worried food would run out before got money to buy more", "Food bought didn't last" and "relied on few kinds of low-cost food to feed children"). Responses of "yes," "often," "sometimes," "almost every month," and "some months but not every month" were coded as affirmative. The households were categorized into four levels of food security based on the sum of affirmative responses using the USDA cut-off values at each year (Raw score of 0 = High food security; raw score of 1-2 = Marginal food security; raw score of 3-7 = Low food security; and raw score 8-18 = Very low food security) [54].

PLOS ONE
To assess stability and change in the food security during Year-3 to 5, at each year households were further categorized as food secure (raw score of 0 to 2 affirmative responses) or food insecure (raw score of 3 or more affirmative responses) [1]. Based on reported food security status at Year-3 and 5, a categorical variable with three levels was created to indicate stability and change in food security (stable food secure = food secure at both time points; unstable food security = food secure at Year-3 and insecure at Year-5 or vice versa; stable food insecure = food insecure at both time points).
Risk indicator variables used in the latent class analysis. Based on prior research [1,32,50,55,56] a range of educational, economic stability, social context, neighborhood and built environment, and maternal health, health care and substance use variables measured at Baseline/Year-1 were identified as potential SDOH risk indicators to be used in the Latent Class Analysis (LCA). All risk factors listed below will be used as independent indicators in the LCA.
Educational variables. Educational risk indicators included self-reported highest level of education of the maternal grandparents (less than high school, high school/equivalent, greater than high school) and maternal highest level of education (less than high school, high school/ equivalent, greater than high school) at Year-1.
Economic stability variables. This included maternal poverty status [<200% of federal poverty level (FPL); low-income], 200-299% FPL (working poor) and �300% of FPL (not poor/ higher income), employment (regular work for pay vs no regular work for pay) and household material hardships at Year-1. Four-items measured hardships related to housing (e.g. not pay the full amount of the rent or mortgage, evicted from your home or apartment for not paying the rent or mortgage), four-items measured hardships related to utility/bill payments (e.g. missed a payment or were late with the gas or electricity bill, gas or electricity ever turned off because the bill was not paid) and a single-item measured medical hardships (needed to see a doctor or go to the hospital but did not go). The response options were yes/no and one or more affirmative responses to each set of items measuring hardships was coded as experiencing the particular hardship.
Incarceration variable. Maternal and paternal reports of child's father ever being in jail by Year-1 was used as an indicator of social/community level risk.
Neighborhood safety variable. At baseline mothers reported on how safe the streets around their home at night are. The response options were 1 = very safe, 2 = safe, 3 = unsafe, and 4 = very unsafe. A dichotomous perceived neighborhood safety variable was created to indicate neighborhood and built environment level risk.
Maternal health and health care characteristics. Maternal health related risk factors included self-reported poor health, health problems limiting work, major depressive disorder (MDD), generalized anxiety disorder (GAD), health insurance (insurance, public insurance, private insurance) at Year-1. Maternal MDD was measured using the Composite International Diagnostic Interview-Short Form (CIDI-SF) and those who answered affirmatively to three or more questions were considered depressed [57]. Maternal GAD was also measured using the CIDI-SF and those who provided affirmative responses to three or more of the seven physiological symptoms were considered as "meeting the anxious criteria" [57].
Maternal substance use. Maternal substance use factors included the self-reported smoking during the past month, heavy drinking in past month and drug use in past month at Year-1. Those who provided affirmative responses to the single item "In past month, did you smoke cigarettes?" were considered as smokers. Heavy drinking was defined as reporting "5 or more drinks in one day" on one or more occasions in the past month. The drug use variable was created using the reposes to the two questions, "In past month did you smoke marijuana or pot?" and "In past month, did you use cocaine/crack/speed/LSD/heroin/other hard drug?". Affirmative responses to either or both questions were coded as drug use in the past month.
Covariates. Prior research findings have indicated that maternal age, race/ethnicity, nativity, marital status, multigenerational households and use of public assistance programs is associated with the risk of food insecurity [1,26,28,32,36,[58][59][60][61][62][63][64]. Thus, we have included maternal age (in years), race/ethnicity (non-Hispanic white, non-Hispanic black [reference], Hispanic and other), nativity (US born [reference] vs. foreign born), marital status (married [reference], cohabitating and single), grandparents living in the household, public food assistance program use (none [reference], use food stamps only, use WIC only, use both programs) and received welfare or TANF benefits during the past year. All covariates were measured at Year-1 except for maternal race/ethnicity and nativity measured at baseline and included in all adjusted models.
Statistical analysis. The descriptive statistics and all regression analysis were conducted using Stata version 16.1 (StataCorp, College Station, Texas). The LCA was conducted using Mplus version 8.4 (Muthen & Muthen) to derive risk profiles among the study sample. The method allows to identify latent classes/subgroups of mothers who are similar to each other based on the selected risk indicators.
The bivariate association between the risk indicator variables to be used in the LCA and food insecurity at Year-3 and -5 were evaluated using chi square test for categorical variables and t-tests for continuous variables. All variables listed were significantly related to food insecurity at both time points except for and maternal drug use (S1 Table). However, based on the consistent research evidence linking drug use [65][66][67] with the risk of food insecurity, we have decided to include all a priori identified risk indicators in the LCA.
A three-step approach was used in the LCA; 1). Identifying the best fitting LCA model to data, 2). Determining the latent class of each participant based on the model posterior probabilities, and 3). Evaluate the association between the assigned class membership and food insecurity. We examined the fit of six models (one-class to six-class) using four fit indices; the Bayesian Information Criteria (BIC), the sample size adjusted Bayesian Information Criteria (adjusted BIC), the Akaike Information Criteria (AIC), and Vuong Lo Mendell Rubin Likelihood Ratio Test (VLMR-LRT). Smaller values for BIC, adjusted BIC and AIC indicated a better model fit. A significant p-value for VLMR LRT indicates the (K-1)-class model has to be rejected in favor of a model with at least K-classes. In addition to the fit indices, we explored the identified class compositions to ensure the subgroups are meaningful. Further, the model entropy values closer to one (range: 0-1) was used as an indicator of a good classification.
After identifying the best fitting model and assigning each participant to the classes with highest membership probability, the class characteristics and food security were compared using chi square test for categorical variables and ANOVA test for continuous variables. Next, multinomial logistic regression analyses were conducted to evaluate the association between the class membership and level of food security at Year-3 (reference outcome = food secure) and the stability and change in food security during Years-3 to 5 (reference outcome = stable food security).

Study participants
The characteristics of the current study sample (n = 2,348) are reported in Table 1. At Year-3 approximately 13% of mothers reported low food security and 4% very low food security. There were 17% mothers reporting unstable food security during Years-3 and -5 and 8% reporting persistent food insecurity. A majority of the sample were non-Hispanic black (50%), US born (88%) with some college or higher education (42%), income <200% of the FPL (45%) and Medicaid/public health insurance (60%).

Latent class analysis
The five-class model was selected as the best fitting model based on the fit statistics and the composition of the classes ( Table 2). The significance tests of the VLMR LRT indicated that a five-class solution was a significantly better fit than a four-class solution but that the six-class solution was not a significantly better solution compared to five-class model. In addition, The BIC value was lowest in the five-class model, and the entropy value was higher for the fiveclass model compared to the six-class model indicating a better fit. The S2 Table shows the average individual posterior probabilities for being assigned to a specific latent class. The values off the diagonal are lower and the diagonal values are between 71% to 91% indicating a good classification.

PLOS ONE
The composition of the identified five risk classes and the bivariate comparisons between classes are reported in Tables 3 and 4.
The classes were labeled based on the highest reported risk factors and the descriptions of the five risk classes were as follows; Class 1: High utility and medical hardship class: This class consists of 57% of mothers with income below the 200% of the FPL, highest rates of utility hardship (83%), second highest rate of medical hardship (12%), and third highest rates of housing hardships (42%) of all other classes. However, in this class a second highest proportion of mothers and maternal grandparents had some college or higher education and highest rate of regular work for pay. Mothers also reported lower rates of health risk factors and substance use, except for moderate rates of depression and smoking.
Class 2: High housing and employment hardship, high substance use, and incarceration class: This class reports highest rates of partner incarceration (63%) and highest rates of mother substance use behaviors (60% smoking, 15% heavy drinking and 8% drug use). This class also consists of highest proportion of mothers with less than high school educated, second highest rate of no regular work for pay, second highest rates of mothers living under 200% of FPL, highest rate of housing hardship, and second highest rate of utility hardships.
Class 3: High housing and medical hardship, poor health, and health care class: This class reports highest rate of mothers with income below 200% of FPL (95%), no regular work for pay (68%) and medial hardships (28%) and second highest rate of housing hardships (49%). This class also consist of the highest rates of poor health (78%), health problems limiting work (44%), depression (63%), GAD (27%) and second highest rate of public health insurance use (84%).
Class 5: Low-risk class: Mothers in this class reported highest rates of protective factors and lower rates of all risk factors. This class consist of a highest proportion of mothers (85%) and maternal grandparents (59%) with some college or higher education, income at or above 300% of the FPL (64%), living in safe neighborhoods (97%), and access to private insurance (89%). Mothers also reported lowest rates of partner incarceration (7%), poor health (2%) and smoking (12%).
The distribution of the levels of food insecurity at Year-3 and -5 and the stability/change in food insecurity during Years-3 to -5 of the five risk classes are reported in Table 4. Mothers in the lowest risk class (Class 5) reported significantly lower rates of low and very low food  security at both time points and lower rates of persistent food insecurity. As there were only two mothers with very low food security (0.4%) at Year-3 and seven with persistent food insecurity during Year-3 and -5 (1.4%) among Class 5, they were excluded from the regression analysis due to inadequate sample size. Table 5 shows the results of the unadjusted and adjusted multinomial logistic regression models predicting level of food security at Year-3. According to the adjusted models, the odds of marginal food security among high utility and medical hardship class (Class 1), high housing and employment hardship, high substance use, and incarceration class (Class 2) and high housing and medical hardship, poor health, and health care class (Class 3) was not significantly different compared to the high employment hardship and lowincome class (Class 4). However, compared to the high employment hardship and lowincome class, all other risk classes had significantly greater odds of reporting low food security and very low food security. The post-estimation tests indicated that the adjusted odds ratios for low food security and very low food security was not significantly different among high utility and medical hardship class (Class 1), high housing and employment hardship, high substance use, and incarceration class (Class 2) and high housing and medical hardship, poor health, and health care class (Class 3). Table 6 shows the results of the unadjusted and adjusted multinomial logistic regression predicting stable/changing food insecurity. According to the adjusted model, compared to the high employment hardship and low-income class (Class 4), high utility and medical hardship class (Class 1), high housing and employment hardship, high substance use, and incarceration class (Class 2) and high housing and medical hardship, poor health, and health care class  PLOS ONE Table 4. The food security status, socio-demographic and household characteristics of the five latent classes identified in the LCA (n = 2,348). The post estimation tests indicated that the adjusted odds ratio of persistent food insecurity and unstable food insecurity was not significantly different among high utility and medical hardship class (Class 1), high housing and employment hardship, high substance use, and incarceration class (Class 2) and high housing and medical hardship, poor health, and health care class (Class 3).

Discussion
The purpose of this study was to identify SDOH risk classes among mothers with young children and to assess whether these risk classes are predictive of the severity and stability/change of food insecurity. This study adds to the food insecurity literature providing evidence of

PLOS ONE
subgroups of mothers with young children that have different levels of risk of food insecurity. Using LCA, we were able to identify five distinctive classes (1. high utility and medical hardship class, 2. high housing and employment hardship, high substance use, and incarceration class, 3. high housing and medical hardship, poor health, and health care class, 4. high employment hardship and low-income class and 5. low-risk class) of mothers varying by their educational, economic stability, social context, neighborhood and built environment, and maternal health, health care and substance use characteristics. O'Reilly, Hager, Harrington, et al. [50], identified three classes among African American urban households; low risk class (with low  probability of all risk factors: unemployment, low education, smoking, stress and depression), high stress/depression class and low education (high school or less) and less stress/depression class. These classes are similar to the classes identified (Class 5 with lower probability of all risks, Class 1 and 3 with higher health risks and two classes with high educational and employment risks, but low in health risks (Class 2 and Class 4) in the current study except, the current study sample was not limited to African American households and used greater number of variables. Thus, it does not warrant direct comparison of the identified classes. The racial/ethnically diverse sample of mothers with young children in the current sample adds to the food insecurity literature. The regression analysis revealed that compared to mothers in the high employment hardship and low-income class (Class 4), mothers in high utility and medical hardship class (Class 1), high housing and employment hardship, high substance use, and incarceration class (Class 2) and high housing and medical hardship, poor health, and health care class (Class 3) were at greater risk of low and very low food security. The models assessing the risk of change/stability of food insecurity showed that compared to mothers in high employment hardship and lowincome class (Class 4), mothers in mothers in high utility and medical hardship class (Class 1), high housing and employment hardship, high substance use, and incarceration class (Class 2) and high housing and medical hardship, poor health, and health care class (Class 3) had greater risk of experiencing unstable food insecurity or persistent food insecurity during years-3 to 5. The differences in risk between mothers in the high employment hardship and low-income class (Class 4) and other classes of mothers could be partially explained by the characteristics of the mothers in these classes. Mothers in Class 4 compared to mothers in Class 2 and 3 were also less likely to report material hardships, partner incarceration, depressive/anxious symptoms and substance use. Related, mothers in Class 1 report greater material hardships compared to mothers in Class 4, even though Class 1 mothers had several protective socioeconomic factors: higher percent of mothers with income 200-299% of the FPL (working poor), educated and employed. Material hardships, partner incarceration, depressive/anxious symptoms, and substance use have been independently linked to lower food insecurity risk [9,29,31,34,37]. Even though the mothers in high employment hardship and low-income class (Class 4) were similar to Classes 1, 2 and 3 in poverty, mothers in Class 4 were more likely to report being married or cohabitating and engaged in a higher use of WIC, which may also

PLOS ONE
explain their lower risk of food insecurity compared to mothers in other high-risk classes [59,68]. Last, a greater proportion of mothers in Class 4 were Hispanic and foreign-born, and also display a better health profile than mother in Class 1, 2, and 3. This finding concurs with the Healthy Immigrant Effect, where immigrants display better health than individuals from host country [69]. Prior research has shown that poor health and substance use behaviors may limit mother's financial stability and the ability to maintain steady employment [70] placing the family at risk of experiencing food insecurity and other material hardships [9,37,71]. The current findings are in line with previous findings [32,50] suggesting households experiencing multiple risks are at greater risk of lowest levels of food security and persistent food insecurity over time. Some families may present fewer risk factors associated with experiencing food insecurity but may have an overall higher risk for experiencing food insecurity compared to families that exhibit numerous risk and protective factors. For example, mothers experiencing other forms of material hardships aside from food insecurity, poorer physical or mental health status and reporting greater substance use are at greater risk of severe food insecurity and to persistency experience food insecurity compared to mothers with some economic disadvantages but less financial hardships and poor health risk factors. However, protective factors are limited in this sample. While greater education (Class 1 and 5), neighborhood safety and access to private insurance (Class 5) and use of public assistance (Class 4) was present, classes were not significantly different in their ability to have financial support. In other words, access to loans, a place to live, and emergency children care was not prevalent in any one of the classes. Thus, the economic instability experienced by mothers may be further exacerbated by the inability to financially depend on others.

Implications for practice
Identifying the characteristic profile of subgroups of families at greater risk of food insecurity and who were more likely to be affected by persistent food insecurity could be useful in providing resources to assist such families to improve access to adequate healthy food. The access to healthy food could be influenced by multiple domains including spatial-temporal, economic, social, service delivery and personal access [72]. The characteristic profiles predicting greater risk of food insecurity also provide an indication of possible domains and dimensions within the domains that could possibly hindering the food access of these families. Thus, the findings provide guidance for policy and environmental efforts needed to improve access to healthy foods among these families at risk. In clinical settings, pediatricians and other health care providers could use the two-question screener for food insecurity [73] to screen for children/families at risk of food insecurity. However, this screener may not be able to identify all potential families at risk as parents/children may not report or feel embarrassed to accept that they are struggling to provide food for the family [74,75]. Specifically, stigma, fear of child being taken away, and shame could prevent caregivers from reporting food insecurity [75]. According to McLeod, Vasudevan, Warnick, et al. [76], patients were more likely to be screened for food insecurity at health maintenance exam or at a new patient visit compared to return visits and may not fully capture food insecurity, particularly among at-risk populations. Thus, having a characteristics profile of families who could be at greater risk of food insecurity could be an additional tool that could aid health care providers to identify families that may need additional resources.

Study strengths and limitations
The current study uses a large ethnically diverse national dataset; however, the FFCWS is oversampled for non-marital births limiting the generalizability of the results. We were able to use multiple variables of educational, economic stability, social context, neighborhood and built environment, health, health care and substance use to develop the various latent profiles. However, this information was collected through self-reported data which are prone to both recall and social desirability bias. Therefore, the study findings need to be interpreted with caution. The longitudinal study design allowed to assess the future risk of food insecurity and make prediction for stability and change in food insecurity over two years' time period. While a SDOH framework was used, there was limited information in the data set about social context and community/neighborhood level factors, in addition to diverse measures of economic instability, such as various forms of debt (e.g., credit card, automobile). Adding such factors to LCA model could provide additional information about the families at risk of food insecurity. Last, this study focused on maternal characteristics and household factors. About 39% of the sample were single mothers and nonresidential father contributions were not considered. Previous research has demonstrated that a combination of nonresidential father formal and inkind support is associated with lower experiences of child food insecurity [77]. To gain a more wholistic view of families at risk for experiencing food hardship, identifying the risk classes (via LCA approach) among a sample of nonresidential fathers is needed. Using a SDOH framework would provide an opportunity to compare findings with this study and prior research [50].
In addition, it is important to acknowledge that the LCA approach categorizes individuals into classes based on probabilities. Therefore, it is hard to determine the exact number of individuals in each class and guarantee proper class assignment. However, the current study results could be used to identify possible classes among mothers with young children with differential risks of food insecurity.

Conclusions
LCA could be used to identify distinctive latent classes of mothers with young children who vary by their educational, economic stability, social context, neighborhood and built environment, and maternal health, health care, and substance use characteristics. Among the identified risk profiles, mothers in high utility and medical hardship class, high housing and employment hardship, high substance use, and incarceration class, and high housing and medical hardship, poor health, and health care class are at greater risk of experiencing low and very low food insecurity compared to mothers in high employment hardship and low-income class. Further, mothers in high utility and medical hardship class, high housing and employment hardship, high substance use, and incarceration class, and high housing and medical hardship, poor health, and health care class are at greater risk of experiencing persistent and unstable food insecurity compared to mothers in high employment hardship and low-income class. The generated risk profiles could be used by health care providers as an additional tool to identify families in need for resources to ensure household food security.
Supporting information S1 Table. Bivariate association between the indicator variables used in the LCA and food insecurity at Year-3 and Year-5 (n = 2,348). Reports the chi 2 test statistics for categorical variables and t-test statistics for continuous variables. a Drug use was not significantly associated with food insecurity in the current sample. However, decided to include it in the LCA due to the literature showing significant association with food insecurity. (DOCX) Hernandez.