Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations.

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

The Human Phenotype Ontology project: linking molecular biology and disease through phenotype data.

The Human Phenotype Ontology (HPO) project, available at http://www.human-phenotype-ontology.org, provides a structured, comprehensive and well-defined set of 10,088 classes (terms) describing human phenotypic abnormalities and 13,326 subclass relations between the HPO classes. In addition we have developed logical definitions for 46% of all HPO classes using terms from ontologies for anatomy, cell types, function, embryology, pathology and other domains. This allows interoperability with several resources, especially those containing phenotype information on model organisms such as mouse and zebrafish. Here we describe the updated HPO database, which provides annotations of 7,278 human hereditary syndromes listed in OMIM, Orphanet and DECIPHER to classes of the HPO. Various meta-attributes such as frequency, references and negations are associated with each annotation. Several large-scale projects worldwide utilize the HPO for describing phenotype information in their datasets. We have therefore generated equivalence mappings to other phenotype vocabularies such as LDDB, Orphanet, MedDRA, UMLS and phenoDB, allowing integration of existing datasets and interoperability with multiple biomedical resources. We have created various ways to access the HPO database content using flat files, a MySQL database, and Web-based tools. All data and documentation on the HPO project can be found online.

A GDF5 point mutation strikes twice--causing BDA1 and SYNS2.

Growth and Differentiation Factor 5 (GDF5) is a secreted growth factor that belongs to the Bone Morphogenetic Protein (BMP) family and plays a pivotal role during limb development. GDF5 is a susceptibility gene for osteoarthritis (OA) and mutations in GDF5 are associated with a wide variety of skeletal malformations ranging from complex syndromes such as acromesomelic chondrodysplasias to isolated forms of brachydactylies or multiple synostoses syndrome 2 (SYNS2). Here, we report on a family with an autosomal dominant inherited combination of SYNS2 and additional brachydactyly type A1 (BDA1) caused by a single point mutation in GDF5 (p.W414R). Functional studies, including chondrogenesis assays with primary mesenchymal cells, luciferase reporter gene assays and Surface Plasmon Resonance analysis, of the GDF5(W414R) variant in comparison to other GDF5 mutations associated with isolated BDA1 (p.R399C) or SYNS2 (p.E491K) revealed a dual pathomechanism characterized by a gain- and loss-of-function at the same time. On the one hand insensitivity to the main GDF5 antagonist NOGGIN (NOG) leads to a GDF5 gain of function and subsequent SYNS2 phenotype. Whereas on the other hand, a reduced signaling activity, specifically via the BMP receptor type IA (BMPR1A), is likely responsible for the BDA1 phenotype. These results demonstrate that one mutation in the overlapping interface of antagonist and receptor binding site in GDF5 can lead to a GDF5 variant with pathophysiological relevance for both, BDA1 and SYNS2 development. Consequently, our study assembles another part of the molecular puzzle of how loss and gain of function mutations in GDF5 affect bone development in hands and feet resulting in specific types of brachydactyly and SYNS2. These novel insights into the biology of GDF5 might also provide further clues on the pathophysiology of OA.

Clinical interpretation of CNVs with cross-species phenotype data.

BACKGROUND: Clinical evaluation of CNVs identified via techniques such as array comparative genome hybridisation (aCGH) involves the inspection of lists of known and unknown duplications and deletions with the goal of distinguishing pathogenic from benign CNVs. A key step in this process is the comparison of the individual's phenotypic abnormalities with those associated with Mendelian disorders of the genes affected by the CNV. However, because often there is not much known about these human genes, an additional source of data that could be used is model organism phenotype data. Currently, almost 6000 genes in mouse and zebrafish are, when knocked out, associated with a phenotype in the model organism, but no disease is known to be caused by mutations in the human ortholog. Yet, searching model organism databases and comparing model organism phenotypes with patient phenotypes for identifying novel disease genes and medical evaluation of CNVs is hindered by the difficulty in integrating phenotype information across species and the lack of appropriate software tools. METHODS: Here, we present an integrated ranking scheme based on phenotypic matching, degree of overlap with known benign or pathogenic CNVs and the haploinsufficiency score for the prioritisation of CNVs responsible for a patient's clinical findings. RESULTS: We show that this scheme leads to significant improvements compared with rankings that do not exploit phenotypic information. We provide a software tool called PhenogramViz, which supports phenotype-driven interpretation of aCGH findings based on multiple data sources, including the integrated cross-species phenotype ontology Uberpheno, in order to visualise gene-to-phenotype relations. CONCLUSIONS: Integrating and visualising cross-species phenotype information on the affected genes may help in routine diagnostics of CNVs.

