Search

BackgroundCommon chromosome 9p21 single nucleotide polymorphisms (SNPs) increase coronary heart disease risk, independent of traditional lipid risk factors. However, lipids comprise large numbers of structurally related molecules not measured in traditional risk measurements, and many have inflammatory bioactivities. Here, we applied lipidomic and genomic approaches to 3 model systems to characterize lipid metabolic changes in common Chr9p21 SNPs, which confer ≈30% elevated coronary heart disease risk associated with altered expression of ANRIL, a long ncRNA.MethodsUntargeted and targeted lipidomics was applied to plasma from NPHSII (Northwick Park Heart Study II) homozygotes for AA or GG in rs10757274, followed by correlation and network analysis. To identify candidate genes, transcriptomic data from shRNA downregulation of ANRIL in HEK-293 cells was mined. Transcriptional data from vascular smooth muscle cells differentiated from induced pluripotent stem cells of individuals with/without Chr9p21 risk, nonrisk alleles, and corresponding knockout isogenic lines were next examined. Last, an in-silico analysis of miRNAs was conducted to identify how ANRIL might control lysoPL (lysophosphospholipid)/lysoPA (lysophosphatidic acid) genes.ResultsElevated risk GG correlated with reduced lysoPLs, lysoPA, and ATX (autotaxin). Five other risk SNPs did not show this phenotype. LysoPL-lysoPA interconversion was uncoupled from ATX in GG plasma, suggesting metabolic dysregulation. Significantly altered expression of several lysoPL/lysoPA metabolizing enzymes was found in HEK cells lacking ANRIL. In the vascular smooth muscle cells data set, the presence of risk alleles associated with altered expression of several lysoPL/lysoPA enzymes. Deletion of the risk locus reversed the expression of several lysoPL/lysoPA genes to nonrisk haplotype levels. Genes that were altered across both cell data sets were DGKA, MBOAT2, PLPP1, and LPL. The in-silico analysis identified 4 ANRIL-regulate

BackgroundCommon chromosome 9p21 single nucleotide polymorphisms (SNPs) increase coronary heart disease risk, independent of traditional lipid risk factors. However, lipids comprise large numbers of structurally related molecules not measured in traditional risk measurements, and many have inflammatory bioactivities. Here, we applied lipidomic and genomic approaches to 3 model systems to characterize lipid metabolic changes in common Chr9p21 SNPs, which confer ≈30% elevated coronary heart disease risk associated with altered expression of ANRIL, a long ncRNA.MethodsUntargeted and targeted lipidomics was applied to plasma from NPHSII (Northwick Park Heart Study II) homozygotes for AA or GG in rs10757274, followed by correlation and network analysis. To identify candidate genes, transcriptomic data from shRNA downregulation of ANRIL in HEK-293 cells was mined. Transcriptional data from vascular smooth muscle cells differentiated from induced pluripotent stem cells of individuals with/without Chr9p21 risk, nonrisk alleles, and corresponding knockout isogenic lines were next examined. Last, an in-silico analysis of miRNAs was conducted to identify how ANRIL might control lysoPL (lysophosphospholipid)/lysoPA (lysophosphatidic acid) genes.ResultsElevated risk GG correlated with reduced lysoPLs, lysoPA, and ATX (autotaxin). Five other risk SNPs did not show this phenotype. LysoPL-lysoPA interconversion was uncoupled from ATX in GG plasma, suggesting metabolic dysregulation. Significantly altered expression of several lysoPL/lysoPA metabolizing enzymes was found in HEK cells lacking ANRIL. In the vascular smooth muscle cells data set, the presence of risk alleles associated with altered expression of several lysoPL/lysoPA enzymes. Deletion of the risk locus reversed the expression of several lysoPL/lysoPA genes to nonrisk haplotype levels. Genes that were altered across both cell data sets were DGKA, MBOAT2, PLPP1, and LPL. The in-silico analysis identified 4 ANRIL-regulate

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

The 9p21.3 cardiovascular disease locus is the most influential common genetic risk factor for coronary artery disease (CAD), accounting for ∼10%-15% of disease in non-African populations. The ∼60 kb risk haplotype is human-specific and lacks coding genes, hindering efforts to decipher its function. Here, we produce induced pluripotent stem cells (iPSCs) from risk and non-risk individuals, delete each haplotype using genome editing, and generate vascular smooth muscle cells (VSMCs). Risk VSMCs exhibit globally altered transcriptional networks that intersect with previously identified CAD risk genes and pathways, concomitant with aberrant adhesion, contraction, and proliferation. Unexpectedly, deleting the risk haplotype rescues VSMC stability, while expressing the 9p21.3-associated long non-coding RNA ANRIL induces risk phenotypes in non-risk VSMCs. This study shows that the risk haplotype selectively predisposes VSMCs to adopt a cell state associated with CAD phenotypes, defines new VSMC-based networks of CAD risk genes, and establishes haplotype-edited iPSCs as powerful tools for functionally annotating the human genome.

The 9p21.3 cardiovascular disease locus is the most influential common genetic risk factor for coronary artery disease (CAD), accounting for ∼10%-15% of disease in non-African populations. The ∼60 kb risk haplotype is human-specific and lacks coding genes, hindering efforts to decipher its function. Here, we produce induced pluripotent stem cells (iPSCs) from risk and non-risk individuals, delete each haplotype using genome editing, and generate vascular smooth muscle cells (VSMCs). Risk VSMCs exhibit globally altered transcriptional networks that intersect with previously identified CAD risk genes and pathways, concomitant with aberrant adhesion, contraction, and proliferation. Unexpectedly, deleting the risk haplotype rescues VSMC stability, while expressing the 9p21.3-associated long non-coding RNA ANRIL induces risk phenotypes in non-risk VSMCs. This study shows that the risk haplotype selectively predisposes VSMCs to adopt a cell state associated with CAD phenotypes, defines new VSMC-based networks of CAD risk genes, and establishes haplotype-edited iPSCs as powerful tools for functionally annotating the human genome.

Background & aimsGenetic factors are believed to affect risk for irritable bowel syndrome (IBS), but there have been no sufficiently powered and adequately sized studies. To identify DNA variants associated with IBS risk, we performed a genome-wide association study (GWAS) of the large UK Biobank population-based cohort, which includes genotype and health data from 500,000 participants.MethodsWe studied 7,287,191 high-quality single nucleotide polymorphisms in individuals who self-reported a doctor's diagnosis of IBS (cases; n = 9576) compared to the remainder of the cohort (controls; n = 336,499) (mean age of study subjects, 40-69 years). Genome-wide significant findings were further investigated in 2045 patients with IBS from tertiary centers and 7955 population controls from Europe and the United States, and a small general population sample from Sweden (n = 249). Functional annotation of GWAS results was carried out by integrating data from multiple biorepositories to obtain biological insights from the observed associations.ResultsWe identified a genome-wide significant association on chromosome 9q31.2 (single nucleotide polymorphism rs10512344; P = 3.57 × 10-8) in a region previously linked to age at menarche, and 13 additional loci of suggestive significance (P < 5.0×10-6). Sex-stratified analyses revealed that the variants at 9q31.2 affect risk of IBS in women only (P = 4.29 × 10-10 in UK Biobank) and also associate with constipation-predominant IBS in women (P = .015 in the tertiary cohort) and harder stools in women (P = .0012 in the population-based sample). Functional annotation of the 9q31.2 locus identified 8 candidate genes, including the elongator complex protein 1 gene (ELP1 or IKBKAP), which is mutated in patients with familial dysautonomia.ConclusionsIn a sufficiently powered GWAS of IBS, we associated variants at the locus 9q31.2 with risk of IBS in women. This observation may provide addit

