Is increasing urbanicity associated with changes in breastfeeding duration in rural India? An analysis of cross-sectional household data from the Andhra Pradesh children and parents study

Objective To investigate whether village-level urbanicity and lower level socioeconomic factors are associated with breastfeeding practices in transitioning rural communities in India. Setting 29 villages in Ranga Reddy district, southern India between 2011 and 2014. Participants 7848 children under 6 years identified via a cross-sectional household survey conducted as part of the Andhra Pradesh Children and Parents Study. Outcome measures Two key indicators of optimal breastfeeding: termination of exclusive breastfeeding before 6 months and discontinuation of breastfeeding by 24 months. Village urbanicity was classified as low, medium or high according to satellite assessed night-light intensity. Results Breastfeeding initiation was almost universal, and approximately two in three children were exclusively breastfed to 6 months and a similar proportion breastfed to 24 months. Using multilevel logistic regression, increasing urbanicity was associated with breastfeeding discontinuation before 24 months (medium urbanicity OR 1.45, 95% CI 0.71 to 2.96; high urbanicity OR 2.96, 95% CI 1.45 to 6.05) but not with early (<6 months) termination of exclusive breastfeeding. Increased maternal education was independently associated with both measures of suboptimal breastfeeding, and higher household socioeconomic position was associated with early termination of exclusive breastfeeding. Conclusion In this transitional Indian rural community, early stage urbanicity was associated with a shorter duration of breastfeeding. Closer surveillance of changes in breastfeeding practices alongside appropriate intervention strategies are recommended for emerging economies.

Objective To investigate whether village-level urbanicity and lower level socioeconomic factors are associated with breastfeeding practices in transitioning rural communities in India. setting 29 villages in Ranga Reddy district, southern India between 2011 and 2014. Participants 7848 children under 6 years identified via a cross-sectional household survey conducted as part of the Andhra Pradesh Children and Parents Study. Outcome measures Two key indicators of optimal breastfeeding: termination of exclusive breastfeeding before 6 months and discontinuation of breastfeeding by 24 months. Village urbanicity was classified as low, medium or high according to satellite assessed night-light intensity. results Breastfeeding initiation was almost universal, and approximately two in three children were exclusively breastfed to 6 months and a similar proportion breastfed to 24 months. Using multilevel logistic regression, increasing urbanicity was associated with breastfeeding discontinuation before 24 months (medium urbanicity OR 1.45, 95% CI 0.71 to 2.96; high urbanicity OR 2.96, 95% CI 1.45 to 6.05) but not with early (<6 months) termination of exclusive breastfeeding. Increased maternal education was independently associated with both measures of suboptimal breastfeeding, and higher household socioeconomic position was associated with early termination of exclusive breastfeeding. Conclusion In this transitional Indian rural community, early stage urbanicity was associated with a shorter duration of breastfeeding. Closer surveillance of changes in breastfeeding practices alongside appropriate intervention strategies are recommended for emerging economies.

IntrOduCtIOn
The promotion of breastfeeding is one of the three interventions identified as having the largest potential impact on global child death. 1 Optimal breastfeeding is defined by the WHO as early breastfeeding initiation, exclusive breastfeeding (EBF) to 6 months and continued breastfeeding to 2 years or beyond alongside appropriate complementary feeding. Many low-income and middle-income countries (LMICs) have a strong tradition of near universal and prolonged breastfeeding, 2 3 though EBF to 6 months (hereafter referred to simply as 'exclusive breastfeeding') is less common. A small increase in the global proportion of children exclusively breastfed between 1995 and 2010 has been reported, but the overall proportion (40%) still falls strikingly short of universal coverage and obscures differences in country-specific trends. 4 strengths and limitations of this study ► Previous studies have investigated the association between urbanisation and breastfeeding using the urban-rural dichotomy. ► We used data from a large rural cohort in southern India that is currently undergoing rapid and uneven urbanisation due to its proximity to a major urban centre. ► The use of night-time light intensity data as an indicator of urbanicity allowed us to examine subtler changes in breastfeeding practices along the urbanrural continuum. ► Sixteen per cent of children were excluded from the analysis due to missing information on breastfeeding practices. ► We relied on maternal retrospective recall of breastfeeding events for our outcome measurement.