Whole exome sequencing identifies FGF16 nonsense mutations as the cause of X-linked recessive metacarpal 4/5 fusion.

BACKGROUND: Metacarpal 4-5 fusion (MF4; MIM %309630) is a rare congenital malformation of the hand characterised by the partial or complete fusion of the fourth and fifth metacarpals. The anomaly occurs as an isolated trait or part of a genetic syndrome. METHODS: To search for disease-causing mutation, whole exome sequencing (WES) was performed on samples from a single trio. Before WES, molecular screening including gene sequencing and array comparative genomic hybridisation was applied. Validation of WES and segregation studies were done using routine Sanger sequencing. RESULTS: Exome sequencing detected a nonsense mutation (c.C535T; p.R179X) in exon 3 of the FGF16 gene, which maps to chromosome Xq21.1. Mutational screening of the FGF16 gene performed in an unrelated proband of different ethnicity showed another nonsense mutation in exon 3 (c.C470A; p.S157X). CONCLUSIONS: This study shows that truncating mutations of FGF16 are associated with X-linked recessive metacarpal 4-5 fusion. The study provides evidence for the involvement of FGF16 in the fine tuning of the human skeleton of the hand.

Getting ready for the Human Phenome Project: the 2012 forum of the Human Variome Project.

A forum of the Human Variome Project (HVP) was held as a satellite to the 2012 Annual Meeting of the American Society of Human Genetics in San Francisco, California. The theme of this meeting was "Getting Ready for the Human Phenome Project." Understanding the genetic contribution to both rare single-gene "Mendelian" disorders and more complex common diseases will require integration of research efforts among many fields and better defined phenotypes. The HVP is dedicated to bringing together researchers and research populations throughout the world to provide the resources to investigate the impact of genetic variation on disease. To this end, there needs to be a greater sharing of phenotype and genotype data. For this to occur, many databases that currently exist will need to become interoperable to allow for the combining of cohorts with similar phenotypes to increase statistical power for studies attempting to identify novel disease genes or causative genetic variants. Improved systems and tools that enhance the collection of phenotype data from clinicians are urgently needed. This meeting begins the HVP's effort toward this important goal.

Proximal and distal 15q25.2 microdeletions-genotype-phenotype delineation of two neurodevelopmental susceptibility loci.

Distal 15q25.2 microdeletions have recently been reported as a copy number variation (CNV) locus for neurodevelopmental and neuropsychiatric disorders with variable outcome. In addition, more proximal microdeletions of 15q25.2 have been described as a susceptibility locus for cognitive deficits, congenital diaphragmatic hernia (CDH), and Diamond-Blackfan anaemia (DBA). We describe two patients with 15q25.2 deletion, one with the more distal deletion and the other with a deletion overlapping both the distal and proximal 15q25.2 deletions and compare them to the 18 so far reported patients with 15q25.2 deletions. We provide a characterization of the 15q25.2 microdeletions and contribute to the genotype-phenotype delineation for these two novel microdeletion syndromes.

Duplications of BHLHA9 are associated with ectrodactyly and tibia hemimelia inherited in non-Mendelian fashion.

