Nikolas Pontikos

The Human Phenotype Ontology in 2017

Deep phenotyping has been defined as the precise and comprehensive analysis of phenotypic abnormalities in which the individual components of the phenotype are observed and described. The three components of the Human Phenotype Ontology (HPO; www.human-phenotype-ontology.org) project are the phenotype vocabulary, disease-phenotype annotations and the algorithms that operate on these. These components are being used for computational deep phenotyping and precision medicine as well as integration of clinical data into translational research. The HPO is being increasingly adopted as a standard for phenotypic abnormalities by diverse groups such as international rare disease organizations, registries, clinical labs, biomedical resources, and clinical software tools and will thereby contribute toward nascent efforts at global data exchange for identifying disease etiologies. This update article reviews the progress of the HPO project since the debut Nucleic Acids Research database article in 2014, including specific areas of expansion such as common (complex) disease, new algorithms for phenotype driven genomic discovery and diagnostics, integration of cross-species mapping efforts with the Mammalian Phenotype Ontology, an improved quality control pipeline, and the addition of patient-friendly terminology.

Standardizing Flow Cytometry Immunophenotyping Analysis from the Human ImmunoPhenotyping Consortium

Standardization of immunophenotyping requires careful attention to reagents, sample handling, instrument setup, and data analysis, and is essential for successful cross-study and cross-center comparison of data. Experts developed five standardized, eight-color panels for identification of major immune cell subsets in peripheral blood. These were produced as pre-configured, lyophilized, reagents in 96-well plates. We present the results of a coordinated analysis of samples across nine laboratories using these panels with standardized operating procedures (SOPs). Manual gating was performed by each site and by a central site. Automated gating algorithms were developed and tested by the FlowCAP consortium. Centralized manual gating can reduce cross-center variability, and we sought to determine whether automated methods could streamline and standardize the analysis. Within-site variability was low in all experiments, but cross-site variability was lower when central analysis was performed in comparison with site-specific analysis. It was also lower for clearly defined cell subsets than those based on dim markers and for rare populations. Automated gating was able to match the performance of central manual analysis for all tested panels, exhibiting little to no bias and comparable variability. Standardized staining, data collection, and automated gating can increase power, reduce variability, and streamline analysis for immunophenotyping.

Mutations in REEP6 Cause Autosomal-Recessive Retinitis Pigmentosa

Retinitis pigmentosa (RP) is the most frequent form of inherited retinal dystrophy. RP is genetically heterogeneous and the genes identified to date encode proteins involved in a wide range of functional pathways, including photoreceptor development, phototransduction, the retinoid cycle, cilia, and outer segment development. Here we report the identification of biallelic mutations in Receptor Expression Enhancer Protein 6 (REEP6) in seven individuals with autosomal-recessive RP from five unrelated families. REEP6 is a member of the REEP/Yop1 family of proteins that influence the structure of the endoplasmic reticulum but is relatively unstudied. The six variants identified include three frameshift variants, two missense variants, and a genomic rearrangement that disrupts exon 1. Human 3D organoid optic cups were used to investigate REEP6 expression and confirmed the expression of a retina-specific isoform REEP6.1, which is specifically affected by one of the frameshift mutations. Expression of the two missense variants (c.383C>T [p.Pro128Leu] and c.404T>C [p.Leu135Pro]) and the REEP6.1 frameshift mutant in cultured cells suggest that these changes destabilize the protein. Furthermore, CRISPR-Cas9-mediated gene editing was used to produce Reep6 knock-in mice with the p.Leu135Pro RP-associated variant identified in one RP-affected individual. The homozygous knock-in mice mimic the clinical phenotypes of RP, including progressive photoreceptor degeneration and dysfunction of the rod photoreceptors. Therefore, our study implicates REEP6 in retinal homeostasis and highlights a pathway previously uncharacterized in retinal dystrophy.

Autosomal-Dominant Corneal Endothelial Dystrophies CHED1 and PPCD1 Are Allelic Disorders Caused by Non-coding Mutations in the Promoter of OVOL2