Open Access
Many LMICs are currently experiencing a rapid increase in the proportion of people living in built-up areas, and the social, cultural and economic changes associated with this process of urbanisation have the potential to impact on traditional breastfeeding practices. Direct threats to optimal breastfeeding include early introduction of other liquids and inappropriate supplementation with solid or semisolid foods. These behaviours may be influenced by changing social norms, for example, increasing numbers of mothers working outside the home. Of all positive health behaviours, breastfeeding is one of the few more prevalent in LMICs compared with high-income countries. 3 Within LMICs, this trend is mirrored by a higher prevalence of suboptimal breastfeeding in urban areas compared with rural areas 5 : a trend also observed in India 6-8 alongside variation by various socioeconomic indicators. [6][7][8][9][10] Although the high-level urbanrural comparison is of interest, there may also be subtler changes in breastfeeding practices along the urbanrural continuum given the periurban effects on villages close to urban centres. These changes can potentially be investigated by using a measure of 'urbanicity' that aims to assess the extent of urbanisation in a given area. A number of different indicators of urbanicity have evolved, including the use of remote light sensing 11 12 and multicomponent scales. 13 14 The early identification of changes in breastfeeding practices accompanying the urbanicity transition-and an understanding of the underlying mechanisms-are necessary for informing appropriate interventions to protect traditionally positive breastfeeding practices in transitioning communities.
The Andhra Pradesh Children and Parents Study (APCAPS) is a rural sociodemographic cohort in southern India that is currently undergoing rapid and uneven urbanisation due to its proximity to a major urban centre (Hyderabad), providing a unique opportunity to examine the association between early stage urbanicity and breastfeeding practices.
MethOds study design APCAPS is an intergenerational cohort originally established to study the long-term effects of early-life undernutrition on risk of cardiovascular disease and subsequently expanded to include transgenerational influences of other environmental and genetic factors on chronic diseases in transitioning rural India.
The original cohort is based on the participants in the Hyderabad Nutrition Trial conducted in 1987-1990 in 29 villages approximately 50-100 km from Hyderabad in Telangana state (formally Andhra Pradesh), southern India. 15 The dataset used in this analysis is based on a cross-sectional household survey conducted between 2011 and 2014 in the study villages. All households (household defined as a group of people living in the same residence and sharing a common kitchen) in the study villages were visited by fieldworkers, and sociodemographic information was collected on each household. In addition, a basic health profile was collected for each child under 6 years of age, comprising information on infant feeding (colostrum intake, total duration of breastfeeding and age of onset of weaning), immunisation and anthropometric measurements. Fieldworkers made repeated visits to households to maximise response and to clarify inconsistencies in collected data. Data were collected from 23 314 households in total, of which 5968 (25.6%) included at least one child under 6 years.
The study received approval from the ethics committees of the National Institute of Nutrition (Hyderabad, India) and London School of Hygiene and Tropical Medicine (London, UK). Approval was also sought from the Indian Council for Medical Research and the village heads and their committees in each of the study villages. Written informed consent (or witnessed thumbprint if illiterate) was obtained from the participants prior to their inclusion in the study.
breastfeeding outcomes and explanatory variables Two breastfeeding outcomes were used in this analysis: termination of EBF before 6 months and discontinuation of breastfeeding before 24 months. These outcomes reflect failure to achieve two of the specific WHO recommendations for optimum feeding practices (EBF to 6 months and continued breastfeeding to 2 years). 16 As part of the basic health profile for children compiled for children under 6 in the household survey, mothers were asked to report the total duration of breastfeeding (in months), and the age (in months) at onset of weaning. 'Weaning' was defined by fieldworkers as the age at which the child was given anything other than mother's milk, that is, age at initiation of complementary feeding. A copy of the questions used in the survey is provided as a supplementary figure (see online supplementary figure S1).
Our primary explanatory factor was urbanicity, measured using remotely sensed village-level night-time light intensity (NTLI) scores, as these are objective, unbiased and easily available over wide areas. Although this analysis represents the first application of NTLI data to the APCAPS population, NTLI data are increasingly being used as an area-based indicator of socioeconomic development. 11 12 The light that is included in the NTLI score include any outside lights, ranging from fires and gas flares to lights related to human settlements. Low-level lights such as from streets and car headlights can be observed if there is a sufficient number of sources, but indoor lights cannot be observed. NTLI scores were calculated for 2012 using the National Oceanic and Atmospheric Administration Stable Lights product that provides yearly average NTLI measures processed and filtered to remove events such as fires and lightening contamination by cloud or moon reflections and background noise, at a 1 km resolution. Scores for each village were calculated by summing the raw NTLI values over each village polygon (digitised using Bing Maps combined with GPS-based surveying by the field teams). The 1 km resolution NTLI data were Open Access upscaled to 100 m resolution to allow more accurate estimation of the NTLI values covered by each polygon, as many villages are small and partially cover multiple 1 km grid cells.
NTLI scores for each village were validated against alternative urbanicity measurements (field worker ranking and a multicomponent urbanicity score based on household-level material assets and village-level availability of infrastructure and services) showing positive correlations (0.65 and 0.53, respectively). Study villages were ranked by their NTLI score and divided into tertiles to represent 'low'' (10 villages), 'medium' (10 villages) and 'high' (nine villages) levels of urbanicity. The NTLI tertile scores matched the field worker ranking in 50% of the villages, and cases of disagreement between NTLI and field worker ranking, the latter was more conservative and ranked villages as medium urbanicity rather than high urbanicity. Only one village had a significant divergence between NTLI and fieldworker ranking.
In addition, we investigated mother-level socioeconomic factors that may be correlated with urbanicity: maternal education (no formal education, primary education, or secondary education and higher), maternal employment (paid work vs no paid work) and a household level standard of living index (SLI). Asset-based SLIs have been established as a valid proxy measure of household wealth. 17 We generated an SLI score for each household, calculated by using information on household assets including house and land ownership, characteristics of the home (electricity, water pump, separate kitchen and separate toilet) and ownership of various assets (tractor, radio, AC, washing machine, bore hole, telephone, TV, fridge, bicycle, two wheeler, four wheeler, bank account, animal cart, sofa, cot/bed, mattress and table). Principal component analysis was used to determine the weights for each component in the index, 18 and households were divided into quintiles according to their weighted score. We also report data on a number of other factors likely to be associated with breastfeeding practices: sex of child, birth order, maternal age (grouped) and household composition (joint/extended or nuclear).

