Search

Objective Next generation sequencing (NGS) has expanded the diagnostic paradigm turning the focus to the growth plate. The aim of the study was to determine the prevalence of variants in genes implicated in skeletal dysplasias in probands with short stature and mild skeletal anomalies. Design Clinical and radiological data were collected from 108 probands with short stature and mild skeletal anomalies. Methods A customized skeletal dysplasia NGS panel was performed. Variants were classified using ACMG recommendations and Sherloc. Anthropometric measurements and skeletal anomalies were subsequently compared in those with or without an identified genetic defect. Results Heterozygous variants were identified in 21/108 probands (19.4%). Variants were most frequently identified in ACAN (n = 10) and IHH (n = 7) whilst one variant was detected in COL2A1, CREBBP, EXT1, and PTPN11. Statistically significant differences (P < 0.05) were observed for sitting height/height (SH/H) ratio, SH/H ratio standard deviation score (SDS), and the SH/H ratio SDS >1 in those with an identified variant compared to those without. Conclusions A molecular defect was elucidated in a fifth of patients. Thus, the prevalence of mild forms of skeletal dysplasias is relatively high in individuals with short stature and mild skeletal anomalies, with variants in ACAN and IHH accounting for 81% of the cases. An elevated SH/H ratio appears to be associated with a greater probability in detecting a variant, but no other clinical or radiological feature has been found determinant to finding a genetic cause. Currently, we cannot perform extensive molecular studies in all short stature individuals so detailed clinical and radiological phenotyping may orientate which are the candidate patients to obtain worthwhile results. In addition, detailed phenotyping of probands and family members will often aid variant classification.

PURPOSE: Germline loss-of-function variants in CTNNB1 cause neurodevelopmental disorder with spastic diplegia and visual defects (NEDSDV; OMIM 615075) and are the most frequent, recurrent monogenic cause of cerebral palsy (CP). We investigated the range of clinical phenotypes owing to disruptions of CTNNB1 to determine the association between NEDSDV and CP.METHODS: Genetic information from 404 individuals with collectively 392 pathogenic CTNNB1 variants were ascertained for the study. From these, detailed phenotypes for 52 previously unpublished individuals were collected and combined with 68 previously published individuals with comparable clinical information. The functional effects of selected CTNNB1 missense variants were assessed using TOPFlash assay.RESULTS: The phenotypes associated with pathogenic CTNNB1 variants were similar. A diagnosis of CP was not significantly associated with any set of traits that defined a specific phenotypic subgroup, indicating that CP is not additional to NEDSDV. Two CTNNB1 missense variants were dominant negative regulators of WNT signaling, highlighting the utility of the TOPFlash assay to functionally assess variants.CONCLUSION: NEDSDV is a clinically homogeneous disorder irrespective of initial clinical diagnoses, including CP, or entry points for genetic testing.

Background Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance. Results In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5–100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels. Conclusion These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

Our purpose was to identify mutations responsible for non-syndromic congenital cataracts through the implementation of next-generation sequencing (NGS) in our center. A sample of peripheral blood was obtained from probands and willing family members and genomic DNA was extracted from leukocytes. DNA was analyzed implementing a panel (OFTv2.1) including 39 known congenital cataracts disease genes. 62 probands from 51 families were recruited. Pathogenic or likely pathogenic variants were identified in 32 patients and 25 families; in 16 families (64%) these were de novo mutations. The mutation detection rate was 49%. Almost all reported mutations were autosomal dominant. Mutations in crystallin genes were found in 30% of the probands. Mutations in membrane proteins were detected in seven families (two in GJA3 and five in GJA8). Mutations in LIM2 and MIP were each found in three families. Other mutations detected affected EPHA2, PAX6, HSF4 and PITX3. Variants classified as of unknown significance were found in 5 families (9.8%), affecting CRYBB3, LIM2, EPHA2, ABCB6 and TDRD7. Mutations lead to different cataract phenotypes within the same family, This research was funded by ONCE grant number 2020/0197782 and by FIS grant number PI18-1234-ISCIII

Early-onset high myopia (EoHM) is a disease that causes a spherical refraction error of ≥−6 diopters before 10 years of age, with potential multiple ocular complications. In this article, we report a clinical and genetic study of 43 families with EoHM recruited in our center. A complete ophthalmological evaluation was performed, and a sample of peripheral blood was obtained from proband and family members. DNA was analyzed using a customized next-generation sequencing panel that included 419 genes related to ophthalmological disorders with a suspected genetic cause, and genes related to EoHM pathogenesis. We detected pathogenic and likely pathogenic variants in 23.9% of the families and detected variants of unknown significance in 76.1%. Of these, 5.7% were found in genes related to non-syndromic EoHM, 48.6% in genes associated with inherited retinal dystrophies that can include a syndromic phenotype, and 45.7% in genes that are not directly related to EoHM or retinal dystrophy. We found no candidate genes in 23% of the patients, which suggests that further studies are needed. We propose a systematic genetic analysis for patients with EoHM because it helps with follow-up, prognosis and genetic counseling.