Congenital hereditary endothelial dystrophy 1 (CHED1) and posterior polymorphous corneal dystrophy 1 (PPCD1) are autosomal-dominant corneal endothelial dystrophies that have been genetically mapped to overlapping loci on the short arm of chromosome 20. We combined genetic and genomic approaches to identify the cause of disease in extensive pedigrees comprising over 100 affected individuals. After exclusion of pathogenic coding, splice-site, and copy-number variations, a parallel approach using targeted and whole-genome sequencing facilitated the identification of pathogenic variants in a conserved region of the OVOL2 proximal promoter sequence in the index families (c.−339_361dup for CHED1 and c.−370T>C for PPCD1). Direct sequencing of the OVOL2 promoter in other unrelated affected individuals identified two additional mutations within the conserved proximal promoter sequence (c.−274T>G and c.−307T>C). OVOL2 encodes ovo-like zinc finger 2, a C2H2 zinc-finger transcription factor that regulates mesenchymal-to-epithelial transition and acts as a direct transcriptional repressor of the established PPCD-associated gene ZEB1. Interestingly, we did not detect OVOL2 expression in the normal corneal endothelium. Our in vitro data demonstrate that all four mutated OVOL2 promoters exhibited more transcriptional activity than the corresponding wild-type promoter, and we postulate that the mutations identified create cryptic cis-acting regulatory sequence binding sites that drive aberrant OVOL2 expression during endothelial cell development. Our data establish CHED1 and PPCD1 as allelic conditions and show that CHED1 represents the extreme of what can be considered a disease spectrum. They also implicate transcriptional dysregulation of OVOL2 as a common cause of dominantly inherited corneal endothelial dystrophies.

Dissection of a complex disease susceptibility region using a Bayesian stochastic search approach to fine mapping

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.

Functional variants in the LRRK2 gene confer shared effects on risk for Crohn’s disease and Parkinson’s disease

Crohn’s disease (CD)–associated variants in the LRRK2 gene for risk (N2081D) and for protection (N551K) mediate shared effects in CD and Parkinson’s disease. A shared history Crohn’s disease (CD), an inflammatory bowel disease, has a relatively high prevalence in Ashkenazi Jewish populations. Hui et al. conducted genome-wide association analysis in 2066 CD patients and 3633 healthy control individuals of Ashkenazi Jewish ancestry and identified two functional variants in the LRRK2 gene. The LRRK2 gene has been previously linked to the development of Parkinson’s disease (PD). The new LRRK2 variants conferred risk for CD (N2081D) or protection from CD (N551K/R1398H). Analysis of other variants within the LRRK2 locus in 24,570 individuals revealed similar genetic effects between CD and PD in both Ashkenazi Jewish and non-Jewish cohorts. The presence of shared LRRK2 alleles in CD and PD provides insight into disease mechanisms and potential treatments. Crohn’s disease (CD), a form of inflammatory bowel disease, has a higher prevalence in Ashkenazi Jewish than in non-Jewish European populations. To define the role of nonsynonymous mutations, we performed exome sequencing of Ashkenazi Jewish patients with CD, followed by array-based genotyping and association analysis in 2066 CD cases and 3633 healthy controls. We detected association signals in the LRRK2 gene that conferred risk for CD (N2081D variant, P = 9.5 × 10−10) or protection from CD (N551K variant, tagging R1398H-associated haplotype, P = 3.3 × 10−8). These variants affected CD age of onset, disease location, LRRK2 activity, and autophagy. Bayesian network analysis of CD patient intestinal tissue further implicated LRRK2 in CD pathogenesis. Analysis of the extended LRRK2 locus in 24,570 CD cases, patients with Parkinson’s disease (PD), and healthy controls revealed extensive pleiotropy, with shared genetic effects between CD and PD in both Ashkenazi Jewish and non-Jewish cohorts. The LRRK2 N2081D CD risk allele is located in the same kinase domain as G2019S, a mutation that is the major genetic cause of familial and sporadic PD. Like the G2019S mutation, the N2081D variant was associated with increased kinase activity, whereas neither N551K nor R1398H variants on the protective haplotype altered kinase activity. We also confirmed that R1398H, but not N551K, increased guanosine triphosphate binding and hydrolyzing enzyme (GTPase) activity, thereby deactivating LRRK2. The presence of shared LRRK2 alleles in CD and PD provides refined insight into disease mechanisms and may have major implications for the treatment of these two seemingly unrelated diseases.

