Causal and Synthetic Associations of Variants in the SERPINA Gene Cluster with Alpha1-antitrypsin Serum Levels

Several infrequent genetic polymorphisms in the SERPINA1 gene are known to substantially reduce concentration of alpha1-antitrypsin (AAT) in the blood. Since low AAT serum levels fail to protect pulmonary tissue from enzymatic degradation, these polymorphisms also increase the risk for early onset chronic obstructive pulmonary disease (COPD). The role of more common SERPINA1 single nucleotide polymorphisms (SNPs) in respiratory health remains poorly understood.We present here an agnostic investigation of genetic determinants of circulating AAT levels in a general population sample by performing a genome-wide association study (GWAS) in 1392 individuals of the SAPALDIA cohort.Five common SNPs, defined by showing minor allele frequencies (MAFs) >5%, reached genome-wide significance, all located in the SERPINA gene cluster at 14q32.13. The top-ranking genotyped SNP rs4905179 was associated with an estimated effect of β = −0.068 g/L per minor allele (P = 1.20*10⁻¹²). But denser SERPINA1 locus genotyping in 5569 participants with subsequent stepwise conditional analysis, as well as exon-sequencing in a subsample (N = 410), suggested that AAT serum level is causally determined at this locus by rare (MAF<1%) and low-frequent (MAF 1–5%) variants only, in particular by the well-documented protein inhibitor S and Z (PI S, PI Z) variants. Replication of the association of rs4905179 with AAT serum levels in the Copenhagen City Heart Study (N = 8273) was successful (P<0.0001), as was the replication of its synthetic nature (the effect disappeared after adjusting for PI S and Z, P = 0.57). Extending the analysis to lung function revealed a more complex situation. Only in individuals with severely compromised pulmonary health (N = 397), associations of common SNPs at this locus with lung function were driven by rarer PI S or Z variants. Overall, our meta-analysis of lung function in ever-smokers does not support a functional role of common SNPs in the SERPINA gene cluster in the general population.     

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.     

Genetic Association Study of TNFAIP3, IFIH1, IRF5 Polymorphisms with Polymyositis/Dermatomyositis in Chinese Han Population

Background

Single-nucleotide polymorphisms (SNPs) in the TNFAIP3, IFIH1, and IRF5 genes have been associated with several auto-inflammation diseases, while the susceptibility between these genes and idiopathic inflammatory myopathies (IIMs) were not reported. This study aimed to investigate whether TNFAIP3, IFIH1, and IRF5 gene polymorphisms confer susceptibility for the IIMs in Chinese Han population.

Methods

A large case–control study of Chinese subjects with polymyositis (PM) (n = 298) and dermatomyositis (DM) (n = 530) was accomplished. 968 healthy and ethnically matched controls were available for comparison. Six SNPs in the TNFAIP3 region (rs2230926 and rs5029939), the IFIH1 gene (rs1990760 and rs3747517) and the IRF5 region (rs4728142 and rs729302) were assessed and genotyped using the Sequenom MassArray iPLEX platform.

Results

Our study indicated a strong allele association was observed in PM/DM and PM patients for rs2230926 (OR: 1.61, 95%CI: 1.20–2.16, P_c = 7.5×10⁻³; OR: 1.88, 95%CI: 1.30–2.74, P_c = 4.0×10⁻³, respectively) and rs5029939 (OR: 1.64, 95%CI: 1.21–2.21, P_c = 6.0×10⁻³; OR: 1.88, 95%CI: 1.28–2.76, P_c = 5.5×10⁻³,respectively). And rs2230926 and rs5029939 were significantly associated with interstitial lung disease (ILD) in PM/DM and PM patients (P_c = 0.04 and P_c = 0.016; P_c = 0.02 and P_c = 0.03, respectively). In addition, rs4728142 allele and genotype had significant association with PM/DM patients (P_c = 0.026 and P_c = 0.048, respectively). Further analysis with three logistic regression genetic models revealed statistically significant difference in the genotypic distribution in the PM/DM, PM or DM patients when the additive and dominant models were used.

Conclusions

This was the first study to reveal TNFAIP3 and IRF5 polymorphisms were associated with PM/DM patients or these patients with ILD, indicating that TNFAIP3 and IRF5 might be the susceptibility gene for PM/DM patients in Chinese Han population.     