Context Heterozygous variants in the Indian hedgehog gene (IHH) have been reported to cause brachydactyly type A1 and mild hand and feet skeletal anomalies with short stature. Genetic screening in individuals with short stature and mild skeletal anomalies has been increasing over recent years, allowing us to broaden the clinical spectrum of skeletal dysplasias. Objective The objective of this article is to describe the genotype and phenotype of 16 probands with heterozygous variants in IHH. Patients and Methods Targeted next-generation sequencing or Sanger sequencing was performed in patients with short stature and/or brachydactyly for which the genetic cause was unknown. Results Fifteen different heterozygous IHH variants were detected, one of which is the first reported complete deletion of IHH. None of the patients showed the classical phenotype of brachydactyly type A1. The most frequently observed clinical characteristics were mild to moderate short stature as well as shortening of the middle phalanx on the fifth finger. The identified IHH variants were demonstrated to cosegregate with the short stature and/or brachydactyly in the 13 probands whose family members were available. However, clinical heterogeneity was observed: Two short-statured probands showed no hand radiological anomalies, whereas another 5 were of normal height but had brachydactyly. Conclusions Short stature and/or mild skeletal hand defects can be caused by IHH variants. Defects in this gene should be considered in individuals with these findings, especially when there is an autosomal dominant pattern of inheritance. Although no genotype-phenotype correlation was observed, cosegregation studies should be performed and where possible functional characterization before concluding that a variant is causative.

Generalized lymphatic anomaly (GLA) is a vascular disorder characterized by diffuse or multifocal lymphatic malformations (LMs). The etiology of GLA is poorly understood. We identified four distinct somatic PIK3CA variants (Glu542Lys, Gln546Lys, His1047Arg, and His1047Leu) in tissue samples from five out of nine patients with GLA. These same PIK3CA variants occur in PIK3CA-related overgrowth spectrum and cause hyperactivation of the PI3K–AKT–mTOR pathway. We found that the mTOR inhibitor, rapamycin, prevented lymphatic hyperplasia and dysfunction in mice that expressed an active form of PIK3CA (His1047Arg) in their lymphatics. We also found that rapamycin reduced pain in patients with GLA. In conclusion, we report that somatic activating PIK3CA mutations can cause GLA, and we provide preclinical and clinical evidence to support the use of rapamycin for the treatment of this disabling and deadly disease.

Transducin β-like 1 X-linked (TBL1X) gene encodes a subunit of the nuclear corepressor-silencing mediator for retinoid and thyroid hormone receptor complex (NCoR-SMRT) involved in repression of thyroid hormone action in the pituitary and hypothalamus. TBL1X defects were recently associated with central hypothyroidism and hearing loss. The current study aims to describe the clinical and genetic characterization of a male diagnosed with central hypothyroidism through thyroid hormone profiling, TRH test, brain MRI, audiometry, and psychological evaluation. Next-generation sequencing of known genes involved in thyroid disorders was implemented. The 6-year-old boy was diagnosed with central hypothyroidism [free T4: 10.42 pmol/L (normal: 12 to 22 pmol/L); TSH: 1.57 mIU/L (normal: 0.7 to 5.7 mIU/L)], with a mildly reduced TSH response to TRH. He was further diagnosed with attention-deficit/hyperactivity disorder (ADHD) at 7 years, alternating episodes of encopresis and constipation, and frequent headaches. MRI showed a normal pituitary but detected a Chiari malformation type I (CMI). At 10 years, audiometry identified poor hearing threshold at high frequencies. Sequencing revealed a nonsense hemizygous mutation in TBL1X [c.1015C>T; p.(Arg339Ter)] largely truncating its WD-40 repeat domain involved in nuclear protein-protein interactions. In conclusion, to our knowledge, we identified the first severely truncating TBL1X mutation in a patient with central hypothyroidism, hypoacusia, and novel clinical features like ADHD, gastrointestinal dysmotility, and CMI. Given the relevance of TBL1X and NCoR-SMRT for the regulation of transcriptional programs at different tissues (pituitary, cochlea, brain, fossa posterior, and cerebellum), severe mutations in TBL1X may lead to a distinct syndrome with a phenotypic spectrum wider than previously reported.

Mosaic Variegated Aneuploidy Syndrome 2 (MVA2; MIM 614114 ) is a rare autosomal recessive disorder, characterized by mosaic aneuploidies involving multiple chromosomes and tissues, caused by biallelic pathogenic variants in the CEP57 gene. Only 10 patients have been reported to date. We report two additional non related cases born to Moroccan consanguineous parents, carrying the previously described c.915_925dup11 CEP57 homozygous variant. Common features of these 12 cases include growth retardation, typically of prenatal onset, distinctive facial features, endocrine, cardiovascular and skeletal, abnormalities while malignancies have not been reported. This report describes the phenotypical spectrum of MVA2.

Craniosynostosis is commonly caused by mutations in fibroblast growth factor receptors (FGFRs), highlighting the essential role of FGF-mediated signaling in skeletal development. We set out to identify the molecular defect in a family referred for craniosynostosis and in whom no mutation was previously detected. Using next-generation sequencing, we identified a novel missense mutation in FGF9. Modeling based upon the crystal structure and functional studies confirmed its pathogenicity showing that it impaired homodimerization and FGFR3 binding. Only one FGF9 mutation has been previously reported in a multigeneration family with multiple synostoses (SYNS3) but no signs of craniosynostosis. In contrast, our family has a greater phenotypic resemblance to that observed in the Fgf9 spontaneous mouse mutant, elbow-knee-synostosis, Eks, with both multiple synostoses and craniosynostosis. We have demonstrated for the first time that mutations in FGF9 cause craniosynostosis in humans and confirm that FGF9 mutations cause multiple synostoses.

