Ayse Bilge Ozel

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10−8). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4–1.88, P = 2.7 × 10−10, and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29–1.68, P = 4.9 × 10−9). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10−4; tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Integrated genomic sequencing reveals mutational landscape of T-cell prolymphocytic leukemia.

The comprehensive genetic alterations underlying the pathogenesis of T-cell prolymphocytic leukemia (T-PLL) are unknown. To address this, we performed whole-genome sequencing (WGS), whole-exome sequencing (WES), high-resolution copy-number analysis, and Sanger resequencing of a large cohort of T-PLL. WGS and WES identified novel mutations in recurrently altered genes not previously implicated in T-PLL including EZH2, FBXW10, and CHEK2. Strikingly, WGS and/or WES showed largely mutually exclusive mutations affecting IL2RG, JAK1, JAK3, or STAT5B in 38 of 50 T-PLL genomes (76.0%). Notably, gain-of-function IL2RG mutations are novel and have not been reported in any form of cancer. Further, high-frequency mutations in STAT5B have not been previously reported in T-PLL. Functionally, IL2RG-JAK1-JAK3-STAT5B mutations led to signal transducer and activator of transcription 5 (STAT5) hyperactivation, transformed Ba/F3 cells resulting in cytokine-independent growth, and/or enhanced colony formation in Jurkat T cells. Importantly, primary T-PLL cells exhibited constitutive activation of STAT5, and targeted pharmacologic inhibition of STAT5 with pimozide induced apoptosis in primary T-PLL cells. These results for the first time provide a portrait of the mutational landscape of T-PLL and implicate deregulation of DNA repair and epigenetic modulators as well as high-frequency mutational activation of the IL2RG-JAK1-JAK3-STAT5B axis in the pathogenesis of T-PLL. These findings offer opportunities for novel targeted therapies in this aggressive leukemia.

Linkage analysis identifies a locus for plasma von Willebrand factor undetected by genome-wide association

The plasma glycoprotein von Willebrand factor (VWF) exhibits fivefold antigen level variation across the normal human population determined by both genetic and environmental factors. Low levels of VWF are associated with bleeding and elevated levels with increased risk for thrombosis, myocardial infarction, and stroke. To identify additional genetic determinants of VWF antigen levels and to minimize the impact of age and illness-related environmental factors, we performed genome-wide association analysis in two young and healthy cohorts (n = 1,152 and n = 2,310) and identified signals at ABO (P < 7.9E-139) and VWF (P < 5.5E-16), consistent with previous reports. Additionally, linkage analysis based on sibling structure within the cohorts, identified significant signals at chromosome 2q12–2p13 (LOD score 5.3) and at the ABO locus on chromosome 9q34 (LOD score 2.9) that explained 19.2% and 24.5% of the variance in VWF levels, respectively. Given its strong effect, the linkage region on chromosome 2 could harbor a potentially important determinant of bleeding and thrombosis risk. The absence of a chromosome 2 association signal in this or previous association studies suggests a causative gene harboring many genetic variants that are individually rare, but in aggregate common. These results raise the possibility that similar loci could explain a significant portion of the “missing heritability” for other complex genetic traits.

Genome-Wide Analysis of Central Corneal Thickness in Primary Open-Angle Glaucoma Cases in the NEIGHBOR and GLAUGEN Consortia

PURPOSE To investigate the effects of central corneal thickness (CCT)-associated variants on primary open-angle glaucoma (POAG) risk using single nucleotide polymorphisms (SNP) data from the Glaucoma Genes and Environment (GLAUGEN) and National Eye Institute (NEI) Glaucoma Human Genetics Collaboration (NEIGHBOR) consortia. METHODS A replication analysis of previously reported CCT SNPs was performed in a CCT dataset (n = 1117) and these SNPs were then tested for association with POAG using a larger POAG dataset (n = 6470). Then a CCT genome-wide association study (GWAS) was performed. Top SNPs from this analysis were selected and tested for association with POAG. cDNA libraries from fetal and adult brain and ocular tissue samples were generated and used for candidate gene expression analysis. RESULTS Association with one of 20 previously published CCT SNPs was replicated: rs12447690, near the ZNF469 gene (P = 0.001; β = -5.08 μm/allele). None of these SNPs were significantly associated with POAG. In the CCT GWAS, no SNPs reached genome-wide significance. After testing 50 candidate SNPs for association with POAG, one SNP was identified, rs7481514 within the neurotrimin (NTM) gene, that was significantly associated with POAG in a low-tension subset (P = 0.00099; Odds Ratio [OR] = 1.28). Additionally, SNPs in the CNTNAP4 gene showed suggestive association with POAG (top SNP = rs1428758; P = 0.018; OR = 0.84). NTM and CNTNAP4 were shown to be expressed in ocular tissues. CONCLUSIONS The results suggest previously reported CCT loci are not significantly associated with POAG susceptibility. By performing a quantitative analysis of CCT and a subsequent analysis of POAG, SNPs in two cell adhesion molecules, NTM and CNTNAP4, were identified and may increase POAG susceptibility in a subset of cases.

