Genome-wide Association Studies Single Nucleotide Polymorphisms Research Articles

A GWAS SNP for Schizophrenia Is Linked to the Internal MIR137 Promoter and Supports Differential Allele-Specific Expression

Single nucleotide polymorphisms (SNPs) within the MIR137 gene locus have been shown to confer risk for schizophrenia through genome-wide association studies (GWAS). The expression levels of microRNA-137 (miR-137) and its validated gene targets have also been shown to be disrupted in several neuropsychiatric conditions, including schizophrenia. Regulation of miR-137 expression is thus imperative for normal neuronal functioning. We previously characterized an internal promoter domain within the MIR137 gene that contained a variable number tandem repeat (VNTR) polymorphism and could alter the in vitro levels of miR-137 in a stimulus-induced and allele-specific manner. We now demonstrate that haplotype tagging-SNP analysis linked the rs1625579 GWAS SNP for schizophrenia to this internal MIR137 promoter through a proxy SNP rs2660304 located at this domain. We postulated that the rs2660304 promoter SNP may act as predisposing factor for schizophrenia through altering the levels of miR-137 expression in a genotype-dependent manner. Reporter gene analysis of the internal MIR137 promoter containing the common VNTR variant demonstrated genotype-dependent differences in promoter activity with respect to rs2660304. In line with previous reports, the major allele of the rs2660304 proxy SNP, which has previously been linked with schizophrenia risk through genetic association, resulted in downregulation of reporter gene expression in a tissue culture model. The genetic influence of the rs2660304 proxy SNP on the transcriptional activity of the internal MIR137 promoter, and thus the levels of miR-137 expression, therefore offers a distinct regulatory mechanism to explain the functional significance of the rs1625579 GWAS SNP for schizophrenia risk.

Genome-wide association study of new-onset atrial fibrillation after coronary artery bypass grafting surgery

Postoperative atrial fibrillation (AF) is a potentially life-threatening complication after coronary artery bypass graft (CABG) surgery. Genetic predisposition may predict risk for developing postoperative AF. Study subjects underwent CABG surgery with cardiopulmonary bypass at Duke University Medical Center. In a discovery cohort of 877 individuals from the Perioperative Genetics and Safety Outcomes Study, we performed a genome-wide association study using a logistic regression model with a covariate adjustment for AF risk index. Single-nucleotide polymorphisms (SNPs) that met a P < 5 × 10(-5) were further tested using a replication dataset of 304 individuals from the CATHeterization GENetics biorepository, followed by meta-analysis. Potential pathways related to postoperative AF were identified through gene enrichment analysis using the top genome-wide association study SNPs (P < 10(-4)). Nineteen SNPs met the a priori defined discovery threshold for replication, but only 3 met nominal significance (P < .05) in the CATHeterization GENetics group, with only one-rs10504554, in the intronic region in lymphocyte antigen 96 (LY96)-showing the same direction of the effect for postoperative AF (odds ratio [OR] 0.48, 95% CI 0.34-0.68, P = 2.9 × 10(-5) vs OR 0.55, 95% CI 0.31-0.99, P = .046) and strong overall association by meta-analysis (meta-P = 4.0 × 10(-6)). Gene enrichment analysis highlighted the role of LY96 in pathways of biologic relevance to activation and modulation of innate immune responses. Our analysis also showed potential association between LY96 and nuclear factor κ-B interaction and postoperative AF through their relevance to inflammatory signaling pathways. In patients undergoing CABG surgery, we found genetic polymorphisms in LY96 associated with decreased risk of postoperative AF.

Genome-wide association study of corticobasal degeneration identifies risk variants shared with progressive supranuclear palsy.

