Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

Genetic Variation in the 3'-Untranslated Region of NBN Gene Is Associated with Gastric Cancer Risk in a Chinese Population

NBN plays a crucial role in carcinogenesis as a core component for both homologous recombination (HR) and non-homologous end-joining (NHEJ) DNA double-strand breaks (DSBs) repair pathways. Genetic variants in the NBN gene have been associated with multiple cancers risk, suggesting pleiotropic effect on cancer. We hypothesized that genetic variants in the NBN gene may modify the risk of gastric cancer. To test this hypothesis, we evaluated the association between four potentially functional single nucleotide polymorphisms in NBN and gastric cancer risk in a case–control study of 1,140 gastric cancer cases and 1,547 controls in a Chinese population. We found that the A allele of rs10464867 (G>A) was significantly associated with a decreased risk of gastric cancer (odds ratio [OR] = 0.81, 95% confidence interval [95% CI] = 0.71–0.94; P = 4.71×10⁻³). Furthermore, the association between A allele of rs10464867 and decreased risk of gastric cancer was more significantly in elder individuals (per-allele OR = 0.72[0.59–0.88], P = 1.07×10⁻³), and male individuals (per-allele OR = 0.73[0.62–0.87], P = 3.68×10⁻⁴). We further conducted a haplotype analysis and identified that the NBN A_rs10464867G_rs14448G_rs1063053 haplotype conferred stronger protective effect on gastric cancer (OR = 0.76[0.65–0.89], P = 6.39×10⁻⁴). In summary, these findings indicate that genetic variants at NBN gene may contribute to gastric cancer susceptibility and may further advance our understanding of NBN gene in cancer development.     

Low-Frequency Coding Variants at 6p21.33 and 20q11.21 Are Associated with Lung Cancer Risk in Chinese Populations

Genome-wide association studies have successfully identified a subset of common variants associated with lung cancer risk. However, these variants explain only a fraction of lung cancer heritability. It has been proposed that low-frequency or rare variants might have strong effects and contribute to the missing heritability. To assess the role of low-frequency or rare variants in lung cancer development, we analyzed exome chips representing 1,348 lung cancer subjects and 1,998 control subjects during the discovery stage and subsequently evaluated promising associations in an additional 4,699 affected subjects and 4,915 control subjects during the replication stages. Single-variant and gene-based analyses were carried out for coding variants with a minor allele frequency less than 0.05. We identified three low-frequency missense variants in BAT2 (rs9469031, c.1544C>T [p.Pro515Leu]; odds ratio [OR] = 0.55, p = 1.28 × 10⁻¹⁰), FKBPL (rs200847762, c.410C>T [p.Pro137Leu]; OR = 0.25, p = 9.79 × 10⁻¹²), and BPIFB1 (rs6141383, c.850G>A [p.Val284Met]; OR = 1.72, p = 1.79 × 10⁻⁷); these variants were associated with lung cancer risk. rs9469031 in BAT2 and rs6141383 in BPIFB1 were also associated with the age of onset of lung cancer (p = 0.001 and 0.006, respectively). BAT2 and FKBPL at 6p21.33 and BPIFB1 at 20q11.21 were differentially expressed in lung tumors and paired normal tissues. Gene-based analysis revealed that FKBPL, in which two independent variants were identified, might account for the association with lung cancer risk at 6p21.33. Our results highlight the important role low-frequency variants play in lung cancer susceptibility and indicate that candidate genes at 6p21.33 and 20q11.21 are potentially biologically relevant to lung carcinogenesis.     

Polymorphisms in a Putative Enhancer at the 10q21.2 Breast Cancer Risk Locus Regulate NRBF2 Expression

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82–0.88]) and ER-negative (OR = 0.87 [0.82–0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91–0.95] and OR = 1.06 [1.03–1.09]) and ER-negative (OR = 0.95 [0.91–0.98] and OR = 1.08 [1.04–1.13]) disease. There was weaker evidence for iCHAV4, located 5′ of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90–0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1–4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.     

Identification of Multiple Genetic Susceptibility Loci in Takayasu Arteritis

