Genetic Variants Associated with Serum Thyroid Stimulating Hormone (TSH) Levels in European Americans and African Americans from the eMERGE Network

Thyroid stimulating hormone (TSH) hormone levels are normally tightly regulated within an individual; thus, relatively small variations may indicate thyroid disease. Genome-wide association studies (GWAS) have identified variants in PDE8B and FOXE1 that are associated with TSH levels. However, prior studies lacked racial/ethnic diversity, limiting the generalization of these findings to individuals of non-European ethnicities. The Electronic Medical Records and Genomics (eMERGE) Network is a collaboration across institutions with biobanks linked to electronic medical records (EMRs). The eMERGE Network uses EMR-derived phenotypes to perform GWAS in diverse populations for a variety of phenotypes. In this report, we identified serum TSH levels from 4,501 European American and 351 African American euthyroid individuals in the eMERGE Network with existing GWAS data. Tests of association were performed using linear regression and adjusted for age, sex, body mass index (BMI), and principal components, assuming an additive genetic model. Our results replicate the known association of PDE8B with serum TSH levels in European Americans (rs2046045 p = 1.85×10⁻¹⁷, β = 0.09). FOXE1 variants, associated with hypothyroidism, were not genome-wide significant (rs10759944: p = 1.08×10⁻⁶, β = −0.05). No SNPs reached genome-wide significance in African Americans. However, multiple known associations with TSH levels in European ancestry were nominally significant in African Americans, including PDE8B (rs2046045 p = 0.03, β = −0.09), VEGFA (rs11755845 p = 0.01, β = −0.13), and NFIA (rs334699 p = 1.50×10⁻³, β = −0.17). We found little evidence that SNPs previously associated with other thyroid-related disorders were associated with serum TSH levels in this study. These results support the previously reported association between PDE8B and serum TSH levels in European Americans and emphasize the need for additional genetic studies in more diverse populations.     

Role of Nicotine Dependence on the Relationship between Variants in the Nicotinic Receptor Genes and Risk of Lung Adenocarcinoma

Several variations in the nicotinic receptor genes have been identified to be associated with both lung cancer risk and smoking in the genome-wide association (GWA) studies. However, the relationships among these three factors (genetic variants, nicotine dependence, and lung cancer) remain unclear. In an attempt to elucidate these relationships, we applied mediation analysis to quantify the impact of nicotine dependence on the association between the nicotinic receptor genetic variants and lung adenocarcinoma risk. We evaluated 23 single nucleotide polymorphisms (SNPs) in the five nicotinic receptor related genes (CHRNB3, CHRNA6, and CHRNA5/A3/B4) previously reported to be associated with lung cancer risk and smoking behavior and 14 SNPs in the four ‘control’ genes (TERT, CLPTM1L, CYP1A1, and TP53), which were not reported in the smoking GWA studies. A total of 661 lung adenocarcinoma cases and 1,347 controls with a smoking history, obtained from the Environment and Genetics in Lung Cancer Etiology case-control study, were included in the study. Results show that nicotine dependence is a mediator of the association between lung adenocarcinoma and gene variations in the regions of CHRNA5/A3/B4 and accounts for approximately 15% of this relationship. The top two CHRNA3 SNPs associated with the risk for lung adenocarcinoma were rs1051730 and rs12914385 (p-value = 1.9×10⁻¹⁰ and 1.1×10⁻¹⁰, respectively). Also, these two SNPs had significant indirect effects on lung adenocarcinoma risk through nicotine dependence (p = 0.003 and 0.007). Gene variations rs2736100 and rs2853676 in TERT and rs401681 and rs31489 in CLPTM1L had significant direct associations on lung adenocarcinoma without indirect effects through nicotine dependence. Our findings suggest that nicotine dependence plays an important role between genetic variants in the CHRNA5/A3/B4 region, especially CHRNA3, and lung adenocarcinoma. This may provide valuable information for understanding the pathogenesis of lung adenocarcinoma and for conducting personalized smoking cessation interventions.     

Contribution of Common PCSK1 Genetic Variants to Obesity in 8,359 Subjects from Multi-Ethnic American Population