Genome-Wide Association Study to Identify the Genetic Determinants of Otitis Media Susceptibility in Childhood

Background

Otitis media (OM) is a common childhood disease characterised by middle ear inflammation and effusion. Susceptibility to recurrent acute OM (rAOM; ≥3 episodes of AOM in 6 months) and chronic OM with effusion (COME; MEE ≥3 months) is 40–70% heritable. Few underlying genes have been identified to date, and no genome-wide association study (GWAS) of OM has been reported.

Methods and Findings

Data for 2,524,817 single nucleotide polymorphisms (SNPs; 535,544 quality-controlled SNPs genotyped by Illumina 660W-Quad; 1,989,273 by imputation) were analysed for association with OM in 416 cases and 1,075 controls from the Western Australian Pregnancy Cohort (Raine) Study. Logistic regression analyses under an additive model undertaken in GenABEL/ProbABEL adjusting for population substructure using principal components identified SNPs at CAPN14 (rs6755194: OR = 1.90; 95%CI 1.47–2.45; P_adj-PCA = 8.3×10⁻⁷) on chromosome 2p23.1 as the top hit, with independent effects (rs1862981: OR = 1.60; 95%CI 1.29–1.99; P_adj-PCA = 2.2×10⁻⁵) observed at the adjacent GALNT14 gene. In a gene-based analysis in VEGAS, BPIFA3 (P_Gene = 2×10⁻⁵) and BPIFA1 (P_Gene = 1.07×10⁻⁴) in the BPIFA gene cluster on chromosome 20q11.21 were the top hits. In all, 32 genomic regions show evidence of association (P_adj-PCA<10⁻⁵) in this GWAS, with pathway analysis showing a connection between top candidates and the TGFβ pathway. However, top and tag-SNP analysis for seven selected candidate genes in this pathway did not replicate in 645 families (793 affected individuals) from the Western Australian Family Study of Otitis Media (WAFSOM). Lack of replication may be explained by sample size, difference in OM disease severity between primary and replication cohorts or due to type I error in the primary GWAS.

Conclusions

This first discovery GWAS for an OM phenotype has identified CAPN14 and GALNT14 on chromosome 2p23.1 and the BPIFA gene cluster on chromosome 20q11.21 as novel candidate genes which warrant further analysis in cohorts matched more precisely for clinical phenotypes.     

The Associations between Immunity-Related Genes and Breast Cancer Prognosis in Korean Women

We investigated the role of common genetic variation in immune-related genes on breast cancer disease-free survival (DFS) in Korean women. 107 breast cancer patients of the Seoul Breast Cancer Study (SEBCS) were selected for this study. A total of 2,432 tag single nucleotide polymorphisms (SNPs) in 283 immune-related genes were genotyped with the GoldenGate Oligonucleotide pool assay (OPA). A multivariate Cox-proportional hazard model and polygenic risk score model were used to estimate the effects of SNPs on breast cancer prognosis. Harrell’s C index was calculated to estimate the predictive accuracy of polygenic risk score model. Subsequently, an extended gene set enrichment analysis (GSEA-SNP) was conducted to approximate the biological pathway. In addition, to confirm our results with current evidence, previous studies were systematically reviewed. Sixty-two SNPs were statistically significant at p-value less than 0.05. The most significant SNPs were rs1952438 in SOCS4 gene (hazard ratio (HR) = 11.99, 95% CI = 3.62–39.72, P = 4.84E-05), rs2289278 in TSLP gene (HR = 4.25, 95% CI = 2.10–8.62, P = 5.99E-05) and rs2074724 in HGF gene (HR = 4.63, 95% CI = 2.18–9.87, P = 7.04E-05). In the polygenic risk score model, the HR of women in the 3^rd tertile was 6.78 (95% CI = 1.48–31.06) compared to patients in the 1^st tertile of polygenic risk score. Harrell’s C index was 0.813 with total patients and 0.924 in 4-fold cross validation. In the pathway analysis, 18 pathways were significantly associated with breast cancer prognosis (P<0.1). The IL-6R, IL-8, IL-10RB, IL-12A, and IL-12B was associated with the prognosis of cancer in data of both our study and a previous study. Therefore, our results suggest that genetic polymorphisms in immune-related genes have relevance to breast cancer prognosis among Korean women.     

Low VHL mRNA Expression is Associated with More Aggressive Tumor Features of Papillary Thyroid Carcinoma