Target concentration dependence of DNA melting temperature on oligonucleotide microarrays

The design of microarrays is currently based on studies focusing on DNA hybridization reaction in bulk solution. However, the presence of a surface to which the probe strand is attached can make the solution‐based approximations invalid, resulting in sub‐optimum hybridization conditions. To determine the effect of surfaces on DNA duplex formation, the authors studied the dependence of DNA melting temperature (Tm) on target concentration. An automated system was developed to capture the melting profiles of a 25‐mer perfect‐match probe–target pair initially hybridized at 23°C. Target concentrations ranged from 0.0165 to 15 nM with different probe amounts (0.03–0.82 pmol on a surface area of 1018 Å2), a constant probe density (5 × 1012 molecules/cm2) and spacer length (15 dT). The authors found that Tm for duplexes anchored to a surface is lower than in‐solution, and this difference increases with increasing target concentration. In a representative set, a target concentration increase from 0.5 to 15 nM with 0.82 pmol of probe on the surface resulted in a Tm decrease of 6°C when compared with a 4°C increase in solution. At very low target concentrations, a multi‐melting process was observed in low temperature domains of the curves. This was attributed to the presence of truncated or mismatch probes.

Genome-wide studies of von Willebrand factor propeptide identify loci contributing to variation in propeptide levels and von Willebrand factor clearance.

Essentials Variants at ABO, von Willebrand Factor (VWF) and 2q12 contribute to the variation in plasma in VWF. We performed a genome‐wide association study of plasma VWF propeptide in 3,238 individuals. ABO, VWF and 2q12 loci had weak or no association or linkage with plasma VWFpp levels. VWF associated variants at ABO, VWF and 2q12 loci primarily affect VWF clearance rates.

Mutations in VPS13D lead to a new recessive ataxia with spasticity and mitochondrial defects

To identify novel causes of recessive ataxias, including spinocerebellar ataxia with saccadic intrusions, spastic ataxias, and spastic paraplegia.

Genetic variants in ADAMTS13 as well as smoking are major determinants of plasma ADAMTS13 levels

The metalloprotease ADAMTS13 cleaves von Willebrand factor (VWF) in circulating blood, limiting the size of VWF multimers and regulating VWF activity. Abnormal regulation of VWF contributes to bleeding and to thrombotic disorders. ADAMTS13 levels in plasma are highly variable among healthy individuals, although the heritability and the genetic determinants of this variation are unclear. We performed genome-wide association studies of plasma ADAMTS13 concentrations in 3244 individuals from 2 independent cohorts of healthy individuals. The heritability of ADAMTS13 levels was between 59.1% (all individuals) and 83.5% (siblings only), whereas tobacco smoking was associated with a decrease in plasma ADAMTS13 levels. Meta-analysis identified common variants near the ADAMTS13 locus on chromosome 9q34.2 that were significantly associated with ADAMTS13 levels and collectively explained 20.0% of the variance. The top single nucleotide polymorphism (SNP), rs28673647, resides in an intron of ADAMTS13 (β, 6.7%; P = 1.3E-52). Conditional analysis revealed 3 additional independent signals represented by rs3739893 (β, -22.3%; P = 1.2E-30) and rs3124762 (β, 3.5%; P = 8.9E-9) close to ADAMTS13 and rs4075970 (β, 2.4%; P = 6.8E-9) on 21q22.3. Linkage analysis also identified the region around ADAMTS13 (9q34.2) as the top signal (LOD 3.5), consistent with our SNP association analyses. Two nonsynonymous ADAMTS13 variants in the top 2 independent linkage disequilibrium blocks (Q448E and A732V) were identified and characterized in vitro. This study uncovered specific common genetic polymorphisms that are key genetic determinants of the variation in plasma ADAMTS13 levels in healthy individuals.

Extended regions of suspected mis-assembly in the rat reference genome

We performed whole-genome sequencing for eight inbred rat strains commonly used in genetic mapping studies, and they are the founders of the NIH heterogeneous stock (HS) outbred colony. We provide their sequences and variant calls to the rat genomics community. When analyzing the variant calls we identified regions with unusually high heterozygosity. We show that these regions are consistent across the eight inbred strains, including the BN strain, which was the basis of the rat reference genome. These regions show significantly higher read depth than other regions in the genome. The evidence suggests that these regions may contain segmental duplications that are incorrectly overlaid in the reference genome. We provide masks for these suspected regions of mis-assembly as a resource for the community to flag potentially false interpretations of mapping results or functional data.