Takayasu arteritis is a rare inflammatory disease of large arteries. The etiology of Takayasu arteritis remains poorly understood, but genetic contribution to the disease pathogenesis is supported by the genetic association with HLA-B^∗52. We genotyped ∼200,000 genetic variants in two ethnically divergent Takayasu arteritis cohorts from Turkey and North America by using a custom-designed genotyping platform (Immunochip). Additional genetic variants and the classical HLA alleles were imputed and analyzed. We identified and confirmed two independent susceptibility loci within the HLA region (r² < 0.2): HLA-B/MICA (rs12524487, OR = 3.29, p = 5.57 × 10⁻¹⁶) and HLA-DQB1/HLA-DRB1 (rs113452171, OR = 2.34, p = 3.74 × 10⁻⁹; and rs189754752, OR = 2.47, p = 4.22 × 10⁻⁹). In addition, we identified and confirmed a genetic association between Takayasu arteritis and the FCGR2A/FCGR3A locus on chromosome 1 (rs10919543, OR = 1.81, p = 5.89 × 10⁻¹²). The risk allele in this locus results in increased mRNA expression of FCGR2A. We also established the genetic association between IL12B and Takayasu arteritis (rs56167332, OR = 1.54, p = 2.18 × 10⁻⁸).     

Rare,Low-Frequency,and Common Variants in the Protein-Coding Sequence of Biological Candidate Genes from GWASs Contribute to Risk of Rheumatoid Arthritis

The extent to which variants in the protein-coding sequence of genes contribute to risk of rheumatoid arthritis (RA) is unknown. In this study, we addressed this issue by deep exon sequencing and large-scale genotyping of 25 biological candidate genes located within RA risk loci discovered by genome-wide association studies (GWASs). First, we assessed the contribution of rare coding variants in the 25 genes to the risk of RA in a pooled sequencing study of 500 RA cases and 650 controls of European ancestry. We observed an accumulation of rare nonsynonymous variants exclusive to RA cases in IL2RA and IL2RB (burden test: p = 0.007 and p = 0.018, respectively). Next, we assessed the aggregate contribution of low-frequency and common coding variants to the risk of RA by dense genotyping of the 25 gene loci in 10,609 RA cases and 35,605 controls. We observed a strong enrichment of coding variants with a nominal signal of association with RA (p < 0.05) after adjusting for the best signal of association at the loci (p_enrichment = 6.4 × 10⁻⁴). For one locus containing CD2, we found that a missense variant, rs699738 (c.798C>A [p.His266Gln]), and a noncoding variant, rs624988, reside on distinct haplotypes and independently contribute to the risk of RA (p = 4.6 × 10⁻⁶). Overall, our results indicate that variants (distributed across the allele-frequency spectrum) within the protein-coding portion of a subset of biological candidate genes identified by GWASs contribute to the risk of RA. Further, we have demonstrated that very large sample sizes will be required for comprehensively identifying the independent alleles contributing to the missing heritability of RA.     

A genome-wide association study reveals that variants within the HLA region are associated with risk for nonobstructive azoospermia

A genome-wide association study of Han Chinese subjects was conducted to identify genetic susceptibility loci for nonobstructive azoospermia (NOA). In the discovery stage, 802 azoospermia cases and 1,863 controls were screened for genetic variants in the genome. Promising SNPs were subsequently confirmed in two independent sets of subjects: 818 azoospermia cases and 1,755 controls from northern China, and 606 azoospermia cases and 958 controls from central and southern China. We detected variants at human leukocyte antigen (HLA) regions that were independently associated with NOA (HLA-DRA, rs3129878, p_combine = 3.70 × 10⁻¹⁶, odds ratio [OR] = 1.37; C6orf10 and BTNL2, rs498422, p_combine = 2.43 × 10⁻¹², OR = 1.42). These findings provide additional insight into the pathogenesis of NOA.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

Genetic Variants at 6p21.1 and 7p15.3 Are Associated with Risk of Multiple Cancers in Han Chinese