Common PCSK1 variants (notably rs6232 and rs6235) have been shown to be associated with obesity in European, Asian and Mexican populations. To determine whether common PCSK1 variants contribute to obesity in American population, we conducted association analyses in 8,359 subjects using two multi-ethnic American studies: the Coronary Artery Risk Development in Young Adults (CARDIA) study and the Multi-Ethnic Study of Atherosclerosis (MESA). By evaluating the contribution of rs6232 and rs6235 in each ethnic group, we found that in European-American subjects from CARDIA, only rs6232 was associated with BMI (P = 0.006) and obesity (P = 0.018) but also increased the obesity incidence during the 20 years of follow-up (HR = 1.53 [1.07–2.19], P = 0.019). Alternatively, in African-American subjects from CARDIA, rs6235 was associated with BMI (P = 0.028) and obesity (P = 0.018). Further, by combining the two case-control ethnic groups from the CARDIA study in a meta-analysis, association between rs6235 and obesity risk remained significant (OR = 1.23 [1.05–1.45], P = 9.5×10⁻³). However, neither rs6232 nor rs6235 was associated with BMI or obesity in the MESA study. Interestingly, rs6232 was associated with BMI (P = 4.2×10⁻³) and obesity (P = 3.4×10⁻³) in the younger European-American group combining samples from the both studies [less than median age (53 years)], but not among the older age group (P = 0.756 and P = 0.935 for BMI and obesity, respectively). By combining all the case-control ethnic groups from CARDIA and MESA in a meta-analysis, we found no significant association for the both variants and obesity risk. Finally, by exploring the full PCSK1 locus, we observed that no variant remained significant after correction for multiple testing. These results indicate that common PCSK1 variants (notably rs6232 and rs6235) contribute modestly to obesity in multi-ethnic American population. Further, these results suggest that the association of rs6232 with obesity may be age-dependent in European-Americans. However, multiple replication studies in multi-ethnic American population are needed to confirm our findings.     

Non-Alcoholic Fatty Liver Disease Is Closely Associated with Sub-Clinical Inflammation: A Case-Control Study on Asian Indians in North India

Objectives

Association between sub-clinical inflammation and non-alcoholic fatty liver disease (NAFLD) has not been studied in Asian Indians. In this case-control study, we aimed to analyse association of NAFLD with the sub-clinical inflammation and metabolic profile in Asian Indians in north India.

Methods

Ultrasound diagnosed 120 cases of NAFLD were compared to 152 healthy controls without NAFLD. Anthropometric profile [body mass index (BMI), waist circumference (WC), hip circumference (HC)], high-sensitivity C-reactive protein (hs-CRP), metabolic profile [fasting blood glucose (FBG), lipid profile] and hepatic function tests [alanine aminotransferase (ALT) and aspartate aminotransferase (AST)] were recorded.

Results

