Urinary and Dietary Analysis of 18,470 Bangladeshis Reveal a Correlation of Rice Consumption with Arsenic Exposure and Toxicity

Background

We utilized data from the Health Effects of Arsenic Longitudinal Study (HEALS) in Araihazar, Bangladesh, to evaluate the association of steamed rice consumption with urinary total arsenic concentration and arsenical skin lesions in the overall study cohort (N=18,470) and in a subset with available urinary arsenic metabolite data (N=4,517).

Methods

General linear models with standardized beta coefficients were used to estimate associations between steamed rice consumption and urinary total arsenic concentration and urinary arsenic metabolites. Logistic regression models were used to estimate prevalence odds ratios (ORs) and their 95% confidence intervals (CIs) for the associations between rice intake and prevalent skin lesions at baseline. Discrete time hazard models were used to estimate discrete time (HRs) ratios and their 95% CIs for the associations between rice intake and incident skin lesions.

Results

Steamed rice consumption was positively associated with creatinine-adjusted urinary total arsenic (β=0.041, 95% CI: 0.032-0.051) and urinary total arsenic with statistical adjustment for creatinine in the model (β=0.043, 95% CI: 0.032-0.053). Additionally, we observed a significant trend in skin lesion prevalence (P-trend=0.007) and a moderate trend in skin lesion incidence (P-trend=0.07) associated with increased intake of steamed rice.

Conclusions

This study suggests that rice intake may be a source of arsenic exposure beyond drinking water.     

Genome-Wide Association Study of Parity in Bangladeshi Women

Human fertility is a complex trait determined by gene-environment interactions in which genetic factors represent a significant component. To better understand inter-individual variability in fertility, we performed one of the first genome-wide association studies (GWAS) of common fertility phenotypes, lifetime number of pregnancies and number of children in a developing country population. The fertility phenotype data and DNA samples were obtained at baseline recruitment from individuals participating in a large prospective cohort study in Bangladesh. GWAS analyses of fertility phenotypes were conducted among 1,686 married women. One SNP on chromosome 4 was non-significantly associated with number of children at P <10^-7 and number of pregnancies at P <10^-6. This SNP is located in a region without a gene within 1 Mb. One SNP on chromosome 6 was non-significantly associated with extreme number of children at P <10^-6. The closest gene to this SNP is HDGFL1, a hepatoma-derived growth factor. When we excluded hormonal contraceptive users, a SNP on chromosome 5 was non-significantly associated at P <10^-5 for number of children and number of pregnancies. This SNP is located near C5orf64, an open reading frame, and ZSWIM6, a zinc ion binding gene. We also estimated the heritability of these phenotypes from our genotype data using GCTA (Genome-wide Complex Trait Analysis) for number of children (h_g ² = 0.149, SE = 0.24, p-value = 0.265) and number of pregnancies (h_g ² = 0.007, SE = 0.22, p-value = 0.487). Our genome-wide association study and heritability estimates of number of pregnancies and number of children in Bangladesh did not confer strong evidence of common variants for parity variation. However, our results suggest that future studies may want to consider the role of 3 notable SNPs in their analysis.