In patients with Philadelphia chromosome-positive (Ph+) acute lymphoblastic leukemia (ALL) treated with chemotherapy plus a tyrosine kinase inhibitor (TKI), the prognostic impact of additional chromosomal abnormalities (ACAs) is not well-established. We evaluated the prognostic impact of individual ACAs in 152 patients with Ph+ ALL receiving first-line intensive chemotherapy plus either imatinib (n = 36), dasatinib (n = 74), or ponatinib (n = 42). ACAs were identified in 118 patients (78%). Compared to outcomes of patients without ACAs, ACAs were not associated with differences in either relapse-free survival (RFS; P = 0.42) or overall survival (OS; P = 0.51). When individual ACAs were evaluated, +der(22)t(9;22) and/or -9/9p in the absence of high hyperdiploidy (HeH) was present in 16% of patients and constituted a poor-risk ACA group. Patients with one or more poor-risk ACAs in the absence of HeH had significantly shorter RFS (5-year RFS rate 33% versus 59%, P = 0.01) and OS (5-year OS rate 24% versus 63%, P = 0.003). Poor-risk ACAs were prognostic in patients who received imatinib and dasatinib but not in those who received ponatinib. By multivariate analysis, this poor-risk ACA group was independently associated with worse RFS (HR 2.03 [95% CI 1.08-3.30], P = 0.03) and OS (HR 2.02 [95% CI 1.10-3.71], P = 0.02). Patients with Ph+ ALL who have +der(22)t(9;22) and/or -9/9p in the absence of HeH have relatively poor outcomes when treated with chemotherapy plus a TKI.

In patients with Philadelphia chromosome-positive (Ph+) acute lymphoblastic leukemia (ALL) treated with chemotherapy plus a tyrosine kinase inhibitor (TKI), the prognostic impact of additional chromosomal abnormalities (ACAs) is not well-established. We evaluated the prognostic impact of individual ACAs in 152 patients with Ph+ ALL receiving first-line intensive chemotherapy plus either imatinib (n = 36), dasatinib (n = 74), or ponatinib (n = 42). ACAs were identified in 118 patients (78%). Compared to outcomes of patients without ACAs, ACAs were not associated with differences in either relapse-free survival (RFS; P = 0.42) or overall survival (OS; P = 0.51). When individual ACAs were evaluated, +der(22)t(9;22) and/or -9/9p in the absence of high hyperdiploidy (HeH) was present in 16% of patients and constituted a poor-risk ACA group. Patients with one or more poor-risk ACAs in the absence of HeH had significantly shorter RFS (5-year RFS rate 33% versus 59%, P = 0.01) and OS (5-year OS rate 24% versus 63%, P = 0.003). Poor-risk ACAs were prognostic in patients who received imatinib and dasatinib but not in those who received ponatinib. By multivariate analysis, this poor-risk ACA group was independently associated with worse RFS (HR 2.03 [95% CI 1.08-3.30], P = 0.03) and OS (HR 2.02 [95% CI 1.10-3.71], P = 0.02). Patients with Ph+ ALL who have +der(22)t(9;22) and/or -9/9p in the absence of HeH have relatively poor outcomes when treated with chemotherapy plus a TKI.

Population-based genome wide association studies have identified a locus at 9p22.2 associated with ovarian cancer risk, which also modifies ovarian cancer risk in BRCA1 and BRCA2 mutation carriers. We conducted fine-scale mapping at 9p22.2 to identify potential causal variants in BRCA1 and BRCA2 mutation carriers. Genotype data were available for 15,252 (2,462 ovarian cancer cases) BRCA1 and 8,211 (631 ovarian cancer cases) BRCA2 mutation carriers. Following genotype imputation, ovarian cancer associations were assessed for 4,873 and 5,020 SNPs in BRCA1 and BRCA 2 mutation carriers respectively, within a retrospective cohort analytical framework. In BRCA1 mutation carriers one set of eight correlated candidate causal variants for ovarian cancer risk modification was identified (top SNP rs10124837, HR: 0.73, 95%CI: 0.68 to 0.79, p-value 2× 10-16). These variants were located up to 20 kb upstream of BNC2. In BRCA2 mutation carriers one region, up to 45 kb upstream of BNC2, and containing 100 correlated SNPs was identified as candidate causal (top SNP rs62543585, HR: 0.69, 95%CI: 0.59 to 0.80, p-value 1.0 × 10-6). The candidate causal in BRCA1 mutation carriers did not include the strongest associated variant at this locus in the general population. In sum, we identified a set of candidate causal variants in a region that encompasses the BNC2 transcription start site. The ovarian cancer association at 9p22.2 may be mediated by different variants in BRCA1 mutation carriers and in the general population. Thus, potentially different mechanisms may underlie ovarian cancer risk for mutation carriers and the general population.

Paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) is the most common cancer of childhood, yet little is known about BCP-ALL predisposition. In this study, in 2,187 cases of European ancestry and 5,543 controls, we discover and replicate a locus indexed by rs77728904 at 9p21.3 associated with BCP-ALL susceptibility (Pcombined=3.32 × 10(-15), OR=1.72) and independent from rs3731217, the previously reported ALL-associated variant in this region. Of correlated SNPs tagged by this locus, only rs662463 is significant in African Americans, suggesting it is a plausible causative variant. Functional analysis shows that rs662463 is a cis-eQTL for CDKN2B, with the risk allele associated with lower expression, and suggests that rs662463 influences BCP-ALL risk by regulating CDKN2B expression through CEBPB signalling. Functional analysis of rs3731217 suggests it is associated with BCP-ALL by acting within a splicing regulatory element determining CDKN2A exon 3 usage (P=0.01). These findings provide new insights into the critical role of the CDKN2 locus in BCP-ALL aetiology.

Paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) is the most common cancer of childhood, yet little is known about BCP-ALL predisposition. In this study, in 2,187 cases of European ancestry and 5,543 controls, we discover and replicate a locus indexed by rs77728904 at 9p21.3 associated with BCP-ALL susceptibility (Pcombined=3.32 × 10(-15), OR=1.72) and independent from rs3731217, the previously reported ALL-associated variant in this region. Of correlated SNPs tagged by this locus, only rs662463 is significant in African Americans, suggesting it is a plausible causative variant. Functional analysis shows that rs662463 is a cis-eQTL for CDKN2B, with the risk allele associated with lower expression, and suggests that rs662463 influences BCP-ALL risk by regulating CDKN2B expression through CEBPB signalling. Functional analysis of rs3731217 suggests it is associated with BCP-ALL by acting within a splicing regulatory element determining CDKN2A exon 3 usage (P=0.01). These findings provide new insights into the critical role of the CDKN2 locus in BCP-ALL aetiology.