Metabolic parameters [FBG, total cholesterol (TC), serum triglycerides (TG),low-density lipoprotein (LDL-c)], hs-CRP and prevalence of the metabolic syndrome were higher in cases as compared to controls (p-value<0.05 for all). The median (range) of hs-CRP (mg/L) for cases [2.6(0.2–13.4)] were significantly higher than in controls [1.4(0.03–11.4), p = 0.01]. Similarly, higher values of hs-CRP were obtained when subgroups of cases with obesity, abdominal obesity and the metabolic syndrome were compared to controls [2.75 (0.03–14.3) vs. 1.52 (0.04–14.3), p = 0.0010; 2.8 (0.03–14.3) vs. 1.5 (0.06–14.3), p = 0.0014 and 2.7 (0.5–14.3) vs. 1.6 (0.06–8.5), p = 0.0013, respectively. On multivariate logistic regression analysis BMI (p = 0.001), WC (p = 0.001), FBG (p = 0.002), TC (p = 0.008), TG (p = 0.002), blood pressure (p = 0.005), metabolic syndrome (p = 0.001) and hs-CRP (p = 0.003) were significantly and independently associated with NAFLD. After adjusting for significant variables, the association between high hs-CRP and NAFLD remained large and statistically significant [adjusted OR = 1.17, 95% confidence interval (CI) = 1.05–1.29]. An increase in 1 mg/dl of hs-CRP level calculated to increase the risk of developing NAFLD by 1.7 times as compared to controls after adjusting for significant variables associated with NAFLD.

Conclusions

In this cohort of Asian Indians in North India, presence of NAFLD showed independent relationships with sub-clinical inflammation.     

Association of LIN28B with Adult Adiposity-Related Traits in Females

Context

Pubertal timing is under strong genetic control and its early onset associates with several adverse health outcomes in adulthood, including obesity, type 2 diabetes and cardiovascular disease. Recent data indicate strong association between pubertal timing and genetic variants near LIN28B, but it is currently unknown whether the gene contributes to the association between puberty and adult disease.

Objective

To elucidate the putative genetic link between early puberty and adult disease risk, we examined the association of two genetic variants near LIN28B with adult body size and metabolic profiles in randomly ascertained adult Finnish males and females.

Methods

Two single nucleotide polymorphisms (SNPs), rs7759938, the lead SNP previously associated with pubertal timing and height, and rs314279, previously also associated with menarcheal age but only partially correlated with rs7759938 (r² = 0.30), were genotyped in 26,636 study subjects participating in the Finnish population survey FINRISK. Marker associations with adult height, weight, body mass index (BMI), hip and waist circumference, blood glucose, serum insulin and lipid/lipoprotein levels were determined by linear regression analyses.

Results

Both rs7759938 and rs314279 associated with adult height in both sexes (p = 2×10⁻⁶ and p = 0.001). Furthermore, rs314279 associated with increased weight in females (p = 0.001). Conditioned analyses including both SNPs in the regression model verified that rs314279 independently associates with adult female weight, BMI and hip circumference (p<0.005). Neither SNP associated with glucose, lipid, or lipoprotein levels.

Conclusion

Genetic variants near the puberty-associated gene LIN28B associate with adult weight and body shape in females, suggesting that the gene may tag molecular pathways influencing adult adiposity-related traits.     

Whole-Genome Sequencing of Individuals from a Founder Population Identifies Candidate Genes for Asthma

Asthma is a complex genetic disease caused by a combination of genetic and environmental risk factors. We sought to test classes of genetic variants largely missed by genome-wide association studies (GWAS), including copy number variants (CNVs) and low-frequency variants, by performing whole-genome sequencing (WGS) on 16 individuals from asthma-enriched and asthma-depleted families. The samples were obtained from an extended 13-generation Hutterite pedigree with reduced genetic heterogeneity due to a small founding gene pool and reduced environmental heterogeneity as a result of a communal lifestyle. We sequenced each individual to an average depth of 13-fold, generated a comprehensive catalog of genetic variants, and tested the most severe mutations for association with asthma. We identified and validated 1960 CNVs, 19 nonsense or splice-site single nucleotide variants (SNVs), and 18 insertions or deletions that were out of frame. As follow-up, we performed targeted sequencing of 16 genes in 837 cases and 540 controls of Puerto Rican ancestry and found that controls carry a significantly higher burden of mutations in IL27RA (2.0% of controls; 0.23% of cases; nominal p = 0.004; Bonferroni p = 0.21). We also genotyped 593 CNVs in 1199 Hutterite individuals. We identified a nominally significant association (p = 0.03; Odds ratio (OR) = 3.13) between a 6 kbp deletion in an intron of NEDD4L and increased risk of asthma. We genotyped this deletion in an additional 4787 non-Hutterite individuals (nominal p = 0.056; OR = 1.69). NEDD4L is expressed in bronchial epithelial cells, and conditional knockout of this gene in the lung in mice leads to severe inflammation and mucus accumulation. Our study represents one of the early instances of applying WGS to complex disease with a large environmental component and demonstrates how WGS can identify risk variants, including CNVs and low-frequency variants, largely untested in GWAS.     

CDKAL1-Related Single Nucleotide Polymorphisms Are Associated with Insulin Resistance in a Cross-Sectional Cohort of Greek Children

Five novel loci recently found to be associated with body mass in two GWAS of East Asian populations were evaluated in two cohorts of Swedish and Greek children and adolescents. These loci are located within, or in the proximity of: CDKAL1, PCSK1, GP2, PAX6 and KLF9. No association with body mass has previously been reported for these loci in GWAS performed on European populations. The single nucleotide polymorphisms (SNPs) with the strongest association at each loci in the East Asian GWAS were genotyped in two cohorts, one obesity case control cohort of Swedish children and adolescents consisting of 496 cases and 520 controls and one cross-sectional cohort of 2293 nine-to-thirteen year old Greek children and adolescents. SNPs were surveyed for association with body mass and other phenotypic traits commonly associated with obesity, including adipose tissue distribution, insulin resistance and daily caloric intake. No association with body mass was found in either cohort. However, among the Greek children, association with insulin resistance could be observed for the two CDKAL1-related SNPs: rs9356744 (β = 0.018, p = 0.014) and rs2206734 (β = 0.024, p = 0.001). CDKAL1-related variants have previously been associated with type 2 diabetes and insulin response. This study reports association of CDKAL1-related SNPs with insulin resistance, a clinical marker related to type 2 diabetes in a cross-sectional cohort of Greek children and adolescents of European descent.     

Population-Specific Haplotype Association of the Postsynaptic Density Gene DLG4 with Schizophrenia,in Family-Based Association Studies

The post-synaptic density (PSD) of glutamatergic synapses harbors a multitude of proteins critical for maintaining synaptic dynamics. Alteration of protein expression levels in this matrix is a marked phenomenon of neuropsychiatric disorders including schizophrenia, where cognitive functions are impaired. To investigate the genetic relationship of genes expressed in the PSD with schizophrenia, a family-based association analysis of genetic variants in PSD genes such as DLG4, DLG1, PICK1 and MDM2, was performed, using Japanese samples (124 pedigrees, n = 376 subjects). Results showed a significant association of the rs17203281 variant from the DLG4 gene, with preferential transmission of the C allele (p = 0.02), although significance disappeared after correction for multiple testing. Replication analysis of this variant, found no association in a Chinese schizophrenia cohort (293 pedigrees, n = 1163 subjects) or in a Japanese case-control sample (n = 4182 subjects). The DLG4 expression levels between postmortem brain samples from schizophrenia patients showed no significant changes from controls. Interestingly, a five marker haplotype in DLG4, involving rs2242449, rs17203281, rs390200, rs222853 and rs222837, was enriched in a population specific manner, where the sequences A-C-C-C-A and G-C-C-C-A accumulated in Japanese (p = 0.0009) and Chinese (p = 0.0007) schizophrenia pedigree samples, respectively. However, this could not be replicated in case-control samples. None of the variants in other examined candidate genes showed any significant association in these samples. The current study highlights a putative role for DLG4 in schizophrenia pathogenesis, evidenced by haplotype association, and warrants further dense screening for variants within these haplotypes.     

Genetic Variation in the Base Excision Repair Pathway,Environmental Risk Factors,and Colorectal Adenoma Risk

Cigarette smoking, high alcohol intake, and low dietary folate levels are risk factors for colorectal adenomas. Oxidative damage caused by these three factors can be repaired through the base excision repair pathway (BER). We hypothesized that genetic variation in BER might modify colorectal adenoma risk. In a sigmoidoscopy-based study, we examined associations between 182 haplotype tagging SNPs in 14 BER genes, and colorectal adenoma risk, and examined their potential role as modifiers of the effect cigarette smoking, alcohol intake, and dietary folate levels. Among all individuals, no statistically significant associations between BER SNPs and adenoma risk persisted after correction for multiple comparisons. However, among Asian-Pacific Islanders we observed two SNPs in FEN1 and one in NTHL1, and among African-Americans one SNP in APEX1 that were associated with colorectal adenoma risk. Significant associations were also observed between SNPs in the NEIL2 gene and rectal adenoma risk. Three SNPS modified the effect of smoking (MUTYH interaction p = 0.002; OGG1 interaction p = 0.013); FEN1 interaction p = 0.013)), one SNP in LIG3 modified the effect of alcohol consumption (interaction p = 0.024) and two SNPs in LIG3 modified the effect of dietary folate (interaction p = 0.001 and p = 0.08) on colorectal adenoma risk. These findings support a role for genetic variants in the BER pathway as potential modifiers of colorectal adenoma risk. Our findings strengthen the role of oxidative damage induced by key lifestyle and dietary risk factors in colorectal adenoma formation.     

