Ernest Turro

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.

Haplotype and isoform specific expression estimation using multi-mapping RNA-seq reads

We present a novel pipeline and methodology for simultaneously estimating isoform expression and allelic imbalance in diploid organisms using RNA-seq data. We achieve this by modeling the expression of haplotype-specific isoforms. If unknown, the two parental isoform sequences can be individually reconstructed. A new statistical method, MMSEQ, deconvolves the mapping of reads to multiple transcripts (isoforms or haplotype-specific isoforms). Our software can take into account non-uniform read generation and works with paired-end reads.

Transcriptional diversity during lineage commitment of human blood progenitors

Introduction Blood production in humans culminates in the daily release of around 1011 cells into the circulation, mainly platelets and red blood cells. All blood cells originate from a minute population of hematopoietic stem cells (HSCs) that expands and differentiates into progenitor cells with increasingly restricted lineage choice. Characterizing alternative splicing events involved in hematopoiesis is critical for interpreting the effects of mutations leading to inherited disorders and blood cancers and for the rational design of strategies to advance transplantation and regenerative medicine. Overview of methodology. RNA-sequencing reads from human blood progenitors [opaque cells in (A)] were mapped to the transcriptome to quantify gene and transcript expression. Reads were also mapped to the genome to identify novel splice junctions and characterize alternative splicing events (B). Rationale To address this, we explored the transcriptional diversity of human blood progenitors by sequencing RNA from six progenitor and two precursor populations representing the classical myeloid commitment stages of hematopoiesis and the main lymphoid stage. Data were aligned to the human reference transcriptome and genome to quantify known transcript isoforms and to identify novel splicing events, respectively. We used Bayesian polytomous model selection to classify transcripts into distinct expression patterns across the three cell types that comprise each differentiation step. Results We identified extensive transcriptional changes involving 6711 genes and 10,724 transcripts and validated a number of these. Many of the changes at the transcript isoform level did not result in significant changes at the gene expression level. Moreover, we identified transcripts unique to each of the progenitor populations, observing enrichment in non–protein-coding elements at the early stages of differentiation. We discovered 7881 novel splice junctions and 2301 differentially used alternative splicing events, enriched in genes involved in regulatory processes and often resulting in the gain or loss of functional domains. Of the alternative splice sites displaying differential usage, 73% resulted in exon-skipping events involving at least one protein domain (38.5%) or introducing a premature stop codon (26%). Enrichment analysis of RNA-binding motifs provided insights into the regulation of cell type–specific splicing events. To demonstrate the importance of specific isoforms in driving lineage fating events, we investigated the role of a transcription factor highlighted by our analyses. Our data show that nuclear factor I/B (NFIB) is highly expressed in megakaryocytes and that it is transcribed from an unannotated transcription start site preceding a novel exon. The novel NFIB isoform lacks the DNA binding/dimerization domain and therefore is unable to interact with its binding partner, NFIC. We further show that NFIB and NFIC are important in megakaryocyte differentiation. Conclusion We produced a quantitative catalog of transcriptional changes and splicing events representing the early progenitors of human blood. Our analyses unveil a previously undetected layer of regulation affecting cell fating, which involves transcriptional isoforms switching without noticeable changes at the gene level and resulting in the gain or loss of protein functions. A BLUEPRINT of immune cell development To determine the epigenetic mechanisms that direct blood cells to develop into the many components of our immune system, the BLUEPRINT consortium examined the regulation of DNA and RNA transcription to dissect the molecular traits that govern blood cell differentiation. By inducing immune responses, Saeed et al. document the epigenetic changes in the genome that underlie immune cell differentiation. Cheng et al. demonstrate that trained monocytes are highly dependent on the breakdown of sugars in the presence of oxygen, which allows cells to produce the energy needed to mount an immune response. Chen et al. examine RNA transcripts and find that specific cell lineages use RNA transcripts of different length and composition (isoforms) to form proteins. Together, the studies reveal how epigenetic effects can drive the development of blood cells involved in the immune system. Science, this issue 10.1126/science.1251086, 10.1126/science.1250684, 10.1126/science.1251033 RNA sequencing identifies how different cell fate decisions are made during blood cell differentiation. Blood cells derive from hematopoietic stem cells through stepwise fating events. To characterize gene expression programs driving lineage choice, we sequenced RNA from eight primary human hematopoietic progenitor populations representing the major myeloid commitment stages and the main lymphoid stage. We identified extensive cell type–specific expression changes: 6711 genes and 10,724 transcripts, enriched in non–protein-coding elements at early stages of differentiation. In addition, we found 7881 novel splice junctions and 2301 differentially used alternative splicing events, enriched in genes involved in regulatory processes. We demonstrated experimentally cell-specific isoform usage, identifying nuclear factor I/B (NFIB) as a regulator of megakaryocyte maturation—the platelet precursor. Our data highlight the complexity of fating events in closely related progenitor populations, the understanding of which is essential for the advancement of transplantation and regenerative medicine.