Population-based genome wide association studies have identified a locus at 9p22.2 associated with ovarian cancer risk, which also modifies ovarian cancer risk in BRCA1 and BRCA2 mutation carriers. We conducted fine-scale mapping at 9p22.2 to identify potential causal variants in BRCA1 and BRCA2 mutation carriers. Genotype data were available for 15,252 (2,462 ovarian cancer cases) BRCA1 and 8,211 (631 ovarian cancer cases) BRCA2 mutation carriers. Following genotype imputation, ovarian cancer associations were assessed for 4,873 and 5,020 SNPs in BRCA1 and BRCA 2 mutation carriers respectively, within a retrospective cohort analytical framework. In BRCA1 mutation carriers one set of eight correlated candidate causal variants for ovarian cancer risk modification was identified (top SNP rs10124837, HR: 0.73, 95%CI: 0.68 to 0.79, p-value 2× 10-16). These variants were located up to 20 kb upstream of BNC2. In BRCA2 mutation carriers one region, up to 45 kb upstream of BNC2, and containing 100 correlated SNPs was identified as candidate causal (top SNP rs62543585, HR: 0.69, 95%CI: 0.59 to 0.80, p-value 1.0 × 10-6). The candidate causal in BRCA1 mutation carriers did not include the strongest associated variant at this locus in the general population. In sum, we identified a set of candidate causal variants in a region that encompasses the BNC2 transcription start site. The ovarian cancer association at 9p22.2 may be mediated by different variants in BRCA1 mutation carriers and in the general population. Thus, potentially different mechanisms may underlie ovarian cancer risk for mutation carriers and the general population.

The presence of AML1-ETO (RUNX1-CBF2T1), a fusion oncoprotein resulting from a t(8;21) chromosomal translocation, has been implicated as a necessary but insufficient event in the development of a subset of acute myeloid leukemias (AML). While AML1-ETO prolongs survival and inhibits differentiation of hematopoietic stem cells (HSC), other contributory events are needed for cell proliferation and leukemogenesis. We have postulated that specific tumor suppressor genes keep the leukemic potential of AML1-ETO in check. In studying del(9q), one of the most common concomitant chromosomal abnormalities with t(8;21), we identified the loss of an apparent tumor suppressor, TLE4, that appears to cooperate with AML1-ETO to confer a leukemic phenotype. This study sought to identify the molecular basis of this cooperation. We show that the loss of TLE4 confers proliferative advantage to leukemic cells, simultaneous with an upregulation of a pro- inflammatory signature mediated through aberrant increases in Wnt signaling activity. We further demonstrate that inhibition of cyclooxygenase (COX) activity partly reverses the pro-leukemic phenotype due to TLE4 knockdown, pointing towards a novel therapeutic approach for myeloid leukemia.

Heritable genetic variants can significantly affect the lifetime risk of developing cancer, including polyposis and colorectal cancer (CRC). Variants in genes currently known to be associated with a high risk for polyposis or CRC, however, explain only a limited number of hereditary cases. The identification of additional genetic causes is, therefore, crucial to improve CRC prevention, detection and treatment. We have performed genome-wide and targeted DNA copy number profiling and resequencing in early-onset and familial polyposis/CRC patients, and show that deletions affecting the open reading frame of the tumour suppressor gene FOCAD are recurrent and significantly enriched in CRC patients compared with unaffected controls. All patients carrying FOCAD deletions exhibited a personal or family history of polyposis. RNA in situ hybridization revealed FOCAD expression in epithelial cells in the colonic crypt, the site of tumour initiation, as well as in colonic tumours and organoids. Our data suggest that monoallelic germline deletions in the tumour suppressor gene FOCAD underlie moderate genetic predisposition to the development of polyposis and CRC.

Obsessive-compulsive disorder (OCD) is a psychiatric condition characterized by intrusive thoughts and urges and repetitive, intentional behaviors that cause significant distress and impair functioning. The OCD Collaborative Genetics Association Study (OCGAS) is comprised of comprehensively assessed OCD patients with an early age of OCD onset. After application of a stringent quality control protocol, a total of 1065 families (containing 1406 patients with OCD), combined with population-based samples (resulting in a total sample of 5061 individuals), were studied. An integrative analyses pipeline was utilized, involving association testing at single-nucleotide polymorphism (SNP) and gene levels (via a hybrid approach that allowed for combined analyses of the family- and population-based data). The smallest P-value was observed for a marker on chromosome 9 (near PTPRD, P=4.13 × 10(-)(7)). Pre-synaptic PTPRD promotes the differentiation of glutamatergic synapses and interacts with SLITRK3. Together, both proteins selectively regulate the development of inhibitory GABAergic synapses. Although no SNPs were identified as associated with OCD at genome-wide significance level, follow-up analyses of genome-wide association study (GWAS) signals from a previously published OCD study identified significant enrichment (P=0.0176). Secondary analyses of high-confidence interaction partners of DLGAP1 and GRIK2 (both showing evidence for association in our follow-up and the original GWAS study) revealed a trend of association (P=0.075) for a set of genes such as NEUROD6, SV2A, GRIA4, SLC1A2 and PTPRD. Analyses at the gene level revealed association of IQCK and C16orf88 (both P<1 × 10(-)(6), experiment-wide significant), as well as OFCC1 (P=6.29 × 10(-)(5)). The suggestive findings in this study await replication in larger samples.

Heritable genetic variants can significantly affect the lifetime risk of developing cancer, including polyposis and colorectal cancer (CRC). Variants in genes currently known to be associated with a high risk for polyposis or CRC, however, explain only a limited number of hereditary cases. The identification of additional genetic causes is, therefore, crucial to improve CRC prevention, detection and treatment. We have performed genome-wide and targeted DNA copy number profiling and resequencing in early-onset and familial polyposis/CRC patients, and show that deletions affecting the open reading frame of the tumour suppressor gene FOCAD are recurrent and significantly enriched in CRC patients compared with unaffected controls. All patients carrying FOCAD deletions exhibited a personal or family history of polyposis. RNA in situ hybridization revealed FOCAD expression in epithelial cells in the colonic crypt, the site of tumour initiation, as well as in colonic tumours and organoids. Our data suggest that monoallelic germline deletions in the tumour suppressor gene FOCAD underlie moderate genetic predisposition to the development of polyposis and CRC.