Corticobasal degeneration (CBD) is a neurodegenerative disorder affecting movement and cognition, definitively diagnosed only at autopsy. Here, we conduct a genome-wide association study (GWAS) in CBD cases (n=152) and 3,311 controls, and 67 CBD cases and 439 controls in a replication stage. Associations with meta-analysis were 17q21 at MAPT (P=1.42 × 10−12), 8p12 at lnc-KIF13B-1, a long non-coding RNA (rs643472; P=3.41 × 10−8), and 2p22 at SOS1 (rs963731; P=1.76 × 10−7). Testing for association of CBD with top progressive supranuclear palsy (PSP) GWAS single-nucleotide polymorphisms (SNPs) identified associations at MOBP (3p22; rs1768208; P=2.07 × 10−7) and MAPT H1c (17q21; rs242557; P=7.91 × 10−6). We previously reported SNP/transcript level associations with rs8070723/MAPT, rs242557/MAPT, and rs1768208/MOBP and herein identified association with rs963731/SOS1. We identify new CBD susceptibility loci and show that CBD and PSP share a genetic risk factor other than MAPT at 3p22 MOBP (myelin-associated oligodendrocyte basic protein).

Analysis of Prostate Cancer Susceptibility Variants in South African Men: Replicating Associations on Chromosomes 8q24 and 10q11.

Genome-wide association studies (GWAS) have implicated single nucleotide polymorphisms (SNPs) on chromosomes 2p15, 6q25, 7p15.2, 7q21, 8q24, 10q11, 10q26, 11q13, 17q12, 17q24, 19q13, and Xp11, with prostate cancer (PCa) susceptibility and/or tumour aggressiveness, in populations of African, European, and Asian ancestry. The objective of this study was to confirm these associations in South African Mixed Ancestry and White men. We evaluated 17 prioritised GWAS SNPs in South African cases (331 Mixed Ancestry and 155 White) and controls (178 Mixed Ancestry and 145 White). The replicated SNP associations for the different South African ethnic groups were rs7008482 (8q24) (p = 2.45 × 10−5), rs6983267 (8q24) (p = 4.48 × 10−7), and rs10993994 (10q11) (p = 1.40 × 10−3) in Mixed Ancestry men and rs10993994 (p = 1.56 × 10−9) in White men. No significant associations were observed for the analyses stratified by disease aggressiveness in the individual and the combined population group analysis. The present study demonstrates that a number of known PCa susceptibility variants may contribute to disease susceptibility in South African men. Larger genetic investigations extended to other South African population groups are warranted to confirm the role of these and other SNPs in disease susceptibility.

Genetic and environmental determinants of 25-hydroxyvitamin D levels in multiple sclerosis

Background: Evidence is accumulating supporting a beneficial effect of vitamin D in multiple sclerosis (MS). Genome-wide association studies (GWAS) have shown significant associations between 25-hydroxyvitamin D (25(OH)D) and single nucleotide polymorphisms (SNPs) in key genes in the vitamin D metabolism. Objective: To examine the association between 25(OH)D and six GWAS SNPs and environmental factors in 1497 MS patients. Methods: Blood samples and lifestyle questionnaires were collected between 2009 and 2012. Genotyping of GC-, CYP2R1- and NADSYN1-SNPs was performed by TaqMan allelic discrimination (Life Technologies). Results: We found significant associations between 25(OH)D and SNPs in GC (rs7041, p = 0.01 and rs2282679, p = 0.03) and CYP2R1 (rs10741657, p =1.8 × 10−4). Season of blood sampling (p = 2.8 × 10−31), sex (p = 1.9 × 10−5), BMI (p = 2.3 × 10−5), vitamin supplements (p = 7.0 × 10−22), and fish intake (p = 0.02) also had significant effects on 25(OH)D. Conclusion: In this cross-sectional study, we found significant effects of environmental factors and SNPs in GC and CYP2R1 on 25(OH)D in MS patients. Since 25(OH)D might have protective effects in MS, and vitamin D supply is a modifiable factor, it may be important to include this in the MS treatment regimen.

Integrative genomics identifies 7p11.2 as a novel locus for fever and clinical stress response in humans.