The STK33-Linked SNP rs4929949 Is Associated with Obesity and BMI in Two Independent Cohorts of Swedish and Greek Children

Recent genome wide association studies (GWAS) have identified a locus on chromosome 11p15.5, closely associated with serine/threonine kinase 33 (STK33), to be associated with body mass. STK33, a relatively understudied protein, has been linked to KRAS mutation-driven cancers and explored as a potential antineoplastic drug target. The strongest association with body mass observed at this loci in GWAS was rs4929949, a single nucleotide polymorphism located within intron 1 of the gene encoding STK33. The functional implications of rs4929949 or related variants have not been explored as of yet. We have genotyped rs4929949 in two cohorts, an obesity case-control cohort of 991 Swedish children, and a cross-sectional cohort of 2308 Greek school children. We found that the minor allele of rs4929949 was associated with obesity in the cohort of Swedish children and adolescents (OR = 1.199 (95%CI: 1.002–1.434), p = 0.047), and with body mass in the cross-sectional cohort of Greek children (β = 0.08147 (95% CI: 0.1345–0.1618), p = 0.021). We observe the effects of rs4929949 on body mass to be detectable already at adolescence. Subsequent analysis did not detect any association of rs4929949 to phenotypic measurements describing body adiposity or to metabolic factors such as insulin levels, triglycerides and insulin resistance (HOMA).     

