Hostname: page-component-7c8c6479df-r7xzm Total loading time: 0 Render date: 2024-03-28T11:01:22.081Z Has data issue: false hasContentIssue false

Comprehensive analysis of the association of seasonal variability with maternal and neonatal nutrition in lowland Nepal

Published online by Cambridge University Press:  23 August 2021

Naomi M Saville*
Affiliation:
Institute for Global Health (IGH), University College London (UCL), London, UK
Mario Cortina-Borja
Affiliation:
Great Ormond Street Institute of Child Health (ICH), University College London (UCL), London, UK
Bianca L De Stavola
Affiliation:
Great Ormond Street Institute of Child Health (ICH), University College London (UCL), London, UK
Emma Pomeroy
Affiliation:
Department of Archaeology, University of Cambridge, Cambridge, UK
Akanksha Marphatia
Affiliation:
Department of Geography, University of Cambridge, Cambridge, UK
Alice Reid
Affiliation:
Department of Geography, University of Cambridge, Cambridge, UK
Dharma S Manandhar
Affiliation:
Mother and Infant Research Activities (MIRA), Kathmandu, Nepal
Jonathan CK Wells
Affiliation:
Great Ormond Street Institute of Child Health (ICH), University College London (UCL), London, UK
*
*Corresponding author: Email n.saville@ucl.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Objective:

To provide a comprehensive seasonal analysis of pregnant mothers’ eating behaviour and maternal/newborn nutritional status in an undernourished population from lowland rural Nepal, where weather patterns, agricultural labour, food availability and disease prevalence vary seasonally.

Design:

Secondary analysis of cluster-randomised Low Birth Weight South Asia Trial data, applying cosinor analysis to predict seasonal patterns.

Outcomes:

Maternal mid-upper arm circumference (MUAC), BMI, dietary diversity, meals per day, eating down and food aversion in pregnancy (≥31 weeks’ gestation) and neonatal z-scores of length-for-age (LAZ), weight-for-age (WAZ) and head circumference-for-age (HCAZ) and weight-for-length (WLZ).

Setting:

Rural areas of Dhanusha and Mahottari districts in plains of Nepal.

Participants:

2831 mothers aged 13–50 and 3330 neonates.

Results:

We found seasonal patterns in newborn anthropometry and pregnant mothers’ anthropometry, meal frequency, dietary diversity, food aversion and eating down. Seasonality in intake varied by food group. Offspring anthropometry broadly tracked mothers’. Annual amplitudes in mothers’ MUAC and BMI were 0·27 kg/m2 and 0·22 cm, with peaks post-harvest and nadirs in October when food insecurity peaked. Annual LAZ, WAZ and WLZ amplitudes were 0·125, 0·159 and 0·411 z-scores, respectively. Neonates were the shortest but least thin (higher WLZ) in winter (December/January). In the hot season, WLZ was the lowest (May/June) while LAZ was the highest (March and August). HCAZ did not vary significantly. Food aversion and eating down peaked pre-monsoon (April/May).

Conclusions:

Our analyses revealed complex seasonal patterns in maternal nutrition and neonatal size. Seasonality should be accounted for when designing and evaluating public heath nutrition interventions.

Type
Research Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press on behalf of The Nutrition Society

Highly seasonal environments impose multiple stresses that may impact various aspects of human nutrition, including dietary intake, physical activity patterns and somatic nutritional status(Reference Ulijaszek and Strickland1Reference Wells3). Seasonal effects during early development, especially in the first 1000 d, may have long-term consequences(Reference Madan, Haas and Menon4). Seasonal effects have previously been reported in maternal nutrition and birth weight(Reference Rush5Reference Moore, Fulford and Streatfield8), neonatal survival(Reference Moore, Fulford and Streatfield8) and physical activity and dietary intakes of pregnant mothers(Reference Rao, Kanade and Yajnik9).

In the plains of South Asia, weather patterns, agricultural labour activities, food availability and disease prevalence all vary greatly by season. During the pre-monsoon/monsoon seasons (April–September), high temperatures and humidity cause heat stress(Reference Pradhan, Shrestha and Shrestha10), associated with lower birth weight(Reference Wells and Cole11) and increased prevalence of infectious diseases(Reference Pradhan, Shrestha and Shrestha10,Reference Moffat12Reference Rohwerder15) . Conversely, low temperatures stimulate energy expenditure (heat production)(Reference Shrestha, Moore and Peel16), but are also associated with increased neonatal mortality from hypothermia(Reference Mullany, Katz and Khatry17), higher prevalence of acute respiratory infections(Reference Mathisen, Strand and Sharma18,Reference Banstola and Banstola19) and higher indoor air pollution from cooking and heating(Reference Devakumar, Semple and Osrin20,Reference Bartington, Bakolis and Devakumar21) , which may impact negatively upon birth weight(Reference Li, Huang and Jiao22).

Agricultural cycles are associated with seasonal food shortages(13) that may particularly impact women in South Asia, due to their large burden of agricultural work(Reference Panter-Brick23Reference Halbrendt, Kimura and Gray25) and their responsibility for preparing family meals(Reference Rao, Gokhale and Kanade26). During periods of peak agricultural work, pregnant women may experience both high energy expenditure and restricted food intake, due to lack of time to manage meals(Reference Rao and Raju27), while seasonal availability affects intake of micronutrient-rich foods(Reference Campbell, Talegawkar and Christian28). Intra-household allocation of food may become more inequitable during seasonal food shortages(Reference Harris-Fry, Shrestha and Costello29), which adversely affect women due to their lower household status(Reference Babu, Thirumaran and Mohanam30).

Whilst seasonal patterns of food security, energy expenditure and disease have been described in Nepal, previous studies have generally included only limited nutritional outcomes, resulting in a fragmentary overview. This study provides a comprehensive seasonal analysis of maternal and newborn nutrition in an undernourished population from the plains (Terai) to better understand the timing of nutritional vulnerability in mothers and their newborns, and inform interventions aimed at tackling undernutrition in pregnancy and infancy. We model seasonal patterns in somatic outcomes, pregnancy diet quality/quantity and eating behaviour using cosinor analysis.

Methods

Study setting and participants

We analysed data from a cluster-randomised controlled trial conducted with Maithili-speaking women in Nepal’s plains, where low maternal BMI (<18·5 kg/m2) and anaemia and child wasting are common, and women’s agency is poor(31). Supplemental Figure S1 illustrates monthly patterns of weather, agricultural activities and the lean season in our study area. There is substantial variability in temperature, humidity and precipitation, exerting contrasting composite ecological stresses by season.

The protocol, methods and results of the cluster-randomised controlled trial have been reported(Reference Saville, Shrestha and Style32Reference Style, Beard and Harris-Fry36). Briefly, the four-arm trial tested the impact of community interventions in pregnancy upon birth weight and infant weight-for-age z-score (WAZ) at the trial endpoint. The unit of randomisation was a ‘Village Development Committee’ geopolitical administrative area. All pregnant women in each cluster could enrol at any gestational age. Interventions tested involved participatory learning and action (PLA) women’s groups focussed upon nutrition in pregnancy, with and without transfers of food or cash to pregnant mothers. Whilst birth weight increased in PLA only, PLA + cash and PLA + food arms compared with the control arm, this was only significant in the PLA + food arm. No anthropometric differences were found at endpoint, when child age averaged 9 months.

The trial enrolled 25 090 women. All data were collected between January 2014 and April 2015 by local interviewers, educated mostly to bachelors’ level and given several weeks of training on interviewing, dietary assessment and anthropometry. Varying numbers were captured at different time points (Fig S2). Infants measured >3 d of birth were excluded, resulting in 3330 measures ≤72 h of birth and 2846 women in pregnancy. Loss to follow-up at different time points was due to data collector overload (more women enrolled than predicted), political unrest affecting data collector performance, women migrating to their parental home in late pregnancy for delivery and halting of funding before all outcomes could be collected. Technical error of measurement and coefficient of reliability for intra- and inter-observer variability was evaluated after each training session, and data collectors were retrained as needed.

Data collection