BACKGROUND: Split-hand/foot malformation (SHFM)-also known as ectrodactyly-is a congenital disorder characterised by severe malformations of the distal limbs affecting the central rays of hands and/or feet. A distinct entity termed SHFLD presents with SHFM and long bone deficiency. Mouse models suggest that a defect of the central apical ectodermal ridge leads to the phenotype. Although six different loci/mutations (SHFM1-6) have been associated with SHFM, the underlying cause in a large number of cases is still unresolved. METHODS: High resolution array comparative genomic hybridisation (CGH) was performed in patients with SHFLD to detect copy number changes. Candidate genes were further evaluated for expression and function during limb development by whole mount in situ hybridisation and morpholino knock-down experiments. RESULTS: Array CGH showed microduplications on chromosome 17p13.3, a locus previously associated with SHFLD. Detailed analysis of 17 families revealed that this copy number variation serves as a susceptibility factor for a highly variable phenotype with reduced penetrance, particularly in females. Compared to other known causes for SHFLD 17p duplications appear to be the most frequent cause of SHFLD. A ~11.8 kb minimal critical region was identified encompassing a single gene, BHLHA9, a putative basic loop helix transcription factor. Whole mount in situ hybridisation showed expression restricted to the limb bud mesenchyme underlying the apical ectodermal ridge in mouse and zebrafish embryos. Knock down of bhlha9 in zebrafish resulted in shortening of the pectoral fins. CONCLUSIONS: Genomic duplications encompassing BHLHA9 are associated with SHFLD and non-Mendelian inheritance characterised by a high degree of non-penetrance with sex bias. Knock-down of bhlha9 in zebrafish causes severe reduction defects of the pectoral fin, indicating a role for this gene in limb development.

Ontological phenotype standards for neurogenetics.

Neurological disorders comprise one of the largest groups of human diseases. Due to the myriad symptoms and the extreme degree of clinical variability characteristic of many neurological diseases, the differential diagnosis process is extremely challenging. Even though most neurogenetic diseases are individually rare, collectively, the subgroup of neurogenetic disorders is large, comprising more than 2,400 different disorders. Recently, increasing efforts have been undertaken to unravel the molecular basis of neurogenetic diseases and to correlate pathogenetic mechanisms with clinical signs and symptoms. In order to enable computer-based analyses, the systematic representation of the neurological phenotype is of major importance. We demonstrate how the Human Phenotype Ontology (HPO) can be incorporated into these efforts by providing a systematic semantic representation of phenotypic abnormalities encountered in human genetic diseases. The combination of the HPO together with the Orphanet disease classification represents a promising resource for automated disease classification, performing computational clustering and analysis of the neurogenetic phenome. Furthermore, standardized representations of neurologic phenotypic abnormalities employing the HPO link neurological phenotypic abnormalities to anatomical and functional entities represented in other biomedical ontologies through the semantic references provided by the HPO.

MouseFinder: Candidate disease genes from mouse phenotype data.

Mouse phenotype data represents a valuable resource for the identification of disease-associated genes, especially where the molecular basis is unknown and there is no clue to the candidate gene's function, pathway involvement or expression pattern. However, until recently these data have not been systematically used due to difficulties in mapping between clinical features observed in humans and mouse phenotype annotations. Here, we describe a semantic approach to solve this problem and demonstrate highly significant recall of known disease-gene associations and orthology relationships. A Web application (MouseFinder; www.mousemodels.org) has been developed to allow users to search the results of our whole-phenome comparison of human and mouse. We demonstrate its use in identifying ARTN as a strong candidate gene within the 1p34.1-p32 mapped locus for a hereditary form of ptosis.

Double heterozygosity for mutations in BRCA1 and BRCA2 in German breast cancer patients: implications on test strategies and clinical management.

Double heterozygosity for disease-causing BRCA1 and BRCA2 mutations is a very rare condition in most populations. Here we describe genetic and clinical data of eight female double heterozygotes (DH) for BRCA1 and BRCA2 mutations found in a cohort of 8162 German breast/ovarian cancer families and compare it with the data of their single heterozygous relatives and of the index patients of the German Consortium for Hereditary Breast and Ovarian Cancer. Furthermore, we analyze the phenotypic features of these patients with respect to age at onset of first cancer, first breast/ovarian cancer and the number of disease manifestations and compare them to that of published Caucasian female DHs and their single heterozygous female relatives. German DHs were not significantly younger at diagnosis of first breast cancer than the single heterozygous index patients of the German Consortium. However, if the data of our study were pooled with that of the literature, DHs were substantially younger at onset of first cancer (mean age 40.4 years, 95 % CI = 36.6-44.1) than their single heterozygous female relatives (mean age 51.9 years, 95 % CI = 46.8-57.0). The two groups also differed concerning the onset of first breast cancer (mean age 40.6 years, 95 % CI = 36.6-44.5 vs. 52.6, 95 % CI = 47.5-57.6). In addition, DHs had a more severe disease than their female relatives carrying a single BRCA mutation (1.4 vs. 0.6 manifestations per person). In contrast to Ashkenazi Jewish females, Caucasian DH females might develop breast cancer at an earlier age and have a more severe disease than single heterozygous BRCA mutation carriers. Therefore, DHs may benefit from more intensive surveillance programs/follow-up care and prophylactic surgery.