Obsessive-compulsive disorder (OCD) is a psychiatric condition characterized by intrusive thoughts and urges and repetitive, intentional behaviors that cause significant distress and impair functioning. The OCD Collaborative Genetics Association Study (OCGAS) is comprised of comprehensively assessed OCD patients with an early age of OCD onset. After application of a stringent quality control protocol, a total of 1065 families (containing 1406 patients with OCD), combined with population-based samples (resulting in a total sample of 5061 individuals), were studied. An integrative analyses pipeline was utilized, involving association testing at single-nucleotide polymorphism (SNP) and gene levels (via a hybrid approach that allowed for combined analyses of the family- and population-based data). The smallest P-value was observed for a marker on chromosome 9 (near PTPRD, P=4.13 × 10(-)(7)). Pre-synaptic PTPRD promotes the differentiation of glutamatergic synapses and interacts with SLITRK3. Together, both proteins selectively regulate the development of inhibitory GABAergic synapses. Although no SNPs were identified as associated with OCD at genome-wide significance level, follow-up analyses of genome-wide association study (GWAS) signals from a previously published OCD study identified significant enrichment (P=0.0176). Secondary analyses of high-confidence interaction partners of DLGAP1 and GRIK2 (both showing evidence for association in our follow-up and the original GWAS study) revealed a trend of association (P=0.075) for a set of genes such as NEUROD6, SV2A, GRIA4, SLC1A2 and PTPRD. Analyses at the gene level revealed association of IQCK and C16orf88 (both P<1 × 10(-)(6), experiment-wide significant), as well as OFCC1 (P=6.29 × 10(-)(5)). The suggestive findings in this study await replication in larger samples.

The concept that frontotemporal dementia (FTD) is a purely cortical dementia has largely been refuted by the recognition of its close association with motor neuron disease, and the identification of transactive response DNA-binding protein 43 (TDP-43) as a major pathological substrate underlying both diseases. Genetic findings have transformed this field and revealed connections between disorders that were previous thought clinically unrelated. The discovery that the C9ORF72 locus is responsible for the majority of hereditary FTD, amyotrophic lateral sclerosis (ALS), and FTD-ALS cases and the understanding that repeat-containing RNA plays a crucial role in pathogenesis of both disorders has paved the way for the development of potential biomarkers and therapeutic targets for these devastating diseases. In this review, we summarize the historical aspects leading up to our current understanding of the genetic, clinical, and neuropathological overlap between FTD and ALS, and include brief discussions on chronic traumatic encephalopathy (CTE), given its association with TDP-43 pathology, its associated increased dementia risk, and reports of ALS in CTE patients. In addition, we describe other genetic associations between dementia and neuromuscular disease, such as inclusion body myositis with Paget's disease and FTD.

Whole-genome oligonucleotide single-nucleotide polymorphism (oligo-SNP) arrays enable simultaneous interrogation of copy number variations (CNVs), copy neutral regions of homozygosity (ROH) and uniparental disomy (UPD). Structural variation in the human genome contributes significantly to genetic variation, and often has deleterious effects leading to disease causation. Co-occurrence of CNV and regions of allelic homozygosity in tandem involving the same chromosomal arm are extremely rare. Replication-based mechanisms such as microhomology-mediated break-induced replication (MMBIR) are recent models predicted to induce structural rearrangements and gene dosage aberrations; however, supportive evidence in humans for one-ended DNA break repair coupled with MMBIR giving rise to interstitial copy number gains and distal loss of heterozygosity has not been documented. We report on the identification and characterization of two cases with interstitial triplication followed by uniparental isodisomy (isoUPD) for remainder of the chromosomal arm. Case 1 has a triplication at 9q21.11-q21.33 and segmental paternal isoUPD for 9q21.33-qter, and presented with citrullinemia with a homozygous mutation in the argininosuccinate synthetase gene (ASS1 at 9q34.1). Case 2 has a triplication at 22q12.1-q12.2 and segmental maternal isoUPD 22q12.2-qter, and presented with hearing loss, mild dysmorphic features and bilateral iris coloboma. Interstitial triplication coupled with distal segmental isoUPD is a novel finding that provides human evidence for one-ended DNA break and replication-mediated repair. Both copy number gains and isoUPD may contribute to the phenotype. Significantly, these cases represent the first detailed genomic analysis that provides support for a MMBIR mechanism inducing copy number gains and segmental isoUPD in tandem.

Array comparative genomic hybridization (aCGH) is now commonly used to identify copy number changes in individuals with developmental delay, intellectual disabilities, autism spectrum disorders, and/or multiple congenital anomalies. We report on an infant with multiple congenital anomalies and a novel 2.6 Mb interstitial deletion within 9q21.32q21.33 detected by aCGH. Her clinical presentation included dysmorphic craniofacial features, cleft palate, atrial septal defect, bicornuate uterus, bilateral hip dislocation, hypotonia, and recurrent pneumonia. Parental aCGH studies were negative for copy loss in this region. To our knowledge, no similar deletions have been reported in available databases or published literature. This deletion encompasses 12 genes, and prediction algorithms as well as experimental data suggest that a subset is likely to be haploinsufficient. Included are a neurotrophin receptor (NKG2D), a gene implicated in cilia function (KIF27), an adaptor protein important for ubiquitin-dependent protein quality control (UBQLN1), a gene important for transcription and signaling (HNRNPK), and a gene involved in maintaining genomic stability (RMI1). Identifying additional patients with similar copy losses and further study of these genes will contribute to a better understanding of the pathophysiology of multiple congenital anomalies.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10(-11)). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10(-10)) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.

IntroductionWe examined whether the variant at chromosome 9p21, rs4977574, was associated with long-term cardiovascular mortality in Han Chinese patients with coronary artery disease (CAD).MethodologySubjects who underwent coronary angiography for chest pain were consecutively enrolled. Fasting blood samples were collected for laboratory and genotype assessments. The information was correlated with data collected from the national death database.ResultsThere were 925 cases with CAD and 634 without CAD enrolled in the present study. The G allele conferred a significant increase in risk of CAD (odds ratio = 1.47, P = 0.003 in the dominant model; odds ratio = 1.36, P = 0.018 in the recessive model). During a median of 11 years (inter-quartile range between 5.2 and 12.5 years) of follow-up, neither the total nor the cardiovascular mortality was different among CAD subjects with different genotypes. Using Cox regression analysis, genotypes of rs4977574 still failed to predict cardiovascular mortality (hazard ratio = 1.25, P = 0.138 in the dominant model; hazard ratio = 1.05, P = 0.729 in the recessive model).ConclusionsThe rs4977574 at chromosome 9p21 is associated with presence of CAD in Han Chinese. However, rs4977574 could not predict cardiovascular mortality in these CAD subjects during the eleven-year period of the study.

Array comparative genomic hybridization (aCGH) is now commonly used to identify copy number changes in individuals with developmental delay, intellectual disabilities, autism spectrum disorders, and/or multiple congenital anomalies. We report on an infant with multiple congenital anomalies and a novel 2.6 Mb interstitial deletion within 9q21.32q21.33 detected by aCGH. Her clinical presentation included dysmorphic craniofacial features, cleft palate, atrial septal defect, bicornuate uterus, bilateral hip dislocation, hypotonia, and recurrent pneumonia. Parental aCGH studies were negative for copy loss in this region. To our knowledge, no similar deletions have been reported in available databases or published literature. This deletion encompasses 12 genes, and prediction algorithms as well as experimental data suggest that a subset is likely to be haploinsufficient. Included are a neurotrophin receptor (NKG2D), a gene implicated in cilia function (KIF27), an adaptor protein important for ubiquitin-dependent protein quality control (UBQLN1), a gene important for transcription and signaling (HNRNPK), and a gene involved in maintaining genomic stability (RMI1). Identifying additional patients with similar copy losses and further study of these genes will contribute to a better understanding of the pathophysiology of multiple congenital anomalies.