Fever predicts clinical outcomes in sepsis, trauma and during cardiovascular stress, yet the genetic determinants are poorly understood. We used an integrative genomics approach to identify novel genomic determinants of the febrile response to experimental endotoxemia. We highlight multiple integrated lines of evidence establishing the clinical relevance of this novel fever locus. Through genome-wide association study (GWAS) of evoked endotoxemia (lipopolysaccharide (LPS) 1 ng/kg IV) in healthy subjects of European ancestry we discovered a locus on chr7p11.2 significantly associated with the peak febrile response to LPS (top single nucleotide polymorphism (SNP) rs7805622, P = 2.4 × 10(-12)), as well as with temperature fluctuation over time. We replicated this association in a smaller independent LPS study (rs7805622, P = 0.03). In clinical translation, this locus was also associated with temperature and mortality in critically ill patients with trauma or severe sepsis. The top GWAS SNPs are not located within protein-coding genes, but have significant cis-expression quantitative trait loci (eQTL) associations with expression of a cluster of genes ∼400 kb upstream, several of which (SUMF2, CCT6A, GBAS) are regulated by LPS in vivo in blood cells. LPS- and cold-treatment of adipose stromal cells in vitro suggest genotype-specific modulation of eQTL candidate genes (PSPH). Several eQTL genes were up-regulated in brown and white adipose following cold exposure in mice, highlighting a potential role in thermogenesis. Thus, through genomic interrogation of experimental endotoxemia, we identified and replicated a novel fever locus on chr7p11.2 that modulates clinical responses in trauma and sepsis, and highlight integrated in vivo and in vitro evidence for possible novel cis candidate genes conserved across human and mouse.

Targeted allelic expression profiling in human islets identifies cis-regulatory effects for multiple variants identified by type 2 diabetes genome-wide association studies.

Genome-wide association studies (GWAS) have identified variation at >65 genomic loci associated with susceptibility to type 2 diabetes, but little progress has been made in elucidating the molecular mechanisms behind most of these associations. Using samples heterozygous for transcribed single nucleotide polymorphisms (SNPs), allelic expression profiling is a powerful technique for identifying cis-regulatory variants controlling gene expression. In this study, exonic SNPs, suitable for measuring mature mRNA levels and in high linkage disequilibrium with 65 lead type 2 diabetes GWAS SNPs, were identified and allelic expression determined by real-time PCR using RNA and DNA isolated from islets of 36 white nondiabetic donors. A significant allelic expression imbalance (AEI) was identified for 7/14 (50%) genes tested (ANPEP, CAMK2B, HMG20A, KCNJ11, NOTCH2, SLC30A8, and WFS1), with significant AEI confirmed for five of these genes using other linked exonic SNPs. Lastly, results of a targeted islet expression quantitative trait loci experiment support the AEI findings for ANPEP, further implicating ANPEP as the causative gene at its locus. The results of this study support the hypothesis that changes to cis-regulation of gene expression are involved in a large proportion of SNP associations with type 2 diabetes susceptibility.

Abstract B14: Population genetics analysis of prostate cancer GWAS SNPs to evaluate West African genetic ancestry as a risk factor