Deletions of exons with regulatory activity at the DYNC1I1 locus are associated with split-hand/split-foot malformation: array CGH screening of 134 unrelated families.

BACKGROUND: A growing number of non-coding regulatory mutations are being identified in congenital disease. Very recently also some exons of protein coding genes have been identified to act as tissue specific enhancer elements and were therefore termed exonic enhancers or "eExons". METHODS: We screened a cohort of 134 unrelated families with split-hand/split-foot malformation (SHFM) with high resolution array CGH for CNVs with regulatory potential. RESULTS: In three families with an autosomal dominant non-syndromic SHFM phenotype we detected microdeletions encompassing the exonic enhancer (eExons) 15 and 17 of DYNC1I1. In a fourth family, who had hearing loss in addition to SHFM, we found a larger deletion of 510 kb including the eExons of DYNC1I1 and, in addition, the human brain enhancer hs1642. Exons 15 and 17 of DYNC1I1 are known to act as tissue specific limb enhancers of DLX5/6, two genes that have been shown to be associated with SHFM in mice. In our cohort of 134 unrelated families with SHFM, deletions of the eExons of DYNC1I1 account for approximately 3% of the cases, while 17p13.3 duplications were identified in 13% of the families, 10q24 duplications in 12%, and TP63 mutations were detected in 4%. CONCLUSIONS: We reduce the minimal critical region for SHFM1 to 78 kb. Hearing loss, however, appears to be associated with deletions of a more telomeric region encompassing the brain enhancer element hs1642. Thus, SHFM1 as well as hearing loss at the same locus are caused by deletion of regulatory elements. Deletions of the exons with regulatory potential of DYNC1I1 are an example of the emerging role of exonic enhancer elements and their implications in congenital malformation syndromes.

Homozygous missense and nonsense mutations in BMPR1B cause acromesomelic chondrodysplasia-type Grebe.

Acromesomelic chondrodysplasias (ACDs) are characterized by disproportionate shortening of the appendicular skeleton, predominantly affecting the middle (forearms and forelegs) and distal segments (hands and feet). Here, we present two consanguineous families with missense (c.157T>C, p.(C53R)) or nonsense (c.657G>A, p.(W219*)) mutations in BMPR1B. Homozygous affected individuals show clinical and radiographic findings consistent with ACD-type Grebe. Functional analysis of the missense mutation C53R revealed that the mutated receptor was partially located at the cell membrane. In contrast to the wild-type receptor, C53R mutation hindered the activation of the receptor by its ligand GDF5, as shown by reporter gene assay. Further, overexpression of the C53R mutation in an in vitro chondrogenesis assay showed no effect on cell differentiation, indicating a loss of function. The nonsense mutation (c.657G>A, p.(W219*)) introduces a premature stop codon, which is predicted to be subject to nonsense-mediated mRNA decay, causing reduced protein translation of the mutant allele. A loss-of-function effect of both mutations causing recessive ACD-type Grebe is further supported by the mild brachydactyly or even non-penetrance of these mutations observed in the heterozygous parents. In contrast, dominant-negative BMPR1B mutations described previously are associated with autosomal-dominant brachydactyly-type A2.

Construction and accessibility of a cross-species phenotype ontology along with gene annotations for biomedical research.