IntroductionWe examined whether the variant at chromosome 9p21, rs4977574, was associated with long-term cardiovascular mortality in Han Chinese patients with coronary artery disease (CAD).MethodologySubjects who underwent coronary angiography for chest pain were consecutively enrolled. Fasting blood samples were collected for laboratory and genotype assessments. The information was correlated with data collected from the national death database.ResultsThere were 925 cases with CAD and 634 without CAD enrolled in the present study. The G allele conferred a significant increase in risk of CAD (odds ratio = 1.47, P = 0.003 in the dominant model; odds ratio = 1.36, P = 0.018 in the recessive model). During a median of 11 years (inter-quartile range between 5.2 and 12.5 years) of follow-up, neither the total nor the cardiovascular mortality was different among CAD subjects with different genotypes. Using Cox regression analysis, genotypes of rs4977574 still failed to predict cardiovascular mortality (hazard ratio = 1.25, P = 0.138 in the dominant model; hazard ratio = 1.05, P = 0.729 in the recessive model).ConclusionsThe rs4977574 at chromosome 9p21 is associated with presence of CAD in Han Chinese. However, rs4977574 could not predict cardiovascular mortality in these CAD subjects during the eleven-year period of the study.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10(-11)). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10(-10)) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.

ObjectiveTo investigate whether previously reported 9p21.3 single nucleotide polymorphisms (SNPs) are associated with risk of brain arteriovenous malformations (BAVM), which often have accompanying arterial aneurysms. Common variants in the 9p21.3 locus have been reported to be associated with multiple cardiovascular phenotypes, including coronary artery disease and intracranial aneurysms (rs10757278 and rs1333040).MethodsWe used data from 338 BAVM cases participating in the University of California, San Francisco (UCSF)-Kaiser Brain AVM Study Project and 504 healthy controls to evaluate genotypes for seven common SNPs (minor allele frequency>0.05) that were imputed using 1000 Genomes Phase 1 European data (R(2)>0.87). Association with BAVM was tested using logistic regression adjusting for age, sex and the top three principal components of ancestry. Subgroup analysis included 205 BAVM cases with aneurysm data: 74 BAVM with aneurysm versus 504 controls and 131 BAVM without aneurysm versus 504 controls.ResultsWe observed suggestive association with BAVM and rs10757278-G (OR=1.23, 95% CI 0.99 to 1.53, p=0.064) and rs1333040-T (OR=1.27, 95% CI 1.01 to 1.58, p=0.04). For rs10757278-G, the association was stronger in BAVM cases with aneurysm (OR=1.52, 95% CI 1.03 to 2.22, p=0.032) than in BAVM without aneurysm (OR=0.98, 95% CI 0.72 to 1.34, p=0.91). Similar patterns of effects were observed for rs1333040 and for other SNPs in linkage disequilibrium (r(2)>0.8) with rs10757278.ConclusionsCommon 9p21.3 variants showed similar effect sizes for association with BAVM as previously reported for aneurysmal disease. The association with BAVM appears to be explained by known associations with aneurysms, suggesting that BAVM-associated aneurysms share similar vascular pathology mechanisms with other aneurysm types.

Genome-wide association studies have identified more than 70 common variants that are associated with breast cancer risk. Most of these variants map to non-protein-coding regions and several map to gene deserts, regions of several hundred kilobases lacking protein-coding genes. We hypothesized that gene deserts harbor long-range regulatory elements that can physically interact with target genes to influence their expression. To test this, we developed Capture Hi-C (CHi-C), which, by incorporating a sequence capture step into a Hi-C protocol, allows high-resolution analysis of targeted regions of the genome. We used CHi-C to investigate long-range interactions at three breast cancer gene deserts mapping to 2q35, 8q24.21, and 9q31.2. We identified interaction peaks between putative regulatory elements ("bait fragments") within the captured regions and "targets" that included both protein-coding genes and long noncoding (lnc) RNAs over distances of 6.6 kb to 2.6 Mb. Target protein-coding genes were IGFBP5, KLF4, NSMCE2, and MYC; and target lncRNAs included DIRC3, PVT1, and CCDC26. For one gene desert, we were able to define two SNPs (rs12613955 and rs4442975) that were highly correlated with the published risk variant and that mapped within the bait end of an interaction peak. In vivo ChIP-qPCR data show that one of these, rs4442975, affects the binding of FOXA1 and implicate this SNP as a putative functional variant.

Rare copy number variants (CNVs) disrupting ASTN2 or both ASTN2 and TRIM32 have been reported at 9q33.1 by genome-wide studies in a few individuals with neurodevelopmental disorders (NDDs). The vertebrate-specific astrotactins, ASTN2 and its paralog ASTN1, have key roles in glial-guided neuronal migration during brain development. To determine the prevalence of astrotactin mutations and delineate their associated phenotypic spectrum, we screened ASTN2/TRIM32 and ASTN1 (1q25.2) for exonic CNVs in clinical microarray data from 89 985 individuals across 10 sites, including 64 114 NDD subjects. In this clinical dataset, we identified 46 deletions and 12 duplications affecting ASTN2. Deletions of ASTN1 were much rarer. Deletions near the 3' terminus of ASTN2, which would disrupt all transcript isoforms (a subset of these deletions also included TRIM32), were significantly enriched in the NDD subjects (P = 0.002) compared with 44 085 population-based controls. Frequent phenotypes observed in individuals with such deletions include autism spectrum disorder (ASD), attention deficit hyperactivity disorder (ADHD), speech delay, anxiety and obsessive compulsive disorder (OCD). The 3'-terminal ASTN2 deletions were significantly enriched compared with controls in males with NDDs, but not in females. Upon quantifying ASTN2 human brain RNA, we observed shorter isoforms expressed from an alternative transcription start site of recent evolutionary origin near the 3' end. Spatiotemporal expression profiling in the human brain revealed consistently high ASTN1 expression while ASTN2 expression peaked in the early embryonic neocortex and postnatal cerebellar cortex. Our findings shed new light on the role of the astrotactins in psychopathology and their interplay in human neurodevelopment.

Tourettes syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10(-8)); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10(-6)). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10(-7) for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

Tourette's syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10(-8)); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10(-6)). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10(-7) for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

Tourette's syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10(-8)); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10(-6)). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10(-7) for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