Association of Common Genetic Variants with Lipid Traits in the Indian Population

Genome-wide association studies (GWAS) have been instrumental in identifying novel genetic variants associated with altered plasma lipid levels. However, these quantitative trait loci have not been tested in the Indian population, where there is a poorly understood and growing burden of cardiometabolic disorders. We present the association of six single nucleotide polymorphisms in 1671 sib pairs (3342 subjects) with four lipid traits: total cholesterol, triglycerides, high density lipoprotein cholesterol (HDL-C) and low density lipoprotein cholesterol (LDL-C). We also investigated the interaction effects of gender, location, fat intake and physical activity. Each copy of the risk allele of rs964184 at APOA1 was associated with 1.06 mmol/l increase in triglycerides (SE = 0.049; p = 0.006), rs3764261 at CETP with 1.02 mmol/l increase in both total cholesterol (SE = 0.042; p = 0.017) and HDL-C (SE = 0.041; p = 0.008), rs646776 at CELSR2-PSRC1-SORT1 with 0.96 mmol/l decrease in cholesterol (SE = 0.043; p = 0.0003) and 0.15 mmol/l decrease in LDL-C levels (SE = 0.043; p = 0.0003) and rs2954029 at TRIB1 with 1.02 mmol/l increase in HDL-C (SE = 0.039; p = 0.047). A combined risk score of APOA1 and CETP loci predicted an increase of 1.25 mmol/l in HDL-C level (SE = 0.312; p = 0.0007). Urban location and sex had strong interaction effects on the genetic association of most of the studied loci with lipid traits. To conclude, we validated four genetic variants (identified by GWAS in western populations) associated with lipid traits in the Indian population. The interaction effects found here may explain the sex-specific differences in lipid levels and their heritability. Urbanization appears to influence the nature of the association with GWAS lipid loci in this population. However, these findings will require replication in other Indian populations.     

Analysis of Factors That Determine Weight Gain during Smoking Cessation Therapy

Cigarette smokers are generally known to gain weight after quitting smoking, and such weight gain is thought to contribute to the worsening of glucose tolerance. While smoking cessation therapy such as nicotine replacement is useful to minimize post-cessation weight gain, substantial gain occurs even during the therapy. The purpose of the present study was to identify factors associated with weight gain during smoking cessation therapy. We evaluated 186 patients(132 males and 54 females)who visited our outpatient clinic for smoking cessation, and successfully achieved smoking abstinence. We performed gender-adjusted regression analysis for the rate of BMI increase from the beginning of cessation to 3 months after initiation. Furthermore, we performed multivariate analysis to investigate factors that determine the BMI increase after smoking cessation. The mean BMI significantly (p<0.0001) increased from 23.5±3.6 kg/m² at the initial consultation to 23.9±3.8 kg/m² at 3 months after the start of therapy. There was no significant difference in the extent of BMI increase between nicotine patch and varenicline therapy groups. Factors significantly correlated with the %BMI increase at 3 months after the start of therapy were triglyceride (p = 0.0006, βa = 0.260), high-density lipoprotein cholesterol (p = 0.0386, βa = −0.168), daily cigarette consumption (p = 0.0385, βa = 0.154), and the Fagerström Test for Nicotine Dependence (FTND) score (p = 0.0060, βa = 0.203). Stepwise multivariate analysis demonstrated that triglyceride and the FTND score were the factors determining the post-cessation BMI increase and that the FTND score was the strongest one. The present study demonstrated that smokers with a high FTND score are more likely to gain weight during smoking cessation therapy. Thus, smokers with a high nicotine dependency may require intervention against weight gain in the cessation clinic.     