Alterations of the von Hippel–Lindau (VHL) tumor suppressor gene can cause different hereditary tumors associated with VHL syndrome, but the potential role of the VHL gene in papillary thyroid carcinoma (PTC) has not been characterized. This study set out to investigate the relationship of VHL expression level with clinicopathological features of PTC in an ethnically and geographically homogenous group of 264 patients from Serbia, for the first time. Multivariate logistic regression analysis showed a strong correlation between low level of VHL expression and advanced clinical stage (OR = 5.78, 95% CI 3.17–10.53, P<0.0001), classical papillary morphology of the tumor (OR = 2.92, 95% CI 1.33–6.44, P = 0.008) and multifocality (OR = 1.96, 95% CI 1.06–3.62, P = 0.031). In disease-free survival analysis, low VHL expression had marginal significance (P = 0.0502 by the log-rank test) but did not appear to be an independent predictor of the risk for chance of faster recurrence in a proportion hazards model. No somatic mutations or evidence of VHL downregulation via promoter hypermethylation in PTC were found. The results indicate that the decrease of VHL expression associates with tumor progression but the mechanism of downregulation remains to be elucidated.     

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.     

A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Effects of APOA5 ?1131T>C (rs662799) on Fasting Plasma Lipids and Risk of Metabolic Syndrome: Evidence from a Case-Control Study in China and a Meta-Analysis

The apolipoprotein A5 (APOA5) gene −1131T>C (rs662799) has been suggested to be involved in the pathway of lipid homeostasis and the development of metabolic syndrome (MetS). However, the findings are not consistent. To systematically evaluate the associations between −1131T>C polymorphism and fasting lipid parameters and the risk of MetS, we conducted a case-control study in a Chinese population and a meta-analysis. The findings from 1840 Chinese participants indicated that the C allele carriers had significantly higher fasting total cholesterol (TC), triglycerides (TG) and lower HDL-cholesterol (HDL-C) than the TT homozygotes carriers. The −1131C allele was also found to be significantly associated with increased risk of MetS (OR  =  1.40, 95% confidence interval (CI)  =  1.15, 1.69) compared to the TT homozygotes. In the meta-analysis of 51,868 participants from 46 East Asian studies, 26 European studies and 19 studies of other ethnic groups, the −1131C allele was associated with higher fasting TC (weighted mean difference (WMD)  =  0.08 mmol/L, 95% CI  =  0.05, 0.10, P = 1.74×10⁻⁹), TG (WMD  =  0.30 mmol/L, 95% CI  =  0.26, 0.33, P =  1.87×10⁻⁵⁵), LDL-cholesterol (LDL-C) (WMD  =  0.04 mmol/L, 95% CI  =  0.02, 0.07, P = 0.002), and lower HDL-C (WMD  =  −0.05 mmol/L, 95% CI  =  −0.06,−0.04, P = 1.88×10⁻²¹), respectively. Based on 12 studies with 5,573 MetS cases and 8,290 controls from 5 East Asian studies, 5 European studies and 2 studies of other ethnic groups, the −1131C allele was associated with increased risk of MetS with an OR (95% CI)  =  1.33 (1.16, 1.53) in the overall population, 1.43 (1.29, 1.58) in East Asian and 1.30 (0.94, 1.78) in European populations. In conclusion, the −1131C allele may be associated with elevated levels of fasting TG, TC, LDL-C and decreased HDL-C, and increased risk of MetS, especially in East Asians.     

Interactions between Urinary 4-tert-Octylphenol Levels and Metabolism Enzyme Gene Variants on Idiopathic Male Infertility