Abstract Background: Prostate cancer (Pca) is the most common cancer among men in the U.S., and African American (AA) men have higher incidence compared to European American men. Environmental and sociocultural factors may play a role for increasing Pca risk among AA men, but the causal relationship of any environmental and social/cultural factors with Pca risk is not well demonstrated. On the other hand, there is strong evidence for genetic factors that explain increased incidence in AAs. Admixture mapping identified loci at 8q24 region showing association between African genetic ancestry and Pca risk, and a genome-wide association study (GWAS) in AAs identified a novel locus at 17q21 associated with Pca. However, replication of variants identified in GWAS among European and Asian populations has been difficult in AAs. Here, we performed population genetics analyses of Pca GWAS identified single nucleotide polymorphisms (SNPs) to evaluate whether these SNPs are informative of genetic ancestry and whether West African genetic Ancestry (WAA) is a Pca risk factor. Samples and Methods: We accessed National Human Genome Research Institute GWAS catalog in December, 2012, and there were 80 SNPs reported in 17 studies. We downloaded genotype data of 76 SNPs in 1000 Genomes (1KG) Project Phase I May 2011 and 62 SNPs in Human Genome Diversity Project (HGDP) from SPSmart. We performed principal component analysis using Golden Helix SVS 7. STRUCTURE was also used to estimate genetic ancestry and examine population clustering patterns. Replication study subjects, 470 AA men from Washington, DC (220 Pca cases and 250 controls) were genotyped as a part of GWAS in AAs with Illumina Infinium 1M. Logistic regression analysis was performed adjusting for age and 5 principal components (PCs). Results: Mean heterozygosity of 76 Pca GWAS SNPs in the 1KG Project dataset was lower in African populations than in European populations, possibly due to ascertainment bias. The majority of Pca GWAS SNPs had FST less than 0.20, but rs9623117 on chromosome 22 was highly differentiated (FST=0.63). In the 1KG Project dataset, PC1 separated African populations from non-African populations, and PC2 reflects variation in non-African populations separating East Asian from European populations. AA (ASW) individuals fell between the African and European clusters. Fewer Pca GWAS SNPs were genotyped in HGDP dataset, but they successfully captured global human genome variation forming African and non-African population clusters. In the STRUCTURE analysis with the 1KG project dataset, K=2 separated African from non-African populations, and K=3 separated European from East Asian populations. WAA estimated using 76 Pca GWAS SNPs was significantly correlated with WAA estimated using ancestry informative markers. Only four Pca GWAS SNPs were replicated and three of them were independently associated (P&amp;lt;0.05) in our replication samples. The number of risk alleles were positively associated with Pca risk (PTREND=0.002), and AAs carrying ≥5 risk alleles had 3.5 times greater odds of Pca (95% C.I.; 1.36-9.17; P=0.009) compared to AAs with ≤ 2 risk alleles. Risk allele frequencies of the three independently associated SNPs were higher in West Africans than in Europeans and was positively correlated with WAA estimates (PANOVA&amp;lt;0.001). Conclusion: Results of our analyses provide evidence suggesting that WAA is a risk factor for Pca. Population genetics analyses show that Pca GWAS SNPs can successfully capture global human genome variation and are informative of genetic ancestry. AAs with high WAA are likely to carry more Pca risk alleles than AAs with low WAA. Citation Format: Ken Batai, Ebony Shah, Rick A. Kittles. Population genetics analysis of prostate cancer GWAS SNPs to evaluate West African genetic ancestry as a risk factor. [abstract]. In: Proceedings of the Sixth AACR Conference: The Science of Cancer Health Disparities; Dec 6–9, 2013; Atlanta, GA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2014;23(11 Suppl):Abstract nr B14. doi:10.1158/1538-7755.DISP13-B14

Genetic variation in HTR4 and lung function: GWAS follow-up in mouse.

Human genome-wide association studies (GWASs) have identified numerous associations between single nucleotide polymorphisms (SNPs) and pulmonary function. Proving that there is a causal relationship between GWAS SNPs, many of which are noncoding and without known functional impact, and these traits has been elusive. Furthermore, noncoding GWAS-identified SNPs may exert trans-regulatory effects rather than impact the proximal gene. Noncoding variants in 5-hydroxytryptamine (serotonin) receptor 4 (HTR4) are associated with pulmonary function in human GWASs. To gain insight into whether this association is causal, we tested whether Htr4-null mice have altered pulmonary function. We found that HTR4-deficient mice have 12% higher baseline lung resistance and also increased methacholine-induced airway hyperresponsiveness (AHR) as measured by lung resistance (27%), tissue resistance (48%), and tissue elastance (30%). Furthermore, Htr4-null mice were more sensitive to serotonin-induced AHR. In models of exposure to bacterial lipopolysaccharide, bleomycin, and allergic airway inflammation induced by house dust mites, pulmonary function and cytokine profiles in Htr4-null mice differed little from their wild-type controls. The findings of altered baseline lung function and increased AHR in Htr4-null mice support a causal relationship between genetic variation in HTR4 and pulmonary function identified in human GWAS.

Functional annotation of colon cancer risk SNPs.