Glucogenosis type IX is caused by pathogenic variants of the PHKA2 gene. Herein, we report a patient with clinical symptoms compatible with Glycogen Storage Disease type IXa. PYGL, PHKA1, PHKA2, PHKB and PHKG2 genes were analyzed by Next Generation Sequencing (NGS). We identified the previously undescribed hemizygous missense variant NM_000292.2(PHKA2):c.1963G > A, p.(Glu655Lys) in PHKA2 exon 18. In silico analyses showed two possible pathogenic consequences: it affects a highly conserved amino acid and disrupts the exon 18 canonical splice donor site. The variant was found as a “de novo” event.

[Purpose] Germline mutations in DNA damage repair (DDR) genes are identified in a significant proportion of patients with metastatic prostate cancer, but the clinical implications of these genes remain unclear. This prospective multicenter cohort study evaluated the prevalence and effect of germline DDR (gDDR) mutations on metastatic castration-resistance prostate cancer (mCRPC) outcomes., [Patients and Methods] Unselected patients were enrolled at diagnosis of mCRPC and were screened for gDDR mutations in 107 genes. The primary aim was to assess the impact of ATM/BRCA1/BRCA2/PALB2 germline mutations on cause-specific survival (CSS) from diagnosis of mCRPC. Secondary aims included the association of gDDR subgroups with response outcomes for mCRPC treatments. Combined progression-free survival from the first systemic therapy (PFS) until progression on the second systemic therapy (PFS2) was also explored., [Results] We identified 68 carriers (16.2%) of 419 eligible patients, including 14 with BRCA2, eight with ATM, four with BRCA1, and none with PALB2 mutations. The study did not reach its primary end point, because the difference in CSS between ATM/BRCA1/BRCA2/PALB2 carriers and noncarriers was not statistically significant (23.3 v 33.2 months; P = .264). CSS was halved in germline BRCA2 (gBRCA2) carriers (17.4 v 33.2 months; P = .027), and gBRCA2 mutations were identified as an independent prognostic factor for CCS (hazard ratio [HR], 2.11; P = .033). Significant interactions between gBRCA2 status and treatment type (androgen signaling inhibitor v taxane therapy) were observed (CSS adjusted P = .014; PFS2 adjusted P = .005). CSS (24.0 v 17.0 months) and PFS2 (18.9 v 8.6 months) were greater in gBRCA2 carriers treated in first line with abiraterone or enzalutamide compared with taxanes. Clinical outcomes did not differ by treatment type in noncarriers., [Conclusion] gBRCA2 mutations have a deleterious impact on mCRPC outcomes that may be affected by the first line of treatment used. Determination of gBRCA2 status may be of assistance for the selection of the initial treatment in mCRPC. Nonetheless, confirmatory studies are required before these results can support a change in clinical practice., Supported by an unrestricted grant from Fundación Cris contra el cancer; three investigator awards from the Prostate Cancer Foundation (C.C.P. [2013], D.O. [2014], and E.C. [2017]); and three grants from Fondo de Investigación Sanitaria, Instituto de Salud Carlos III (No. PI13/01287 and PI16/01565 to D.O. and No. PI15/01471 to P.L.). During the conduct of this study, E.C., D.O., P.N., and L.M.-P. were supported by grants from Ministerio de Economía, Industria y Competitividad (No. JCI-2014-19129 [E.C.], No. RYC-2015-18625 [D.O.], No. SVP-2013-067937 [P.N.], No. SVP-2014-068895 [L.M.]); D.O. was also funded by a Return fellowship from Fundación Científica de la Asociación Española Contra el Cancer, 2012-2015; N.R.L. and R.L., by grants from Instituto de Salud Carlos III (No. CM14-00200 to N.R.L. and No. CM17-00221 [R.L.]); and Y.C., by a grant from Ministerio de Educación, Cultura y Deportes (No. FPU15/05126). C.C.P. was supported by a congressional-designated medical research program award (No. CMRP-PC131820).

Background: Congenital aniridia is a complex ocular disorder, usually associated with severe visual impairment, generally caused by mutations on the PAX6 gene. The clinical phenotype of PAX6 mutations is highly variable, making the genotype–phenotype correlations difficult to establish. Methods: we describe the phenotype of eight patients from seven unrelated families with confirmed mutations in PAX6, and very different clinical manifestations. Results: Only two patients had the classical aniridia phenotype while the other two presented with aniridia-related manifestations, such as aniridia-related keratopathy or partial aniridia. Congenital cataracts were the main manifestation in three of the patients in this series. All the patients had nystagmus and low visual acuity. Conclusions: The diagnosis of mild forms of aniridia is challenging, but these patients have a potentially blinding hereditary disease that might present with a more severe phenotype in future generations. Clinicians should be aware of the mild aniridia phenotype and request genetic testing to perform an accurate diagnosis.

