Effect of integration of supplemental nutrition with public health programmes in pregnancy and early childhood on cardiovascular risk in rural Indian adolescents: long term follow-up of Hyderabad nutrition trial

Initial trial design (1987-90)

Using the opportunity afforded by the stepwise expansion of this programme during the 1980s and 1990s, the National Institute of Nutrition in India ran a trial to assess (among other things) the impact of food supplementation in pregnancy on the birth weight of offspring. A cluster of villages with a total population of 30 000 was chosen from each of the two adjacent administrative areas (called “blocks”), one of which already had the Integrated Child Development Services programme in place (intervention arm), whereas the other was awaiting implementation (control arm). As the 100 or more villages in each of the two blocks were spread over an unfeasibly large area, villages were chosen geographically for random selection: contiguous villages falling within a 10 km radius of a prominent central village (deemed to be the approximate centre of the block on visual inspection of the area map) in each block were selected. This process resulted in 15 villages from the intervention arm and 14 villages from the control arm being recruited to the study.

A 12 member team of investigators resided full time in the field for the duration of the study. Following a lead-in period of six months (involving training, pilot data collection, and household enumeration to identify women in the reproductive age group—13-45 years), the trial included all births in the area between 1 January 1987 and 31 December 1990. The “at risk” (of pregnancy) women were monitored monthly to detect missing of menstruation, and those identified as pregnant were followed closely during pregnancy until delivery (with clinical examinations in each of the three trimesters). The field team attempted to visit the home as soon as possible after delivery to collect data on the outcome of the pregnancy and to weigh the newborn within 48 hours. An infant beam balance with a 20 g accuracy (John Chatillon & Sons, NY, USA) was used to measure weight.

The supplement given in this trial was “upma,” a local preparation made from corn-soya blend (120 g) and soybean oil (16 g), providing 2.51 MJ and 20 g protein to the women and half this amount to the children. No other nutrients (such as micronutrients, iodised salt) were added to the supplement. The other universal programmes—immunisation, anaemia control in pregnancy through distribution of iron (60 mg elemental) and folic acid (500 μg) tablets, and the provision of basic health care—existed to a similar extent in both the intervention and control areas, although their uptake may be presumed to have been higher in the intervention area. Crossover of supplements from intervention to control villages was not a problem because of their separation by a large number of villages not involved in the study; in addition, supplements were offered only to the named residents of the village. A preliminary abstract reported the ongoing study,10 but full findings were unreported.

Follow-up survey (2003-5)

We designed this follow-up study to establish the status of women and their offspring who took part in the trial and to clinically examine offspring who were still resident in the area. The parents and children gave written informed consent.

We did two parallel surveys on a village to village basis—the first to identify the study participants and the second to do clinical examinations. The identification survey preceded the clinical survey by three months, and we interchanged the survey area between intervention and control villages every couple of months to even out any variations arising from seasons and experience of the study team. In the identification survey, we first identified women who took part in the baseline trial by using their own and their husband’s name. We interviewed each woman thus identified to collect details on children she had previously borne. All children born in these villages during the initial trial period (1987-90) were potential study participants. However, owing to limited resources and to make maximal use of the baseline data, we invited only those children who could be successfully matched to the previous records to have clinical examinations. We used date of birth and sex for matching, as infants’ names were not recorded in the baseline study (traditionally, naming is delayed and not done at the time of birth in India). Where multiple children of the same woman were eligible, we invited them all and made appropriate adjustments in the analyses.

We used a single clinic site at each village (generally a central point such as the village hall or health centre) to reduce measurement error arising from differences in clinic conditions. We held clinics in the morning, and invited 10-15 children to come fasting each day. We completed a full interviewer administered questionnaire and clinical examination on the children and a brief questionnaire on the mothers (who were asked to accompany the child).

Measurements

We measured height with a portable stadiometer (Leicester height measure; Chasmors, London). The participant stood erect with his or her head in the Frankfort plane, and a gentle upward pressure was applied under the mastoid. We used a digital weighing machine (HD 305; Tanita, Tokyo, Japan) to measure weight. We measured skinfold thickness at four sites (biceps, triceps, subscapular, and suprailiac) in triplicate with the Holtain calliper (supplied by Chasmors, London). We used a validated oscillometric device (HEM 705; Omron, Matsusaka Co, Japan) to measure blood pressure in the supine position with appropriate cuff sizes; we took two measurements and averaged them for analyses. We measured ambient room temperature with a digital thermometer.

We assessed the augmentation index, a measure of global arterial stiffness, by using an applanation tonometry technique with the Sphygmocor apparatus (Vx system; Atcor (PWV) Medical, Sydney, Australia).14 The characteristic pressure waveform produced by blood flow in the arteries changes as the arteries get stiffer with age or under conditions that lead to their premature stiffening (such as atherosclerosis).15 Pulse wave analysis can therefore be used to assess the functioning of the vascular tree. Pressure waveforms are obtained non-invasively by applying a pressure sensitive probe over a peripheral artery and transformed into the corresponding central arterial waveform by using a generalised transfer function validated against invasive pressure recordings.14 The augmentation index (difference between the first and second peaks of the central arterial waveform, expressed as a percentage of pulse pressure) is negative in healthy young adults but becomes increasingly positive as arteries stiffen.15 We measured the augmentation index over the radial artery in the supine position, taking an average of two high quality recordings (quality index >90%).