Extensive compensatory cis-trans regulation in the evolution of mouse gene expression

Gene expression levels are thought to diverge primarily via regulatory mutations in trans within species, and in cis between species. To test this hypothesis in mammals we used RNA-sequencing to measure gene expression divergence between C57BL/6J and CAST/EiJ mouse strains and allele-specific expression in their F1 progeny. We identified 535 genes with parent-of-origin specific expression patterns, although few of these showed full allelic silencing. This suggests that the number of imprinted genes in a typical mouse somatic tissue is relatively small. In the set of nonimprinted genes, 32% showed evidence of divergent expression between the two strains. Of these, 2% could be attributed purely to variants acting in trans, while 43% were attributable only to variants acting in cis. The genes with expression divergence driven by changes in trans showed significantly higher sequence constraint than genes where the divergence was explained by variants acting in cis. The remaining genes with divergent patterns of expression (55%) were regulated by a combination of variants acting in cis and variants acting in trans. Intriguingly, the changes in expression induced by the cis and trans variants were in opposite directions more frequently than expected by chance, implying that compensatory regulation to stabilize gene expression levels is widespread. We propose that expression levels of genes regulated by this mechanism are fine-tuned by cis variants that arise following regulatory changes in trans, suggesting that many cis variants are not the primary targets of natural selection.

Extensive Co-Operation between the Epstein-Barr Virus EBNA3 Proteins in the Manipulation of Host Gene Expression and Epigenetic Chromatin Modification

Epstein-Barr virus (EBV) is able to drive the transformation of B-cells, resulting in the generation of lymphoblastoid cell lines (LCLs) in vitro. EBV nuclear proteins EBNA3A and EBNA3C are necessary for efficient transformation, while EBNA3B is dispensable. We describe a transcriptome analysis of BL31 cells infected with a series of EBNA3-knockout EBVs, including one deleted for all three EBNA3 genes. Using Affymetrix Exon 1.0 ST microarrays analysed with the MMBGX algorithm, we have identified over 1000 genes whose regulation by EBV requires one of the EBNA3s. Remarkably, a third of the genes identified require more than one EBNA3 for their regulation, predominantly EBNA3C co-operating with either EBNA3B, EBNA3A or both. The microarray was validated by real-time PCR, while ChIP analysis of a selection of co-operatively repressed promoters indicates a role for polycomb group complexes. Targets include genes involved in apoptosis, cell migration and B-cell differentiation, and show a highly significant but subtle alteration in genes involved in mitosis. In order to assess the relevance of the BL31 system to LCLs, we analysed the transcriptome of a set of EBNA3B knockout (3BKO) LCLs. Around a third of the genes whose expression level in LCLs was altered in the absence of EBNA3B were also altered in 3BKO-BL31 cell lines. Among these are TERT and TCL1A, implying that EBV-induced changes in the expression of these genes are not required for B-cell transformation. We also identify 26 genes that require both EBNA3A and EBNA3B for their regulation in LCLs. Together, this shows the complexity of the interaction between EBV and its host, whereby multiple EBNA3 proteins co-operate to modulate the behaviour of the host cell.