Maternal weight, height and mid-upper arm circumference (MUAC) after 31 weeks’ gestation and neonatal weight, length and head circumference (HC) were measured using: Tanita BD590 scales (infants, precision 10 g); Tanita solar scales (mothers, precision 100 g); ‘Shorr boards’ for length/height (precision 1 mm) and Seca 212 circumference tapes for MUAC and head circumference. We calculated anthropometric z-scores for WAZ, length-for-age (LAZ), and head circumference-for-age (HCAZ) and weight-for-length (WLZ) using WHO 2006 growth standards and cleaning criteria and maternal BMI(Reference Crowe, Seal and Grijalva-Eternod37). Measures were obtained across the year but with a gap from 20 May to 7 July 2014 inclusive, when training prevented data collection (Supplemental Fig S3). Whilst this gap spanned the late pre-monsoon and early monsoon, there were sufficient measures before and after the gap for these seasons to be well represented in the data.

During pregnancy (≥31 weeks), we asked mothers to recall intake of ten food groups in the previous 24 h using the standard FAO Women’s Dietary Diversity Score (WDDS) questionnaire(38), which was newly released at the time of the study. This tool has been shown to be a good proxy for mean probability of adequacy of diets in Zambia(Reference Caswell, Talegawkar and Siamusantu39) and amongst pregnant women in Bangladesh(Reference Nguyen, Huybregts and Sanghvi40). We also asked at ≥31 weeks: (i) if they had been eating less, more or the same amount in the last 2 weeks compared to pre-pregnancy; (ii) if they suffered severe loss of appetite or food aversion after the 4th month; (iii) if they had vomited in the last 2 weeks; (iv) how many main meals and snacks they had consumed in the preceding 24 h. At the early pregnancy interview (<31 weeks), we used the Months of Adequate Household Food Provisioning questions to estimate food security preceding and during early pregnancy(Reference Bilinsky and Swindale41). For those who reported food insecurity, we asked which Nepali months over the last year had shortage of food.

Statistical methods – modelling of seasonality

We fitted cosinor regression models to detect seasonality in annual and half-yearly (semestral) cycles and to describe the peak and nadir of selected outcomes(Reference Fernandez, Hermida and Mojon42Reference Padhye and Hanneman45). To avoid confusion arising from discrepancies in the definition of Nepalese seasons and use of Gregorian/Nepalese calendars, we analysed day in the year as a continuous variable.

Cosinor linear regression models for each maternal and newborn outcome were specified with annual and semestral cycles, to fully capture their seasonality. To maximise comparability across outcomes, we fitted both terms in all models.

We first create a day in the year variable, $0 \le d \le \;364$ for date of birth or measurement and transform it to angles:

$$t = 2 \times \pi \times d/365.$$

For simplicity, we define below the model for an annual cycle only. The full model specified with annual and semestral cycles and the model specification for binary outcomes are provided in Annexe S5.

An annual cycle model for a continuous outcome ${Y_{ij}}\left( t \right),$ measured in mother (or child) i belonging to cluster j for angle t (in radians, corresponding to transformed day of the year), is specified as a mixed-effects model:

$${Y_{ij}}\left( t \right) = {M_j} + {\beta _1}\;{\rm{sin}}\left( {{{2\pi t} \over D}} \right) + \;{\beta _2}\;{\rm{cos}}\left( {{{2\pi t} \over D}} \right) + {\gamma ^T}\;X + \;{\varepsilon _{ijt}},$$

where $\;{M_j}$ = $\;\left( {M + {u_j}} \right)$ is the random coefficient representing the average outcome (mesor) for cluster j, with $M$ the overall average while ${u_j}\;$ is a random effects term assumed to be normally distributed with mean 0 and variance $\sigma _u^2\;$ ; sin(.) and cos(.) are trigonometric functions of t (with $\pi = 3.14159$ and D = 365, the total number of days in a year); ${\beta _1}$ and ${\beta _{2\;}}$ are the regression coefficients for the annual cycle and ${\varepsilon _{ijt}}$ is a residual term, assumed to be normally distributed independently of ${u_j}$ with mean 0 and variance $\sigma _e^2$ . $X$ is a design matrix of covariates for parameters $\gamma $ .

This model is equivalent to the following linear regression representation:

$${Y_{ij}}\left( t \right) = {M_j} + {\rm{AMP}} \times {\rm{cos}}\left( {{{2\pi t} \over D} + \phi } \right) + \;{\varepsilon _{ijt}},$$

where:

AMP = Amplitude = $\sqrt {{{\left( {{\beta _1}\;{\rm{sin}}\left( {{{2\pi t} \over D}} \right)} \right)}^2}\; + \;{{\left( {{\beta _2}\;{\rm{cos}}\left( {{{2\pi t} \over D}} \right)} \right)}^2}} $

and

$$\phi = \;{\rm{Acrophase}} = {\rm{arctan}}\left( {{{ - \;{\beta _1}\;{\rm{sin}}\left( {{{2\pi t} \over D}} \right)} \over {\;{\beta _2}\;{\rm{cos}}\left( {{{2\pi t} \over D}} \right)}}} \right)$$

The acrophase $\phi $ can be converted into days (from the start of the cycle) as follows:

$${\rm{Acrocalendar}}\;{\rm{days}}\;{\rm{for}}\;{\rm{annual}}\;{\rm{seasonality}} = \left( {{{\left( { - \;\phi *365} \right)} \over {\;2\pi }}} \right)$$

The numerical day in the year can then be converted back into a calendar day (e.g. 100 is 10th April in a non-leap year).

Figure 1 illustrates the role played by mesor, amplitude and acrophase in the cosinor model for a particular cluster and shows how with this model there is one peak and one nadir in the values of ${Y_{ij}}\left( t \right)\;$ per year. Semestral cosinor terms (specifically, ${\rm{sin}}\left( {{{4\pi t} \over D}} \right)$ and ${\rm{cos}}\left( {{{4\pi t} \over D}} \right)$ where D = 182·5) allow for seasonal cycles with two peaks and two nadirs per year. Models were adjusted for age of the infant, gestational age in weeks, primigravid status of the mother and trial study arm.

Fig. 1 Schematic diagram of cosinor analysis output

These models provided: (i) coefficients for the sine and cosine functions that were then used to derive predicted seasonal patterns; (ii) yearly and semestral amplitudes describing the magnitude of seasonal variability measured from the mid-point in the cluster-specific mesor (presented in the units of measurement for continuous outcomes and in probabilities for binary outcomes (with CI derived from bootstrapped standard errors with 1000 samples)); (iii) yearly and semestral acrophases measuring distance from the start of the cycle (1 Jan) to its peak/nadir, which we converted into calendar days of the high/low points (with bootstrapped CI).

The P-values for pairs of annual, and pairs of semestral, sine and cosine functions are presented separately to enable the relative importance of annual and semestral cycles to be compared across outcomes.

We predicted seasonal patterns and their 95 % CI for each outcome by combining the relevant estimated parameters and their variance–covariance matrices and plotted them using the ggplot2 package in R(Reference Wickham46). For continuous variables, plots are presented in the units of measurement whereas plots for binary variables are presented as predicted probabilities of a positive recall for that outcome against day in the year. The scale of the y-axes is expanded to enable easy detection of patterns of predicted seasonal change and comparisons between plots.

Statistical methods – of covariates and outcomes

Our primary objective was to test for seasonality in dietary indicators in pregnancy and anthropometric outcomes in pregnancy and at birth, hence we did not build causal models. Before fitting cosinor regression models on our outcomes of interest, we tested for seasonality in the following variables to investigate whether they might confound the association between seasonality and the outcomes, by fitting regression models with a cosinor term: maternal height, maternal age, gestational age, primigravida status and infant age. We found significant seasonality in primigravida, infant age and gestational age, but not in maternal age or height. In Nepal, where childbearing is very rare outside of marriage, seasonality in primigravidae may be due to marriage seasonality whereby certain months are considered auspicious for marriage. Seasonality in infant and gestational ages may reflect differing timing of interviews. To reduce noise from these factors, we adjusted neonatal outcomes for log-transformed infant age (days), maternal outcomes for gestational age (weeks) and all outcomes for primigravid status. Model intercepts varied by the unit of randomisation to account for clustering of the outcomes within each unit (M j in the equations). Since trial interventions could potentially have influenced seasonal patterns, we also adjusted for study arm. Hence, each plot of predicted seasonal patterns for each outcome represents the reference values of multigravida and control study arm, and average gestational age (in maternal models) or child age (in newborn models).