Cancer susceptibility loci identified in reported genome-wide association studies (GWAS) are often tumor-specific; however, evidence of pleiotropy of some genes/loci has also been observed and biologically plausible. We hypothesized that there are important regions in the genome harboring genetic variants associated with risk of multiple types of cancer. In the current study, we attempted to map genetic variants that have consistent effects on risk of multiple cancers using our existing genome-wide scan data of lung cancer, noncardia gastric cancer, and esophageal squamous-cell carcinoma with overall 5,368 cases and 4,006 controls (GWAS stage), followed by a further evaluation in additional 9,001 cases with one of these cancer types and 11,436 controls (replication stage). Five variants satisfying the criteria of pleiotropy with p values from 1.10 × 10⁻⁸ to 8.96 × 10⁻⁶ for genome-wide scans of three cancer types were further evaluated in the replication stage. We found consistent associations of rs2494938 at 6p21.1 and rs2285947 at 7p15.3 with these three cancers in both GWAS and replication stages. In combined samples of GWAS and replication stages, the minor alleles of rs2494938 and rs2285947 were significantly associated with an increased risk of the cancers (odds ratio [OR] = 1.15, 95% confidence interval [CI], 1.10–1.19 and OR = 1.17, 95% CI, 1.12–1.21), with the p values being 1.20 × 10⁻¹² and 1.26 × 10⁻¹⁶, respectively, which are at a genome-wide significance level. Our findings highlight the potential importance of variants at 6p21.1 and 7p15.3 in the susceptibility to multiple cancers.     

A common variant highly associated with plasma VEGFA levels also contributes to the variation of both LDL-C and HDL-C

Vascular endothelial growth factor A (VEGFA) is among the most-significant stimulators of angiogenesis. Its effect on cardiovascular diseases and on the variation of related risk factors such as lipid parameters is considered important, although as yet unclear. Recently, we identified four common variants (rs6921438, rs4416670, rs6993770, and rs10738760) that explain up to 50% of the heritability of plasma VEGFA levels. In the present study, we aimed at assessing the contribution of these variants to the variation of blood lipid levels (including apoE, triglycerides, total cholesterol, low- and high-density lipoprotein cholesterol levels (LDL-C and HDL-C)] in healthy subjects. The effect of these single-nucleotide polymorphisms (SNPs) on lipid levels was assessed using linear regression in discovery and replication samples (n = 1,006 and n = 1,145; respectively), followed by a meta-analysis. Their gene×gene and gene×environment interactions were also assessed. SNP rs6921438 was associated with HDL-C (β = −0.08 mmol/l, P_overall = 1.2 × 10⁻⁷) and LDL-C (β = 0.13 mmol/l, P_overall = 1.5 × 10⁻⁴). We also identified a significant association between the interaction rs4416670×hypertension and apoE variation (P_overall = 1.7 × 10⁻⁵). Therefore, our present study shows a common genetic regulation between VEGFA and cholesterol homeostasis molecules. The SNP rs6921438 is in linkage disequilibrium with variants located in an enhancer- and promoter-associated histone mark region and could have a regulatory effect in the expression of surrounding genes, including VEGFA.     

Genome-wide Association Study Identifies TNFSF15 and POU2AF1 as Susceptibility Loci for Primary Biliary Cirrhosis in the Japanese Population

For the identification of susceptibility loci for primary biliary cirrhosis (PBC), a genome-wide association study (GWAS) was performed in 963 Japanese individuals (487 PBC cases and 476 healthy controls) and in a subsequent replication study that included 1,402 other Japanese individuals (787 cases and 615 controls). In addition to the most significant susceptibility region, human leukocyte antigen (HLA), we identified two significant susceptibility loci, TNFSF15 (rs4979462) and POU2AF1 (rs4938534) (combined odds ratio [OR] = 1.56, p = 2.84 × 10⁻¹⁴ for rs4979462, and combined OR = 1.39, p = 2.38 × 10⁻⁸ for rs4938534). Among 21 non-HLA susceptibility loci for PBC identified in GWASs of individuals of European descent, three loci (IL7R, IKZF3, and CD80) showed significant associations (combined p = 3.66 × 10⁻⁸, 3.66 × 10⁻⁹, and 3.04 × 10⁻⁹, respectively) and STAT4 and NFKB1 loci showed suggestive association with PBC (combined p = 1.11 × 10⁻⁶ and 1.42 × 10⁻⁷, respectively) in the Japanese population. These observations indicated the existence of ethnic differences in genetic susceptibility loci to PBC and the importance of TNF signaling and B cell differentiation for the development of PBC in individuals of European descent and Japanese individuals.     

