Genome-wide association study of ulcerative colitis identifies three new susceptibility loci, including the HNF4A region

Ulcerative colitis is a common form of inflammatory bowel disease with a complex etiology. As part of the Wellcome Trust Case Control Consortium 2, we performed a genome-wide association scan for ulcerative colitis in 2,361 cases and 5,417 controls. Loci showing evidence of association at P 1 × 10 5 were followed up by genotyping in an independent set of 2,321 cases and 4,818 controls. We find genome-wide significant evidence of association at three new loci, each containing at least one biologically relevant candidate gene, on chromosomes 20q13 (HNF4A; P = 3.2 × 10 17), 16q22 (CDH1 and CDH3; P = 2.8 × 10 8) and 7q31 (LAMB1; P = 3.0 × 10 8). Of note, CDH1 has recently been associated with susceptibility to colorectal cancer, an established complication of longstanding ulcerative colitis. The new associations suggest that changes in the integrity of the intestinal epithelial barrier may contribute to the pathogenesis of ulcerative colitis. © 2009 Nature America, Inc. All rights reserved.

Susceptibility variants for male-pattern baldness on chromosome 20p11

We carried out a genome-wide association study in 296 individuals with male-pattern baldness (androgenetic alopecia) and 347 controls. We then investigated the 30 best SNPs in an independent replication sample and found highly significant association for five SNPs on chromosome 20p11 (rs2180439 combined P = 2.7 x 10(-15)). No interaction was detected with the X-chromosomal androgen receptor locus, suggesting that the 20p11 locus has a role in a yet-to-be-identified androgen-independent pathway.

AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma

Invasive urothelial cell carcinoma (UCC) is characterized by increased chromosomal instability and follows an aggressive clinical course in contrast to non-invasive disease. To identify molecular processes that confer and maintain an aggressive malignant phenotype, we used a high-throughput genome-wide approach to interrogate a cohort of high and low clinical risk UCC tumors. Differential expression analyses highlighted cohesive dysregulation of critical genes involved in the G(2)/M checkpoint in aggressive UCC. Hierarchical clustering based on DNA Damage Response (DDR) genes separated tumors according to a pre-defined clinical risk phenotype. Using array-comparative genomic hybridization, we confirmed that the DDR was disrupted in tumors displaying high genomic instability. We identified DNA copy number gains at 20q13.2-q13.3 (AURKA locus) and determined that overexpression of AURKA accompanied dysregulation of DDR genes in high risk tumors. We postulated that DDR-deficient UCC tumors are advantaged by a selective pressure for AURKA associated override of M phase barriers and confirmed this in an independent tissue microarray series. This mechanism that enables cancer cells to maintain an aggressive phenotype forms a rationale for targeting AURKA as a therapeutic strategy in advanced stage UCC.

Male-pattern baldness susceptibility locus at 20p11.

We conducted a genome-wide association study for androgenic alopecia in 1,125 men and identified a newly associated locus at chromosome 20p11.22, confirmed in three independent cohorts (n = 1,650; OR = 1.60, P = 1.1 x 10(-14) for rs1160312). The one man in seven who harbors risk alleles at both 20p11.22 and AR (encoding the androgen receptor) has a sevenfold-increased odds of androgenic alopecia (OR = 7.12, P = 3.7 x 10(-15)).

Characterization of IGH locus breakpoints in multiple myeloma indicates a subset of translocations appear to occur in pregerminal center B cells.

Translocations in myeloma are thought to occur solely in mature B cells in the germinal center through class switch recombination (CSR). We used a targeted captured technique followed by massively parallel sequencing to determine the exact breakpoints in both the immunoglobulin heavy chain (IGH) locus and the partner chromosome in 61 presentation multiple myeloma samples. The majority of samples (62%) have a breakpoint within the switch regions upstream of the IGH constant genes and are generated through CSR in a mature B cell. However, the proportion of CSR translocations is not consistent between cytogenetic subgroups. We find that 100% of t(4;14) are CSR-mediated; however, 21% of t(11;14) and 25% of t(14;20) are generated through DH-JH recombination activation gene-mediated mechanisms, indicating they occur earlier in B-cell development at the pro-B-cell stage in the bone marrow. These 2 groups also generate translocations through receptor revision, as determined by the breakpoints and mutation status of the segments used in 10% and 50% of t(11;14) and t(14;20) samples, respectively. The study indicates that in a significant number of cases the translocation-based etiological events underlying myeloma may arise at the pro-B-cell hematological progenitor cell level, much earlier in B-cell development than was previously thought.