X-linked myotubular myopathy (XLMTM; OMIM 310400) is a centronuclear congenital muscular disorder of X-linked recessive inheritance. Although female carriers are typically asymptomatic, affected heterozygous females have been described. Here, we describe the case of a sporadic female patient with suspicion of centronuclear myopathy and a heterozygous large deletion at Xq28 encompassing the MAMLD1, MTM1, MTMR1, CD99L2, and HMGB3 genes. The deletion was first detected using a custom next generation sequencing (NGS)-based multigene panel and finally characterized by comparative genomic hybridization array and multiplex ligation probe assay techniques. In this patient we have confirmed, by MTM1 mRNA quantification, a MTM1 gene expression less than the expected 50 percent in patient muscle. The significant 20% reduction in MTM1 mRNA expression in muscle, precludes low level of the normal myotubularin protein as the cause of the phenotype in this heterozygous female. We have also found that BIN1 expression in patient muscle biopsy was significantly increased, and postulate that BIN1 expression will be increased in XLMTM patient muscle as an attempt to maintain muscle function.

Progressive pseudorheumatoid dysplasia (PPD) is a rare autosomal recessive disorder characterized by spondyloepiphyseal dysplasia associated with pain and stiffness of multiple joints, enlargement of the interphalangeal joints, normal inflammatory parameters, and absence of extra-skeletal manifestations. Homozygous or compound heterozygous WISP3 mutations cause PPD. We report two siblings from a non-consanguineous Ecuadorian family with a late-onset spondyloepiphyseal dysplasia. Mutation screening was undertaken in the two affected siblings using a customized skeletal dysplasia next generation sequencing (NGS) panel and confirmed by Sanger sequencing. Two compound heterozygous mutations were identified in WISP3 exon 2, c.[190G>A];[197G>A] (p.[(Gly64Arg)];[(Ser66Asn)]) in the two siblings, both of which had been inherited. The p. (Gly64Arg) mutation has not been previously described whilst the p. (Ser66Asn) mutation has been reported in two PPD families. The two siblings presented with atypical PPD, as they presented during late childhood, yet the severity was different between them. The progression was particularly aggressive in the male sibling who suffered severe scoliosis by the age of 13 years. This case reaffirms the clinical heterogeneity of this disorder and the clinical utility of NGS to genetically diagnose skeletal dysplasias, enabling adequate management, monitorization, and genetic counseling. © 2016 Wiley Periodicals, Inc.

Pitfalls of trio-based exome sequencing: imprinted genes and parental mosaicism— MAGEL2 as an example

Despite BMP4 signaling being critical to Rathke's pouch induction and maintenance during early stages of pituitary development, its implication in the etiology of combined pituitary hormone deficiency (CPHD) and other clinical presentations of congenital hypopituitarism has not yet been definitely demonstrated. We report here the first CPHD patient with a de novo pathogenic loss-of-function variant in BMP4. A 6-year-old boy, with macrocephaly, myopia/astigmatism, mild psychomotor retardation, anterior pituitary hypoplasia and ectopic posterior pituitary, clinically diagnosed with growth hormone deficiency, and central hypothyroidism, was referred for genetic analysis of CPHD. Targeted NGS analysis with a custom panel (n = 310 genes) identified a novel heterozygous de novo nonsense variant, NM_001202.5:c.794G > A, p.(Trp265*) in BMP4, which introduces a premature stop codon in the BMP4 pro-domain, impairing the transcription of the TGF-β mature peptide domain. Additional relevant variants in other genes implicated in pituitary development signaling pathways such as SMAD4 and E2F4 (BMP/TGF-pathway), ALMS1 (NOTCH-pathway), and TSHZ1 (Prokineticin-pathway), were also identified. Our results support the implication of the BMP/TGF-β signaling pathway in the etiology of CPHD and suggest that oligogenic contribution of additional inherited variants may modify the phenotypic expressivity of BMP4 pathogenic variants.

The authors have carried out an investigation of the two draft maps of the human proteome published in 2014 in Nature. The findings include an abundance of poor spectra, low-scoring peptide-spectrum matches and incorrectly identified proteins in both these studies, highlighting clear issues with the application of false discovery rates. This noise means that the claims made by the two papers - the identification of high numbers of protein coding genes, the detection of novel coding regions and the draft tissue maps themselves - should be treated with considerable caution. The authors recommend that clinicians and researchers do not use the unfiltered data from these studies. Despite this these studies will inspire further investigation into tissue-based proteomics. As long as this future work has proper quality controls, it could help produce a consensus map of the human proteome and improve our understanding of the processes that underlie health and disease. This paper is supported by a National Institutes of Health grant [U41 HG007234] and by the Spanish Ministry of Economics and Competitiveness [BIO2012-40205, BIO2012-37926 PRB2-ProteoRed-PT13/0001/0017 RD07-0067-0014-COMBIOMED, RETICS-RD12-0042-0056]. The authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed. Sí

