Paul Leo

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis–associated haplotypes at 11 loci. Two ankylosing spondylitis–associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.

Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1

We report a genome-wide association study for open-angle glaucoma (OAG) blindness using a discovery cohort of 590 individuals with severe visual field loss (cases) and 3,956 controls. We identified associated loci at TMCO1 (rs4656461[G] odds ratio (OR) = 1.68, P = 6.1 × 10−10) and CDKN2B-AS1 (rs4977756[A] OR = 1.50, P = 4.7 × 10−9). We replicated these associations in an independent cohort of cases with advanced OAG (rs4656461 P = 0.010; rs4977756 P = 0.042) and two additional cohorts of less severe OAG (rs4656461 combined discovery and replication P = 6.00 × 10−14, OR = 1.51, 95% CI 1.35–1.68; rs4977756 combined P = 1.35 × 10−14, OR = 1.39, 95% CI 1.28–1.51). We show retinal expression of genes at both loci in human ocular tissues. We also show that CDKN2A and CDKN2B are upregulated in the retina of a rat model of glaucoma.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55–85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or −4.0 to −1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD–associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of −0.11 standard deviations [SD] per C allele, P = 6.2×10−9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (−0.14 SD per C allele, P = 2.3×10−12, and −0.16 SD per G allele, P = 1.2×10−15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10−9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10−6 and rs2707466: OR = 1.22, P = 7.2×10−6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16−/− mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%–61% (6.5×10−13<P<5.9×10−4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture.

Defects in the IFT-B Component IFT172 Cause Jeune and Mainzer-Saldino Syndromes in Humans

Intraflagellar transport (IFT) depends on two evolutionarily conserved modules, subcomplexes A (IFT-A) and B (IFT-B), to drive ciliary assembly and maintenance. All six IFT-A components and their motor protein, DYNC2H1, have been linked to human skeletal ciliopathies, including asphyxiating thoracic dystrophy (ATD; also known as Jeune syndrome), Sensenbrenner syndrome, and Mainzer-Saldino syndrome (MZSDS). Conversely, the 14 subunits in the IFT-B module, with the exception of IFT80, have unknown roles in human disease. To identify additional IFT-B components defective in ciliopathies, we independently performed different mutation analyses: candidate-based sequencing of all IFT-B-encoding genes in 1,467 individuals with a nephronophthisis-related ciliopathy or whole-exome resequencing in 63 individuals with ATD. We thereby detected biallelic mutations in the IFT-B-encoding gene IFT172 in 12 families. All affected individuals displayed abnormalities of the thorax and/or long bones, as well as renal, hepatic, or retinal involvement, consistent with the diagnosis of ATD or MZSDS. Additionally, cerebellar aplasia or hypoplasia characteristic of Joubert syndrome was present in 2 out of 12 families. Fibroblasts from affected individuals showed disturbed ciliary composition, suggesting alteration of ciliary transport and signaling. Knockdown of ift172 in zebrafish recapitulated the human phenotype and demonstrated a genetic interaction between ift172 and ift80. In summary, we have identified defects in IFT172 as a cause of complex ATD and MZSDS. Our findings link the group of skeletal ciliopathies to an additional IFT-B component, IFT172, similar to what has been shown for IFT-A.

Functional Analyses of Glycyl-tRNA Synthetase Mutations Suggest a Key Role for tRNA-Charging Enzymes in Peripheral Axons

Charcot–Marie–Tooth disease type 2D (CMT2D) and distal spinal muscular atrophy type V (dSMA-V) are axonal neuropathies characterized by a phenotype that is more severe in the upper extremities. We previously implicated mutations in the gene encoding glycyl-tRNA synthetase (GARS) as the cause of CMT2D and dSMA-V. GARS is a member of the family of aminoacyl-tRNA synthetases responsible for charging tRNA with cognate amino acids; GARS ligates glycine to tRNAGly. Here, we present functional analyses of disease-associated GARS mutations and show that there are not any significant mutation-associated changes in GARS expression levels; that the majority of identified GARS mutations modeled in yeast severely impair viability; and that, in most cases, mutant GARS protein mislocalizes in neuronal cells. Indeed, four of the five mutations studied show loss-of-function features in at least one assay, suggesting that tRNA-charging deficits play a role in disease pathogenesis. Finally, we detected endogenous GARS-associated granules in the neurite projections of cultured neurons and in the peripheral nerve axons of normal human tissue. These data are particularly important in light of the recent identification of CMT-associated mutations in another tRNA synthetase gene [YARS (tyrosyl-tRNA synthetase gene)]. Together, these findings suggest that tRNA-charging enzymes play a key role in maintaining peripheral axons.

Short-Rib Polydactyly and Jeune Syndromes Are Caused by Mutations in WDR60

Short-rib polydactyly syndromes (SRPS I-V) are a group of lethal congenital disorders characterized by shortening of the ribs and long bones, polydactyly, and a range of extraskeletal phenotypes. A number of other disorders in this grouping, including Jeune and Ellis-van Creveld syndromes, have an overlapping but generally milder phenotype. Collectively, these short-rib dysplasias (with or without polydactyly) share a common underlying defect in primary cilium function and form a subset of the ciliopathy disease spectrum. By using whole-exome capture and massive parallel sequencing of DNA from an affected Australian individual with SRPS type III, we detected two novel heterozygous mutations in WDR60, a relatively uncharacterized gene. These mutations segregated appropriately in the unaffected parents and another affected family member, confirming compound heterozygosity, and both were predicted to have a damaging effect on the protein. Analysis of an additional 54 skeletal ciliopathy exomes identified compound heterozygous mutations in WDR60 in a Spanish individual with Jeune syndrome of relatively mild presentation. Of note, these two families share one novel WDR60 missense mutation, although haplotype analysis suggested no shared ancestry. We further show that WDR60 localizes at the base of the primary cilium in wild-type human chondrocytes, and analysis of fibroblasts from affected individuals revealed a defect in ciliogenesis and aberrant accumulation of the GLI2 transcription factor at the centrosome or basal body in the absence of an obvious axoneme. These findings show that WDR60 mutations can cause skeletal ciliopathies and suggest a role for WDR60 in ciliogenesis.

A potent Chk1 inhibitor is selectively cytotoxic in melanomas with high levels of replicative stress

There are few effective treatments for metastatic melanoma. Checkpoint kinase 1 (Chk1) inhibitors are being trialled for their efficacy in enhancing conventional chemotherapeutic agents, but their effectiveness as single agents is not known. We have examined the effectiveness of two novel Chk1 selective inhibitors, AR323 and AR678, in a panel of melanoma cell lines and normal cell types. We demonstrate that these drugs display single-agent activity, with IC50s in the low nanomolar range. The drugs produce cytotoxic effects in cell lines that are most sensitive to these drugs, whereas normal cells are only sensitive to these drugs at the higher concentrations where they have cytostatic activity. The cytotoxic effect is the consequence of inhibition of S-phase Chk1, which drives cells prematurely from late S phase into an aberrant mitosis and results in either failure of cytokinesis or cell death through an apoptotic mechanism. The sensitivity to the Chk1 inhibitors was correlated with the level of endogenous DNA damage indicating replicative stress. Chk1 inhibitors are viable single-agent therapies that target melanoma cells with high levels of endogenous DNA damage. This sensitivity suggests that Chk1 is a critical component of an adaptation to replicative stress in these cells. It also suggests that markers of DNA damage may be useful in identifying the melanomas and potentially other tumour types that are more likely to be sensitive to Chk1 inhibitors as single agents.

Major histocompatibility complex associations of ankylosing spondylitis are complex and involve further epistasis with ERAP1.

Ankylosing spondylitis (AS) is a common, highly heritable, inflammatory arthritis for which HLA-B*27 is the major genetic risk factor, although its role in the aetiology of AS remains elusive. To better understand the genetic basis of the MHC susceptibility loci, we genotyped 7,264 MHC SNPs in 22,647 AS cases and controls of European descent. We impute SNPs, classical HLA alleles and amino-acid residues within HLA proteins, and tested these for association to AS status. Here we show that in addition to effects due to HLA-B*27 alleles, several other HLA-B alleles also affect susceptibility. After controlling for the associated haplotypes in HLA-B, we observe independent associations with variants in the HLA-A, HLA-DPB1 and HLA-DRB1 loci. We also demonstrate that the ERAP1 SNP rs30187 association is not restricted only to carriers of HLA-B*27 but also found in HLA-B*40:01 carriers independently of HLA-B*27 genotype.

Integrated genome-wide chromatin occupancy and expression analyses identify key myeloid pro-differentiation transcription factors repressed by Myb

To gain insight into the mechanisms by which the Myb transcription factor controls normal hematopoiesis and particularly, how it contributes to leukemogenesis, we mapped the genome-wide occupancy of Myb by chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-Seq) in ERMYB myeloid progenitor cells. By integrating the genome occupancy data with whole genome expression profiling data, we identified a Myb-regulated transcriptional program. Gene signatures for leukemia stem cells, normal hematopoietic stem/progenitor cells and myeloid development were overrepresented in 2368 Myb regulated genes. Of these, Myb bound directly near or within 793 genes. Myb directly activates some genes known critical in maintaining hematopoietic stem cells, such as Gfi1 and Cited2. Importantly, we also show that, despite being usually considered as a transactivator, Myb also functions to repress approximately half of its direct targets, including several key regulators of myeloid differentiation, such as Sfpi1 (also known as Pu.1), Runx1, Junb and Cebpb. Furthermore, our results demonstrate that interaction with p300, an established coactivator for Myb, is unexpectedly required for Myb-mediated transcriptional repression. We propose that the repression of the above mentioned key pro-differentiation factors may contribute essentially to Myb’s ability to suppress differentiation and promote self-renewal, thus maintaining progenitor cells in an undifferentiated state and promoting leukemic transformation.