Association of Aminoacyl-tRNA Synthetases Gene Polymorphisms with the Risk of Congenital Heart Disease in the Chinese Han Population

Aminoacyl-tRNA synthetases (ARSs) are in charge of cellular protein synthesis and have additional domains that function in a versatile manner beyond translation. Eight core ARSs (EPRS, MRS, QRS, RRS, IRS, LRS, KRS, DRS) combined with three nonenzymatic components form a complex known as multisynthetase complex (MSC).We hypothesize that the single-nucleotide polymorphisms (SNPs) of the eight core ARS coding genes might influence the susceptibility of sporadic congenital heart disease (CHD). Thus, we conducted a case-control study of 984 CHD cases and 2953 non-CHD controls in the Chinese Han population to evaluate the associations of 16 potentially functional SNPs within the eight ARS coding genes with the risk of CHD. We observed significant associations with the risk of CHD for rs1061248 [G/A; odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.81–0.99; P = 3.81×10⁻²], rs2230301 [A/C; OR = 0.73, 95%CI = 0.60–0.90, P = 3.81×10⁻²], rs1061160 [G/A; OR = 1.18, 95%CI = 1.06–1.31; P = 3.53×10⁻³] and rs5030754 [G/A; OR = 1.39, 95%CI = 1.11–1.75; P = 4.47×10⁻³] of EPRS gene. After multiple comparisons, rs1061248 conferred no predisposition to CHD. Additionally, a combined analysis showed a significant dosage-response effect of CHD risk among individuals carrying the different number of risk alleles (P _trend = 5.00×10⁻⁴). Compared with individuals with “0–2” risk allele, those carrying “3”, “4” or “5 or more” risk alleles had a 0.97-, 1.25- or 1.38-fold increased risk of CHD, respectively. These findings indicate that genetic variants of the EPRS gene may influence the individual susceptibility to CHD in the Chinese Han population.     

Polymorphism in miR-31 and miR-584 binding site in the angiotensinogen gene differentially influences body fat distribution in both sexes

Angiotensinogen (AGT), its active fragments and microRNA-31 (miR-31) play an important role in adipocyte differentiation. AGT contains a miR-31 polymorphic binding site. We hypothesize that the rs7079 polymorphism in the miR-31/584 binding site of the AGT gene could influence body fat distribution. A total of 751 subjects (195 men, 556 women) were enrolled in the study. The rs7079 genotypes were determined by qRT-PCR. Anthropometric measurements were taken on all subjects, who were subsequently divided into two groups: obese (>30 kg m⁻²) and non-obese (<30 kg m⁻²). Linear regression models were created to determine the contributions of sex, obesity status and rs7079 to all measured parameters. Adding the rs7079 genotype significantly contributed to the linear regression model for waist circumference (p = 0.013), hip circumference (p = 0.018) and supraspinal skin-fold thickness (p = 1 × 10⁻³). Differences between sexes and between the obese and non-obese groups were observed. Waist circumference was lower in men carrying the A allele (p = 0.022); hip circumference was higher only in obese women carrying the A allele (p = 0.015). While men carrying the A allele had lower supraspinal skin-fold thickness (p = 0.022), this parameter was found to be higher in A allele carrying women (p = 3 × 10⁻³). The higher total sum of skin-fold thickness in A allele carrying women was restricted to obese individuals (p = 0.028). The presence of the A allele was associated with both lower tricipital skin-fold thickness in non-obese women (p = 0.023) and a trend of higher thickness in non-obese men (p = 0.065). Significant associations of rs7079 in the AGT gene and body fat distribution were observed. The distribution followed opposing patterns in both sexes.     

Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p_meta-Euro = 2.08 × 10⁻¹⁰), transmembrane protein 39A (TMEM39A; rs1132200; p_meta-all = 8.62 × 10⁻⁹), and 17q21 (rs1453560; p_meta-all = 3.48 × 10⁻¹⁰) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10⁻⁸ < p_meta-Euro < 9.99 × 10⁻⁵) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.     