Colorectal cancer (CRC) is a leading cause of cancer-related deaths in the United States. Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with increased risk for CRC. A molecular understanding of the functional consequences of this genetic variation has been complicated because each GWAS SNP is a surrogate for hundreds of other SNPs, most of which are located in non-coding regions. Here we use genomic and epigenomic information to test the hypothesis that the GWAS SNPs and/or correlated SNPs are in elements that regulate gene expression, and identify 23 promoters and 28 enhancers. Using gene expression data from normal and tumour cells, we identify 66 putative target genes of the risk-associated enhancers (10 of which were also identified by promoter SNPs). Employing CRISPR nucleases, we delete one risk-associated enhancer and identify genes showing altered expression. We suggest that similar studies be performed to characterize all CRC risk-associated enhancers.

Replication study for the association of seven genome- GWAS-identified Loci with susceptibility to ovarian cancer in the Polish population.

We investigated the previously-demonstrated association of seven genome-wide association studies (GWAS) single nucleotide polymorphisms (SNPs), including rs2072590 (HOXD-AS1), rs2665390 (TIPARP), rs10088218 and rs10098821 (8q24), rs3814113 (9p22), rs9303542 (SKAP1) and rs2363956 (ANKLE1), as risk factors of epithelial ovarian tumors (EOTs). These SNPs were genotyped in two hundred seventy three patients with EOTs and four hundred sixty four unrelated healthy females from the Polish population. We observed the lowest p values of the trend test for the 9p22 rs3814113 and 8q24 rs10098821 SNPs in patients with all subtypes of ovarian cancer (ptrend = 0.010 and ptrend = 0.014, respectively). There were also significant p values for the trend of the 9p22 rs3814113 and the 8q24 rs10098821 SNPs for serous histological subtypes of ovarian cancer (ptrend = 0.006, ptrend = 0.033, respectively). Moreover, stratification of the patients based on their histological type of cancer demonstrated, in the dominant hereditary model, a significant association of the 9p22 rs3814113 SNP with serous ovarian carcinoma OR = 0.532 (95 % CI = 0.342 - 0.827, p = 0.005, pcorr = 0.035). Despite the relatively small sample size of cases and controls, our studies confirmed some of the previously-demonstrated GWAS SNPs as genetic risk factors for EOTs.Electronic supplementary materialThe online version of this article (doi:10.1007/s12253-014-9822-6) contains supplementary material, which is available to authorized users.

Aromatase inhibitor-associated bone fractures: a case-cohort GWAS and functional genomics.

Bone fractures are a major consequence of osteoporosis. There is a direct relationship between serum estrogen concentrations and osteoporosis risk. Aromatase inhibitors (AIs) greatly decrease serum estrogen levels in postmenopausal women, and increased incidence of fractures is a side effect of AI therapy. We performed a discovery case-cohort genome-wide association study (GWAS) using samples from 1071 patients, 231 cases and 840 controls, enrolled in the MA.27 breast cancer AI trial to identify genetic factors involved in AI-related fractures, followed by functional genomic validation. Association analyses identified 20 GWAS single nucleotide polymorphism (SNP) signals with P < 5E-06. After removal of signals in gene deserts and those composed entirely of imputed SNPs, we applied a functional validation "decision cascade" that resulted in validation of the CTSZ-SLMO2-ATP5E, TRAM2-TMEM14A, and MAP4K4 genes. These genes all displayed estradiol (E2)-dependent induction in human fetal osteoblasts transfected with estrogen receptor-α, and their knockdown altered the expression of known osteoporosis-related genes. These same genes also displayed SNP-dependent variation in E2 induction that paralleled the SNP-dependent induction of known osteoporosis genes, such as osteoprotegerin. In summary, our case-cohort GWAS identified SNPs in or near CTSZ-SLMO2-ATP5E, TRAM2-TMEM14A, and MAP4K4 that were associated with risk for bone fracture in estrogen receptor-positive breast cancer patients treated with AIs. These genes displayed E2-dependent induction, their knockdown altered the expression of genes related to osteoporosis, and they displayed SNP genotype-dependent variation in E2 induction. These observations may lead to the identification of novel mechanisms associated with fracture risk in postmenopausal women treated with AIs.