Octylphenol (OP) and Trichlorophenol (TCP) act as endocrine disruptors and have effects on male reproductive function. We studied the interactions between 4-tert-Octylphenol (4-t-OP), 4-n- Octylphenol (4-n-OP), 2,3,4-Trichlorophenol (2,3,4-TCP), 2,4,5-Trichlorophenol (2,4,5-TCP) urinary exposure levels and polymorphisms in selected xenobiotic metabolism enzyme genes among 589 idiopathic male infertile patients and 396 controls in a Han-Chinese population. Ultra high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) was used to measure alkylphenols and chlorophenols in urine. Polymorphisms were genotyped using the SNPstream platform and the Taqman method. Among four phenols that were detected, we found that only exposure to 4-t-OP increased the risk of male infertility (P _trend = 1.70×10⁻⁷). The strongest interaction was between 4-t-OP and rs4918758 in CYP2C9 (P _inter = 6.05×10⁻⁷). It presented a significant monotonic increase in risk estimates for male infertility with increasing 4-t-OP exposure levels among men with TC/CC genotype (low level compared with non-exposed, odds ratio (OR) = 2.26, 95% confidence intervals (CI) = 1.06, 4.83; high level compared with non-exposed, OR = 9.22, 95% CI = 2.78, 30.59), but no associations observed among men with TT genotype. We also found interactions between 4-t-OP and rs4986894 in CYP2C19, and between rs1048943 in CYP1A1, on male infertile risk (P _inter = 8.09×10⁻⁷, P _inter = 3.73×10⁻⁴, respectively).We observed notable interactions between 4-t-OP exposure and metabolism enzyme gene polymorphisms on idiopathic infertility in Han-Chinese men.     

Association between the rs2910164 Polymorphism in Pre-Mir-146a Sequence and Thyroid Carcinogenesis

Background

Rs2910164, a Single nucleotide polymorphism (SNP) located in the precursor microRNA sequence of miR-146a, is the only MicroRNA sequence SNP studied in papillary thyroid cancer (PTC). Association studies had been performed in US and UK-Northern European populations, but results were inconsistence. This study evaluated the association between rs2910164 and the risk of PTC as well as benign thyroid tumor (BN), and examined the clinicopathological characteristics of PTC and BN for different genotypes.

Methods

This case-control study genotyped rs2910164 in 753 PTCs, 484 BNs and 760 controls in a Chinese Han population. Clinicopathological and genetic data were collected and compared. Multivariate logistic regression was performed to calculate adjusted odds ratios (ORs).

Results

There were no differences in rs2910164 genotype distributions between the three groups. PTC cases with three genotypes (CC, CG, GG) had similar clinicopathological characteristics except the existence of “para-cancer” BN (PTC/BN, P = 0.006). PTC/BN patients were older (P = 0.009), and had smaller cancer lesions (P<0.001), lower serum thyrotropin levels (1.82±1.42 vs. 2.21±1.74, P = 0.04), and lower rates of level VI lymph node metastasis (20.8% vs. 52.7%, P<0.001) and lateral neck lymph node metastasis (11.5% vs. 23.0%, P = 0.011) compared with PTC only. Then we supposed a possible progression from BN to PTC which may involve rs2910164 in and performed a multivariate logistic regression analysis of PTC/BN and BN cases to determine risk factors of this progression. Results showed that the rs2910164 GG homozygote (OR = 2.25, 95% CI 1.22–4.14, P = 0.01) was the only risk factor in this study.

Conclusion

Rs2910164 was not associated with increased risk of PTC and BN in Chinese patients, but may play a latent role in the transformation from BN to PTC.     

Pathway-Based Analysis Using Genome-wide Association Data from a Korean Non-Small Cell Lung Cancer Study

Pathway-based analysis, used in conjunction with genome-wide association study (GWAS) techniques, is a powerful tool to detect subtle but systematic patterns in genome that can help elucidate complex diseases, like cancers. Here, we stepped back from genetic polymorphisms at a single locus and examined how multiple association signals can be orchestrated to find pathways related to lung cancer susceptibility. We used single-nucleotide polymorphism (SNP) array data from 869 non-small cell lung cancer (NSCLC) cases from a previous GWAS at the National Cancer Center and 1,533 controls from the Korean Association Resource project for the pathway-based analysis. After mapping single-nucleotide polymorphisms to genes, considering their coding region and regulatory elements (±20 kbp), multivariate logistic regression of additive and dominant genetic models were fitted against disease status, with adjustments for age, gender, and smoking status. Pathway statistics were evaluated using Gene Set Enrichment Analysis (GSEA) and Adaptive Rank Truncated Product (ARTP) methods. Among 880 pathways, 11 showed relatively significant statistics compared to our positive controls (P_GSEA≤0.025, false discovery rate≤0.25). Candidate pathways were validated using the ARTP method and similarities between pathways were computed against each other. The top-ranked pathways were ABC Transporters (P_GSEA<0.001, P_ARTP = 0.001), VEGF Signaling Pathway (P_GSEA<0.001, P_ARTP = 0.008), G1/S Check Point (P_GSEA = 0.004, P_ARTP = 0.013), and NRAGE Signals Death through JNK (P_GSEA = 0.006, P_ARTP = 0.001). Our results demonstrate that pathway analysis can shed light on post-GWAS research and help identify potential targets for cancer susceptibility.     