Continuous outcomes included mothers’ MUAC, BMI, WDDS and number of main meals and snacks per day (meals) at ≥31 weeks’ gestation, and neonatal HAZ, WAZ, WLZ and HCAZ (≤72 h from birth) by day in the year. To get an indication of the peak and nadir values of continuous variables, it is possible to add or subtract the amplitude from the mesor. Binary outcomes, analysed with mixed-effects logistic regression cosinor models, were eating less than when not pregnant (eating down), vomiting, diarrhoea in the last 2 weeks, aversion to food/lack of appetite and consumption of each of nine individual food groups (meat/fish, dairy, eggs, pulses, nuts and seeds, green leafy vegetables (GLV), vitamin A-rich fruits and vegetables, other vegetables and other fruits) in the preceding 24 h, all measured ≥31weeks’ pregnancy.

Analyses were performed in Stata 15·1 (StataCorp LLC, College Station, Texas) and R 4·0·2 (The R Foundation for Statistical Computing, Vienna, Austria).

Results

Respondent characteristics

Table 1 contrasts respondent characteristics for women sampled at pregnancy and/or birth v. women not followed-up at these time points. Due to implementation challenges, such as political unrest and cessation of funding, only a small proportion of the 25 090 pregnancies enrolled were captured at ≥31 weeks’ gestation and/or within 3 d of birth. However, differences between captured and not captured groups were small (0–3 percentage points), so we deem the pregnancy and delivery samples representative of those enrolled in the study.

Table 1 Respondent characteristics

* Minimum age was 12 years in the sample analysed but there were only 1–3 individuals in the 12- and 13-year age categories. In the group that was not sampled at either time point, there was one mother each in the age categories of 9, 10 and 11 years old. It was not possible to verify the validity of these young maternal ages since respondents usually are unable to recall their date or birth.

The population comprises 84 % Hindu and 16 % Muslim, Maithili-speaking Terai castes, including 21 % marginalised Dalit groups and very few women of hills ethnicity. Three quarters of households had no toilet and ˜95 % cooked on a mixture of dried dung, firewood and agricultural residues that elevate indoor pollution. Most (86 %) belonged to farming families that sourced staple foods via production from their own or share-cropped land or through agricultural labour exchange. Most households purchase locally grown food (especially fruits and vegetables) from the local ‘haat bazaar’ market.

Landholding sizes were small, half a hectare on average. Over one-third of the pregnant mothers were primigravid. Most were adolescent (40 %) or in their early twenties (mean age 21·7 years, sd 4·7) with 8 % aged 15–16 years and just 1–2 and 2–3 individual mothers aged 12 and 13 years, respectively, in the analysed samples. The majority (89 %) had married below 18 years (mean 15·2 years, sd 2·1). Two-thirds of women had no schooling and 20 % had primary education only, while husbands had only slightly better educational opportunities. Mean gestational age at sampling in pregnancy was 35·6 weeks. After exclusion of newborns measured >72 h, mean age at sampling was 1·2 d. Boys comprised 53 % of the sample.

Supplemental Figure S4 shows recalled Nepalese month of food insecurity amongst 217 mothers who experienced months of food shortage in the year before the early pregnancy interview (recalling the period 27 Jan 2013–31 March 2014). The main ‘lean period’ falls in the post-monsoon months of mid-September to mid-November, just preceding rice harvest. Food insecurity falls to the lowest prevalence during the rice harvest in mid-November to mid-December and remains uncommon while rice stocks last until mid-February. From then (beginning of the hot pre-monsoon) until August (monsoon), food insecurity remains moderate with a slight improvement around pulse and wheat crop harvest in March/April. These food stocks remain up to the end of the monsoon at which point food insecurity peaks sharply.

Cosinor analyses

All fitted cosinor models showed evidence of significant clustering in the data, which required the use of random intercepts to resolve geographical clustering. Figs. 25 show predicted seasonal patterns for all relevant outcomes. Table 2 complements these plots, providing the relative size (amplitude), timing (acrophase) and statistical significance of the annual and semestral cycles of seasonality after adjustment for study arm, primigravida, gestational age and age of the infant. Table 2 also provides summaries over the whole sampling period including: mean newborn and pregnant mothers’ anthropometry, WDDS and meals; and probabilities of eating down, vomiting, food aversion, diarrhoea and food group consumption.

Fig. 2 Predicted seasonal patterns of neonatal anthropometric z-scores for length-for-age (LAZ), weight-for-age (WAZ), weight-for-length (WLZ) and head circumference-for-age (HCAZ) at birth plotted using cosinor analysis outputs analysed in relation to day of birth. a = significant annual cycle; s = significant semestral cycle; LAZ n 3291, WAZ n 3330; WLZ n 3029, HCAZ n 3200

Fig. 3 Predicted seasonal patterns of mothers’ BMI and MUAC in pregnancy plotted using cosinor analysis outputs in relation to day of measurement. a = significant annual cycle; s = significant semestral cycle; MUAC n 2816 and BMI n 2768; Note that the y axis range spans only 4 kg/m2 and 3 cm for BMI and MUAC, reflecting small yet statistically significant seasonal changes in maternal nutritional status

Fig. 4 Predicted seasonal patterns of mothers’ number of eating occasions and Women’s Dietary Diversity Score (WDDS), eating less than when not pregnant and food aversion in pregnancy plotted using cosinor analysis outputs in relation to day of measurement. a = significant annual cycle; s = significant semestral cycle; n 2831 for each outcome

Fig. 5 Predicted seasonal patterns of eight food groups included in the Women’s Dietary Diversity Score (WDDS) in pregnancy plotted using cosinor analysis outputs in relation to day of measurement. a = significant annual cycle; s = significant semestral cycle; n 2831

Table 2 Cosinor analysis results for annual and semestral (half-yearly) seasonality of newborn anthropometry and continuous and binary anthropometric and dietary outcomes in pregnancy

GLV, green leafy vegetables; HCAZ, head circumference-for-age z-score; LAZ, length-for-age z-score; meals/day, includes meals and snacks per day; MUAC, mid-upper arm circumference; WAZ, weight-for-age z-score; WDDS, Women’s Dietary Diversity Score; WLZ, weight-for-length z-score; VitA fruit/veg, vitamin A-rich fruits and vegetables: eating down: eating less than when not pregnant; diarrhoea recall in last 2 weeks; food aversion includes loss of appetite.

Predicted plots of seasonal variability combine both the annual and semestral patterns into one curve. They refer to the timing of the peaks and nadirs by outcome, whilst P-values and amplitudes are taken from Table 2.

Analysis of seasonality of neonatal anthropometry by month of birth

Table 2 shows the mean z-scores for neonatal LAZ, WAZ, WLZ and HCAZ. The percentages of newborns that were moderately and severely undernourished were: 12·3 % <−2 LAZ and 3·8 % <−3 LAZ (moderately and severely stunted, respectively); 11·8 % <−2 WAZ and 3·7 % <−3 WAZ (moderately and severely underweight, respectively); and 9·2 % <−2 WLZ and 3·2 % <−3 WLZ (moderately and severely wasted, respectively).

Figure 2 and Table 2 indicate significant annual and semestral cycles of seasonality for LAZ, WAZ and WLZ but not HCAZ at birth (adjusted for primigravida, age of the infant and study arm). The amplitude of the annual seasonal changes was larger (0·125–0·411 z-scores) than the semestral variation, which was close to 0·1 z-score but both are large enough to be of clinical significance. The annual nadirs and peaks calculated by adding or subtracting the amplitude to the mesor are: LAZ −1·019 to −0·770; WAZ −1·238 to −0·920; WLZ −0·778 to 0·045.

Children were born shorter but heavier in the winter (December/January), so had higher WLZ. Conversely, in the hot pre-/early-monsoon season, children had the lowest WLZ and WAZ (May/June), giving relatively high LAZ in July.

WLZ had the highest annual seasonal variation of 0·411 (95 % CI 0·283,0·540) z-scores (3–4 times greater than variability in WAZ or LAZ), with newborns least thin (higher WLZ) if born in the winter (January) and thinnest (lower WLZ) if born in the hot pre-monsoon season (May/June). The small yet significant semestral seasonality in WLZ was difficult to detect in plots. WAZ, like WLZ, was highest amongst those born in the winter (February) and lowest in those born in pre-monsoon (May/June) but the amplitude of seasonal WAZ variation was less than the variation in WLZ at 0·159 (95 % CI 0·085,0·233) z-scores. LAZ was the lowest (i.e. newborns were born shortest) in December, and 0·125 (95 % CI 0·022,0·227) higher in hot pre-monsoon and monsoon (March/April and July/August). A small dip in LAZ of 0·087 (95 % CI 0·012, 0·0162) amplitude was found during June, at the end of the hot pre-monsoon period (detected as a significant semestral cosinor term). HCAZ showed no significant seasonality in cosinor analysis, but in plots appeared very slightly lower in newborns born in the hot pre-monsoon (May) and higher in the winter (February) and monsoon (August).