ObjectiveHow 9p21 variation affects risk of cardiovascular disease is unclear, so we assessed whether 9p21 variants are associated with arterial elasticity or retinal microvascular findings.MethodsIn the prospective Multi-Ethnic Study of Atherosclerosis (MESA) we assessed 378 SNPs in the 9p21 locus. Within four ethnic groups, we used an additive genetic model to relate the 9p21 SNPs to five vascular phenotypes: small and large elasticity derived from radial diastolic pulse contour analysis; Young's elastic modulus from carotid artery ultrasound measurements; and the diameter of the central retinal arteries and veins.ResultsIn neither ethnic-specific nor pooled data was there any statistically significant association between any of the 9p21 SNPs and any of the five vascular phenotypes.ConclusionOur study does not support an association of 9p21 variation with arterial elasticity or retinal microvascular abnormalities.

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Five single nucleotide polymorphisms (SNPs) associated with thyroid cancer (TC) risk have been reported: rs2910164 (5q24); rs6983267 (8q24); rs965513 and rs1867277 (9q22); and rs944289 (14q13). Most of these associations have not been replicated in independent populations and the combined effects of the SNPs on risk have not been examined. This study genotyped the five TC SNPs in 781 patients recruited through the TCUKIN study. Genotype data from 6122 controls were obtained from the CORGI and Wellcome Trust Case-Control Consortium studies. Significant associations were detected between TC and rs965513A (p=6.35×10(-34)), rs1867277A (p=5.90×10(-24)), rs944289T (p=6.95×10(-7)), and rs6983267G (p=0.016). rs6983267 was most strongly associated under a recessive model (P(GG vs GT + TT)=0.004), in contrast to the association of this SNP with other cancer types. However, no evidence was found of an association between rs2910164 and disease under any risk model (p>0.7). The rs1867277 association remained significant (p=0.008) after accounting for genotypes at the nearby rs965513 (p=2.3×10(-13)) and these SNPs did not tag a single high risk haplotype. The four validated TC SNPs accounted for a relatively large proportion (∼11%) of the sibling relative risk of TC, principally owing to the large effect size of rs965513 (OR 1.74).

Diaphyseal medullary stenosis with malignant fibrous histiocytoma (DMS-MFH) is an autosomal-dominant syndrome characterized by bone dysplasia, myopathy, and bone cancer. We previously mapped the DMS-MFH tumor-suppressing-gene locus to chromosomal region 9p21-22 but failed to identify mutations in known genes in this region. We now demonstrate that DMS-MFH results from mutations in the most proximal of three previously uncharacterized terminal exons of the gene encoding methylthioadenosine phosphorylase, MTAP. Intriguingly, two of these MTAP exons arose from early and independent retroviral-integration events in primate genomes at least 40 million years ago, and since then, their genomic integration has gained a functional role. MTAP is a ubiquitously expressed homotrimeric-subunit enzyme critical to polyamine metabolism and adenine and methionine salvage pathways and was believed to be encoded as a single transcript from the eight previously described exons. Six distinct retroviral-sequence-containing MTAP isoforms, each of which can physically interact with archetype MTAP, have been identified. The disease-causing mutations occur within one of these retroviral-derived exons and result in exon skipping and dysregulated alternative splicing of all MTAP isoforms. Our results identify a gene involved in the development of bone sarcoma, provide evidence of the primate-specific evolution of certain parts of an existing gene, and demonstrate that mutations in parts of this gene can result in human disease despite its relatively recent origin.

Diaphyseal medullary stenosis with malignant fibrous histiocytoma (DMS-MFH) is an autosomal-dominant syndrome characterized by bone dysplasia, myopathy, and bone cancer. We previously mapped the DMS-MFH tumor-suppressing-gene locus to chromosomal region 9p21-22 but failed to identify mutations in known genes in this region. We now demonstrate that DMS-MFH results from mutations in the most proximal of three previously uncharacterized terminal exons of the gene encoding methylthioadenosine phosphorylase, MTAP. Intriguingly, two of these MTAP exons arose from early and independent retroviral-integration events in primate genomes at least 40 million years ago, and since then, their genomic integration has gained a functional role. MTAP is a ubiquitously expressed homotrimeric-subunit enzyme critical to polyamine metabolism and adenine and methionine salvage pathways and was believed to be encoded as a single transcript from the eight previously described exons. Six distinct retroviral-sequence-containing MTAP isoforms, each of which can physically interact with archetype MTAP, have been identified. The disease-causing mutations occur within one of these retroviral-derived exons and result in exon skipping and dysregulated alternative splicing of all MTAP isoforms. Our results identify a gene involved in the development of bone sarcoma, provide evidence of the primate-specific evolution of certain parts of an existing gene, and demonstrate that mutations in parts of this gene can result in human disease despite its relatively recent origin.

Five single nucleotide polymorphisms (SNPs) associated with thyroid cancer (TC) risk have been reported: rs2910164 (5q24); rs6983267 (8q24); rs965513 and rs1867277 (9q22); and rs944289 (14q13). Most of these associations have not been replicated in independent populations and the combined effects of the SNPs on risk have not been examined. This study genotyped the five TC SNPs in 781 patients recruited through the TCUKIN study. Genotype data from 6122 controls were obtained from the CORGI and Wellcome Trust Case-Control Consortium studies. Significant associations were detected between TC and rs965513A (p=6.35×10(-34)), rs1867277A (p=5.90×10(-24)), rs944289T (p=6.95×10(-7)), and rs6983267G (p=0.016). rs6983267 was most strongly associated under a recessive model (P(GG vs GT + TT)=0.004), in contrast to the association of this SNP with other cancer types. However, no evidence was found of an association between rs2910164 and disease under any risk model (p>0.7). The rs1867277 association remained significant (p=0.008) after accounting for genotypes at the nearby rs965513 (p=2.3×10(-13)) and these SNPs did not tag a single high risk haplotype. The four validated TC SNPs accounted for a relatively large proportion (∼11%) of the sibling relative risk of TC, principally owing to the large effect size of rs965513 (OR 1.74).

BackgroundOur recent genome-wide association study identified a novel breast cancer susceptibility locus at 9q31.2 (rs865686).MethodsTo further investigate the rs865686-breast cancer association, we conducted a replication study within the Breast Cancer Association Consortium, which comprises 37 case-control studies (48,394 cases, 50,836 controls).ResultsThis replication study provides additional strong evidence of an inverse association between rs865686 and breast cancer risk [study-adjusted per G-allele OR, 0.90; 95% confidence interval (CI), 0.88; 0.91, P = 2.01 × 10(-29)] among women of European ancestry. There were ethnic differences in the estimated minor (G)-allele frequency among controls [0.09, 0.30, and 0.38 among, respectively, Asians, Eastern Europeans, and other Europeans; P for heterogeneity (P(het)) = 1.3 × 10(-143)], but no evidence of ethnic differences in per allele OR (P(het) = 0.43). rs865686 was associated with estrogen receptor-positive (ER(+)) disease (per G-allele OR, 0.89; 95% CI, 0.86-0.91; P = 3.13 × 10(-22)) but less strongly, if at all, with ER-negative (ER(-)) disease (OR, 0.98; 95% CI, 0.94-1.02; P = 0.26; P(het) = 1.16 × 10(-6)), with no evidence of independent heterogeneity by progesterone receptor or HER2 status. The strength of the breast cancer association decreased with increasing age at diagnosis, with case-only analysis showing a trend in the number of copies of the G allele with increasing age at diagnosis (P for linear trend = 0.0095), but only among women with ER(+) tumors.ConclusionsThis study is the first to show that rs865686 is a susceptibility marker for ER(+) breast cancer.ImpactThe findings further support the view that genetic susceptibility varies according to tumor subtype.

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