Association of Variants at UMOD with Chronic Kidney Disease and Kidney Stones—Role of Age and Comorbid Diseases

Chronic kidney disease (CKD) is a worldwide public health problem that is associated with substantial morbidity and mortality. To search for sequence variants that associate with CKD, we conducted a genome-wide association study (GWAS) that included a total of 3,203 Icelandic cases and 38,782 controls. We observed an association between CKD and a variant with 80% population frequency, rs4293393-T, positioned next to the UMOD gene (GeneID: 7369) on chromosome 16p12 (OR = 1.25, P = 4.1×10⁻¹⁰). This gene encodes uromodulin (Tamm-Horsfall protein), the most abundant protein in mammalian urine. The variant also associates significantly with serum creatinine concentration (SCr) in Icelandic subjects (N = 24,635, P = 1.3×10⁻²³) but not in a smaller set of healthy Dutch controls (N = 1,819, P = 0.39). Our findings validate the association between the UMOD variant and both CKD and SCr recently discovered in a large GWAS. In the Icelandic dataset, we demonstrate that the effect on SCr increases substantially with both age (P = 3.0×10⁻¹⁷) and number of comorbid diseases (P = 0.008). The association with CKD is also stronger in the older age groups. These results suggest that the UMOD variant may influence the adaptation of the kidney to age-related risk factors of kidney disease such as hypertension and diabetes. The variant also associates with serum urea (P = 1.0×10⁻⁶), uric acid (P = 0.0064), and suggestively with gout. In contrast to CKD, the UMOD variant confers protection against kidney stones when studied in 3,617 Icelandic and Dutch kidney stone cases and 43,201 controls (OR = 0.88, P = 5.7×10⁻⁵).     

Hybrid Models Identified a 12-Gene Signature for Lung Cancer Prognosis and Chemoresponse Prediction

Background

Lung cancer remains the leading cause of cancer-related deaths worldwide. The recurrence rate ranges from 35–50% among early stage non-small cell lung cancer patients. To date, there is no fully-validated and clinically applied prognostic gene signature for personalized treatment.

Methodology/Principal Findings

From genome-wide mRNA expression profiles generated on 256 lung adenocarcinoma patients, a 12-gene signature was identified using combinatorial gene selection methods, and a risk score algorithm was developed with Naïve Bayes. The 12-gene model generates significant patient stratification in the training cohort HLM & UM (n = 256; log-rank P = 6.96e-7) and two independent validation sets, MSK (n = 104; log-rank P = 9.88e-4) and DFCI (n = 82; log-rank P = 2.57e-4), using Kaplan-Meier analyses. This gene signature also stratifies stage I and IB lung adenocarcinoma patients into two distinct survival groups (log-rank P<0.04). The 12-gene risk score is more significant (hazard ratio = 4.19, 95% CI: [2.08, 8.46]) than other commonly used clinical factors except tumor stage (III vs. I) in multivariate Cox analyses. The 12-gene model is more accurate than previously published lung cancer gene signatures on the same datasets. Furthermore, this signature accurately predicts chemoresistance/chemosensitivity to Cisplatin, Carboplatin, Paclitaxel, Etoposide, Erlotinib, and Gefitinib in NCI-60 cancer cell lines (P<0.017). The identified 12 genes exhibit curated interactions with major lung cancer signaling hallmarks in functional pathway analysis. The expression patterns of the signature genes have been confirmed in RT-PCR analyses of independent tumor samples.

Conclusions/Significance

The results demonstrate the clinical utility of the identified gene signature in prognostic categorization. With this 12-gene risk score algorithm, early stage patients at high risk for tumor recurrence could be identified for adjuvant chemotherapy; whereas stage I and II patients at low risk could be spared the toxic side effects of chemotherapeutic drugs.     

Ancestry-Shift Refinement Mapping of the C6orf97-ESR1 Breast Cancer Susceptibility Locus