Seasonal patterns in pregnant mothers’ anthropometry

Results of adjusted cosinor analysis of mothers’ anthropometry by month of measurement are given in Fig. 3 and Table 2. Twenty-eight percent and 8·3 % were thin (MUAC 21–22·9 cm) and very thin (<21 cm), respectively. Mothers’ MUAC and BMI showed significant annual seasonality with amplitudes of 22 mm and 0·27 kg/m2, respectively, dipping to their lowest points of 22·9 cm and 18·78 kg/m2 in October during the pre-harvest period of highest food insecurity (Fig S4). BMI peaked in winter a few months after rice harvest (February/March) at 19·32 kg/m2, while MUAC gradually increased to a peak of 23·33 cm between March and June (pre-/early-monsoon) before falling over the monsoon months (July–September).

Seasonal patterns in pregnancy diet

Figure 4 shows the combined annual and semestral predicted WDDS and meals and snacks (meals) per day and predicted probabilities of eating down and food aversion. WDDS showed larger and differing seasonal variation from other dietary outcomes. Meals per day varied by <0·5 between a trough in November during the busy food insecure harvest time, increasing over the winter to March and remaining high to peak in July/August (monsoon) when food security was moderate and agricultural (rice planting) workload heavy. Both eating down and food aversion were common (36 % and 29 %, respectively) and showed similar patterns to meals per day. They peaked in the hot/dry pre-monsoon season (April/May) and declined in the winter (November/December). Vomiting (9 %) and diarrhoea (13 %) in pregnancy showed no significant seasonality.

On average, mothers consumed 4·4 (SD 1·5) out of the 10 WDDS food groups in the preceding 24 h, one of which was invariably starchy staples including rice, wheat and potatoes. Cosinor results for WDDS showed significant seasonal variation with similar amplitude (of 0·22–0·25 food groups) in annual and semestral cycles. The cosinor plot showed the lowest WDDS in April during the extremely hot dry season, and the highest dietary diversity was in the late monsoon in August with a second peak in January and a smaller trough in October/November.

Consumption of food groups varied widely, except for starchy staples, consistently consumed daily year-round. The most commonly consumed groups in the preceding 24 h (Table 2 and Fig. 5) were other vegetables (82 %), pulses (71 %), dairy (59 %) and GLV (36 %). Nuts/seeds, meat/fish and other fruits were less frequently consumed in 27, 29 and 25 % of cases, respectively. Only 6 % of pregnant mothers ate eggs, while 8 % ate orange/yellow-coloured vitamin A-rich fruits and vegetables. Cosinor analysis of nine of the ten individual food groups that make up the WDDS (excluding starchy staples) (Table 2, Fig. 5) showed that seasonality varied widely by food group. Annual and semestral associations were found (in order of decreasing amplitude) for pulses, dairy, nuts/seeds and meat/fish. We found annual (but not semestral) seasonality for other vegetables and fruits, and semestral but not annual seasonality for GLV, vitamin A-rich fruits/vegetables and egg.

Plots in Fig. 5 demonstrate how seasonality varies by food group. Consumption of GLV had peaks in winter (January) and early monsoon (June) with troughs in April and September. Vitamin A-rich fruits and vegetables peaked in mango season in July with troughs in March and October and a much smaller peak in January. Meat/fish consumption was lowest in the hot season of May and highest at the time of the Dashain/Tihar festivals in September/October. Pulse consumption peaked in May/June after the April harvest, with the lowest consumption in the winter (February). Nuts/seeds and other vegetable consumption both dropped dramatically in May. Nuts/seeds consumption was the highest in February and September, whereas or other vegetables consumption stayed high between October and February. Conversely, other fruit consumption was highest in monsoon (July) and lowest in winter (December/January). Dairy consumption had a large peak in July and a smaller peak in January with troughs in April and October.

Discussion

We found strong and differing seasonal patterns in maternal and newborn anthropometry and maternal food aversion, eating down, dietary diversity and consumption of individual micronutrient-rich foods in pregnancy, in particular, vegetables, dairy, pulses, fruits and nuts/seeds. Whilst our study was not designed to attribute causal associations between maternal and neonatal traits, any seasonality in foetal phenotype must be mediated by maternal biology. Hence, below we highlight potential factors that might explain these associations. Fig. 6 summarises the timing of major peaks and nadirs for outcomes with significant seasonality and seasonal patterns of weather, agriculture and food security to assist with the visualisation of these potential associations. In general, we observed that offspring seasonal associations track maternal associations with some time lag, which appeared to vary by neonatal outcome. These patterns are likely complicated by the fact that both ecological and cultural factors can contribute to seasonal changes in eating behaviour, food security and workload.

Fig. 6 Schematic diagram of combined seasonal associations for pregnant mothers and neonates. Legend: lowest = major nadirs on plots; highest = major peaks on plots; ns = NS; green shading = positive effect on or indicator of nutrition; red shading = negative effect on or indicator of nutrition; W = winter; S = hot pre-monsoon; M = monsoon; PM =post-monsoon; H = harvest; P = rice planting; L = lean season

Seasonal variation in birth weight has been attributed to seasonality in food insecurity(Reference Hort47), food intake, eating habits, workload and disease exposure (e.g. malaria)(Reference Chodick, Flash and Deoitch48) although malaria is uncommon in our setting(Reference Dhimal, O’Hara and Karki49). In our study, food insecurity peaks before rice harvest (mid-September to mid-November), and pregnant mothers’ diets in October/November are low in vitamin A-rich fruits and vegetables, GLV and dairy. Conversely, the major festivals fall in these months with increased consumption of meat/fish, other vegetables and nuts/seeds. The long period of increasing food insecurity over the monsoon (July–September) which peaks in October, may lead to shortage of nutrients for foetal growth in pregnancy. Maximal foetal growth rates occur earlier in pregnancy for length than for weight(Reference Nyaruhucha, Msuya and Ngowi50Reference Ziegler, O’Donnell and Nelson52). Hence, being born short in the winter (December) may be associated with low MUAC and BMI of the mother some months earlier in October/November. Similarly, the LAZ peak in March occurs some months after the post-harvest improvements in food security and workloads. Poor maternal nutrition in the last trimester of pregnancy appears to be the likely cause of low birth weight in low-income settings(Reference Strand, Barnett and Tong53). Hughes et al found small for gestational age was the highest in Nepal in November, which coincides with the period of pre-harvest food insecurity(Reference Hughes, Katz and Mullany54).

In winter, post-rice harvest, a combination of increased food security, lack of food aversion or heat stress, lower workloads and increased consumption of GLV, dairy and vitamin-A rich vegetables/fruits may be associated with rapid improvements in mothers’ energy and nutrient intake and faster weight gain of the infant in utero. This is consistent with BMI and MUAC peaks after the end of winter in February/March (and the sustained higher MUAC until June) and with newborns being least underweight and thin in winter (January/February). Foetal weight gain in late pregnancy tends to be due to fat accretion(Reference Ziegler, O’Donnell and Nelson52,Reference Muthayya55) , which might explain the winter WLZ peak. The February/March maternal peak in BMI and MUAC coincides with, or slightly lags behind, the peaks in WLZ and WAZ of the newborn (February). We speculate that additional maternal intake may be channeled to the newborn before maternal weight gain, or that gestational age at birth may be higher in winter, but further studies are needed to investigate this.

In the hot, early monsoon (June), mean WLZ and WAZ fell markedly (i.e. increased thinness) together with a small dip in LAZ. The pre-monsoon (March–May) is a time of heat stress(Reference Pradhan, Shrestha and Shrestha10). High temperatures in pregnancy have been associated with low birth weight(Reference Wells and Cole11,Reference Basu, Rau and Pearson56Reference de Bernabé, Soriano and Albaladejo58) and preterm delivery(Reference Basu, Malig and Ostro59,Reference Carolan-Olah and Frankowska60) , hence the low WLZ and WAZ we detected in June could reflect more preterm deliveries. Other studies in South Asia found low birth weight and neonatal morbidity are higher(Reference Bang, Reddy and Baitule61), newborns smaller(Reference Hughes, Katz and Mullany54) and neonatal mortality higher(Reference Moore, Fulford and Streatfield8,Reference Hughes, Katz and Mullany54) in the summer/monsoon. In our study, despite relatively good food security in April/May, food aversion and eating down peaked, and dietary diversity was lowest, with troughs in consumption of meat/fish, GLV, vitamin A-rich fruits and vegetables, nuts and seeds, dairy and other vegetables. Only pulse consumption increased post-harvest in May. Pre-monsoon low availability of certain fruits and vegetables, short shelf-life and heat stress may contribute to decreased intakes, while harvesting of winter crops in April may increase energy expenditure(Reference Panter-Brick23). These factors may limit foetal nutrient availability and lead to babies being born thinnest in June, consistent with falling pregnancy BMI (but not MUAC) over this time.