[Purpose]: CLAPO syndrome is a rare vascular disorder characterized by capillary malformation of the lower lip, lymphatic malformation predominant on the face and neck, asymmetry, and partial/generalized overgrowth. Here we tested the hypothesis that, although the genetic cause is not known, the tissue distribution of the clinical manifestations in CLAPO seems to follow a pattern of somatic mosaicism., [Methods]: We clinically evaluated a cohort of 13 patients with CLAPO and screened 20 DNA blood/tissue samples from 9 patients using high-throughput, deep sequencing., [Results]: We identified five activating mutations in the PIK3CA gene in affected tissues from 6 of the 9 patients studied; one of the variants (NM_006218.2:c.248T>C; p.Phe83Ser) has not been previously described in developmental disorders., [Conclusion]: We describe for the first time the presence of somatic activating PIK3CA mutations in patients with CLAPO. We also report an update of the phenotype and natural history of the syndrome., This research was supported by the project “Genetics of vascular and lymphatic malformations” financed with funds donated by Asociación Ultrafondo and Villareal FC, cofinanced by project IP-17 from the funding call “Todos Somos Raros” (Telemaraton TVE promoted by Fundación Isabel Gemio, Federación ASEM, and Federación Española de Enfermedades Raras), cofinanced by the Instituto de Salud Carlos III, FEDER FUNDS FIS PI15/01481, and IIS-Fundación Jiménez Díaz UAM Genome Medicine Chair.

[Background]: Osteogenesis imperfecta (OI) is a heterogeneous bone disorder characterized by recurrent fractures. Although most cases of OI have heterozygous mutations in COL1A1 or COL1A2 and show autosomal dominant inheritance, during the last years there has been an explosion in the number of genes responsible for both recessive and dominant forms of this condition. Herein, we have analyzed a cohort of patients with OI, all offspring of unaffected parents, to determine the spectrum of variants accounting for these cases. Twenty patients had nonrelated parents and were sporadic, and 21 were born to consanguineous relationships., [Methods]: Mutation analysis was performed using a next‐generation sequencing gene panel, homozygosity mapping, and whole exome sequencing (WES)., [Results]: Patients offspring of nonconsanguineous parents were mostly identified with COL1A1 or COL1A2 heterozygous changes, although there were also a few cases with IFITM5 and WNT1 heterozygous mutations. Only one sporadic patient was a compound heterozygote for two recessive mutations. Patients offspring of consanguineous parents showed homozygous changes in a variety of genes including CRTAP,FKBP10,LEPRE1,PLOD2,PPIB,SERPINF1,TMEM38B, and WNT1. In addition, two patients born to consanguineous parents were found to have de novo COL1A1 heterozygous mutations demonstrating that causative variants in the collagen I structural genes cannot be overlooked in affected children from consanguineous couples. Further to this, WES analysis in probands lacking mutations in OI genes revealed deleterious variants in SCN9A,NTRK1, and SLC2A2, which are associated with congenital indifference to pain (CIP) and Fanconi–Bickel syndrome (FBS)., [Conclusion]: This work provides useful information for clinical and genetic diagnosis of OI patients with no positive family history of this disease. Our data also indicate that CIP and FBS are conditions to be considered in the differential diagnosis of OI and suggest a positive role of SCN9A and NTRK1 in bone development., Spanish Ministry of Economy and Competitiveness (Grant/Award Number: ‘SAF2013-43365-R’), CIBERER (Grant/Award Number: ‘ER14-PR04-ACCI13-760’).