We used an approach that we term ancestry-shift refinement mapping to investigate an association, originally discovered in a GWAS of a Chinese population, between rs2046210[T] and breast cancer susceptibility. The locus is on 6q25.1 in proximity to the C6orf97 and estrogen receptor α (ESR1) genes. We identified a panel of SNPs that are correlated with rs2046210 in Chinese, but not necessarily so in other ancestral populations, and genotyped them in breast cancer case∶control samples of Asian, European, and African origin, a total of 10,176 cases and 13,286 controls. We found that rs2046210[T] does not confer substantial risk of breast cancer in Europeans and Africans (OR = 1.04, P = 0.099, and OR = 0.98, P = 0.77, respectively). Rather, in those ancestries, an association signal arises from a group of less common SNPs typified by rs9397435. The rs9397435[G] allele was found to confer risk of breast cancer in European (OR = 1.15, P = 1.2×10⁻³), African (OR = 1.35, P = 0.014), and Asian (OR = 1.23, P = 2.9×10⁻⁴) population samples. Combined over all ancestries, the OR was 1.19 (P = 3.9×10⁻⁷), was without significant heterogeneity between ancestries (P_het = 0.36) and the SNP fully accounted for the association signal in each ancestry. Haplotypes bearing rs9397435[G] are well tagged by rs2046210[T] only in Asians. The rs9397435[G] allele showed associations with both estrogen receptor positive and estrogen receptor negative breast cancer. Using early-draft data from the 1,000 Genomes project, we found that the risk allele of a novel SNP (rs77275268), which is closely correlated with rs9397435, disrupts a partially methylated CpG sequence within a known CTCF binding site. These studies demonstrate that shifting the analysis among ancestral populations can provide valuable resolution in association mapping.     

Fine Mapping and Functional Analysis of the Multiple Sclerosis Risk Gene CD6

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2^nd SRCR domain with susceptibility to MS (P _{max(T) permutation} = 1×10⁻⁴). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. – CD4⁺ naïve cells, P = 0.0001; CD8⁺ naïve cells, P<0.0001; CD4⁺ and CD8⁺ central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4⁺ and CD8⁺ T cells.     

The Pro12Ala Polymorphism in the Peroxisome Proliferator-Activated Receptor Gamma-2 Gene (PPARγ2) Is Associated with Increased Risk of Coronary Artery Disease: A Meta-Analysis

Background

Contradictory results have been reported regarding the association between Pro12Ala polymorphism of PPARγ2 and coronary artery disease (CAD). We sought to estimate the inconsistent results by performing a comprehensive meta-analysis.

Methods

Studies in English or Chinese publications were identified by screening MEDLINE, Embase, CNKI, Wanfang and CBM. 22 studies including 8948 cases and 14427 controls were selected. A random-effects model was applied to combine the divergent outcomes of the individual studies, while addressing between-study heterogeneity and publication bias.

Results

The Pro12Ala polymorphism of control population followed Hardy-Weinberg equilibrium for all studies (P>0.05). Overall, a marginal increased risk of CAD under the recessive genetic model (AlaAla vs ProAla+ProPro: P = 0.04, OR = 1.31, 95%CI 1.01–1.69, P_{heterogeneity} = 0.67, I² = 0%) and the homozygote comparison (AlaAla vs ProPro: P = 0.04,OR = 1.30, 95%CI 1.01–1.68, P_{heterogeneity} = 0.68, I² = 0%) was observed. In the subgroup analysis by ethnicity, carriers of AlaAla homozygotes had a significant increased risk for CAD among Caucasians (AlaAla vs ProAla+ProPro: P = 0.01, OR = 1.45, 95%CI 1.08–1.96, P_{heterogeneity} = 0.48, I² = 0%; AlaAla vs ProPro: P = 0.02,OR = 1.44, 95%CI 1.07–1.93, P_{heterogeneity} = 0.46, I² = 0%). After dividing into population source, the CAD risk magnitude of hospital-based studies was distinctly strengthened under the recessive model (P = 0.03,OR = 1.85,95%CI 1.07–3.19, P_{heterogeneity} = 0.87,I² = 0%) and the homozygote comparison (P = 0.03,OR = 1.83, 95%CI 1.06–3.16, P_{heterogeneity} = 0.88, I² = 0%). There was no observable publication bias as reflected by funnel plot and Egger’s linear regression test (t = -0.12, P = 0.91).

Conclusion:

Our results demonstrated that the PPARγ2 Pro12Ala polymorphism might be risk-conferring locus for the progression of CAD among Caucasians, but not among Asians.