Effects of GWAS-associated genetic variants on lncRNAs within IBD and T1D candidate loci.

Long non-coding RNAs are a new class of non-coding RNAs that are at the crosshairs in many human diseases such as cancers, cardiovascular disorders, inflammatory and autoimmune disease like Inflammatory Bowel Disease (IBD) and Type 1 Diabetes (T1D). Nearly 90% of the phenotype-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) lie outside of the protein coding regions, and map to the non-coding intervals. However, the relationship between phenotype-associated loci and the non-coding regions including the long non-coding RNAs (lncRNAs) is poorly understood. Here, we systemically identified all annotated IBD and T1D loci-associated lncRNAs, and mapped nominally significant GWAS/ImmunoChip SNPs for IBD and T1D within these lncRNAs. Additionally, we identified tissue-specific cis-eQTLs, and strong linkage disequilibrium (LD) signals associated with these SNPs. We explored sequence and structure based attributes of these lncRNAs, and also predicted the structural effects of mapped SNPs within them. We also identified lncRNAs in IBD and T1D that are under recent positive selection. Our analysis identified putative lncRNA secondary structure-disruptive SNPs within and in close proximity (+/−5 kb flanking regions) of IBD and T1D loci-associated candidate genes, suggesting that these RNA conformation-altering polymorphisms might be associated with diseased-phenotype. Disruption of lncRNA secondary structure due to presence of GWAS SNPs provides valuable information that could be potentially useful for future structure-function studies on lncRNAs.

Identification of genetic variants associated with alternative splicing using sQTLseekeR.

Identification of genetic variants affecting splicing in RNA sequencing population studies is still in its infancy. Splicing phenotype is more complex than gene expression and ought to be treated as a multivariate phenotype to be recapitulated completely. Here we represent the splicing pattern of a gene as the distribution of the relative abundances of a gene’s alternative transcript isoforms. We develop a statistical framework that uses a distance-based approach to compute the variability of splicing ratios across observations, and a non-parametric analogue to multivariate analysis of variance. We implement this approach in the R package sQTLseekeR and use it to analyze RNA-Seq data from the Geuvadis project in 465 individuals. We identify hundreds of single nucleotide polymorphisms (SNPs) as splicing QTLs (sQTLs), including some falling in genome-wide association study SNPs. By developing the appropriate metrics, we show that sQTLseekeR compares favorably with existing methods that rely on univariate approaches, predicting variants that behave as expected from mutations affecting splicing.

Common vitamin D pathway gene variants reveal contrasting effects on serum vitamin D levels in African Americans and European Americans.

Vitamin D deficiency is more common among African Americans (AAs) than among European Americans (EAs), and epidemiologic evidence links vitamin D status to many health outcomes. Two genome-wide association studies (GWAS) in European populations identified vitamin D pathway gene single-nucleotide polymorphisms (SNPs) associated with serum vitamin D [25(OH)D] levels, but a few of these SNPs have been replicated in AAs. Here, we investigated the associations of 39 SNPs in vitamin D pathway genes, including 19 GWAS-identified SNPs, with serum 25(OH)D concentrations in 652 AAs and 405 EAs. Linear and logistic regression analyses were performed adjusting for relevant environmental and biological factors. The pattern of SNP associations was distinct between AAs and EAs. In AAs, six GWAS-identified SNPs in GC, CYP2R1, and DHCR7/NADSYN1 were replicated, while nine GWAS SNPs in GC and CYP2R1 were replicated in EAs. A CYP2R1 SNP, rs12794714, exhibited the strongest signal of association in AAs. In EAs, however, a different CYP2R1 SNP, rs1993116, was the most strongly associated. Our models, which take into account genetic and environmental variables, accounted for 20 and 28 % of the variance in serum vitamin D levels in AAs and EAs, respectively.Electronic supplementary materialThe online version of this article (doi:10.1007/s00439-014-1472-y) contains supplementary material, which is available to authorized users.

Top associated SNPs in prostate cancer are significantly enriched incis-expression quantitative trait loci and at transcription factor binding sites