Functional Characterization Improves Associations between Rare Non-Synonymous Variants in CHRNB4 and Smoking Behavior

Smoking is the leading cause of preventable death worldwide. Accordingly, effort has been devoted to determining the genetic variants that contribute to smoking risk. Genome-wide association studies have identified several variants in nicotinic acetylcholine receptor genes that contribute to nicotine dependence risk. We previously undertook pooled sequencing of the coding regions and flanking sequence of the CHRNA5, CHRNA3, CHRNB4, CHRNA6 and CHRNB3 genes and found that rare missense variants at conserved residues in CHRNB4 are associated with reduced risk of nicotine dependence among African Americans. We identified 10 low frequency (<5%) non-synonymous variants in CHRNB4 and investigated functional effects by co-expression with normal α3 or α4 subunits in human embryonic kidney cells. Voltage-clamp was used to obtain acetylcholine and nicotine concentration–response curves and qRT-PCR, western blots and cell-surface ELISAs were performed to assess expression levels. These results were used to functionally weight genetic variants in a gene-based association test. We find that there is a highly significant correlation between carrier status weighted by either acetylcholine EC₅₀ (β = −0.67, r² = 0.017, P = 2×10⁻⁴) or by response to low nicotine (β = −0.29, r² = 0.02, P = 6×10⁻⁵) when variants are expressed with the α3 subunit. In contrast, there is no significant association when carrier status is unweighted (β = −0.04, r² = 0.0009, P = 0.54). These results highlight the value of functional analysis of variants and the advantages to integrating such data into genetic studies. They also suggest that an increased sensitivity to low concentrations of nicotine is protective from the risk of developing nicotine dependence.     

Association of the IL-10 Gene Family Locus on Chromosome 1 with Juvenile Idiopathic Arthritis (JIA)

Background

The cytokine IL-10 and its family members have been implicated in autoimmune diseases and we have previously reported that genetic variants in IL-10 were associated with a rare group of diseases called juvenile idiopathic arthritis (JIA). The aim of this study was to fine map genetic variants within the IL-10 cytokine family cluster on chromosome 1 using linkage disequilibrium (LD)-tagging single nucleotide polymorphisms (tSNPs) approach with imputation and conditional analysis to test for disease associations.

Methodology/Principal Findings

Fifty-three tSNPs were tested for association between Caucasian paediatric cohorts [219 systemic JIA (sJIA), 187 persistent oligoarticular JIA (pOJIA), and 139 extended OJIA (eOJIA) patients], and controls (Wellcome Trust control cohort, WTCCC2). Significant association with sJIA was detected at rs1400986 in the promoter of IL-20 (odds ratio 1.53; 95% CI 1.21–1.93; p = 0.0004), but in no other subtypes. Imputation analysis identified additional associated SNPs for pOJIA at IL-20 and IL-24, including a rare, functional, missense variant at IL-24 with a p = 0.0002. Penalised logistic regression analysis with HyperLasso and conditional analysis identified several further associations with JIA subtypes. In particular, haplotype analysis refined the sJIA association, with a joint effect at rs1400986 and rs4129024 in intron 1 of MAPKAPK2 (p = 3.2E−5). For pOJIA, a 3-SNP haplotype including rs1878672 in intron 3 of IL-10 showed evidence for association (p = 0.0018). In eOJIA, rs10863962 (3′UTR of FCAMR) and rs12409577 (intron of IL-19) haplotype showed some evidence of association (p = 0.0003).

Conclusions

This study supports previous association of IL-20 with sJIA. Haplotype analyses provided stronger association signals than single point analyses, while a penalised logistic regression approach also suggested multiple independent association signals. Replication studies are required to confirm or refute these findings. The results indicate that combined effects with unknown/rare variants remain to be characterised in JIA, and represent a possible example of synthetic association in this region.     

Association between Milk and Milk Product Consumption and Anthropometric Measures in Adult Men and Women in India: A Cross-Sectional Study

Background

The nutritional aetiology of obesity remains unclear, especially with regard to the role of dairy products in developing countries.

Objective

To examine whether milk/milk product consumption is associated with obesity and high waist circumference among adult Indians.

Methods