BackgroundColorectal cancer (CRC) is the second cause of cancer-related death in the Western world. Much of the CRC genetic risk remains unidentified and may be attributable to a large number of common, low-penetrance genetic variants. Genetic linkage studies in CRC families have reported additional association with regions 9q22-31, 3q21-24, 7q31, 11q, 14q and 22q. There are several plausible candidate genes for CRC susceptibility within the aforementioned linkage regions including PTCH1, XPA and TGFBR1 in 9q22-31, and EPHB1 and MRAS in 3q21-q24.MethodsCRC cases and matched controls were from EPICOLON, a prospective, multicentre, nationwide Spanish initiative, composed of two independent phases. Phase 1 corresponded to 515 CRC cases and 515 controls, whereas phase 2 consisted of 901 CRC cases and 909 controls. Genotyping was performed for 172 single-nucleotide polymorphisms (SNPs) in 84 genes located within regions 9q22-31 and 3q21-q24.ResultsNone of the 172 SNPs analysed in our study could be formally associated with CRC risk. However, rs1444601 (TOPBP1) and rs13088006 (CDV3) in region 3q22 showed interesting results and may have an effect on CRC risk.ConclusionsTOPBP1 and CDV3 genetic variants on region 3q22 may modulate CRC risk. Further validation and meta-analysis should be undertaken in larger CRC cohorts.

BackgroundColorectal cancer (CRC) is the second cause of cancer-related death in the Western world. Much of the CRC genetic risk remains unidentified and may be attributable to a large number of common, low-penetrance genetic variants. Genetic linkage studies in CRC families have reported additional association with regions 9q22-31, 3q21-24, 7q31, 11q, 14q and 22q. There are several plausible candidate genes for CRC susceptibility within the aforementioned linkage regions including PTCH1, XPA and TGFBR1 in 9q22-31, and EPHB1 and MRAS in 3q21-q24.MethodsCRC cases and matched controls were from EPICOLON, a prospective, multicentre, nationwide Spanish initiative, composed of two independent phases. Phase 1 corresponded to 515 CRC cases and 515 controls, whereas phase 2 consisted of 901 CRC cases and 909 controls. Genotyping was performed for 172 single-nucleotide polymorphisms (SNPs) in 84 genes located within regions 9q22-31 and 3q21-q24.ResultsNone of the 172 SNPs analysed in our study could be formally associated with CRC risk. However, rs1444601 (TOPBP1) and rs13088006 (CDV3) in region 3q22 showed interesting results and may have an effect on CRC risk.ConclusionsTOPBP1 and CDV3 genetic variants on region 3q22 may modulate CRC risk. Further validation and meta-analysis should be undertaken in larger CRC cohorts.

We conducted a two-stage genome-wide association study of pancreatic cancer, a cancer with one of the lowest survival rates worldwide. We genotyped 558,542 SNPs in 1,896 individuals with pancreatic cancer and 1,939 controls drawn from 12 prospective cohorts plus one hospital-based case-control study. We conducted a combined analysis of these groups plus an additional 2,457 affected individuals and 2,654 controls from eight case-control studies, adjusting for study, sex, ancestry and five principal components. We identified an association between a locus on 9q34 and pancreatic cancer marked by the SNP rs505922 (combined P = 5.37 x 10(-8); multiplicative per-allele odds ratio 1.20; 95% confidence interval 1.12-1.28). This SNP maps to the first intron of the ABO blood group gene. Our results are consistent with earlier epidemiologic evidence suggesting that people with blood group O may have a lower risk of pancreatic cancer than those with groups A or B.

Phosphotidylinositol (PtdIns) signaling is tightly regulated both spatially and temporally by subcellularly localized PtdIns kinases and phosphatases that dynamically alter downstream signaling events. Joubert syndrome is characterized by a specific midbrain-hindbrain malformation ('molar tooth sign'), variably associated retinal dystrophy, nephronophthisis, liver fibrosis and polydactyly and is included in the newly emerging group of 'ciliopathies'. In individuals with Joubert disease genetically linked to JBTS1, we identified mutations in the INPP5E gene, encoding inositol polyphosphate-5-phosphatase E, which hydrolyzes the 5-phosphate of PtdIns(3,4,5)P3 and PtdIns(4,5)P2. Mutations clustered in the phosphatase domain and impaired 5-phosphatase activity, resulting in altered cellular PtdIns ratios. INPP5E localized to cilia in major organs affected by Joubert syndrome, and mutations promoted premature destabilization of cilia in response to stimulation. These data link PtdIns signaling to the primary cilium, a cellular structure that is becoming increasingly recognized for its role in mediating cell signals and neuronal function.

Many candidate genes have been studied for asthma, but replication has varied. Novel candidate genes have been identified for various complex diseases using genome-wide association studies (GWASs). We conducted a GWAS in 492 Mexican children with asthma, predominantly atopic by skin prick test, and their parents using the Illumina HumanHap 550 K BeadChip to identify novel genetic variation for childhood asthma. The 520,767 autosomal single nucleotide polymorphisms (SNPs) passing quality control were tested for association with childhood asthma using log-linear regression with a log-additive risk model. Eleven of the most significantly associated GWAS SNPs were tested for replication in an independent study of 177 Mexican case-parent trios with childhood-onset asthma and atopy using log-linear analysis. The chromosome 9q21.31 SNP rs2378383 (p = 7.10x10(-6) in the GWAS), located upstream of transducin-like enhancer of split 4 (TLE4), gave a p-value of 0.03 and the same direction and magnitude of association in the replication study (combined p = 6.79x10(-7)). Ancestry analysis on chromosome 9q supported an inverse association between the rs2378383 minor allele (G) and childhood asthma. This work identifies chromosome 9q21.31 as a novel susceptibility locus for childhood asthma in Mexicans. Further, analysis of genome-wide expression data in 51 human tissues from the Novartis Research Foundation showed that median GWAS significance levels for SNPs in genes expressed in the lung differed most significantly from genes not expressed in the lung when compared to 50 other tissues, supporting the biological plausibility of our overall GWAS findings and the multigenic etiology of childhood asthma.

We conducted a two-stage genome-wide association study of pancreatic cancer, a cancer with one of the lowest survival rates worldwide. We genotyped 558,542 SNPs in 1,896 individuals with pancreatic cancer and 1,939 controls drawn from 12 prospective cohorts plus one hospital-based case-control study. We conducted a combined analysis of these groups plus an additional 2,457 affected individuals and 2,654 controls from eight case-control studies, adjusting for study, sex, ancestry and five principal components. We identified an association between a locus on 9q34 and pancreatic cancer marked by the SNP rs505922 (combined P = 5.37 x 10(-8); multiplicative per-allele odds ratio 1.20; 95% confidence interval 1.12-1.28). This SNP maps to the first intron of the ABO blood group gene. Our results are consistent with earlier epidemiologic evidence suggesting that people with blood group O may have a lower risk of pancreatic cancer than those with groups A or B.