Genetic Complexity of Crohn's Disease in Two Large Ashkenazi Jewish Families.

Background & Aims Crohn’s disease (CD) is a highly heritable disease that is particularly common in the Ashkenazi Jewish population. We studied 2 large Ashkenazi Jewish families with a high prevalence of CD in an attempt to identify novel genetic risk variants. Methods Ashkenazi Jewish patients with CD and a positive family history were recruited from the University College London Hospital. We used genome-wide, single-nucleotide polymorphism data to assess the burden of common CD-associated risk variants and for linkage analysis. Exome sequencing was performed and rare variants that were predicted to be deleterious and were observed at a high frequency in cases were prioritized. We undertook within-family association analysis after imputation and assessed candidate variants for evidence of association with CD in an independent cohort of Ashkenazi Jewish individuals. We examined the effects of a variant in DUOX2 on hydrogen peroxide production in HEK293 cells. Results We identified 2 families (1 with >800 members and 1 with >200 members) containing 54 and 26 cases of CD or colitis, respectively. Both families had a significant enrichment of previously described common CD-associated risk variants. No genome-wide significant linkage was observed. Exome sequencing identified candidate variants, including a missense mutation in DUOX2 that impaired its function and a frameshift mutation in CSF2RB that was associated with CD in an independent cohort of Ashkenazi Jewish individuals. Conclusions In a study of 2 large Ashkenazi Jewish with multiple cases of CD, we found the genetic basis of the disease to be complex, with a role for common and rare genetic variants. We identified a frameshift mutation in CSF2RB that was replicated in an independent cohort. These findings show the value of family studies and the importance of the innate immune system in the pathogenesis of CD.

A hybrid qPCR/SNP array approach allows cost efficient assessment of KIR gene copy numbers in large samples

BackgroundKiller Immunoglobulin-like Receptors (KIRs) are surface receptors of natural killer cells that bind to their corresponding Human Leukocyte Antigen (HLA) class I ligands, making them interesting candidate genes for HLA-associated autoimmune diseases, including type 1 diabetes (T1D). However, allelic and copy number variation in the KIR region effectively mask it from standard genome-wide association studies: single nucleotide polymorphism (SNP) probes targeting the region are often discarded by standard genotype callers since they exhibit variable cluster numbers. Quantitative Polymerase Chain Reaction (qPCR) assays address this issue. However, their cost is prohibitive at the sample sizes required for detecting effects typically observed in complex genetic diseases.ResultsWe propose a more powerful and cost-effective alternative, which combines signals from SNPs with more than three clusters found in existing datasets, with qPCR on a subset of samples. First, we showed that noise and batch effects in multiplexed qPCR assays are addressed through normalisation and simultaneous copy number calling of multiple genes. Then, we used supervised classification to impute copy numbers of specific KIR genes from SNP signals. We applied this method to assess copy number variation in two KIR genes, KIR3DL1 and KIR3DS1, which are suitable candidates for T1D susceptibility since they encode the only KIR molecules known to bind with HLA-Bw4 epitopes. We find no association between KIR3DL1/3DS1 copy number and T1D in 6744 cases and 5362 controls; a sample size twenty-fold larger than in any previous KIR association study. Due to our sample size, we can exclude odds ratios larger than 1.1 for the common KIR3DL1/3DS1 copy number groups at the 5% significance level.ConclusionWe found no evidence of association of KIR3DL1/3DS1 copy number with T1D, either overall or dependent on HLA-Bw4 epitope. Five other KIR genes, KIR2DS4, KIR2DL3, KIR2DL5, KIR2DS5 and KIR2DS1, in high linkage disequilibrium with KIR3DL1 and KIR3DS1, are also unlikely to be significantly associated. Our approach could potentially be applied to other KIR genes to allow cost effective assaying of gene copy number in large samples.