Birth size has been associated with micronutrient-rich food intake(Reference Rao, Yajnik and Kanade62Reference Godfrey, Robinson and Barker64), as has gestational weight gain(Reference Wei, He and Lin65). In Bangladesh(Reference Stevens, Watt and Brimbecombe66) and Nepal(Reference Campbell, Talegawkar and Christian28,Reference Dulal, Mundy and Sawal67) , studies found seasonal variation in pregnant mothers’ dietary diversity and food group consumption. Consistent with our findings, Campbell et al found the lowest meat/fish consumption in the hot season (May), and peaks of consumption of GLV in winter and mango in summer(Reference Campbell, Talegawkar and Christian28). In Nepal, Jiang et al(Reference Jiang, Christian and Khatry68) found seasonality in pregnant women’s intake of specific micronutrients which is likely to affect foetal growth and development considerably. Fasting in Ramadan might mimic the effect of seasonal changes in food intake upon pregnant women and has been hypothesised to be associated with lower birth weight. However, findings on this are inconsistent(Reference Savitri, Yadegari and Bakker69,Reference Glazier, Hayes and Hussain70) with both higher(Reference Sakar, Balsak and Verit71) and lower(Reference Alwasel, Abotalib and Aljarallah72) placental weight reported in fasting women.

In Nepal, inequitable intra-household food allocation(Reference Sudo, Sekiyama and Maharjan73,Reference Gittelsohn74) reduces the adequacy of pregnant mothers’ diets compared to their male household heads(Reference Harris-Fry, Paudel and Shrestha75). These inequalities may increase during seasons of shortage, especially for micronutrient-rich foods groups(Reference Harris-Fry, Shrestha and Costello29) and during periods of peak agricultural activity(Reference Babu, Thirumaran and Mohanam30). The low status of young reproductive age women in our setting(Reference Gram76Reference Gram, Morrison and Sharma78) may magnify the impact of seasonal stresses, since ‘junior’ daughters-in-law usually cook and eat last, feeding micronutrient-rich foods to others while leaving little or none for themselves(Reference Morrison, Dulal and Harris-Fry79).

Although we did not detect seasonality in recall of diarrhoea in pregnant women, other studies in Nepal have detected increased diarrhoea in the monsoon(Reference Moffat12,13) and in association with increased temperature(Reference Bhandari, Bi and Sherchand80). Reduced appetite and increased demands on the body during illness may limit foetal growth and increase the risk of preterm delivery(Reference Vogel, Lee and Souza81,Reference Ashorn, Hallamaa and Allen82) . In rural India, higher maternal energy intakes, lower activity levels and exposure to winter (harvest) season predicted larger birth weight(Reference Rao, Kanade and Yajnik9). In Nepal, energy expenditure is highest during rice planting (pre-monsoon/monsoon) and lower in the winter(Reference Panter-Brick23,Reference Panter-Brick24,Reference Panter-Brick and Ellison83,Reference Panter-Brick84) and women gain weight when workloads decrease(Reference Panter-Brick84). This is consistent with our finding that MUAC and BMI are low in October at the end of the monsoon and higher after the rice harvest (from March). A study in Bangladesh found, like us, that women were thinnest in October(Reference Moore, Fulford and Streatfield8).

The seasonal effects we have detected upon neonatal z-scores for LAZ, WAZ and WLZ are of clinical significance and are equivalent to effect sizes of public health interventions designed to improve newborn size such as lipid nutrient(Reference Mridha, Matias and Chaparro85) and iron folic acid(Reference Nisar, Aguayo and Billah86) supplementation of pregnant women. Lower LAZ at birth predicts linear growth retardation(Reference Berngard, Berngard and Krebs87). The seasonal changes in MUAC and BMI amount to 0·45 cm and 0·54 kg/m2 from nadir to peak, respectively. Though small, these differences in nutritional status of pregnant mothers are large enough to make a difference to newborn size at a population level.

Our findings may be helpful in the design of targeted pregnancy interventions. For example, the increased risk of newborn thinness and underweight in the hot season may call for pre-monsoon protein-energy supplements to be provided to pregnant women together with interventions to reduce heat stress and illness in this season. Similarly, supplements for pregnant women over the monsoon until rice harvest might be helpful in preventing maternal thinness in October and newborn stunting in the winter.

With increased frequency of extreme heat waves, flooding and drought in Nepal associated with climate change(Reference Tiwari, Awasthi and Balla88,Reference Singh, Bassignana-Khadka and Singh Karky89) , seasonal stresses and related diseases may increase over time(Reference Bhandari, Bi and Sherchand80).

In general, our study confirms that Nepal has seasonal variation in the proximal determinants of undernutrition in the first 1000 d, as found in Madan et al’s review(Reference Madan, Haas and Menon4). However, we provide information on seasonality of eating behaviours such as food aversion and ‘eating down’ in pregnancy, which have not been described elsewhere, and provide a more comprehensive view of seasonal variability in pregnancy diet and both maternal and newborn anthropometry.

Our large, detailed dataset on maternal and newborn anthropometry, maternal diets, eating behaviours and food security enables seasonal assessment of many outcomes within the same population, and provides insights into associated ecological and cultural factors.

Limitations include insufficient cases of linked newborn and maternal anthropometry/diet to conduct analyses of seasonal associations on both child and maternal outcomes within the same models, so we cannot infer causality. However, since all foetal environmental stimuli and stresses are transduced by maternal phenotype(Reference Wells90), maternal exposures must in some way explain seasonality in newborn size. Low capture rates of enrolled women that resulted from premature stopping of data collection due to lack of funding, political unrest and migration of pregnant women to deliver at their parental home may have introduced biases, though these appear to be small(Reference Saville, Manandhar and Wells91). Since our sample was from Province 2, findings may have limited external validity to the hills or mountains of Nepal but are likely to be generalisable to Maithili-speaking populations across Nepal and North India. There was a gap in data collection between 20 May and 7 July 2014 and some months had higher intensity of data collection than others. Nevertheless, we detected seasonal patterns in maternal behaviour, diet and nutritional status and neonatal anthropometry.

Conclusions

Numerous components of maternal nutrition show seasonality in lowland Nepal, and this pattern extends to newborn size. Seasonality should therefore be accounted for when designing, implementing and evaluating public heath interventions to prevent maternal and child undernutrition in South Asia. Studies should collect data across different seasons and adjust for seasonality in analyses. Interventions could also deliberately target windows of seasonal vulnerability to wasting and anaemia. For longitudinal studies where data collection across the seasons cannot be arranged, data should be collected in the same season at each round, to avoid misinterpretation of apparent time trends that might be attributable to differing season of measurement.

Acknowledgements

Acknowledgements: We thank those involved in the Low Birth Weight South Asia Trial (LBWSAT) who made it possible for these data to be collected, especially the employees Mother and Infant Research Activities, in particular Bhim Prasad Shrestha, Aman Sen, Sonali Jha, Anjana Rai, Puskar Paudel, Bishnu Bhandari as well as 66 data collectors their 16 supervisors and 720 volunteer enumerators; and team members from University College London Institute for Global Health including Sarah Style, Helen Harris-Fry, B. James Beard, Andrew Copas, Joanna Morrison, Lu Gram, Jayne Harthan, Jolene Skordis-Worrall, David Osrin and Anthony Costello.