statistical analysis
We included in the analysis all children under 6 years who were breastfed at least once and for whom information was available on feeding history. The analysis investigating termination of EBF before 6 months was restricted to children 6 months or older at the time of survey, and correspondingly only those children aged 24 months or older were included in the analysis of discontinuation of breastfeeding before 24 months. A small proportion of children (3%) had missing information on one or more variables of interest and were excluded.
We hypothesised that urbanicity would be associated with less favourable breastfeeding practices. This could operate through at least two different indirect pathways (see figure 1): through increasing individual-level employment, education or assets so that households are less likely to maintain breastfeeding, or due to more urbanised villages have a different 'collective' attitude to breastfeeding. We investigated these hypotheses by using multilevel logistic regression modelling with children (level 1) nested within mothers (level 2, max n=5477) nested within villages (level 3, n=29). This approach allowed us to model the variation in breastfeeding outcomes at each level (random effects) and to estimate the effect of specific mother and village-level factors on breastfeeding practices (fixed effects). We initially fitted a null model (model 1) for each of the two outcomes with random intercepts only in order to estimate the baseline between-mother and between-village variance. We then fitted a series of models for each breastfeeding outcome, adding covariates as fixed effects to the included random effects, where fixed effects were interpreted as the average effect on the specified breastfeeding outcome across all mothers and villages. These models included individual demographic factors and mother-level socioeconomic indicators (model 2), individual demographic factors and village-level urbanicity (model 3) and all variables (model 4). Due to the correlation between socioeconomic indicators and urbanicity, we considered estimates from model three our main results. Proportional change in variance was calculated as a measure of change in mother-level (level 2) and village-level (level 3) variance between the null model and subsequent models, and (for village-level variance only) the measure of change between a model with (model 4) and without (model 2) the village-level urbanicity variable included.
Estimates of the association between mother-level socioeconomic variables and breastfeeding outcomes were derived from model 2 (adjusted for individual-level demographic variables, but not urbanicity).
We hypothesised a priori that the association between village-level urbanicity and breastfeeding may vary by household SLI and maternal education. We investigated these cross-level interactions in further models (for SLI, comparing the richest two quintiles to the three poorest quintiles; for education, comparing secondary education vs no or primary education).
All statistical analyses were conducted using Stata version14. Open Access results Characteristics of the sample Information on breastfeeding was available on a total of 7848 children (5390 households), 99% (n=7839) of whom were breastfed at least once (figure 2).
The characteristics of ever-breastfed children by urbanicity of village are presented in table 1. There was little variation in infant sex, birth order or maternal age by urbanicity of village. Children residing in villages classified as more urbanised had mothers that were more likely to have been educated to secondary level, less likely to have mothers in paid employment and a higher standard of living. Joint/extended families were slightly less prevalent in high urbanicity villages. termination of exclusive breastfeeding by 6 months Among the 7142 children no longer exclusively breastfed (88 children ≥6 months were still exclusively breastfed at the time of survey), the mean age at termination of EBF was 6.1 months (SD 1.8), median 6.0 months and IQR 5-6 months. One-third of children (33.5%, n=2420) were EBF for a period of less than 6 months (table 1).