Inherited genetic variants associated with glucocorticoid sensitivity in leukaemia cells

Identification of genetic variants associated with glucocorticoids (GC) sensitivity of leukaemia cells may provide insight into potential drug targets and tailored therapy. In the present study, within 72 leukaemic cell lines derived from Japanese patients with B‐cell precursor acute lymphoblastic leukaemia (ALL), we conducted genome‐wide genotyping of single nucleotide polymorphisms (SNP) and attempted to identify genetic variants associated with GC sensitivity and NR3C1 (GC receptor) gene expression. IC50 measures for prednisolone (Pred) and dexamethasone (Dex) were available using an alamarBlue cell viability assay. IC50 values of Pred showed the strongest association with rs904419 (P = 4.34 × 10⁻⁸), located between the FRMD4B and MITF genes. The median IC50 values of prednisolone for cell lines with rs904419 AA (n = 13), AG (n = 31) and GG (n = 28) genotypes were 0.089, 0.139 and 297 µmol/L, respectively. For dexamethasone sensitivity, suggestive association was observed for SNP rs2306888 (P = 1.43 × 10⁻⁶), a synonymous SNP of the TGFBR3 gene. For NR3C1 gene expression, suggestive association was observed for SNP rs11982167 (P = 6.44 × 10⁻⁸), located in the PLEKHA8 gene. These genetic variants may affect GC sensitivity of ALL cells and may give rise to opportunities in personalized medicine for effective and safe chemotherapy in ALL patients.     

Common Variants of FTO Are Associated with Childhood Obesity in a Cross-Sectional Study of 3,126 Urban Indian Children

Background

FTO variants are robustly associated with obesity and related traits in many population and shown to have variable impact during life course. Although studies have shown association of FTO variants with adiposity in adult Indian, its association in Indian children is yet to be confirmed.

Methods

Here we examined association of FTO variants (rs9939609 and rs8050136) with obesity and related anthropometric and biochemical traits in 3,126 Indian children (aged 11–17 years) including 2,230 normal-weight and 896 over-weight/obese children. We also compared effects observed in the present study with that observed in previous studies on South Asian adults and children of other ethnic groups.

Results

The variant rs9939609 showed significant association with risk of obesity [OR = 1.21, P = 2.5×10⁻³] and its measures BMI, weight, waist circumference and hip circumference [β range = 0.11 to 0.14 Z-score units; P range = 1.3×10⁻⁴ to 1.6×10⁻⁷] in children. The observed effect sizes in Indian children were similar to those reported for European children. Variant rs9939609 explained 0.88% of BMI variance in Indian children. The effect sizes of rs9939609 on BMI and WC were ∼2 fold higher in children than adults. Interestingly rs9939609 was also associated with serum levels of thyroid stimulating hormone (TSH) [β = 0.10 Z-score, P = 5.8×10⁻³]. The other variant rs8050136 was in strong linkage disequilibrium with rs9939609 (r² = 0.97) and provided similar association results.

Conclusion

The study provides first report of association of FTO variants with obesity and related anthropometric traits in Indian children with higher impact in children compared to adults. We also demonstrated association of FTO variant with serum levels of TSH, indicating putative influence of FTO in hypothalamic-pituitary-thyroid axis.     

Systematic Evaluation of Genetic Variants for Polycystic Ovary Syndrome in a Chinese Population