Phenotype analyses, e.g. investigating metabolic processes, tissue formation, or organism behavior, are an important element of most biological and medical research activities. Biomedical researchers are making increased use of ontological standards and methods to capture the results of such analyses, with one focus being the comparison and analysis of phenotype information between species. We have generated a cross-species phenotype ontology for human, mouse and zebrafish that contains classes from the Human Phenotype Ontology, Mammalian Phenotype Ontology, and generated classes for zebrafish phenotypes. We also provide up-to-date annotation data connecting human genes to phenotype classes from the generated ontology. We have included the data generation pipeline into our continuous integration system ensuring stable and up-to-date releases. This article describes the data generation process and is intended to help interested researchers access both the phenotype annotation data and the associated cross-species phenotype ontology. The resource described here can be used in sophisticated semantic similarity and gene set enrichment analyses for phenotype data across species. The stable releases of this resource can be obtained from http://purl.obolibrary.org/obo/hp/uberpheno/.

Phenotypic overlap in the contribution of individual genes to CNV pathogenicity revealed by cross-species computational analysis of single-gene mutations in humans, mice and zebrafish.

Numerous disease syndromes are associated with regions of copy number variation (CNV) in the human genome and, in most cases, the pathogenicity of the CNV is thought to be related to altered dosage of the genes contained within the affected segment. However, establishing the contribution of individual genes to the overall pathogenicity of CNV syndromes is difficult and often relies on the identification of potential candidates through manual searches of the literature and online resources. We describe here the development of a computational framework to comprehensively search phenotypic information from model organisms and single-gene human hereditary disorders, and thus speed the interpretation of the complex phenotypes of CNV disorders. There are currently more than 5000 human genes about which nothing is known phenotypically but for which detailed phenotypic information for the mouse and/or zebrafish orthologs is available. Here, we present an ontology-based approach to identify similarities between human disease manifestations and the mutational phenotypes in characterized model organism genes; this approach can therefore be used even in cases where there is little or no information about the function of the human genes. We applied this algorithm to detect candidate genes for 27 recurrent CNV disorders and identified 802 gene-phenotype associations, approximately half of which involved genes that were previously reported to be associated with individual phenotypic features and half of which were novel candidates. A total of 431 associations were made solely on the basis of model organism phenotype data. Additionally, we observed a striking, statistically significant tendency for individual disease phenotypes to be associated with multiple genes located within a single CNV region, a phenomenon that we denote as pheno-clustering. Many of the clusters also display statistically significant similarities in protein function or vicinity within the protein-protein interaction network. Our results provide a basis for understanding previously un-interpretable genotype-phenotype correlations in pathogenic CNVs and for mobilizing the large amount of model organism phenotype data to provide insights into human genetic disorders.

Functional analysis of alleged NOGGIN mutation G92E disproves its pathogenic relevance.

We identified an amino acid change (p.G92E) in the Bone Morphogenetic Protein antagonist NOGGIN in a 22-month-old boy who presented with a unilateral brachydactyly type B phenotype. Brachydactyly type B is a skeletal malformation that has been associated with increased Bone Morphogenetic Protein pathway activation in other patients. Previously, the amino acid change p.G92E in NOGGIN was described as causing fibrodysplasia ossificans progressiva, a rare genetic disorder characterized by limb malformations and progressive heterotopic bone formation in soft tissues that, like Brachydactyly type B, is caused by increased activation of Bone Morphogenetic Protein signaling. To determine whether G92E-NOGGIN shows impaired antagonism that could lead to increased Bone Morphogenetic Protein signaling, we performed functional assays to evaluate inhibition of BMP signaling. Interestingly, wt-NOGGIN shows different inhibition efficacies towards various Bone Morphogenetic Proteins that are known to be essential in limb development. However, comparing the biological activity of G92E-NOGGIN with wt-NOGGIN, we observed that G92E-NOGGIN inhibits activation of bone morphogenetic protein signaling with equal efficiency as wt-NOGGIN, supporting that G92E-NOGGIN does not cause pathological effects. Genetic testing of the child's parents revealed the same amino acid change in the healthy father, further supporting that p.G92E is a neutral amino acid substitution in NOGGIN. We conclude that p.G92E represents a rare polymorphism of the NOGGIN gene-- causing neither brachydactyly nor fibrodysplasia ossificans progressiva. This study highlights that a given genetic variation should not be considered pathogenic unless supported by functional analyses.