The Human Salivary Microbiome Is Shaped by Shared Environment Rather than Genetics: Evidence from a Large Family of Closely Related Individuals

ABSTRACT The human microbiome is affected by multiple factors, including the environment and host genetics. In this study, we analyzed the salivary microbiomes of an extended family of Ashkenazi Jewish individuals living in several cities and investigated associations with both shared household and host genetic similarities. We found that environmental effects dominated over genetic effects. While there was weak evidence of geographical structuring at the level of cities, we observed a large and significant effect of shared household on microbiome composition, supporting the role of the immediate shared environment in dictating the presence or absence of taxa. This effect was also seen when including adults who had grown up in the same household but moved out prior to the time of sampling, suggesting that the establishment of the salivary microbiome earlier in life may affect its long-term composition. We found weak associations between host genetic relatedness and microbiome dissimilarity when using family pedigrees as proxies for genetic similarity. However, this association disappeared when using more-accurate measures of kinship based on genome-wide genetic markers, indicating that the environment rather than host genetics is the dominant factor affecting the composition of the salivary microbiome in closely related individuals. Our results support the concept that there is a consistent core microbiome conserved across global scales but that small-scale effects due to a shared living environment significantly affect microbial community composition. IMPORTANCE Previous research shows that the salivary microbiomes of relatives are more similar than those of nonrelatives, but it remains difficult to distinguish the effects of relatedness and shared household environment. Furthermore, pedigree measures may not accurately measure host genetic similarity. In this study, we include genetic relatedness based on genome-wide single nucleotide polymorphisms (SNPs) (rather than pedigree measures) and shared environment in the same analysis. We quantify the relative importance of these factors by studying the salivary microbiomes in members of a large extended Ashkenazi Jewish family living in different locations. We find that host genetics plays no significant role and that the dominant factor is the shared environment at the household level. We also find that this effect appears to persist in individuals who have moved out of the parental household, suggesting that aspects of salivary microbiome composition established during upbringing can persist over a time scale of years. IMPORTANCE Previous research shows that the salivary microbiomes of relatives are more similar than those of nonrelatives, but it remains difficult to distinguish the effects of relatedness and shared household environment. Furthermore, pedigree measures may not accurately measure host genetic similarity. In this study, we include genetic relatedness based on genome-wide single nucleotide polymorphisms (SNPs) (rather than pedigree measures) and shared environment in the same analysis. We quantify the relative importance of these factors by studying the salivary microbiomes in members of a large extended Ashkenazi Jewish family living in different locations. We find that host genetics plays no significant role and that the dominant factor is the shared environment at the household level. We also find that this effect appears to persist in individuals who have moved out of the parental household, suggesting that aspects of salivary microbiome composition established during upbringing can persist over a time scale of years.

Mutations in the Spliceosome Component CWC27 Cause Retinal Degeneration with or without Additional Developmental Anomalies

Pre-mRNA splicing factors play a fundamental role in regulating transcript diversity both temporally and spatially. Genetic defects in several spliceosome components have been linked to a set of non-overlapping spliceosomopathy phenotypes in humans, among which skeletal developmental defects and non-syndromic retinitis pigmentosa (RP) are frequent findings. Here we report that defects in spliceosome-associated protein CWC27 are associated with a spectrum of disease phenotypes ranging from isolated RP to severe syndromic forms. By whole-exome sequencing, recessive protein-truncating mutations in CWC27 were found in seven unrelated families that show a range of clinical phenotypes, including retinal degeneration, brachydactyly, craniofacial abnormalities, short stature, and neurological defects. Remarkably, variable expressivity of the human phenotype can be recapitulated in Cwc27 mutant mouse models, with significant embryonic lethality and severe phenotypes in the complete knockout mice while mice with a partial loss-of-function allele mimic the isolated retinal degeneration phenotype. Our study describes a retinal dystrophy-related phenotype spectrum as well as its genetic etiology and highlights the complexity of the spliceosomal gene network.