The South Asian genome.

The genetic sequence variation of people from the Indian subcontinent who comprise one-quarter of the worlds population, is not well described. We carried out whole genome sequencing of 168 South Asians, along with whole-exome sequencing of 147 South Asians to provide deeper characterisation of coding regions. We identify 12,962,155 autosomal sequence variants, including 2,946,861 new SNPs and 312,738 novel indels. This catalogue of SNPs and indels amongst South Asians provides the first comprehensive map of genetic variation in this major human population, and reveals evidence for selective pressures on genes involved in skin biology, metabolism, infection and immunity. Our results will accelerate the search for the genetic variants underlying susceptibility to disorders such as type-2 diabetes and cardiovascular disease which are highly prevalent amongst South Asians.

Human phenotype ontology annotation and cluster analysis to unravel genetic defects in 707 cases with unexplained bleeding and platelet disorders.

BackgroundHeritable bleeding and platelet disorders (BPD) are heterogeneous and frequently have an unknown genetic basis. The BRIDGE-BPD study aims to discover new causal genes for BPD by high throughput sequencing using cluster analyses based on improved and standardised deep, multi-system phenotyping of cases.MethodsWe report a new approach in which the clinical and laboratory characteristics of BPD cases are annotated with adapted Human Phenotype Ontology (HPO) terms. Cluster analyses are then used to characterise groups of cases with similar HPO terms and variants in the same genes.ResultsWe show that 60% of index cases with heritable BPD enrolled at 10 European or US centres were annotated with HPO terms indicating abnormalities in organ systems other than blood or blood-forming tissues, particularly the nervous system. Cases within pedigrees clustered closely together on the bases of their HPO-coded phenotypes, as did cases sharing several clinically suspected syndromic disorders. Cases subsequently found to harbour variants in ACTN1 also clustered closely, even though diagnosis of this recently described disorder was not possible using only the clinical and laboratory data available to the enrolling clinician.ConclusionsThese findings validate our novel HPO-based phenotype clustering methodology for known BPD, thus providing a new discovery tool for BPD of unknown genetic basis. This approach will also be relevant for other rare diseases with significant genetic heterogeneity.

A high-throughput sequencing test for diagnosing inherited bleeding, thrombotic, and platelet disorders

Inherited bleeding, thrombotic, and platelet disorders (BPDs) are diseases that affect ∼300 individuals per million births. With the exception of hemophilia and von Willebrand disease patients, a molecular analysis for patients with a BPD is often unavailable. Many specialized tests are usually required to reach a putative diagnosis and they are typically performed in a step-wise manner to control costs. This approach causes delays and a conclusive molecular diagnosis is often never reached, which can compromise treatment and impede rapid identification of affected relatives. To address this unmet diagnostic need, we designed a high-throughput sequencing platform targeting 63 genes relevant for BPDs. The platform can call single nucleotide variants, short insertions/deletions, and large copy number variants (though not inversions) which are subjected to automated filtering for diagnostic prioritization, resulting in an average of 5.34 candidate variants per individual. We sequenced 159 and 137 samples, respectively, from cases with and without previously known causal variants. Among the latter group, 61 cases had clinical and laboratory phenotypes indicative of a particular molecular etiology, whereas the remainder had an a priori highly uncertain etiology. All previously detected variants were recapitulated and, when the etiology was suspected but unknown or uncertain, a molecular diagnosis was reached in 56 of 61 and only 8 of 76 cases, respectively. The latter category highlights the need for further research into novel causes of BPDs. The ThromboGenomics platform thus provides an affordable DNA-based test to diagnose patients suspected of having a known inherited BPD.