Characterising the TP53-deleted subgroup of chronic lymphocytic leukemia: an analysis of additional cytogenetic abnormalities detected by interphase fluorescence in situ hybridisation and array-based comparative genomic hybridisation.

Deletion of the TP53 gene on chromosome 17p13.1 is the prognostic factor associated with the shortest survival in CLL. We used array-based comparative genomic hybridisation (arrayCGH) to identify additional DNA copy number changes in peripheral blood samples from 74 LRF CLL4 trial patients, 37 with >or=5% and 37 without TP53-deleted cells. ArrayCGH reliably detected deletions on 17p, including the TP53 locus, in cases with >or=50%TP53-deleted cells detected by fluorescence in situ hybridisation, plus seven additional cases with deleted regions on 17p excluding TP53. Losses on chromosomal regions 18p and/or 20p were found exclusively in cases with >or=5%TP53-deleted cells (por=5%TP53-deleted cases (p=0.02). In particular, amplification of 2p and deletion of 6q were both more frequent. Cases with >20%TP53-deleted cells had the worst prognosis in the LRF CLL4 trial.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Compilation of somatic mutations of the CDKN2 gene in human cancers : non-random distribution of base substitutions

The CDKN2 gene, encoding the cyclin-dependent kinase inhibitor p16, is a tumour suppressor gene that maps to chromosome band 9p21-p22. The most common mechanism of inactivation of this gene in human cancers is through homozygous deletion; however, in a smaller proportion of tumours and tumour cell lines intragenic mutations occur. In this study we have compiled a database of over 120 published point mutations in the CDKN2 gene from a wide variety of tumour types. A further 50 deletions, insertions, and splice mutations in CDKN2 have also been compiled. Furthermore, we have standardised the numbering of all mutations according to the full-length 156 amino acid form of p16. From this study we are able to define several hot spots, some of which occur at conserved residues within the ankyrin domains of p16. While many of the hotspots are shared by a number of cancers, the relative importance of each position varies, possibly reflecting the role of different carcinogens in the development of certain tumours. As reported previously, the mutational spectrum of CDKN2 in melanomas differs from that of internal malignancies and supports the involvement of UV in melanoma tumorigenesis. Notably, 52% of all substitutions in melanoma-derived samples occurred at just six nucleotide positions. Nonsense mutations comprise a comparatively high proportion of mutations present in the CDKN2 gene, and possible explanations for this are discussed.

Evidence for three tumor suppressor loci on chromosome 9p involved in melanoma development

Cytogenetic and loss of heterozygosity (LOH) studies have long indicated the presence of a tumor suppressor gene (TSG) on 9p involved in the development of melanoma. Although LOH at 9p has been reported in approximately 60% of melanoma tumors, only 5-10% of these tumors have been shown to carry CDKN2A mutations, raising the possibility that another TSG involved in melanoma maps to chromosome 9p. To investigate this possibility, a panel of 37 melanomas derived from 35 individuals was analyzed for CDKN2A mutations by single-strand conformation polymorphism analysis and sequencing. The melanoma samples were then typed for 15 markers that map to 9p13-24 to investigate LOH trends in this region. In those tumors demonstrating retention of heterozygosity at markers flanking CDKN2A and LOH on one or both sides of the gene, multiplex microsatellite PCR was performed to rule out homozygous deletion of the region encompassing CDKN2A. CDKN2A mutations were found in tumors from 5 patients [5 (14%) of 35], 4 of which demonstrated LOH across the entire region examined. The remaining tumor with no observed LOH carried two point mutations, one on each allele. Although LOH was identified at one or more markers in 22 (59%) of 37 melanoma tumors corresponding to 20 (57%) of 35 individuals, only 11 tumors from 9 individuals [9 (26%) of 35] demonstrated LOH at D9S942 and D9S1748 the markers closest to CDKN2A. Of the remaining 11 tumors with LOH 9 demonstrated LOH at two or more contiguous markers either centromeric and/or telomeric to CDKN2A while retaining heterozygosity at several markers adjacent to CDKN2A. Multiplex PCR revealed one tumor carried a homozygous deletion extending from D9S1748 to the IFN-alpha locus. In the remaining eight tumors, multiplex PCR demonstrated that the observed heterozygosity was not attributable to homozygous deletion and stromal contamination at D9S1748, D9S942, or D9S974, as measured by comparative amplification strengths, which indicates that retention of heterozygosity with flanking LOH does not always indicate a homozygous deletion. This report supports the conclusions of previous studies that a least two TSGs involved in melanoma development in addition to CDKN2A may reside on chromosome 9p.

PTEN inactivation is rare in melanoma tumours but occurs frequently in melanoma cell lines

Deletions detected in cytogenetic and loss of heterozygosity (LOH) studies indicate that at least one tumour suppressor gene maps to the long arm of chromosome 10. Previous deletion mapping studies have observed LOH on 10q in about 30% of melanomas analysed. The PTEN gene, mapping to chromosome band 10q23.3, encodes a protein with both lipid and protein phosphatase activity. Somatic mutations and deletions in have been detected in a variety of cell lines and tumours, including melanoma samples. We performed mutation analyses and extensive allelic loss studies to investigate the role this gene plays in melanoma pathogenesis. We found that a total of 34 out of 57 (60%) melanoma cell lines carried hemizygous deletions of chromosome 10q encompassing the PTEN locus. A further three cell lines carried smaller deletions excluding PTEN. Inactivation of both PTEN alleles by exon-specific homozygous deletion or mutation was observed in 13 out of 57 (23%) melanoma cell lines. The mutation spectrum observed does not indicate an important role for ultraviolet radiation in the genesis of these mutations, and evidence from three cell lines supports the acquisition of PTEN aberrations in culture. Ten out of 49 (20%) matched melanoma tumour/normal samples harboured hemizygous deletions of either the whole chromosome or most of the long arm. Mutations within were detected in only one of the 10 tumours demonstrating LOH at 10q23 that were analysed. These results suggest that PTEN inactivation may be important for the propagation of melanoma cells in culture, and that another chromosome 10 tumour suppressor gene may be important for melanoma pathogenesis.

Genome-wide association study identifies new multiple sclerosis susceptibility loci on chromosomes 12 and 20

To identify multiple sclerosis (MS) susceptibility loci, we conducted a genome-wide association study (GWAS) in 1,618 cases and used shared data for 3,413 controls. We performed replication in an independent set of 2,256 cases and 2,310 controls, for a total of 3,874 cases and 5,723 controls. We identified risk-associated SNPs on chromosome 12q13-14 (rs703842, P = 5.4 x 10(-11); rs10876994, P = 2.7 x 10(-10); rs12368653, P = 1.0 x 10(-7)) and upstream of CD40 on chromosome 20q13 (rs6074022, P = 1.3 x 10(-7); rs1569723, P = 2.9 x 10(-7)). Both loci are also associated with other autoimmune diseases. We also replicated several known MS associations (HLA-DR15, P = 7.0 x 10(-184); CD58, P = 9.6 x 10(-8); EVI5-RPL5, P = 2.5 x 10(-6); IL2RA, P = 7.4 x 10(-6); CLEC16A, P = 1.1 x 10(-4); IL7R, P = 1.3 x 10(-3); TYK2, P = 3.5 x 10(-3)) and observed a statistical interaction between SNPs in EVI5-RPL5 and HLA-DR15 (P = 0.001).

The B subunit of coagulation factor XIII is linked to renin and the Duffy blood group to α-spectrin on human chromosome 1

Family linkage studies were used to detect two linkage relationships on human chromosome 1. The B subunit of coagulation factor XIII showed significant linkage to renin with a maximum lod score of 5.071 at a distance of 10 cM. Significant linkage was also shown between the Duffy blood group and α-spectrin with linkage results giving a combined lod score of 3.194 at 5 cM.

Regional chromosomal assignment of human renin gene to 1q12→qter and use in linkage studies in Charcot-Marie-Tooth disease

The gene for renin, previously mapped to human chromosome 1, was further localized to 1q12 → qter using human-mouse somatic cell hybrid DNAs. The renin DNA probe used (λ HR5) could detect a HindIII restriction fragment length polymorphism. When used in studies of 12 informative families, no linkage could be found between the renin and Charcot-Marie-Tooth disease. Furthermore, an association of any renin allele with hypertension was not apparent.

Isolation and use of chromosome 1 probes for linkage studies on Charcot-Marie-Tooth disease

Nine probes were isolated from a human chromosome 1 enriched library and mapped to regions of chromosome 1 using somatic cell hybrid lines. One clone, LR67, which mapped 1q12→q23 detected a BglI RFLP. This probe, as well as 4 other known chromosome 1 markers, α-spectrin, Factor XIIIB, DR10 and DR78, were used for linkage studies in 15 Charcot-Marie-Tooth disease (CMT1) families. Close linking of CMT1 to any of the 5 markers was not indicated. Total lod scores excluded linkage of CMT1 to LR67 and to DR10 at 5 cM or less, to DR78 and 10 cM or less, α-spectrin at 15 cM or less and Factor XIIIB at 20 cM or less. Possible linkage, however, was shown between LR67 and CMT1 at a distance of 30 cM. Also linkage at a distance of 5 cM was detected between this probe and α-spectrin.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.