Faculty Bibliography

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating approximately 10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Background:Securing a diagnosis of ovarian cancer and establishing means to predict outcomes to therapeutics remain formidable clinical challenges. Early diagnosis is particularly important since survival rates are markedly improved if tumour is detected early.Methods:Comprehensive miRNA profiles were generated on presurgical plasma samples from 42 women with confirmed serous epithelial ovarian cancer, 36 women diagnosed with a benign neoplasm, and 23 comparably age-matched women with no known pelvic mass.Results:Twenty-two miRNAs were differentially expressed between healthy controls and the ovarian cancer group (P<0.05), while a six miRNA profile subset distinguished presurgical plasma from benign and ovarian cancer patients. There were also significant differences in miRNA profiles in presurgical plasma from women diagnosed with ovarian cancer who had short overall survival when compared to women with long overall survival (P<0.05).Conclusion:Our preliminary data support the utility of circulating plasma miRNAs to distinguish women with ovarian cancer from those with a benign mass and identify women likely to benefit from currently available treatment for serous epithelial ovarian cancer from those who may not.British Journal of Cancer advance online publication, 24 December 2013; doi:10.1038/bjc.2013.795 www.bjcancer.com.

We describe a copy-number variant (CNV) for which deletion alleles confer a protective affect against rheumatoid arthritis (RA). This CNV reflects net unit deletions and expansions to a normal two-unit tandem duplication located on human chr12p13.31, a region with conserved synteny to the rat RA susceptibility quantitative trait loci Oia2. Genotyping, using the paralogue ratio test and SNP intensity data, in Swedish samples (2,403 cases, 1,269 controls) showed that the frequency of deletion variants is significantly lower in cases (P = 0.0012, OR = 0.442 [95%CI 0.258-0.755]). Reduced frequencies of deletion variants were also seen in replication materials comprising 9,201 UK samples (1,846 cases, 7,355 controls) and 2,963 US samples (906 controls, 1,967 cases) (Mantel-Haenszel P = 0.036, OR = 0.559 [95%CI 0.323-0.966]). Combining the three datasets produces a Mantel-Haenszel OR of 0.497 (P < 0.0002). The deletion variant lacks 129-kb of DNA containing SLC2A3, NANOGP1, and SLC2A14. SLC2A3 encodes a high-affinity glucose transporter important in the immune response and chondrocyte metabolism, both key aspects of RA pathogenesis. The large effect size of this association, its potential relevance to other diseases in which SLC2A3 is implicated, and the possibility of targeting drugs to inhibit SLC2A3, argue for further examination of the genetics and the biology of this CNV.

In this study, 1,833 systemic sclerosis (SSc) cases and 3,466 controls were genotyped with the Immunochip array. Classical alleles, amino acid residues, and SNPs across the human leukocyte antigen (HLA) region were imputed and tested. These analyses resulted in a model composed of six polymorphic amino acid positions and seven SNPs that explained the observed significant associations in the region. In addition, a replication step comprising 4,017 SSc cases and 5,935 controls was carried out for several selected non-HLA variants, reaching a total of 5,850 cases and 9,401 controls of European ancestry. Following this strategy, we identified and validated three SSc risk loci, including DNASE1L3 at 3p14, the SCHIP1-IL12A locus at 3q25, and ATG5 at 6q21, as well as a suggested association of the TREH-DDX6 locus at 11q23. The associations of several previously reported SSc risk loci were validated and further refined, and the observed peak of association in PXK was related to DNASE1L3. Our study has increased the number of known genetic associations with SSc, provided further insight into the pleiotropic effects of shared autoimmune risk factors, and highlighted the power of dense mapping for detecting previously overlooked susceptibility loci.