Many candidate genes have been studied for asthma, but replication has varied. Novel candidate genes have been identified for various complex diseases using genome-wide association studies (GWASs). We conducted a GWAS in 492 Mexican children with asthma, predominantly atopic by skin prick test, and their parents using the Illumina HumanHap 550 K BeadChip to identify novel genetic variation for childhood asthma. The 520,767 autosomal single nucleotide polymorphisms (SNPs) passing quality control were tested for association with childhood asthma using log-linear regression with a log-additive risk model. Eleven of the most significantly associated GWAS SNPs were tested for replication in an independent study of 177 Mexican case-parent trios with childhood-onset asthma and atopy using log-linear analysis. The chromosome 9q21.31 SNP rs2378383 (p = 7.10x10(-6) in the GWAS), located upstream of transducin-like enhancer of split 4 (TLE4), gave a p-value of 0.03 and the same direction and magnitude of association in the replication study (combined p = 6.79x10(-7)). Ancestry analysis on chromosome 9q supported an inverse association between the rs2378383 minor allele (G) and childhood asthma. This work identifies chromosome 9q21.31 as a novel susceptibility locus for childhood asthma in Mexicans. Further, analysis of genome-wide expression data in 51 human tissues from the Novartis Research Foundation showed that median GWAS significance levels for SNPs in genes expressed in the lung differed most significantly from genes not expressed in the lung when compared to 50 other tissues, supporting the biological plausibility of our overall GWAS findings and the multigenic etiology of childhood asthma.

Phosphotidylinositol (PtdIns) signaling is tightly regulated both spatially and temporally by subcellularly localized PtdIns kinases and phosphatases that dynamically alter downstream signaling events. Joubert syndrome is characterized by a specific midbrain-hindbrain malformation ('molar tooth sign'), variably associated retinal dystrophy, nephronophthisis, liver fibrosis and polydactyly and is included in the newly emerging group of 'ciliopathies'. In individuals with Joubert disease genetically linked to JBTS1, we identified mutations in the INPP5E gene, encoding inositol polyphosphate-5-phosphatase E, which hydrolyzes the 5-phosphate of PtdIns(3,4,5)P3 and PtdIns(4,5)P2. Mutations clustered in the phosphatase domain and impaired 5-phosphatase activity, resulting in altered cellular PtdIns ratios. INPP5E localized to cilia in major organs affected by Joubert syndrome, and mutations promoted premature destabilization of cilia in response to stimulation. These data link PtdIns signaling to the primary cilium, a cellular structure that is becoming increasingly recognized for its role in mediating cell signals and neuronal function.

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease affecting both joints and extra-articular tissues. Although some genetic risk factors for RA are well-established, most notably HLA-DRB1 and PTPN22, these markers do not fully account for the observed heritability. To identify additional susceptibility loci, we carried out a multi-tiered, case-control association study, genotyping 25,966 putative functional SNPs in 475 white North American RA patients and 475 matched controls. Significant markers were genotyped in two additional, independent, white case-control sample sets (661 cases/1322 controls from North America and 596 cases/705 controls from The Netherlands) identifying a SNP, rs1953126, on chromosome 9q33.2 that was significantly associated with RA (OR(common) = 1.28, trend P(comb) = 1.45E-06). Through a comprehensive fine-scale-mapping SNP-selection procedure, 137 additional SNPs in a 668 kb region from MEGF9 to STOM on 9q33.2 were chosen for follow-up genotyping in a staged-approach. Significant single marker results (P(comb)<0.01) spanned a large 525 kb region from FBXW2 to GSN. However, a variety of analyses identified SNPs in a 70 kb region extending from the third intron of PHF19 across TRAF1 into the TRAF1-C5 intergenic region, but excluding the C5 coding region, as the most interesting (trend P(comb): 1.45E-06 --> 5.41E-09). The observed association patterns for these SNPs had heightened statistical significance and a higher degree of consistency across sample sets. In addition, the allele frequencies for these SNPs displayed reduced variability between control groups when compared to other SNPs. Lastly, in combination with the other two known genetic risk factors, HLA-DRB1 and PTPN22, the variants reported here generate more than a 45-fold RA-risk differential.

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease affecting both joints and extra-articular tissues. Although some genetic risk factors for RA are well-established, most notably HLA-DRB1 and PTPN22, these markers do not fully account for the observed heritability. To identify additional susceptibility loci, we carried out a multi-tiered, case-control association study, genotyping 25,966 putative functional SNPs in 475 white North American RA patients and 475 matched controls. Significant markers were genotyped in two additional, independent, white case-control sample sets (661 cases/1322 controls from North America and 596 cases/705 controls from The Netherlands) identifying a SNP, rs1953126, on chromosome 9q33.2 that was significantly associated with RA (OR(common) = 1.28, trend P(comb) = 1.45E-06). Through a comprehensive fine-scale-mapping SNP-selection procedure, 137 additional SNPs in a 668 kb region from MEGF9 to STOM on 9q33.2 were chosen for follow-up genotyping in a staged-approach. Significant single marker results (P(comb)<0.01) spanned a large 525 kb region from FBXW2 to GSN. However, a variety of analyses identified SNPs in a 70 kb region extending from the third intron of PHF19 across TRAF1 into the TRAF1-C5 intergenic region, but excluding the C5 coding region, as the most interesting (trend P(comb): 1.45E-06 --> 5.41E-09). The observed association patterns for these SNPs had heightened statistical significance and a higher degree of consistency across sample sets. In addition, the allele frequencies for these SNPs displayed reduced variability between control groups when compared to other SNPs. Lastly, in combination with the other two known genetic risk factors, HLA-DRB1 and PTPN22, the variants reported here generate more than a 45-fold RA-risk differential.

PurposeUsing fluorescence in situ hybridization (FISH) and high-density single nucleotide polymorphism (SNP) mapping genome array, we comparatively evaluated chromosome 3 status and other chromosomal aberrations within a series of choroidal melanomas biopsied by fine needle aspiration (FNAB).MethodsTransscleral FNAB was performed in 59 patients (59 eyes) who had a clinical diagnosis of choroidal melanoma. Biopsies were processed for chromosome 3 status by centromeric interphase FISH, cytopathology, cell culture, and simultaneous genomic DNA and RNA mapping array analysis.ResultsFISH yielded chromosome 3 status in 38 of 59 (64%) eyes, while high-density SNP mapping array yielded chromosome 3 status in 43 of 59 (73%) eyes. Monosomy 3 was detected by FISH in 15 of 38 (39%) cases, and high-density SNP mapping array data confirmed the finding in 13 of the 15 cases. Furthermore, high-density SNP mapping array revealed five additional cases of significant chromosome 3 aberration not detected by FISH. High-density genomic mapping also provided detailed patterns of chromosomal gain and loss on chromosomes 1, 6, 8, and 9 which segregated into two groups characterized by either monosomy 3 or chromosome 6p gain.ConclusionsHigh-density SNP mapping array was better than FISH in detecting chromosome 3 aberrations and monosomy in our melanoma samples. More importantly, the mapping arrays detected additional patterns of chromosomal aberration, which suggest specific pathways for cytogenetic rearrangements in choroidal melanoma and may improve prognostic testing.