BACKGROUND: A major bottleneck in our understanding of the molecular underpinnings of life is the assignment of function to proteins. While molecular experiments provide the most reliable annotation of proteins, their relatively low throughput and restricted purview have led to an increasing role for computational function prediction. However, assessing methods for protein function prediction and tracking progress in the field remain challenging. RESULTS: We conducted the second critical assessment of functional annotation (CAFA), a timed challenge to assess computational methods that automatically assign protein function. We evaluated 126 methods from 56 research groups for their ability to predict biological functions using Gene Ontology and gene-disease associations using Human Phenotype Ontology on a set of 3681 proteins from 18 species. CAFA2 featured expanded analysis compared with CAFA1, with regards to data set size, variety, and assessment metrics. To review progress in the field, the analysis compared the best methods from CAFA1 to those of CAFA2. CONCLUSIONS: The top-performing methods in CAFA2 outperformed those from CAFA1. This increased accuracy can be attributed to a combination of the growing number of experimental annotations and improved methods for function prediction. The assessment also revealed that the definition of top-performing algorithms is ontology specific, that different performance metrics can be used to probe the nature of accurate predictions, and the relative diversity of predictions in the biological process and human phenotype ontologies. While there was methodological improvement between CAFA1 and CAFA2, the interpretation of results and usefulness of individual methods remain context-dependent., We acknowledge the contributions of Maximilian Hecht, Alexander Grün, Julia Krumhoff, My Nguyen Ly, Jonathan Boidol, Rene Schoeffel, Yann Spöri, Jessika Binder, Christoph Hamm and Karolina Worf. This work was partially supported by the following grants: National Science Foundation grants DBI-1458477 (PR), DBI-1458443 (SDM), DBI-1458390 (CSG), DBI-1458359 (IF), IIS-1319551 (DK), DBI-1262189 (DK), and DBI-1149224 (JC); National Institutes of Health grants R01GM093123 (JC), R01GM097528 (DK), R01GM076990 (PP), R01GM071749 (SEB), R01LM009722 (SDM), and UL1TR000423 (SDM); the National Natural Science Foundation of China grants 3147124 (WT) and 91231116 (WT); the National Basic Research Program of China grant 2012CB316505 (WT); NSERC grant RGPIN 371348-11 (PP); FP7 infrastructure project TransPLANT Award 283496 (ADJvD); Microsoft Research/FAPESP grant 2009/53161-6 and FAPESP fellowship 2010/50491-1 (DCAeS); Biotechnology and Biological Sciences Research Council grants BB/L020505/1 (DTJ), BB/F020481/1 (MJES), BB/K004131/1 (AP), BB/F00964X/1 (AP), and BB/L018241/1 (CD); the Spanish Ministry of Economics and Competitiveness grant BIO2012-40205 (MT); KU Leuven CoE PFV/10/016 SymBioSys (YM); the Newton International Fellowship Scheme of the Royal Society grant NF080750 (TN). CSG was supported in part by the Gordon and Betty Moore Foundation’s Data-Driven Discovery Initiative grant GBMF4552. Computational resources were provided by CSC – IT Center for Science Ltd., Espoo, Finland (TS). This work was supported by the Academy of Finland (TS). RCL and ANM were supported by British Heart Foundation grant RG/13/5/30112. PD, RCL, and REF were supported by Parkinson’s UK grant G-1307, the Alexander von Humboldt Foundation through the German Federal Ministry for Education and Research, Ernst Ludwig Ehrlich Studienwerk, and the Ministry of Education, Science and Technological Development of the Republic of Serbia grant 173001. This work was a Technology Development effort for ENIGMA – Ecosystems and Networks Integrated with Genes and Molecular Assemblies (http://enigma.lbl.gov), a Scientific Focus Area Program at Lawrence Berkeley National Laboratory, which is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Biological & Environmental Research grant DE-AC02-05CH11231. ENIGMA only covers the application of this work to microbial proteins. NSF DBI-0965616 and Australian Research Council grant DP150101550 (KMV). NSF DBI-0965768 (ABH). NIH T15 LM00945102 (training grant for CSF). FP7 FET grant MAESTRA ICT-2013-612944 and FP7 REGPOT grant InnoMol (FS). NIH R01 GM60595 (PCB). University of Padova grants CPDA138081/13 (ST) and GRIC13AAI9 (EL). Swiss National Science Foundation grant 150654 and UK BBSRC grant BB/M015009/1 (COD). PRB2 IPT13/0001 - ISCIII-SGEFI / FEDER (JMF)., This is the final version of the article. It first appeared from BioMed Central at http://dx.doi.org/10.1186/s13059-016-1037-6.

Advances in high-throughput mass spectrometry are making proteomics an increasingly important tool in genome annotation projects. Peptides detected in mass spectrometry experiments can be used to validate gene models and verify the translation of putative coding sequences (CDSs). Here, we have identified peptides that cover 35% of the genes annotated by the GENCODE consortium for the human genome as part of a comprehensive analysis of experimental spectra from two large publicly available mass spectrometry databases. We detected the translation to protein of “novel” and “putative” protein-coding transcripts as well as transcripts annotated as pseudogenes and nonsense-mediated decay targets. We provide a detailed overview of the population of alternatively spliced protein isoforms that are detectable by peptide identification methods. We found that 150 genes expressed multiple alternative protein isoforms. This constitutes the largest set of reliably confirmed alternatively spliced proteins yet discovered. Three groups of genes were highly overrepresented. We detected alternative isoforms for 10 of the 25 possible heterogeneous nuclear ribonucleoproteins, proteins with a key role in the splicing process. Alternative isoforms generated from interchangeable homologous exons and from short indels were also significantly enriched, both in human experiments and in parallel analyses of mouse and Drosophila proteomics experiments. Our results show that a surprisingly high proportion (almost 25%) of the detected alternative isoforms are only subtly different from their constitutive counterparts. Many of the alternative splicing events that give rise to these alternative isoforms are conserved in mouse. It was striking that very few of these conserved splicing events broke Pfam functional domains or would damage globular protein structures. This evidence of a strong bias toward subtle differences in CDS and likely conserved cellular function and structure is remarkable and strongly suggests that the translation of alternative transcripts may be subject to selective constraints.

5071Background: Germline mutations in BRCA2 have been identified in 3-5% of mCRPC patients. PROREPAIR-B (Castro et al ESMO 2017) is the first prospective study to report a worse survival from mCRPC...

164 Background: Germline mutations in DNA repair genes have been associated with poor prostate cancer outcomes in retrospectives studies. Such defects have been identified in 12% of mCRPC patients. Several studies are ongoing to assess the benefit of these patients from platinum-based chemotherapy and PARP inhibitors, but no conclusive data are available with regards to currently approved therapies for mCRPC, as Abiraterone or Enzalutamide. Methods: PROREPAIR-B (NCT03075735) is a prospective multicentre observational cohort study. Patients diagnosed with mCRPC, with unknown mutational status at study entry and who were going to start a first-line treatment for mCRPC were eligible. For this sub-analysis patients who received Abiraterone or Enzalutamide as first androgen receptor targeted therapy (ART) were selected. The endpoints of this sub-analysis included to assess the impact of BRCA1, BRCA2, ATM, PALB2 and other germline mutations in DNA repair genes on cause-specific survival (CSS), progression-free survival (PFS), time to PSA progression (bPFS) and response to the first ART received as 1st or 2nd line therapy. Results: 337 patients were eligible for this analysis. CSS from mCRPC was not significantly different between gDDR carriers and non-carriers. However, CSS from mCRPC in BRCA2 carriers was significantly shorter than in non-carriers (23.3 Vs 34.6 months, p = 0.02). CSS from first ART, PFS and response-rates were not significantly different between both groups. However, the bPFS was significantly shorter in patients harbouring gDDR mutations (7.3 Vs 3.8 months, p = 0.04), especially in BRCA2 carriers (7.3 Vs 3.0 months, p = 0.03). Conclusions: This is the first study to prospectively follow-up DNA repair germline mutations to determine the outcome on standard treatment for mCRPC. The results suggest that different gDDR defects may have different impact on mCRPC outcomes. Clinical trial information: NCT03075735.

Primordial dwarfism encompasses rare conditions characterized by severe intrauterine growth retardation and growth deficiency throughout life. Recently, three POC1A mutations have been reported in six families with the primordial dwarfism, SOFT syndrome (Short stature, Onychodysplasia, Facial dysmorphism, and hypoTrichosis). Using a custom-designed Next-generation sequencing skeletal dysplasia panel, we have identified two novel homozygous POC1A mutations in two individuals with primordial dwarfism. The severe growth retardation and the facial profiles are strikingly similar between our patients and those described previously. However, one of our patients was diagnosed with severe foramen magnum stenosis and subglottic tracheal stenosis, malformations not previously associated with this syndrome. Our findings confirm that POC1A mutations cause SOFT syndrome and that mutations in this gene should be considered in patients with severe pre- and postnatal short stature, symmetric shortening of long bones, triangular facies, sparse hair and short, thickened distal phalanges.

Alternative splicing of messenger RNA can generate a wide variety of mature RNA transcripts, and these transcripts may produce protein isoforms with diverse cellular functions. While there is much supporting evidence for the expression of alternative transcripts, the same is not true for the alternatively spliced protein products. Large-scale mass spectroscopy experiments have identified evidence of alternative splicing at the protein level, but with conflicting results. Here we carried out a rigorous analysis of the peptide evidence from eight large-scale proteomics experiments to assess the scale of alternative splicing that is detectable by high-resolution mass spectroscopy. We find fewer splice events than would be expected: we identified peptides for almost 64% of human protein coding genes, but detected just 282 splice events. This data suggests that most genes have a single dominant isoform at the protein level. Many of the alternative isoforms that we could identify were only subtly different from the main splice isoform. Very few of the splice events identified at the protein level disrupted functional domains, in stark contrast to the two thirds of splice events annotated in the human genome that would lead to the loss or damage of functional domains. The most striking result was that more than 20% of the splice isoforms we identified were generated by substituting one homologous exon for another. This is significantly more than would be expected from the frequency of these events in the genome. These homologous exon substitution events were remarkably conserved—all the homologous exons we identified evolved over 460 million years ago—and eight of the fourteen tissue-specific splice isoforms we identified were generated from homologous exons. The combination of proteomics evidence, ancient origin and tissue-specific splicing indicates that isoforms generated from homologous exons may have important cellular roles., Author Summary Alternative splicing is thought to be one means for generating the protein diversity necessary for the whole range of cellular functions. While the presence of alternatively spliced transcripts in the cell has been amply demonstrated, the same cannot be said for alternatively spliced proteins. The quest for alternative protein isoforms has focused primarily on the analysis of peptides from large-scale mass spectroscopy experiments, but evidence for alternative isoforms has been patchy and contradictory. A careful analysis of the peptide evidence is needed to fully understand the scale of alternative splicing detectable at the protein level. Here we analysed peptides from eight large-scale data sets, identifying just 282 splice events among 12,716 genes. This suggests that most genes have a single dominant isoform. Many of the alternative isoforms that we identified were only subtly different from the main splice variant, and one in five was generated by substitution of homologous exons by swapping one related exon for another. Remarkably, the alternative isoforms generated from homologous exons were highly conserved, first appearing 460 million years ago, and several appear to have tissue-specific roles in the brain and heart. Our results suggest that these particular isoforms are likely to have important cellular roles.

Protein kinases are involved in relevant physiological functions and a broad number of mutations in this superfamily have been reported in the literature to affect protein function and stability. Unfortunately, the exploration of the consequences on the phenotypes of each individual mutation remains a considerable challenge. The wKinMut web-server offers direct prediction of the potential pathogenicity of the mutations from a number of methods, including our recently developed prediction method based on the combination of information from a range of diverse sources, including physicochemical properties and functional annotations from FireDB and Swissprot and kinase-specific characteristics such as the membership to specific kinase groups, the annotation with disease-associated GO terms or the occurrence of the mutation in PFAM domains, and the relevance of the residues in determining kinase subfamily specificity from S3Det. This predictor yields interesting results that compare favourably with other methods in the field when applied to protein kinases. Together with the predictions, wKinMut offers a number of integrated services for the analysis of mutations. These include: the classification of the kinase, information about associations of the kinase with other proteins extracted from iHop, the mapping of the mutations onto PDB structures, pathogenicity records from a number of databases and the classification of mutations in large-scale cancer studies. Importantly, wKinMut is connected with the SNP2L system that extracts mentions of mutations directly from the literature, and therefore increases the possibilities of finding interesting functional information associated to the studied mutations. wKinMut facilitates the exploration of the information available about individual mutations by integrating prediction approaches with the automatic extraction of information from the literature (text mining) and several state-of-the-art databases. wKinMut has been used during the last year for the analysis of the consequences of mutations in the context of a number of cancer genome projects, including the recent analysis of Chronic Lymphocytic Leukemia cases and is publicly available at http://wkinmut.bioinfo.cnio.es .

e20068Background: This is an innovative project that aims to implement NGS to detect new genetic alterations in the early stages of NSCLC and their potential role as biomarkers of the tumoral evolu...

Chimeric RNAs comprise exons from two or more different genes and have the potential to encode novel proteins that alter cellular phenotypes. To date, numerous putative chimeric transcripts have been identified among the ESTs isolated from several organisms and using high throughput RNA sequencing. The few corresponding protein products that have been characterized mostly result from chromosomal translocations and are associated with cancer. Here, we systematically establish that some of the putative chimeric transcripts are genuinely expressed in human cells. Using high throughput RNA sequencing, mass spectrometry experimental data, and functional annotation, we studied 7424 putative human chimeric RNAs. We confirmed the expression of 175 chimeric RNAs in 16 human tissues, with an abundance varying from 0.06 to 17 RPKM (Reads Per Kilobase per Million mapped reads). We show that these chimeric RNAs are significantly more tissue-specific than non-chimeric transcripts. Moreover, we present evidence that chimeras tend to incorporate highly expressed genes. Despite the low expression level of most chimeric RNAs, we show that 12 novel chimeras are translated into proteins detectable in multiple shotgun mass spectrometry experiments. Furthermore, we confirm the expression of three novel chimeric proteins using targeted mass spectrometry. Finally, based on our functional annotation of exon organization and preserved domains, we discuss the potential features of chimeric proteins with illustrative examples and suggest that chimeras significantly exploit signal peptides and transmembrane domains, which can alter the cellular localization of cognate proteins. Taken together, these findings establish that some chimeric RNAs are translated into potentially functional proteins in humans. The work of M.F-M. is supported by the CNIO (Caja Navarra International Postdoctoral Program) and the Miguel Servet (FIS) grant. The work of V.L. is supported by the French ANR MIRI BLAN08-1335497 Project and the ERC Advanced Grant Sisyphe. This study is supported by the Spanish National Bioinformatics Institute (INB-ISCIII) and by grants from the Spanish Government (CONSOLIDER CSD2007-00050, BIO2007-666855, and BIO2011-26205), European Commission FP7 project ASSET (HEALTH-F4-2010-259348), and the NHGRI-NIH ENCODE grants (U54 HG00455-04 and U54 HG004557)

Most of the many mutations described in human protein kinases are tolerated without significant disruption of the corresponding structures or molecular functions, while some of them have been associated to a variety of human diseases, including cancer. In the last decade, a plethora of computational methods to predict the effect of missense single-nucleotide variants (SNVs) have been developed. Still, current high-throughput sequencing efforts and the concomitant need for massive interpretation of protein sequence variants will demand for more efficient and/or accurate computational methods in the forthcoming years. We present KinMut, a support vector machine (SVM) approach, to identify pathogenic mutations in the protein kinase superfamily. KinMut relays on a combination of sequence-derived features that describe mutations at different levels: (1) Gene level: membership to a specific group in Kinbase and the annotation with GO terms; (2) Domain level: annotated PFAM domains; and (3) Residue level: physicochemical features of amino acids, specificity determining positions, and functional annotations from SwissProt and FireDB. The system has been trained with the set of 3492 human kinase mutations in UniProt for which experimental validation of their pathogenic or neutral character exists. In addition, we discuss the relative importance of these independent properties and their combination for the development of a kinase-specific predictor. Finally, we compare KinMut with other state-of-the-art prediction methods. Family-specific features appear among the most discriminative information sources, which allow us to produce accurate results in a reliable and very simple way with minimal supervision. Our study aims to broaden the knowledge on the mechanisms by which mutations in the human kinome contribute to disease with a particular focus in cancer. The classifier as well as further documentation is available at http://kinmut.bioinfo.cnio.es/ .

Stuve-Wiedemann syndrome (SWS; MIM 601559) is a rare autosomal recessive disease caused by mutations in the leukemia inhibitor factor receptor gene (LIFR). Common clinical and radiological findings are often observed, and high neonatal mortality occurs due to respiratory distress and hyperthermic episodes. Despite initially considered as a lethal disorder during the newborn period, in recent years, several SWS childhood survivors have been reported. We report a detailed clinical and radiological characterization of four unrelated childhood SWS molecularly confirmed patients and review 22 previously reported childhood surviving cases. We contribute to the definition of the childhood survival phenotype of SWS, emphasizing the evolving phenotype, characterized by skeletal abnormalities with typical radiological findings, distinctive dysmorphic features, and dysautonomia. Based on the typical features and clinical course, early diagnosis is possible and crucial to plan appropriate management and prevent potential complications. Genetic confirmation is advisable in order to improve genetic counseling to the patients and their families.