While genome-wide association studies (GWAS) have revealed thousands of disease risk single nucleotide polymorphisms (SNPs), their functions remain largely unknown. Recent studies have suggested the regulatory roles of GWAS risk variants in several common diseases; however, the complex regulatory structure in prostate cancer is unclear. We investigated the potential regulatory roles of risk variants in two prostate cancer GWAS datasets by their interactions with expression quantitative trait loci (eQTL) and/or transcription factor binding sites (TFBSs) in three populations. Our results indicated that the moderately associated GWAS SNPs were significantly enriched with cis-eQTLs and TFBSs in Caucasians (CEU), but not in African Americans (AA) or Japanese (JPT); this was also observed in an independent pan-cancer related SNPs from the GWAS Catalog. We found that the eQTL enrichment in the CEU population was tissue-specific to eQTLs from CEU lymphoblastoid cell lines. Importantly, we pinpointed two SNPs, rs2861405 and rs4766642, by overlapping results from cis-eQTL and TFBS as applied to the CEU data. These results suggested that prostate cancer associated SNPs and pan-cancer associated SNPs are likely to play regulatory roles in CEU. However, the negative enrichment results in AA or JPT and the potential mechanisms remain to be elucidated in additional samples.

Sequence variation in TMEM18 in association with body mass index: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study.

Genome-wide association studies for body mass index (BMI) previously identified a locus near TMEM18. We conducted targeted sequencing of this region to investigate the role of common, low-frequency, and rare variants influencing BMI. We sequenced TMEM18 and regions downstream of TMEM18 on chromosome 2 in 3976 individuals of European ancestry from 3 community-based cohorts (Atherosclerosis Risk in Communities, Cardiovascular Health Study, and Framingham Heart Study), including 200 adults selected for high BMI. We examined the association between BMI and variants identified in the region from nucleotide position 586 432 to 677 539 (hg18). Rare variants (minor allele frequency, <1%) were analyzed using a burden test and the sequence kernel association test. Results from the 3 cohort studies were meta-analyzed. We estimate that mean BMI is 0.43 kg/m(2) higher for each copy of the G allele of single-nucleotide polymorphism rs7596758 (minor allele frequency, 29%; P=3.46×10(-4)) using a Bonferroni threshold of P<4.6×10(-4). Analyses conditional on previous genome-wide association study single-nucleotide polymorphisms associated with BMI in the region led to attenuation of this signal and uncovered another independent (r(2)<0.2), statistically significant association, rs186019316 (P=2.11×10(-4)). Both rs186019316 and rs7596758 or proxies are located in transcription factor binding regions. No significant association with rare variants was found in either the exons of TMEM18 or the 3' genome-wide association study region. Targeted sequencing around TMEM18 identified 2 novel BMI variants with possible regulatory function.

ADAM19 and HTR4 variants and pulmonary function: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study.

The pulmonary function measures of forced expiratory volume in 1 second (FEV1) and its ratio to forced vital capacity (FVC) are used in the diagnosis and monitoring of lung diseases and predict cardiovascular mortality in the general population. Genome-wide association studies (GWASs) have identified numerous loci associated with FEV1 and FEV1/FVC, but the causal variants remain uncertain. We hypothesized that novel or rare variants poorly tagged by GWASs may explain the significant associations between FEV1/FVC and 2 genes: ADAM19 and HTR4. We sequenced ADAM19 and its promoter region along with the ≈21-kb portion of HTR4 harboring GWAS single-nucleotide polymorphisms for pulmonary function and analyzed associations with FEV1/FVC among 3983 participants of European ancestry from Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Meta-analysis of common variants in each region identified statistically significant associations (316 tests; P<1.58×10(-4)) with FEV1/FVC for 14 ADAM19 single-nucleotide polymorphisms and 24 HTR4 single-nucleotide polymorphisms. After conditioning on the sentinel GWASs hit in each gene (ADAM19 rs1422795, minor allele frequency=0.33 and HTR4 rs11168048, minor allele frequency=0.40], 1 single-nucleotide polymorphism remained statistically significant (ADAM19 rs13155908, minor allele frequency=0.12; P=1.56×10(-4)). Analysis of rare variants (minor allele frequency <1%) using sequence kernel association test did not identify associations with either region. Sequencing identified 1 common variant associated with FEV1/FVC independent of the sentinel ADAM19 GWAS hit and supports the original HTR4 GWAS findings. Rare variants do not seem to underlie GWAS associations with pulmonary function for common variants in ADAM19 and HTR4.