A gain-of-function variant in DIAPH1 causes dominant macrothrombocytopenia and hearing loss

Macrothrombocytopenia (MTP) is a heterogeneous group of disorders characterized by enlarged and reduced numbers of circulating platelets, sometimes resulting in abnormal bleeding. In most MTP, this phenotype arises because of altered regulation of platelet formation from megakaryocytes (MKs). We report the identification of DIAPH1, which encodes the Rho-effector diaphanous-related formin 1 (DIAPH1), as a candidate gene for MTP using exome sequencing, ontological phenotyping, and similarity regression. We describe 2 unrelated pedigrees with MTP and sensorineural hearing loss that segregate with a DIAPH1 R1213* variant predicting partial truncation of the DIAPH1 diaphanous autoregulatory domain. The R1213* variant was linked to reduced proplatelet formation from cultured MKs, cell clustering, and abnormal cortical filamentous actin. Similarly, in platelets, there was increased filamentous actin and stable microtubules, indicating constitutive activation of DIAPH1. Overexpression of DIAPH1 R1213* in cells reproduced the cytoskeletal alterations found in platelets. Our description of a novel disorder of platelet formation and hearing loss extends the repertoire of DIAPH1-related disease and provides new insight into the autoregulation of DIAPH1 activity.

A dominant gain-of-function mutation in universal tyrosine kinase SRC causes thrombocytopenia, myelofibrosis, bleeding, and bone pathologies

E527K hyperactive SRC results in megakaryocytes with increased podosome formation, thrombocytopenia, myelofibrosis, bleeding, and bone pathologies. SRC shows its stripes The nonreceptor tyrosine kinase SRC is a proto-oncogene that has been associated with cancer progression. Now, Turro et al. find a gain-of-function mutation in SRC in nine patients with myelofibrosis, bleeding, and bone disorders. This mutation prevented SRC from inhibiting itself, and the overactive SRC resulted in enhanced tyrosine phosphorylation in a zebrafish model as well as in patient-derived cells. In patients with myelofibrosis, this SRC mutation was associated with increased outgrowth of myeloid and megakaryocyte colonies, with abnormal platelet production, which could be rescued by SRC kinase inhibition. These findings may be important for understanding the severe bleeding in cancer patients treated with Src family kinase inhibitors. The Src family kinase (SFK) member SRC is a major target in drug development because it is activated in many human cancers, yet deleterious SRC germline mutations have not been reported. We used genome sequencing and Human Phenotype Ontology patient coding to identify a gain-of-function mutation in SRC causing thrombocytopenia, myelofibrosis, bleeding, and bone pathologies in nine cases. Modeling of the E527K substitution predicts loss of SRC’s self-inhibitory capacity, which we confirmed with in vitro studies showing increased SRC kinase activity and enhanced Tyr419 phosphorylation in COS-7 cells overexpressing E527K SRC. The active form of SRC predominates in patients’ platelets, resulting in enhanced overall tyrosine phosphorylation. Patients with myelofibrosis have hypercellular bone marrow with trilineage dysplasia, and their stem cells grown in vitro form more myeloid and megakaryocyte (MK) colonies than control cells. These MKs generate platelets that are dysmorphic, low in number, highly variable in size, and have a paucity of α-granules. Overactive SRC in patient-derived MKs causes a reduction in proplatelet formation, which can be rescued by SRC kinase inhibition. Stem cells transduced with lentiviral E527K SRC form MKs with a similar defect and enhanced tyrosine phosphorylation levels. Patient-derived and E527K-transduced MKs show Y419 SRC–positive stained podosomes that induce altered actin organization. Expression of mutated src in zebrafish recapitulates patients’ blood and bone phenotypes. Similar studies of platelets and MKs may reveal the mechanism underlying the severe bleeding frequently observed in cancer patients treated with next-generation SFK inhibitors.