Fixed effects
There was no statistically significant trend regarding early termination of EBF and village level urbanicity.
The prevalence of early termination of EBF was lowest in medium urbanicity villages (27.2%), higher in lower urbanicity villages (33.6%) and highest in high urbanicity villages (36.5%). In multivariable analysis, there was no evidence that urbanicity was associated with termination of EBF by 6 months (model 3, table 2), with little change in estimates after the addition of demographic and socioeconomic covariates to the model.
After adjustment for other individual-level and mother-level covariates, both children of mothers with primary education and children of mothers with secondary education were more likely to be EBF for less than 6 months when compared with children of mothers with no formal education (primary education: OR 3.37, 95% CI 2.13 to 5.31; secondary education: OR 1.69, 95% CI 1.12 to 2.54; model 2, table 2). Increasing SLI quintile was associated with up to twice the odds of early termination of EBF compared with children from the poorest households (richest quintile OR 2.11, 95% CI 1.22 to 3.63; p value for trend=0.003). There was some evidence that maternal employment was also associated with early termination of EBF (OR 1.43, 95% CI 1.00 to 2.03). The estimates for socioeconomic variables did not change with the addition of urbanicity to the model.

Interaction effects
There was no evidence of interaction between urbanicity and either household SLI or maternal education.

Random effects
There was statistically significant unexplained variance estimates at both the mother and village level (model 1, table 1). Unexplained variability was consistently higher at the mother level compared with the village level. The addition of individual level parameters resulted in a slight decline in community-level variation (variance 1.50 and 1.44 in models 1 and 2, respectively). There was a further decline in village-level variance when the urbanicity variable was added to the model (variance 1.32 in model 4). Comparing the village-level variance between model 2 and model 4 suggests that 8.5% of the observed villagelevel variation can be explained by urbanicity.

discontinuation of breastfeeding by 24 months
At the time of survey, 784 children aged ≥24 months were still being breastfed. Among those children no longer breastfed, the mean and median age at discontinuation of breastfeeding was 21.4 (SD 8.5) and 24 months, respectively, and the IQR was 15-24 months. Nearly 4 in 10 children (37.8%, n=2037) were breastfed for less than 24 months in total.