Mitochondrial DNA variants in obesity.

Heritability estimates for body mass index (BMI) variation are high. For mothers and their offspring higher BMI correlations have been described than for fathers. Variation(s) in the exclusively maternally inherited mitochondrial DNA (mtDNA) might contribute to this parental effect. Thirty-two to 40 mtDNA single nucleotide polymorphisms (SNPs) were available from genome-wide association study SNP arrays (Affymetrix 6.0). For discovery, we analyzed association in a case-control (CC) sample of 1,158 extremely obese children and adolescents and 435 lean adult controls. For independent confirmation, 7,014 population-based adults were analyzed as CC sample of n = 1,697 obese cases (BMI≥30 kg/m2) and n = 2,373 normal weight and lean controls (BMI<25 kg/m2). SNPs were analyzed as single SNPs and haplogroups determined by HaploGrep. Fisher's two-sided exact test was used for association testing. Moreover, the D-loop was re-sequenced (Sanger) in 192 extremely obese children and adolescents and 192 lean adult controls. Association testing of detected variants was performed using Fisher's two-sided exact test. For discovery, nominal association with obesity was found for the frequent allele G of m.8994G/A (rs28358887, p = 0.002) located in ATP6. Haplogroup W was nominally overrepresented in the controls (p = 0.039). These findings could not be confirmed independently. For two of the 252 identified D-loop variants nominal association was detected (m.16292C/T, p = 0.007, m.16189T/C, p = 0.048). Only eight controls carried the m.16292T allele, five of whom belonged to haplogroup W that was initially enriched among these controls. m.16189T/C might create an uninterrupted poly-C tract located near a regulatory element involved in replication of mtDNA. Though follow-up of some D-loop variants still is conceivable, our hypothesis of a contribution of variation in the exclusively maternally inherited mtDNA to the observed larger correlations for BMI between mothers and their offspring could not be substantiated by the findings of the present study.

Identification of candidate susceptibility genes for colorectal cancer through eQTL analysis

In this study, we aim to identify the genes responsible for colorectal cancer risk behind the loci identified in genome-wide association studies (GWAS). These genes may be candidate targets for developing new strategies for prevention or therapy. We analyzed the association of genotypes for 26 GWAS single nucleotide polymorphisms (SNPs) with the expression of genes within a 2 Mb region (cis-eQTLs). Affymetrix Human Genome U219 expression arrays were used to assess gene expression in two series of samples, one of healthy colonic mucosa (n = 47) and other of normal mucosa adjacent to colon cancer (n = 97, total 144). Paired tumor tissues (n = 97) were also analyzed but did not provide additional findings. Partial Pearson correlation (r), adjusted for sample type, was used for the analysis. We have found Bonferroni-significant cis-eQTLs in three loci: rs3802842 in 11q23.1 associated to C11orf53, COLCA1 (C11orf92) and COLCA2 (C11orf93; r = 0.60); rs7136702 in 12q13.12 associated to DIP2B (r = 0.63) and rs5934683 in Xp22.3 associated to SHROOM2 and GPR143 (r = 0.47). For loci in chromosomes 11 and 12, we have found other SNPs in linkage disequilibrium that are more strongly associated with the expression of the identified genes and are better functional candidates: rs7130173 for 11q23.1 (r = 0.66) and rs61927768 for 12q13.12 (r = 0.86). These SNPs are located in DNA regions that may harbor enhancers or transcription factor binding sites. The analysis of trans-eQTLs has identified additional genes in these loci that may have common regulatory mechanisms as shown by the analysis of protein-protein interaction networks.