Integrating genetic data from families with highly penetrant forms of disease together with genetic data from outbred populations represents a promising strategy to uncover the complete frequency spectrum of risk alleles for complex traits such as rheumatoid arthritis (RA). Here, we demonstrate that rare, low-frequency and common alleles at one gene locus, phospholipase B1 (PLB1), might contribute to risk of RA in a 4-generation consanguineous pedigree (Middle Eastern ancestry) and also in unrelated individuals from the general population (European ancestry). Through identity-by-descent (IBD) mapping and whole-exome sequencing, we identified a non-synonymous c.2263G>C (p.G755R) mutation at the PLB1 gene on 2q23, which significantly co-segregated with RA in family members with a dominant mode of inheritance (P = 0.009). We further evaluated PLB1 variants and risk of RA using a GWAS meta-analysis of 8,875 RA cases and 29,367 controls of European ancestry. We identified significant contributions of two independent non-coding variants near PLB1 with risk of RA (rs116018341 [MAF = 0.042] and rs116541814 [MAF = 0.021], combined P = 3.2x10(-6)). Finally, we performed deep exon sequencing of PLB1 in 1,088 RA cases and 1,088 controls (European ancestry), and identified suggestive dispersion of rare protein-coding variant frequencies between cases and controls (P = 0.049 for C-alpha test and P = 0.055 for SKAT). Together, these data suggest that PLB1 is a candidate risk gene for RA. Future studies to characterize the full spectrum of genetic risk in the PLB1 genetic locus are warranted.

Few studies had investigated genome-wide methylation in glioblastoma multiforme (GBM). Our goals were to study differential methylation across the genome in gene promoters using an array-based method, as well as repetitive elements using surrogate global methylation markers. The discovery sample set for this study consisted of 54 GBM from Columbia University and Case Western Reserve University, and 24 brain controls from the New York Brain Bank. We assembled a validation dataset using methylation data of 162 TCGA GBM and 140 brain controls from dbGAP. HumanMethylation27 Analysis Bead-Chips (Illumina) were used to interrogate 26,486 informative CpG sites in both the discovery and validation datasets. Global methylation levels were assessed by analysis of L1 retrotransposon (LINE1), 5 methyl-deoxycytidine (5m-dC) and 5 hydroxylmethyl-deoxycytidine (5hm-dC) in the discovery dataset. We validated a total of 1548 CpG sites (1307 genes) that were differentially methylated in GBM compared to controls. There were more than twice as many hypomethylated genes as hypermethylated ones. Both the discovery and validation datasets found 5 tumor methylation classes. Pathway analyses showed that the top ten pathways in hypomethylated genes were all related to functions of innate and acquired immunities. Among hypermethylated pathways, transcriptional regulatory network in embryonic stem cells was the most significant. In the study of global methylation markers, 5m-dC level was the best discriminant among methylation classes, whereas in survival analyses, high level of LINE1 methylation was an independent, favorable prognostic factor in the discovery dataset. Based on a pathway approach, hypermethylation in genes that control stem cell differentiation were significant, poor prognostic factors of overall survival in both the discovery and validation datasets. Approaches that targeted these methylated genes may be a future therapeutic goal.

OBJECTIVE: A recent high-density fine-mapping (ImmunoChip) study of genetic associations in rheumatoid arthritis (RA) identified 14 risk loci with validated genome-wide significance, as well as a number of loci showing associations suggestive of significance (P = 5 x 10(-5) < 5 x 10(-8) ), but these have yet to be replicated. The aim of this study was to determine whether these potentially significant loci are involved in the pathogenesis of RA, and to explore whether any of the loci are associated with a specific RA serotype. METHODS: A total of 16 single-nucleotide polymorphisms (SNPs) were selected for genotyping and association analyses in 2 independent validation cohorts, comprising 6,106 RA cases and 4,290 controls. A meta-analysis of the data from the original ImmunoChip discovery cohort and from both validation cohorts was carried out, for a combined total of 17,581 RA cases and 20,160 controls. In addition, stratified analysis of patient subsets, defined according to their anti-cyclic citrullinated peptide (anti-CCP) antibody status, was performed. RESULTS: A significant association with RA risk (P < 0.05) was replicated for 6 of the SNPs assessed in the validation cohorts. All SNPs in the validation study had odds ratios (ORs) for RA susceptibility in the same direction as those in the ImmunoChip discovery study. One SNP, rs72928038, mapping to an intron of BACH2, achieved genome-wide significance in the meta-analysis (P = 1.2 x 10(-8) , OR 1.12), and a second SNP, rs911263, mapping to an intron of RAD51B, was significantly associated in the anti-CCP-positive RA subgroup (P = 4 x 10(-8) , OR 0.89), confirming that both are RA susceptibility loci. CONCLUSION: This study provides robust evidence for an association of RA susceptibility with genes involved in B cell differentiation (BACH2) and DNA repair (RAD51B). The finding that the RAD51B gene exhibited different associations based on serologic subtype adds to the expanding knowledge base in defining subgroups of RA.