Fixed effects
Discontinuation of breastfeeding by 24 months was more common in high urbanicity villages (42.1%) and least common in low urbanicity villages (29.6%). After adjustment for individual-level demographic factors, high urbanicity was associated with increased odds of breastfeeding discontinuation before 24 months (OR 2.64, 95% CI 1.29  When compared with children of mothers with no formal schooling, children of mothers with secondary education were at significantly higher odds of breastfeeding discontinuation after adjustment for all demographic and socioeconomic factors (OR 1.63, 95% CI 1.23 to 2.16; model 2, table 3). Maternal employment was associated with a slight reduction in the odds of breastfeeding discontinuation before 24 months (OR 0.77, 95% CI 0.60 to 0.99). There was no evidence that SLI quintile was independently associated with breastfeeding discontinuation by 24 months. The inclusion of urbanicity in the model did not alter the socioeconomic estimates of effect.

Interaction effects
There was no evidence of interaction between urbanicity and mother-level socioeconomic factors (household SLI and maternal education).

Random effects
The random effects parameters for models investigating discontinuation of breastfeeding before 24 months are presented in table 3. In the null model (model 1), the proportion of residual variance attributable to mothers (level 2, 53.7%) was much higher than the variance attributable to villages (level 3, 8.5%). The addition of urbanicity to a model including individual and mother-level factors resulted in a decrease of 8.5% in village-level variance.

dIsCussIOn summary of main findings
In this study, approximately two in three children were exclusively breastfed to 6 months and a similar proportion breastfed to 24 months. At the village level, high urbanicity was associated with breastfeeding discontinuation before 24 months, but there was no evidence that urbanicity was associated with early termination of EBF. At the mother level, increased maternal education was independently associated with both indicators of suboptimal breastfeeding and high SLI associated with an increased odds of EBF for less than 6 months. Maternal employment showed a variable association with breastfeeding. The residual variation in breastfeeding outcomes suggested greater heterogeneity within villages than between villages.

Consistency with previous studies
Our estimates of breastfeeding prevalence are largely consistent with those derived from other population-based †Restricted to ever-breastfed infants aged at least 24 months still breastfed/with age at cessation of breastfeeding (figure 3) (n=5382). Table 1 Continued Table 2 Results of multilevel logistic regression models for the association between urbanicity or individual/household socioeconomic factors and termination of exclusive breastfeeding <6 months, among ever-breastfed children under six in the APCAPS household survey*.   19 and a study of 600 mother-child pairs in Andhra Pradesh reports that 75% of infants aged 3-5 months were exclusively breastfed. 20 Some of the younger infants included in these two study samples will have ceased breastfeeding by 6 months, suggesting that our study sample has a slightly higher proportion of EBF to 6 months. The overall proportion of children breastfed until at least 24 months in our study was almost identical to an analysis of all-India NFHS-2 data: (62.2% vs 63%). 6 Very few existing studies have investigated the association between urbanicity and breastfeeding. In one study based in the Philippines, Dahly and Adair reported that length of breastfeeding was negatively correlated with increasing urbanicity (using a multicomponent measure). 13 The persisting association between high urbanicity and increased odds of breastfeeding discontinuation <24 months-after adjustment for lower level socioeconomic circumstances-reported in our study support the findings from Dahly and Adair.
Increasing urbanicity is associated with positive socioeconomic changes such as improved education for women and increased income and household wealth. A number of other studies from India and other LMICs have demonstrated a negative association between improved socioeconomic position and breastfeeding practices. 6-8 21-24 We found similar results with regard to household SLI and increased maternal education, and early termination of EBF. One explanation for this trend could be the greater affordability and/or social desirability of commercial breastmilk substitutes. The association between education and early termination of EBF was strongest for primary education (primary education: OR 3.37, 95% CI 2.13 to 5.31; secondary education: OR 1.69, 95% CI 1.12 to 2.54). This suggests that while education in general is associated with a reduction in the length of EBF, higher levels of education partially ameliorate this effect. Interestingly, there was some evidence that maternal employment had a protective effect on breastfeeding discontinuation by 24 months, though the opposite trend was observed with regard to early cessation of EBF. There is some evidence of a U-shaped association between education and women's employment in India, with paid employment outside the home common among women with little or no formal education, lower among women with moderate levels of education and rising again with high levels of education. 25 Mothers in employment are likely to be a heterogeneous group, making it difficult to draw any firm conclusions about the association between paid employment and breastfeeding practices in this sample.

strengths and limitations
The APCAPS cohort provides a unique opportunity to investigate current health behaviour and outcomes set against the backdrop of rapid urbanisation and economic transition in rural India. While studies based on the high level urban-rural comparison may help to predict the impact of 'total' urbanisation on breastfeeding practices, they obscure the temporal emergence of subtler changes in the urban environment that may be amenable to intervention. The use of multilevel models enabled us to explore the role of factors at different levels: individual, mother and village.
Although the vast majority of all under 6s in the study villages were included in our analysis, 15.7% (n=1464) were excluded due to missing information on feeding history due to the mother living elsewhere, travelling or deceased. In a comparison of included and excluded children, there was no evidence that infant sex, infant age, number of under 6s living in the household or household SLI differed by missing status (see online supplementary table 1). A slightly higher proportion of excluded children resided in high urbanicity villages (p 0.09).
We relied on maternal recall of breastfeeding events for our outcome measurement. For the analysis of EBF at 6 months, the recall period ranged from 0 to 5.5 years, and for breastfeeding continuation at 24 months the recall period was 0-4 years. A review of 11 studies assessing the validity and reliability of maternal recall of breastfeeding concluded that maternal recall of breastfeeding duration is good, especially when the recall period is short (<3 years). 26 A more recent study, conducted in a population where breastfeeding initiation was near universal and duration long, found that even after 20 years, 64% of women recalled duration correctly to within 1 month (90% within 3 months). 27 However, there is some evidence that recall of age at introduction of complementary foods or non-breastmilk fluids is less accurate. 26 It is unclear whether any misclassification of breastfeeding behaviour is independent of other characteristics, but where differential misclassification has been suggested, more highly educated or wealthier mothers have tended to over-report breastfeeding. 28 Given that these characteristics were associated with suboptimal breastfeeding practices in this study, we may have underestimated any true difference in breastfeeding by sociodemographic characteristics.
Our measure of urbanicity was derived from NLTI data, information that is objective, regularly updated and free to use. Additionally, data on NLTI are available over a number of years and could be used in future studies to investigate trends in urbanicity over time. However, it must be noted that urbanicity is an ecological indicator and as such may not accurately reflect individual environment, particularly given that many women may travel regularly outside their home village for work or family reasons.

Implications
Nearly a quarter (24%) of all global under-five deaths occur in India. 29 In light of the failure to achieve the Millennium Development Goal infant mortality rate (IMR) target reduction, 30 a new target of reducing the Open Access IMR to 20 per 1000 live births by 2020 has recently been proposed. 31 Early results from the latest NFHS-4 data collected in Telegana state report a current IMR of 28 (20 in urban areas, 35 in rural areas). 19 An increase in optimal breastfeeding practices will help to achieve improvements in infant survival, in addition to reducing the considerable burden of infant morbidity. 3 32 India faces an ever-increasing epidemic of chronic disease in common with many other LMICs. Several studies have suggested that breastfeeding has a protective effect on long-term outcomes such as obesity and diabetes in adulthood, 33 though residual confounding is difficult to exclude, 34 and the most recent data from the PROBIT RCT do not support an association between breastfeeding and adiposity in late childhood. 35 A substantial proportion of infants in India are exclusively breastfed for less than the 6 months recommended by WHO, 7 9 10 and a recent study reported that there was little change in the prevalence of EBF in India between 1992-1993 and 2005-2006 8 . The lack of country-specific holistic and coordinated policy programmes supporting breastfeeding has also been highlighted. 36 Therefore, research to further understand the determinants of suboptimal breastfeeding practices in India is timely.
Our findings suggest that in LMICs with a strong tradition of breastfeeding, negative changes in breastfeeding behaviour may be observed during early stages of the urbanicity transition. Reduced duration of breastfeeding among more educated mothers may be one of the earliest markers of this change. India is currently undergoing rapid urbanisation, with the proportion of the population living in towns and cities is set to increase from an estimated 28% in 2011 to 38% by 2026 37 . Many more individuals live in areas which though traditionally described as rural are increasingly displaying many of the characteristics of urban areas. There is good evidence that breastfeeding behaviours are amenable to change through interventions delivered at the household and community level, as well as those targeting health systems. 38 39 Intervention programmes to protect and promote breastfeeding should be considered in transitioning communities to counteract changes in breastfeeding practices, preferably targeted at those mothers identified as most at risk of suboptimal breastfeeding practices.