We also thank all those mothers and their families who gave up their time to participate in the study and permitted us to measure their newborns, and the District Public Health Offices of Dhanusha and Mahottari who supported the trial on the ground. Financial Support: Primary data collection for this study was funded by UKaid from Department for International Development (DFID) South Asia Research Hub (PO 5675) (https://www.gov.uk/dfid-research-outputs/protocol-of-the-low-birth-weight-south-asia-trial-lbwsat). Analyses conducted here were funded by The Leverhulme Trust as part of a project entitled ‘Age at marriage and human capital outcomes in mothers and offspring in Nepal’ (Leverhulme project number RPG-2017-264 https://www.leverhulme.ac.uk/listings?field_grant_scheme_target_id=13). This research was partly supported by the National Institute for Health Research (NIHR) Great Ormond Street Hospital Biomedical Research Centre (https://www.gosh.nhs.uk/our-research/nihr-great-ormond-street-hospital-brc). The views expressed are those of the authors and not necessarily those of the National Health Service (NHS), the NIHR or the UK Department of Health. The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript. Conflict of Interest: None.Authorship: JW and NS conceived the paper and the analyses. NS lead data collection on the ground, conducted all the analyses and wrote the first draft of the paper for publication. BDS advised on analysis in stata and wrote programmes for extracting results and drawing graphs. MCB advised on statistical analysis and programmed graphs in R. EP advised on analysis and interpretation providing initial excel spreadsheets for cosinor analyses. DSM managed the field data collection team in Nepal. All authors reviewed the manuscript, provided input and approved the paper for publication. Ethics of Human Subject Participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving research study participants were approved by the Nepal Health Research Council (292/2018) and the UCL Ethical Review Committee (0326/015). Ethical approval for the original study was obtained from the Nepal Health Research Council (108/2012) and the UCL Ethical Review Committee (4198/001). The Government of Nepal Ministry of Health and the Social Welfare Council granted permission for the primary data collection to be conducted. Written informed consent was obtained from all subjects

Supplementary material

For supplementary material accompanying this paper visit https://doi.org/10.1017/S1368980021003633

References

Ulijaszek, S & Strickland, S (editors) (1993) Seasonality and human ecology. In 35th Symposium of the Society for the Study of Human Biology, pp. 516. Cambridge, UK: Cambridge University Press.Google Scholar
Ferro-Luzzi, A & Branca, F (1993) Nutritional seasonality: the dimensions of the problem. In Seasonality and Human Ecology, pp. 149165 [Ulijaszek, S & Strickland, S, editors]. Cambridge, UK: Cambridge University Press.CrossRefGoogle Scholar
Wells, JC (2010) The Evolutionary Biology of Human Body Fatness: Thrift and Control. Cambridge Studies in Biological and Evolutionary Anthropology, vol. 58. Cambridge, UK: Cambridge University Press.Google Scholar
Madan, EM, Haas, JD, Menon, P et al. (2018) Seasonal variation in the proximal determinants of undernutrition during the first 1000 days of life in rural South Asia: a comprehensive review. Global Food Secur 19, 1123.CrossRefGoogle Scholar
Rush, D (2000) Nutrition and maternal mortality in the developing world. Am J Clin Nutr 72, 212S240S.Google ScholarPubMed
Prentice, AM, Cole, TJ, Foord, FA et al. (1987) Increased birthweight after prenatal dietary supplementation of rural African women. Am J Clin Nutr 46, 912925.CrossRefGoogle ScholarPubMed
Shaheen, R, de Francisco, A, Arifeen, S et al. (2006) Effect of prenatal food supplementation on birth weight: an observational study from Bangladesh. Am J Clin Nutr 83, 13551361.CrossRefGoogle ScholarPubMed
Moore, SE, Fulford, AJ, Streatfield, PK et al. (2004) Comparative analysis of patterns of survival by season of birth in rural Bangladeshi and Gambian populations. Int J Epidemiol 33, 137143.CrossRefGoogle ScholarPubMed
Rao, S, Kanade, AN, Yajnik, CS et al. (2009) Seasonality in maternal intake and activity influence offspring’s birth size among rural Indian mothers – Pune Maternal Nutrition Study. Int J Epidemiol 38, 10941103.CrossRefGoogle ScholarPubMed
Pradhan, B, Shrestha, S, Shrestha, R et al. (2013) Assessing climate change and heat stress responses in the Tarai region of Nepal. Ind Health 51, 101112.CrossRefGoogle ScholarPubMed
Wells, JC & Cole, TJ (2002) Birth weight and environmental heat load: a between-population analysis. Am J Phys Anthropol 119, 276282.CrossRefGoogle ScholarPubMed
Moffat, T (2003) Diarrhea, respiratory infections, protozoan gastrointestinal parasites, and child growth in Kathmandu, Nepal. Am J Phys Anthropol 122, 8597.CrossRefGoogle ScholarPubMed
WFP Nepal (2010) The Food Security Atlas of Nepal. Kathmandu: WFP.Google Scholar
Dahal, S (2008) Climatic determinants of malaria and kala-azar in Nepal. Regional Health Forum WHO South-East Asia Region 12, Special Issue on Protecting Health from Climate Change, 3237.Google Scholar
Rohwerder, B (2016) Seasonal Vulnerability and Risk Calendar in Nepal (GSDRC Helpdesk Research Report 1358). Birmingham, UK: GSDRC, University of Birmingham.Google Scholar
Shrestha, S, Moore, GA & Peel, MC (2018) Trends in winter fog events in the Terai region of Nepal. Agric For Meteorol 259, 118130.CrossRefGoogle Scholar
Mullany, LC, Katz, J, Khatry, SK et al. (2010) Incidence and seasonality of hypothermia among newborns in southern Nepal. Arch Pediatr Adolesc Med 164, 7177.CrossRefGoogle ScholarPubMed
Mathisen, M, Strand, TA, Sharma, BN et al. (2009) RNA viruses in community-acquired childhood pneumonia in semi-urban Nepal; a cross-sectional study. BMC Med 7, 35.CrossRefGoogle ScholarPubMed
Banstola, A & Banstola, A (2013) The epidemiology of hospitalization for pneumonia in children under five in the rural western region of Nepal: a descriptive study. PLoS One 8, e71311.CrossRefGoogle ScholarPubMed
Devakumar, D, Semple, S, Osrin, D et al. (2014) Biomass fuel use and the exposure of children to particulate air pollution in southern Nepal. Environ Int 66, 7987.CrossRefGoogle ScholarPubMed
Bartington, SE, Bakolis, I, Devakumar, D et al. (2017) Patterns of domestic exposure to carbon monoxide and particulate matter in households using biomass fuel in Janakpur, Nepal. Environ Pollut 220, 3845.CrossRefGoogle ScholarPubMed
Li, X, Huang, S, Jiao, A et al. (2017) Association between ambient fine particulate matter and preterm birth or term low birth weight: an updated systematic review and meta-analysis. Environmental Pollut 227, 596605.CrossRefGoogle ScholarPubMed
Panter-Brick, C (1993) Seasonality of energy expenditure during pregnancy and lactation for rural Nepali women. Am J Clin Nutr 57, 620628.CrossRefGoogle ScholarPubMed
Panter-Brick, C (1996) Seasonal and sex variation in physical activity levels among agro-pastoralists in Nepal. Am J Phys Anthropol 100, 721.3.0.CO;2-Y>CrossRefGoogle ScholarPubMed
Halbrendt, J, Kimura, AH, Gray, SA et al. (2014) Implications of conservation agriculture for men’s and women’s workloads among marginalized farmers in the central Middle Hills of Nepal. Mt Res Dev 34, 214222.CrossRefGoogle Scholar
Rao, S, Gokhale, M & Kanade, A (2008) Energy costs of daily activities for women in rural India. Public Health Nutr 11, 142150.CrossRefGoogle ScholarPubMed
Rao, N & Raju, S (2020) Gendered time, seasonality, and nutrition: insights from two Indian districts. Feminist Econ 26(2), 95125.CrossRefGoogle Scholar
Campbell, RK, Talegawkar, SA, Christian, P et al. (2014) Seasonal dietary intakes and socioeconomic status among women in the Terai of Nepal. J Health Popul Nutr 32, 198216.Google ScholarPubMed
Harris-Fry, H, Shrestha, N, Costello, A et al. (2017) Determinants of intra-household food allocation between adults in South Asia – a systematic review. Int J Equity Health 16, 107.CrossRefGoogle ScholarPubMed
Babu, SC, Thirumaran, S & Mohanam, TC (1993) Agricultural productivity, seasonality and gender bias in rural nutrition: empirical evidence from south India. Soc Sci Med 37, 13131319.CrossRefGoogle ScholarPubMed
Ministry of Health Nepal & New ERA ICF (2017) Nepal Demographic and Health Survey 2016. Nepal: Ministry of Health.Google Scholar
Saville, NM, Shrestha, BP, Style, S et al. (2016) Protocol of the low birth weight South Asia Trial (LBWSAT), a cluster-randomised controlled trial testing impact on birth weight and infant nutrition of Participatory Learning and Action through women’s groups, with and without unconditional transfers of fortified food or cash during pregnancy in Nepal. BMC Pregnancy Childbirth 16, 320.CrossRefGoogle ScholarPubMed
Saville, NM, Shrestha, BP, Style, S et al. (2018) Impact on birth weight and child growth of Participatory Learning and Action women’s groups with and without transfers of food or cash during pregnancy: findings of the low birth weight South Asia cluster-randomised controlled trial (LBWSAT) in Nepal. PLoS One 13, e0194064.CrossRefGoogle Scholar
Harris-Fry, HA, Paudel, P, Harrisson, T et al. (2018) Participatory women’s groups with cash transfers can increase dietary diversity and micronutrient adequacy during pregnancy, whereas women’s groups with food transfers can increase equity in intrahousehold energy allocation. J Nutr 148, 14721483.CrossRefGoogle ScholarPubMed
Gram, L, Morrison, J, Saville, N et al. (2019) Do participatory learning and action women’s groups alone or combined with cash or food transfers expand women’s agency in rural Nepal? J Dev Stud 55, 16701686.CrossRefGoogle ScholarPubMed
Style, S, Beard, BJ, Harris-Fry, H et al. (2017) Experiences in running a complex electronic data capture system using mobile phones in a large-scale population trial in southern Nepal. Global Health Action 10, 1330858.CrossRefGoogle Scholar
Crowe, S, Seal, A, Grijalva-Eternod, C et al. (2014) Effect of nutrition survey ‘cleaning criteria’ on estimates of malnutrition prevalence and disease burden: secondary data analysis. PeerJ 2, e380.CrossRefGoogle ScholarPubMed
FAO & FHI 360 (2016) Minimum Dietary Diversity for Women: A Guide for Measurement. Rome: FAO.Google Scholar
Caswell, BL, Talegawkar, SA, Siamusantu, W et al. (2018) A 10-food group dietary diversity score outperforms a 7-food group score in characterizing seasonal variability and micronutrient adequacy in rural Zambian Children. J Nutr 148, 131139.Google ScholarPubMed
Nguyen, PH, Huybregts, L, Sanghvi, TG et al. (2018) Dietary diversity predicts the adequacy of micronutrient intake in pregnant adolescent girls and women in Bangladesh, but use of the 5-Group Cutoff Poorly identifies individuals with inadequate intake. J Nutr 148, 790797.CrossRefGoogle ScholarPubMed
Bilinsky, P & Swindale, A (2007) Months of Adequate Household Food Provisioning (MAHFP) for Measurement of Household Food Access: Indicator Guide. Washington DC: Food and Nutrition Technical Assistance (FANTA) Project: Academy for Educational Development.Google Scholar
Fernandez, JR, Hermida, RC & Mojon, A (2009) Chronobiological analysis techniques. Application to blood pressure. Philos Trans A Math Phys Eng Sci 367, 431445.Google ScholarPubMed
Pomeroy, E, Wells, JC, Stanojevic, S et al. (2014) Birth month associations with height, head circumference, and limb lengths among Peruvian children. Am J Phys Anthropol 154, 115124.CrossRefGoogle ScholarPubMed
Lewis, KM, De Stavola, B & Hardelid, P (2020) Geospatial and seasonal variation of bronchiolitis in England: a cohort study using hospital episode statistics. Thorax 75, 262268.CrossRefGoogle ScholarPubMed
Padhye, NS & Hanneman, SK (2007) Cosinor analysis for temperature time series data of long duration. Biol Res Nurs 9, 3041.CrossRefGoogle ScholarPubMed
Wickham, H (2016) ggplot2: Elegant Graphics for Data Analysis. Cham, Switzerland: Springer International Publishing.Google Scholar
Hort, KP (1987) Seasonal variation of birthweight in Bangladesh. Ann Trop Paediatr 7, 6671.CrossRefGoogle ScholarPubMed
Chodick, G, Flash, S, Deoitch, Y et al. (2009) Seasonality in birth weight: review of global patterns and potential causes. Hum Biol 81(4), 463477.CrossRefGoogle ScholarPubMed
Dhimal, M, O’Hara, RB, Karki, R et al. (2014) Spatio-temporal distribution of malaria and its association with climatic factors and vector-control interventions in two high-risk districts of Nepal. Malar J 13, 457.CrossRefGoogle ScholarPubMed
Nyaruhucha, CN, Msuya, JM, Ngowi, B et al. (2006) Maternal weight gain in second and third trimesters and their relationship with birth weights in Morogoro Municipality, Tanzania. Tanzan Health Res Bull 8, 4144.Google ScholarPubMed
Valero De Bernabe, J, Soriano, T, Albaladejo, R et al. (2004) Risk factors for low birth weight: a review. Eur J Obstet Gynecol Reprod Biol 116, 315.CrossRefGoogle ScholarPubMed
Ziegler, EE, O’Donnell, AM, Nelson, SE et al. (1976) Body composition of the reference fetus. Growth 40, 329341.Google ScholarPubMed
Strand, LB, Barnett, AG & Tong, S (2011) The influence of season and ambient temperature on birth outcomes: a review of the epidemiological literature. Environ Res 111, 451462.CrossRefGoogle ScholarPubMed
Hughes, MM, Katz, J, Mullany, LC et al. (2014) Seasonality of birth outcomes in rural Sarlahi District, Nepal: a population-based prospective cohort. BMC Pregnancy Childbirth 14, 310.CrossRefGoogle ScholarPubMed
Muthayya, S (2009) Maternal nutrition & low birth weight – what is really important? Indian J Med Res 130, 600608.Google ScholarPubMed
Basu, R, Rau, R, Pearson, D et al. (2018) Temperature and term low birth weight in California. Am J Epidemiol 187, 23062314.CrossRefGoogle ScholarPubMed
Wells, JC (2002) Thermal environment and human birth weight. J Theor Biol 214, 413425.CrossRefGoogle ScholarPubMed
de Bernabé, JV, Soriano, T, Albaladejo, R et al. (2004) Risk factors for low birth weight: a review. Eur J Obstet Gynecol Reprod Biol 116, 315.CrossRefGoogle Scholar
Basu, R, Malig, B & Ostro, B (2010) High ambient temperature and the risk of preterm delivery. Am J Epidemiol 172, 11081117.CrossRefGoogle ScholarPubMed
Carolan-Olah, M & Frankowska, D (2014) High environmental temperature and preterm birth: a review of the evidence. Midwifery 30, 5059.CrossRefGoogle Scholar
Bang, AT, Reddy, HM, Baitule, SB et al. (2005) The incidence of morbidities in a cohort of neonates in rural Gadchiroli, India: seasonal and temporal variation and aa hypothesis about prevention. J Perinatol 25, S18S28.CrossRefGoogle Scholar
Rao, S, Yajnik, CS, Kanade, A et al. (2001) Intake of micronutrient-rich foods in rural Indian mothers is associated with the size of their babies at birth: Pune Maternal Nutrition Study. J Nutr 131, 12171224.CrossRefGoogle ScholarPubMed
Malhotra, N, Upadhyay, RP, Bhilwar, M et al. (2014) The role of maternal diet and iron-folic acid supplements in influencing birth weight: evidence from India’s National Family Health Survey. J Trop Pediatr 60, 454460.CrossRefGoogle ScholarPubMed
Godfrey, K, Robinson, S, Barker, DJ et al. (1996) Maternal nutrition in early and late pregnancy in relation to placental and fetal growth. BMJ 312, 410414.CrossRefGoogle ScholarPubMed
Wei, X, He, JR, Lin, Y et al. (2019) The influence of maternal dietary patterns on gestational weight gain: a large prospective cohort study in China. Nutrition 59, 9095.CrossRefGoogle ScholarPubMed
Stevens, B, Watt, K, Brimbecombe, J et al. (2017) The role of seasonality on the diet and household food security of pregnant women living in rural Bangladesh: a cross-sectional study. Public Health Nutr 20, 121129.CrossRefGoogle ScholarPubMed
Dulal, B, Mundy, G, Sawal, R et al. (2017) Homestead food production and maternal and child dietary diversity in Nepal: variations in association by season and agroecological zone. Food Nutr Bull 38, 338353.CrossRefGoogle ScholarPubMed
Jiang, TN, Christian, P, Khatry, SK et al. (2005) Micronutrient deficiencies in early pregnancy are common, concurrent, and vary by season among rural nepali pregnant women. J Nutr 135, 11061112.CrossRefGoogle ScholarPubMed
Savitri, AI, Yadegari, N, Bakker, J et al. (2014) Ramadan fasting and newborn’s birth weight in pregnant Muslim women in The Netherlands. Br J Nutr 112, 15031509.CrossRefGoogle ScholarPubMed
Glazier, JD, Hayes, DJL, Hussain, S et al. (2018) The effect of Ramadan fasting during pregnancy on perinatal outcomes: a systematic review and meta-analysis. BMC Pregnancy Childbirth 18, 421.CrossRefGoogle ScholarPubMed
Sakar, MN, Balsak, D, Verit, FF et al. (2016) The effect of Ramadan fasting and maternal hypoalbuminaemia on neonatal anthropometric parameters and placental weight. J Obstet Gynaecol 36, 483486.CrossRefGoogle ScholarPubMed
Alwasel, SH, Abotalib, Z, Aljarallah, JS et al. (2010) Changes in placental size during Ramadan. Placenta 31, 607610.CrossRefGoogle ScholarPubMed
Sudo, N, Sekiyama, M, Maharjan, M et al. (2006) Gender differences in dietary intake among adults of Hindu communities in lowland Nepal: assessment of portion sizes and food consumption frequencies. Eur J Clin Nutr 60, 469477.CrossRefGoogle ScholarPubMed
Gittelsohn, J (1991) Opening the box: intrahousehold food allocation in rural Nepal. Soc Sci Med 33, 11411154.CrossRefGoogle ScholarPubMed
Harris-Fry, HA, Paudel, P, Shrestha, N et al. (2018) Status and determinants of intra-household food allocation in rural Nepal. Eur J Clin Nutr 72, 15241536.CrossRefGoogle ScholarPubMed
Gram, L (2018) Women’s empowerment in a complex public health intervention in rural Nepal. PhD Thesis, University College London.Google Scholar
Gram, L, Skordis-Worrall, J, Saville, N et al. (2019) ‘There is no point giving cash to women who don’t spend it the way they are told to spend it’ – exploring women’s agency over cash in a combined participatory women’s groups and cash transfer programme to improve low birthweight in rural Nepal. Soc Sci Med 221, 918.CrossRefGoogle Scholar
Gram, L, Morrison, J, Sharma, N et al. (2017) Validating an agency-based tool for measuring women’s empowerment in a complex public health trial in rural Nepal. J Human Dev Capabil 18, 107135.CrossRefGoogle Scholar
Morrison, J, Dulal, S, Harris-Fry, H et al. (2018) Formative qualitative research to develop community-based interventions addressing low birth weight in the plains of Nepal. Public Health Nutr 21, 377384.CrossRefGoogle ScholarPubMed
Bhandari, D, Bi, P, Sherchand, JB et al. (2020) Assessing the effect of climate factors on childhood diarrhoea burden in Kathmandu, Nepal. Int J Hyg Environ Health 223, 199206.CrossRefGoogle ScholarPubMed
Vogel, JP, Lee, AC & Souza, JP (2014) Maternal morbidity and preterm birth in 22 low- and middle-income countries: a secondary analysis of the WHO Global Survey dataset. BMC Pregnancy Childbirth 14, 56.CrossRefGoogle ScholarPubMed
Ashorn, P, Hallamaa, L, Allen, LH et al. (2018) Co-causation of reduced newborn size by maternal undernutrition, infections, and inflammation. Matern Child Nutr 14, e12585.CrossRefGoogle Scholar
Panter-Brick, C & Ellison, PT (1994) Seasonality of workloads and ovarian function in Nepali women. Ann N Y Acad Sci 709, 234235.CrossRefGoogle ScholarPubMed
Panter-Brick, C (1995) Inter-individual and seasonal weight variation in rural Nepali women. J Biosoc Sci 27, 215233.CrossRefGoogle ScholarPubMed
Mridha, MK, Matias, SL, Chaparro, CM et al. (2016) Lipid-based nutrient supplements for pregnant women reduce newborn stunting in a cluster-randomized controlled effectiveness trial in Bangladesh. Am J Clin Nutr 103, 236249.CrossRefGoogle Scholar
Nisar, YB, Aguayo, VM, Billah, SM et al. (2020) Antenatal iron-folic acid supplementation is associated with improved linear growth and reduced risk of stunting or severe stunting in south asian children less than two years of age: a pooled analysis from seven countries. Nutrients 12(9), 26322651.CrossRefGoogle ScholarPubMed
Berngard, SC, Berngard, JB, Krebs, NF et al. (2013) Newborn length predicts early infant linear growth retardation and disproportionately high weight gain in a low-income population. Early Hum Dev 89, 967972.CrossRefGoogle Scholar
Tiwari, KR, Awasthi, KD, Balla, MK et al. (2010) Local people’s perception on Climate Change, its impact and adaptation practices in Himalaya to Terai regions of Nepal. Himalayan Res Pap Arch. https://digitalrepository.unm.edu/nsc_research/48/ (accessed January 2020).Google Scholar
Singh, SP, Bassignana-Khadka, I, Singh Karky, B et al. (2011) Climate Change in the Hindu Kush-Himalayas: The State of Current Knowledge. Kathmandu, Nepal: International Centre for Integrated Mountain Development (ICIMOD).Google Scholar
Wells, JC (2010) Maternal capital and the metabolic ghetto: an evolutionary perspective on the transgenerational basis of health inequalities. Am J Hum Biol 22, 117.CrossRefGoogle ScholarPubMed
Saville, NM, Manandhar, DS & Wells, JCK (2021) Trial experience and data capture in the Low Birth Weight South Asia Trial, a large cluster-randomised controlled trial in lowland Nepal. UCL Discovery. https://discovery.ucl.ac.uk/id/eprint/10120616/ (accessed August 2020).Google Scholar
Figure 0

Fig. 1 Schematic diagram of cosinor analysis output

Figure 1

Table 1 Respondent characteristics

Figure 2

Fig. 2 Predicted seasonal patterns of neonatal anthropometric z-scores for length-for-age (LAZ), weight-for-age (WAZ), weight-for-length (WLZ) and head circumference-for-age (HCAZ) at birth plotted using cosinor analysis outputs analysed in relation to day of birth. a = significant annual cycle; s = significant semestral cycle; LAZ n 3291, WAZ n 3330; WLZ n 3029, HCAZ n 3200

Figure 3

Fig. 3 Predicted seasonal patterns of mothers’ BMI and MUAC in pregnancy plotted using cosinor analysis outputs in relation to day of measurement. a = significant annual cycle; s = significant semestral cycle; MUAC n 2816 and BMI n 2768; Note that the y axis range spans only 4 kg/m2 and 3 cm for BMI and MUAC, reflecting small yet statistically significant seasonal changes in maternal nutritional status

Figure 4

Fig. 4 Predicted seasonal patterns of mothers’ number of eating occasions and Women’s Dietary Diversity Score (WDDS), eating less than when not pregnant and food aversion in pregnancy plotted using cosinor analysis outputs in relation to day of measurement. a = significant annual cycle; s = significant semestral cycle; n 2831 for each outcome

Figure 5

Fig. 5 Predicted seasonal patterns of eight food groups included in the Women’s Dietary Diversity Score (WDDS) in pregnancy plotted using cosinor analysis outputs in relation to day of measurement. a = significant annual cycle; s = significant semestral cycle; n 2831

Figure 6

Table 2 Cosinor analysis results for annual and semestral (half-yearly) seasonality of newborn anthropometry and continuous and binary anthropometric and dietary outcomes in pregnancy

Figure 7

Fig. 6 Schematic diagram of combined seasonal associations for pregnant mothers and neonates. Legend: lowest = major nadirs on plots; highest = major peaks on plots; ns = NS; green shading = positive effect on or indicator of nutrition; red shading = negative effect on or indicator of nutrition; W = winter; S = hot pre-monsoon; M = monsoon; PM =post-monsoon; H = harvest; P = rice planting; L = lean season

Supplementary material: File

Saville et al. supplementary material

Saville et al. supplementary material

Download Saville et al. supplementary material(File)
File 691 KB