Information on plain milk, tea, curd and buttermilk/lassi consumption assessed using a Food Frequency Questionnaire was obtained from the cross-sectional sib-pair designed Indian Migration Study (3698 men and 2659 women), conducted at four factory locations across north, central and south India. The anthropometric measures included were Body Mass Index (BMI) and Waist Circumference (WC). Mixed-effect logistic regression models were conducted to accommodate sib-pair design and adjust for potential confounders.

Results

After controlling for potential confounders, the risk of being obese (BMI≥25 kg/m²) was lower among women (OR = 0.57;95%CI:0.43−0.76;p≤0.0001) and men (OR = 0.67;95%CI: 0.51−0.87;p = 0.005), and the risk of a high WC (men: >90 cm; women: >80 cm) was lower among men (OR = 0.71;95%CI:0.54−0.93;p = 0.005) and women (OR = 0.79;95%CI:0.59−1.05;p>0.05) who consume ≥1 portions of plain milk daily than those who do not consume any milk. The inverse association between daily plain milk consumption and obesity was also confirmed in sibling-pair analyses. Daily tea consumption of ≥1 portion was associated with obesity (OR = 1.51;95%CI:1.00−2.25;p>0.050) and high WC (OR = 1.65;95%CI:1.08−2.51;p>0.019) among men but not among women but there was no strong evidence of association of curd and buttermilk/lassi consumption with obesity and high waist circumference among both men and women.

Conclusions

The independent, inverse association of daily plain milk consumption with the risk of being obese suggests that high plain milk intake may lower the risk of obesity in adult Indians. However, this is an observational finding and uncontrolled confounding cannot be excluded as an explanation for the association. Therefore, confirmatory studies are needed to clarify this relationship.     

Genetic Variation in DNA Repair Pathways and Risk of Non-Hodgkin's Lymphoma

Molecular and genetic evidence suggests that DNA repair pathways may contribute to lymphoma susceptibility. Several studies have examined the association of DNA repair genes with lymphoma risk, but the findings from these reports have been inconsistent. Here we provide the results of a focused analysis of genetic variation in DNA repair genes and their association with the risk of non-Hodgkin''s lymphoma (NHL). With a population of 1,297 NHL cases and 1,946 controls, we have performed a two-stage case/control association analysis of 446 single nucleotide polymorphisms (SNPs) tagging the genetic variation in 81 DNA repair genes. We found the most significant association with NHL risk in the ATM locus for rs227060 (OR = 1.27, 95% CI: 1.13–1.43, p = 6.77×10⁻⁵), which remained significant after adjustment for multiple testing. In a subtype-specific analysis, associations were also observed for the ATM locus among both diffuse large B-cell lymphomas (DLBCL) and small lymphocytic lymphomas (SLL), however there was no association observed among follicular lymphomas (FL). In addition, our study provides suggestive evidence of an interaction between SNPs in MRE11A and NBS1 associated with NHL risk (OR = 0.51, 95% CI: 0.34–0.77, p = 0.0002). Finally, an imputation analysis using the 1,000 Genomes Project data combined with a functional prediction analysis revealed the presence of biologically relevant variants that correlate with the observed association signals. While the findings generated here warrant independent validation, the results of our large study suggest that ATM may be a novel locus associated with the risk of multiple subtypes of NHL.     

Investigation of Cigarette Smoking among Male Schizophrenia Patients

Male schizophrenia patients are known to have a heavier smoking pattern compared with the general population. However, the mechanism for this association is not known, though hypothesis that smoking could alleviate symptomatology of schizophrenia and reduce side effects of antipsychotics has been suggested. The aims of this study were to validate the heavier smoking pattern among male schizophrenia patients and to investigate the possible mechanisms for the association. To enhance the reliability of the study, we recruited two large independent samples with 604 and 535 male Chinese schizophrenia patients, and compared their smoking pattern with that of 535 healthy male controls recruited from general population. Validated multiple indicators and multiple causes structure equation model and regression models were used to investigate the association of smoking with factors of schizophrenia symptomatology and with the usage of antipsychotics and their extra-pyramidal side effects (EPS). Schizophrenia patients had significantly heavier smoking pattern compared with healthy controls in our sample (42.4% vs. 16.8%, p<0.001 for current smoking prevalence; 23.5% vs. 43.3%, p<0.001 for smoking cessation rate; 24.5% vs. 3.0%, p<0.001 for heavy smoker proportion). Their smoking status was also found to be consistently and significantly associated with reduced negative factor scores for schizophrenia symptomatology (β = −0.123, p = 0.051 for sample-A; β = −0.103, p = 0.035 for sample-B; β = −0.082, p = 0.017 for the combined sample). However, no significant association was found between smoking and antipsychotics usage or risk of EPS. These results support that smoking is associated with improved negative symptoms, which could account for the heavier smoking pattern among schizophrenia patients.     