We classified sexual maturation into four stages (corresponding to Tanner’s early, middle, and late puberty and post-puberty) on the basis of time since the onset of menstruation (girls) and testicular volume (boys).16 The boys self assessed testicular volume in private, using Prader’s orchidometer (a chain of 12 wooden testes with volumes ranging from 1 ml to 25 ml). We measured socioeconomic position with the standard of living index, which is a household level, asset based scale devised for use in India.17 We classified the participants as having low (0-14), medium (15-24), or high (25-67) socioeconomic position. We assessed urbanisation of the villages by population size (people) and classified them into low (<2000), medium (2000-5000), and high (>5000).

We collected fasting blood samples (at least eight hours) in appropriate vacutainers, transferred them within one to two hours (in icebox at 4-8° C), and processed them within four hours. We did assays for glucose, triglycerides, total cholesterol, and high density lipoprotein cholesterol on the same day with an autoanalyser (ACE Clinical System; Schiapparelli Biosystems, NJ, USA) and the recommended kits (Alfa Wasserman, NJ, USA). We estimated insulin concentrations by radioimmunoassay in batches within four to six weeks.18

Quality control

We produced detailed protocols and used them regularly to standardise the work of the fieldwork team. “Blinding” of the fieldworkers to the group assignment was not an option; however, all measurements except anthropometry were to a large extent automated (for example, biochemical assays, blood pressure, and arterial stiffness), thus reducing the possibility of bias. Only one observer made each measurement (including biochemical assays), to eliminate interobserver bias. We assessed reproducibility of clinic measurements by repeating the measurements on a random subsample (5%) of participants after one to three weeks and found it to be consistently high (intraclass correlation coefficients of >0.98 for anthropometric measures and >0.85 for blood pressures and augmentation index). We put in place internal and external quality control arrangements (with Cardiac Biochemistry laboratory at the All India Institute of Medical Sciences, Delhi, which is part of UKNEQAS system coordinated from Newcastle, UK) for biochemical assays, and split assays were done on a 5% subsample (intraclass correlation coefficients >0.94). We validated testicular self assessment technique against a trained observer in a separate substudy and found it to be highly accurate (mean difference in model ranks (self reported minus directly observed)—0.07 (95% confidence interval −0.11 to 0.25)—with no evidence of systematic bias on a Bland-Altman plot).

Statistical analyses

We used the log of the sum of four skinfolds to calculate the percentage of body fat19 20; we then converted this into fat and fat-free mass by using body weight and expressed them as the corresponding indices after dividing them by the square of height in metres.21 We assessed central adiposity by the ratio of central (subscapular plus suprailiac) to peripheral (biceps plus triceps) skinfolds. We estimated low density lipoprotein cholesterol from triglycerides, total cholesterol, and high density lipoprotein cholesterol by using the Friedewald-Fredrickson equation,22 and we calculated insulin resistance by homoeostasis model assessment (HOMA), excluding those with fasting glucose ≥7 mmol/l.23

We applied suitable transformations to outcome variables that deviated markedly from a normal distribution (triglycerides, insulin, and HOMA scores). We used linear regression models to investigate association of supplemental nutrition with cardiovascular disease risk factors. Analysis was on an intention to treat basis, using area of birth as proxy (irrespective of whether the participant took the supplement or not). We fitted four predefined models to adjust incrementally for the main domains of potential confounding or intermediary variables: model 1 (physiological variables—age, sex, pubertal stage); model 2 (current socioeconomic status—household standard of living index, village population); model 3 (stature—height); and model 4 (body composition—fat mass index, fat-free mass index, and central-peripheral skinfold ratio). We additionally adjusted blood pressure for ambient room temperature and augmentation index for heart rate, as these factors can artefactually affect their values.24 25 We assessed association of supplemental nutrition with sexual maturation (by ordinal and binary logistic regression), as improved nutrition may influence risk of cardiovascular disease through earlier onset of puberty.26

Data could be clustered at the level of the village and the household (in case of multiple children). To take account of village level clustering, we used robust standard errors in all the models, with village as the level of cluster. This technique uses the cluster level residuals to derive the standard errors, so that the resulting standard errors are larger and valid in the presence of clustering but the parameter estimate remains unchanged. This technique was inappropriate for household level clustering, as very few households had multiple children taking part in the study. We therefore examined the impact of household level clustering on the results by excluding the second of the two children from the same household (no household contained more than two children). We similarly estimated the impact of a few migrant children (who usually lived outside the area but happened to be at home during the study) on the results by excluding these children from the final models. We examined interaction between the intervention and the sex of the participant, as evidence exists of preferential feeding of male children in this setting,27 as well as possible sex differential effects of early undernutrition on risk of cardiovascular disease.28 To retain maximum study power, we examined only one interaction. We used Stata version 9 for statistical analyses.

Sample size calculations done before the start of the study suggested that the anticipated sample (estimated as 1268 overall) was going to be adequate to detect important differences in most outcomes (for example, at 80% power, 5% significance level, and 0.01 intraclass correlation coefficient for village level clustering, the mean detectable differences were 1.8-2.9 mm Hg for systolic blood pressure, 1.5% for augmentation index, and 1.7-3.2 mU/l for insulin).