INTRODUCTION: Studies suggest that respiratory exposures including smoking, proximity to traffic and air pollution might be associated with development of rheumatoid arthritis (RA). RA-related autoantibodies are predictive of the development of RA. OBJECTIVE: We evaluated the relationship between RA-related autoantibodies and exposure to particulate matter (PM), a measure of air pollution of interest to health, in individuals without RA. METHODS: The Studies of the Etiology of Rheumatoid Arthritis (SERA) is a multicentre study following first-degree relatives (FDRs) of a proband with RA. FDRs are without the 1987 ACR (American College of Rheumatology) classifiable RA at enrolment and are followed for the development of RA-related autoimmunity. RA-related autoantibody outcomes as well as tender and swollen joint outcomes were assessed. Exposure to PM was assigned using ambient air pollution monitoring data and interpolated with inverse distance weighting spatial analyses using Geographic Information Systems. PM exposures were linked to FDR's residential zip codes. RESULTS: RA-related autoantibodies as well as tender or swollen joints are not associated with ambient PM concentrations. DISCUSSION: While other respiratory exposures may be associated with increased risk of RA, our data suggest that ambient PM is not associated with autoantibodies and joint signs among individuals without RA, but at increased risk of developing RA.

OBJECTIVE: To identify new genetic associations with juvenile and adult dermatomyositis (DM). METHODS: We performed a genome-wide association study (GWAS) of adult and juvenile DM patients of European ancestry (n = 1,178) and controls (n = 4,724). To assess genetic overlap with other autoimmune disorders, we examined whether 141 single-nucleotide polymorphisms (SNPs) outside the major histocompatibility complex (MHC) locus, and previously associated with autoimmune diseases, predispose to DM. RESULTS: Compared to controls, patients with DM had a strong signal in the MHC region consisting of GWAS-level significance (P < 5 x 10(-8)) at 80 genotyped SNPs. An analysis of 141 non-MHC SNPs previously associated with autoimmune diseases showed that 3 SNPs linked with 3 genes were associated with DM, with a false discovery rate (FDR) of <0.05. These genes were phospholipase C-like 1 (PLCL1; rs6738825, FDR = 0.00089), B lymphoid tyrosine kinase (BLK; rs2736340, FDR = 0.0031), and chemokine (C-C motif) ligand 21 (CCL21; rs951005, FDR = 0.0076). None of these genes was previously reported to be associated with DM. CONCLUSION: Our findings confirm the MHC as the major genetic region associated with DM and indicate that DM shares non-MHC genetic features with other autoimmune diseases, suggesting the presence of additional novel risk loci. This first identification of autoimmune disease genetic predispositions shared with DM may lead to enhanced understanding of pathogenesis and novel diagnostic and therapeutic approaches.

BACKGROUND: Synaesthesia is a neurodevelopmental condition in which a sensation in one modality triggers a perception in a second modality. Autism (shorthand for Autism Spectrum Conditions) is a neurodevelopmental condition involving social-communication disability alongside resistance to change and unusually narrow interests or activities. Whilst on the surface they appear distinct, they have been suggested to share common atypical neural connectivity. METHODS: In the present study, we carried out the first prevalence study of synaesthesia in autism to formally test whether these conditions are independent. After exclusions, 164 adults with autism and 97 controls completed a synaesthesia questionnaire, Autism Spectrum Quotient, and Test of Genuineness-Revised (ToG-R) online. RESULTS: The rate of synaesthesia in adults with autism was 18.9% (31 out of 164), almost three times greater than in controls (7.22%, 7 out of 97, P <0.05). ToG-R proved unsuitable for synaesthetes with autism. CONCLUSIONS: The significant increase in synaesthesia prevalence in autism suggests that the two conditions may share some common underlying mechanisms. Future research is needed to develop more feasible validation methods of synaesthesia in autism.