To date, eleven genome-wide significant (GWS) loci (P < 5×10⁻⁸) for polycystic ovary syndrome (PCOS) have been identified through genome-wide association studies (GWAS). Some of the risk loci have been selected for replications and validated in multiple ethnicities, however, few previous studies investigated all loci. Scanning all the GWAS variants would demonstrate a more informative profile of variance they explained. Thus, we analyzed all the 17 single nucleotide polymorphisms (SNPs) mapping to the 11 GWAS loci in an independent sample set of 800 Chinese subjects with PCOS and 1110 healthy controls systematically. Variants of rs3802457 in C9orf3 locus (P = 5.99×10⁻⁴) and rs13405728 in LHCGR locus (P = 3.73×10⁻⁴) were significantly associated with PCOS after the strict Bonferroni correction in our data set. The further haplotype analysis indicated that in the block of C9orf3 gene (rs4385527 and rs3802457), GA haplotype played a protective role in PCOS (8.7 vs 5.0, P = 9.85×10⁻⁶, OR = 0.548, 95%CI = 0.418–0.717), while GG haplotype was found suffering from an extraordinarily increased risk of PCOS (73.6% vs79.2%, P = 3.41×10⁻⁵, OR = 1.394, 95%CI = 1.191–1.632). Moreover, the directions of effects for all SNPs were consistent with previous GWAS reports (P = 1.53×10⁻⁵). Polygenic score analysis demonstrated that these 17 SNPs have a significant capacity on predicting case-control status in our samples (P = 7.17×10⁻⁹), meanwhile all these gathered 17 SNPs explained about 2.40% of variance. Our findings supported that C9orf3 and LHCGR loci variants were vital susceptibility of PCOS.     

Real-life use of vitamin D3-fortified bread and milk during a winter season: the effects of CYP2R1 and GC genes on 25-hydroxyvitamin D concentrations in Danish families,the VitmaD study

Common genetic variants rs10741657 and rs10766197 in CYP2R1 and rs4588 and rs842999 in GC and a combined genetic risk score (GRS) of these four variants influence late summer 25-hydroxyvitamin D (25(OH)D) concentrations. The objectives were to identify those who are most at risk of developing low vitamin D status during winter and to assess whether vitamin D₃-fortified bread and milk will increase 25(OH)D concentrations in those with genetically determined low 25(OH)D concentrations at late summer. We used data from the VitmaD study. Participants were allocated to either vitamin D₃-fortified bread and milk or non-fortified bread and milk during winter. In the fortification group, CYP2R1 (rs10741657) and GC (rs4588 and rs842999) were statistically significantly associated with winter 25(OH)D concentrations and CYP2R1 (rs10766197) was borderline significant. There was a negative linear trend between 25(OH)D concentrations and carriage of 0–8 risk alleles (p < 0.0001). No association was found for the control group (p = 0.1428). There was a significant positive linear relationship between different quintiles of total vitamin D intake and the increase in 25(OH)D concentrations among carriers of 0–2 (p = 0.0012), 3 (p = 0.0001), 4 (p = 0.0118) or 5 (p = 0.0029) risk alleles, but not among carriers of 6–8 risk alleles (p = 0.1051). Carriers of a high GRS were more prone to be vitamin D deficient compared to carriers of a low GRS. Furthermore, rs4588-AA carriers have a low but very stable 25(OH)D concentration, and interestingly, also low PTH level.

Electronic supplementary material

The online version of this article (doi:10.1007/s12263-014-0413-7) contains supplementary material, which is available to authorized users.     

Genetic loci associated with circulating levels of very long-chain saturated fatty acids

Very long-chain saturated fatty acids (VLSFAs) are saturated fatty acids with 20 or more carbons. In contrast to the more abundant saturated fatty acids, such as palmitic acid, there is growing evidence that circulating VLSFAs may have beneficial biological properties. Whether genetic factors influence circulating levels of VLSFAs is not known. We investigated the association of common genetic variation with plasma phospholipid/erythrocyte levels of three VLSFAs by performing genome-wide association studies in seven population-based cohorts comprising 10,129 subjects of European ancestry. We observed associations of circulating VLSFA concentrations with common variants in two genes, serine palmitoyl-transferase long-chain base subunit 3 (SPTLC3), a gene involved in the rate-limiting step of de novo sphingolipid synthesis, and ceramide synthase 4 (CERS4). The SPTLC3 variant at rs680379 was associated with higher arachidic acid (20:0 , P = 5.81 × 10⁻¹³). The CERS4 variant at rs2100944 was associated with higher levels of 20:0 (P = 2.65 × 10⁻⁴⁰) and in analyses that adjusted for 20:0, with lower levels of behenic acid (P = 4.22 × 10⁻²⁶) and lignoceric acid (P = 3.20 × 10⁻²¹). These novel associations suggest an inter-relationship of circulating VLSFAs and sphingolipid synthesis.