REV-ERB ALPHA Polymorphism Is Associated with Obesity in the Spanish Obese Male Population

REV-ERB ALPHA has been shown to link metabolism with circadian rhythms. We aimed to identify new polymorphisms in the promoter of REV-ERB ALPHA and tested whether these polymorphisms could be associated with obesity in the Spanish population. Of the 1197 subjects included in our study, 779 were obese (BMI 34.38±3.1 kg/m²) and 418 lean (BMI 23.27±1.5 kg/m²). In the obese group, 469 of the 779 had type 2 diabetes. Genomic DNA from all the subjects was obtained from peripheral blood cells and the genotyping in the REV-ERB ALPHA promoter was analyzed by High Resolution Melting. We found six polymorphisms in the REV-ERB ALPHA promoter and identified rs939347 as a SNP with the highest frequency in the total population. We did not find any association between rs939347 and type 2 diabetes (p = 0.101), but rs939347 was associated with obesity (p = 0.036) with the genotype AA exhibiting higher frequency in the obese (5.2% in total obese vs 2.4% in lean). This association was found only in men (p = 0.031; 6.5% AA-carriers in obese men vs 1.9% AA-carriers in lean men), with no association found in the female population (p = 0.505; 4.4% AA-carriers in obese women vs 2.7% AA-carriers in lean women). Our results suggest that the REV-ERB ALPHA rs939347 polymorphism could modulate body fat mass in men. The present work supports the role of REV-ERB ALPHA in the development of obesity as well as a potential target for the treatment of obesity.     

Prospective Validation of Candidate SNPs of VEGF/VEGFR Pathway in Metastatic Colorectal Cancer Patients Treated with First-Line FOLFIRI Plus Bevacizumab

Purpose

The potential impact of different SNPs of VEGF/VEGFR pathway on the clinical outcome of mCRC patients receiving bev-containing regimens has been investigated in retrospective experiences with contrasting results. We previously reported the association of VEGFA rs833061 C/T variants with PFS in metastatic colorectal cancer patients treated with first-line FOLFIRI plus bevacizumab. The primary objective of this work was to prospectively validate that retrospective finding. A confirmatory analysis of other SNPs of VEGF/VEGFR pathway genes was included.

Experimental design

To detect a HR for PFS of 1.7 for VEGFA rs833061 T/T compared to C- variants in metastatic colorectal cancer patients treated with first-line FOLFIRI plus bevacizumab, setting two-sided α = 0.05 and β = 0.20, 199 events were required. VEGFA rs699946 A/G, rs699947 A/C, VEGFR1 rs9582036 A/C and rs7993418 A/G, VEGFR2 rs11133360 C/T, rs12505758 C/T and rs2305948 C/T and EPAS1 rs4145836 A/G were also tested. Germ-line DNA was extracted from peripheral blood. SNPs were analyzed by PCR and sequencing.

Results

Four-hundred-twenty-four pts were included. At the univariate analysis, no differences according to VEGFA rs833061 C/T variants were observed in PFS (p = 0.38) or OS (p = 0.95). Among analyzed SNPs, only VEGFR2 rs12505758 C- variants, compared to T/T, were associated to shorter PFS (HR: 1.36 [1.05–1.75], p = 0.015, dominant genetic model) and OS, with a trend toward significance (HR: 1.34 [0.95–1.88], p = 0.088). In the multivariate model, this association retained significance (HR: 1.405 [1.082–1.825], p = 0.012) in PFS, that was lost by applying multiple testing correction (p = 0.14).

Conclusion

This prospective experience failed to validate the hypothesized predictive impact of VEGFA rs833061 variants. Retrospective findings on different candidate SNPs were not confirmed. Only VEGFR2 rs12505758 variants, whose prognostic and not predictive impact was previously reported, correlated with PFS. Given the complexity of angiogenesis, it is rather unlike that a single germ-line SNP might be a good predictor of benefit from bevacizumab.