Diagnostic utility of next-generation sequencing-based panel testing in 543 patients with suspected skeletal dysplasia.

BACKGROUND: Skeletal dysplasia is typically diagnosed using a combination of radiographic imaging, clinical examinations, and molecular testing. Identifying a molecular diagnosis for an individual with a skeletal dysplasia can lead to improved clinical care, guide future medical management and treatment, and inform assessment of risk for familial recurrence. The molecular diagnostic utility of multi-gene panel testing using next-generation sequencing (NGS) has not yet been characterized for an unselected population of individuals with suspected skeletal dysplasia. In this study, we retrospectively reviewed patient reports to assess the diagnostic yield, reported variant characteristics, impact of copy number variation, and performance in prenatal diagnostics of panel tests for variants in genes associated with skeletal dysplasia and growth disorders. RESULTS: Clinical reports of consecutive patients with a clinical indication of suspected skeletal dysplasia who underwent panel testing were examined. The 543 patients included in the study submitted samples for diagnostic genetic testing with an indication of suspected skeletal dysplasia or growth disorder and received one of three nested panel tests. A molecular diagnosis was established in 42.0% of patients (n = 228/543). Diagnostic variants were identified in 71 genes, nearly half of which (n = 35, 49.3%) contributed uniquely to a molecular diagnosis for a single patient in this cohort. Diagnostic yield was significantly higher among fetal samples (58.0%, n = 51/88) than postnatal samples (38.9%, n = 177/455; z = 3.32, p < 0.0009). Diagnostic variants in fetal cases were identified across 18 genes. Thirteen diagnostic CNVs were reported, representing 5.7% of diagnostic findings and ranging in size from 241-bp to whole chromosome aneuploidy. Additionally, 11.4% (36/315) of non-diagnostic patient reports had suspicious variants of unknown significance (VUS), in which additional family studies that provide segregation data and/or functional characterization may result in reclassification to likely pathogenic. CONCLUSIONS: These findings demonstrate the utility of panel testing for individuals with a suspected skeletal dysplasia or growth disorder, with a particularly high diagnostic yield seen in prenatal cases. Pursuing comprehensive panel testing with high-resolution CNV analysis can provide a diagnostic benefit, given the considerable phenotype overlap amongst skeletal dysplasia conditions.

Diagnostic yield of genetic testing in a heterogeneous cohort of 1376 HCM patients.

BACKGROUND: Genetic testing in hypertrophic cardiomyopathy (HCM) is a published guideline-based recommendation. The diagnostic yield of genetic testing and corresponding HCM-associated genes have been largely documented by single center studies and carefully selected patient cohorts. Our goal was to evaluate the diagnostic yield of genetic testing in a heterogeneous cohort of patients with a clinical suspicion of HCM, referred for genetic testing from multiple centers around the world. METHODS: A retrospective review of patients with a suspected clinical diagnosis of HCM referred for genetic testing at Blueprint Genetics was undertaken. The analysis included syndromic, myopathic and metabolic etiologies. Genetic test results and variant classifications were extracted from the database. Variants classified as pathogenic (P) or likely pathogenic (LP) were considered diagnostic. RESULTS: A total of 1376 samples were analyzed. Three hundred and sixty-nine tests were diagnostic (26.8%); 373 P or LP variants were identified. Only one copy number variant was identified. The majority of diagnostic variants involved genes encoding the sarcomere (85.0%) followed by 4.3% of diagnostic variants identified in the RASopathy genes. Two percent of diagnostic variants were in genes associated with a cardiomyopathy other than HCM or an inherited arrhythmia. Clinical variables that increased the likelihood of identifying a diagnostic variant included: an earlier age at diagnosis (p < 0.0001), a higher maximum wall thickness (MWT) (p < 0.0001), a positive family history (p < 0.0001), the absence of hypertension (p = 0.0002), and the presence of an implantable cardioverter-defibrillator (ICD) (p = 0.0004). CONCLUSION: The diagnostic yield of genetic testing in this heterogeneous cohort of patients with a clinical suspicion of HCM is lower than what has been reported in well-characterized patient cohorts. We report the highest yield of diagnostic variants in the RASopathy genes identified in a laboratory cohort of HCM patients to date. The spectrum of genes implicated in this unselected cohort highlights the importance of pre-and post-test counseling when offering genetic testing to the broad HCM population.

Discovery of a new mutation in the desmin gene in a young patient with cardiomyopathy and muscular weakness.

A 25-year-old woman with a five years history of syncope, mild left ventricular hypertrophy and moderately enlarged atria, was diagnosed with third degree atrioventricular heart block alternating with atrioventricular heart block 2:1, and received a dual chamber pacemaker. After three years of evolution, she developed atrial fibrillation, marked biatrial enlargement, severely depressed longitudinal myocardial velocities, associated with mild girdle weakness and slight increase in creatine kinase level. The diagnosis of restrictive cardiomyopathy with mild skeletal myopathy imposed the screening for a common etiology. Skeletal muscle biopsy revealed the morphological picture of myofibrillar myopathy with sarcoplasmic aggregates, immunoreactive for desmin and other ectopic proteins on immunohistochemistry, appearing as granulofilamentous material at ultrastructural level. Western blot analysis confirmed the desmin overexpression. Genetic testing identified a heterozygous missense variant DES rs869025381, c.1297C>A, p.(Pro433Thr), not previously reported. This is not only the first confirmed Romanian patient with myofibrillar myopathy with clinical features of severe restrictive cardiomyopathy associated with mild skeletal myopathy, but also a case which adds up to the known mutational spectrum in desminopathy.

The rare Costello variant HRAS c.173C>T (p.T58I) with severe neonatal hypertrophic cardiomyopathy.

We report a 10-year-old girl presenting with severe neonatal hypertrophic cardiomyopathy (HCM), feeding difficulties, mildly abnormal facial features, and progressive skeletal muscle symptoms but with normal cognitive development. Targeted oligonucleotide-selective sequencing of 101 cardiomyopathy genes revealed the genetic diagnosis, and the mutation was verified by Sanger sequencing in the patient and her parents. To offer insights into the potential mechanism of patient mutation, protein structural analysis was performed using the resolved structure of human activated HRAS protein with bound GTP analogue (PDB id 5P21) in Discovery Studio 4.5 (Dassault Systèmes Biovia, San Diego, CA). The patient with hypertrophic cardiomyopathy and normal cognitive development was diagnosed with an HRAS mutation c.173C>T (p.T58I), a milder variant of Costello syndrome affecting a highly conserved amino acid, threonine 58. Our analysis suggests that the p.G12 mutations slow GTP hydrolysis rendering HRAS unresponsive to GTPase activating proteins, and resulting in permanently active state. The p.T58I mutation likely affects binding of guanidine-nucleotide-exchange factors, thereby promoting the active state but also allowing for slow inactivation. Patients with the HRAS mutation c.173C>T (p.T58I) might go undiagnosed because of the milder phenotype compared with other mutations causing Costello syndrome. We expand the clinical and molecular picture of the rare HRAS mutation by reporting the first case in Europe and the fourth case in the literature. Our protein structure analysis offers insights into the mechanism of the mildly activating p.T58I mutation. © 2016 Wiley Periodicals, Inc.

Cardiovascular magnetic resonance findings in patients with PRKAG2 gene mutations.

BACKGROUND: Autosomal dominantly inherited PRKAG2 cardiac syndrome is due to a unique defect of the cardiac cell metabolism and has a distinctive histopathology with excess intracellular glycogen, and prognosis different from sarcomeric hypertrophic cardiomyopathy. We aimed to define the distinct characteristics of PRKAG2 using cardiovascular magnetic resonance (CMR). METHODS: CMR (1.5 T) and genetic testing were performed in two families harboring PRKAG2 mutations. On CMR, segmental analysis of left ventricular (LV) hypertrophy (LVH), function, native T1 mapping, and late gadolinium enhancement (LGE) were performed. RESULTS: Six individuals (median age 23 years, range 16-48; two females) had a PRKAG2 mutation: five with an R302Q mutation (family 1), and one with a novel H344P mutation (family 2). Three of six mutation carriers had LV mass above age and gender limits (203 g/m2, 157 g/m2 and 68 g/m2) and others (with R302Q mutation) normal LV masses. All mutation carriers had LVH in at least one segment, with the median maximal wall thickness of 13 mm (range 11-37 mm). Two R302Q mutation carriers with markedly increased LV mass (203 g/m2 and 157 g/m2) showed a diffuse pattern of hypertrophy but predominantly in the interventricular septum, while other mutation carriers exhibited a non-symmetric mid-infero-lateral pattern of hypertrophy. In family 1, the mutation negative male had a mean T1 value of 963 ms, three males with the R302Q mutation, LVH and no LGE a mean value of 918 ± 11 ms, and the oldest male with the R302Q mutation, extensive hypertrophy and LGE a mean value of 973 ms. Of six mutations carriers, two with advanced disease had LGE with 11 and 22 % enhancement of total LV volume. CONCLUSIONS: PRKAG2 cardiac syndrome may present with eccentric distribution of LVH, involving focal mid-infero-lateral pattern in the early disease stage, and more diffuse pattern but focusing on interventricular septum in advanced cases. In patients at earlier stages of disease, without LGE, T1 values may be reduced, while in the advanced disease stage T1 mapping may result in higher values caused by fibrosis. CMR is a valuable tool in detecting diffuse and focal myocardial abnormalities in PRKAG2 cardiomyopathy.

Accurate genetic diagnosis of Finnish pulmonary arterial hypertension patients using oligonucleotide-selective sequencing.

The genetic basis of pulmonary arterial hypertension (PAH) among Finnish PAH patients is poorly understood. We adopted a novel-targeted next-generation sequencing (NGS) approach called Oligonucleotide-Selective Sequencing (OS-Seq) and developed a custom data analysis and interpretation pipeline to identify pathogenic base substitutions, insertions, and deletions in seven genes associated with PAH (BMPR2, BMPR1B, ACVRL1, ENG, SMAD9, CAV1, and KCNK3) from Finnish PAH patients. This study represents the first clinical study with OS-Seq technology on patients suffering from a rare genetic disorder. We analyzed DNA samples from 21 Finnish PAH patients, whose BMPR2 and ACVRL1 mutation status had been previously studied using Sanger sequencing. Our sequencing panel covered 100% of the targeted base pairs with >15× sequencing depth. Pathogenic base substitutions were identified in the BMPR2 gene in 29% of the Finnish PAH cases. Two of the pathogenic variant-positive patients had been previously tested negative using Sanger sequencing. No clinically significant variants were identified in the six other PAH genes. Our study validates the use of targeted OS-Seq for genetic diagnostics of PAH and revealed pathogenic variants that had been previously missed using Sanger sequencing.

Genetics and genotype-phenotype correlations in Finnish patients with dilated cardiomyopathy.

AIMS: Despite our increased understanding of the genetic basis of dilated cardiomyopathy (DCM), the clinical utility and yield of clinically meaningful findings of comprehensive next-generation sequencing (NGS)-based genetic diagnostics in DCM has been poorly described. We utilized a high-quality oligonucleotide-selective sequencing (OS-Seq)-based targeted sequencing panel to investigate the genetic landscape of DCM in Finnish population and to evaluate the utility of OS-Seq technology as a novel comprehensive diagnostic tool. METHODS AND RESULTS: Using OS-Seq, we targeted and sequenced the coding regions and splice junctions of 101 genes associated with cardiomyopathies in 145 unrelated Finnish patients with DCM. We developed effective bioinformatic variant filtering strategy and implemented strict variant classification scheme to reveal diagnostic yield and genotype-phenotype correlations. Implemented OS-Seq technology provided high coverage of the target region (median coverage 410× and 99.42% of the nucleotides were sequenced at least 15× read depth). Diagnostic yield was 35.2% (familial 47.6% and sporadic 25.6%, P = 0.004) when both pathogenic and likely pathogenic variants are considered as disease causing. Of these, 20 (53%) were titin (TTN) truncations (non-sense and frameshift) affecting all TTN transcripts. TTN truncations accounted for 20.6% and 14.6% of the familial and sporadic DCM cases, respectively. CONCLUSION: Panel-based, high-quality NGS enables high diagnostic yield especially in the familial form of DCM, and bioinformatic variant filtering is a reliable step in the process of interpretation of genomic data in a clinical setting.

Expanding the phenotype of Timothy syndrome type 2: an adolescent with ventricular fibrillation but normal development.

Timothy syndrome is a rare multiorgan disorder with prolonged QTc interval, congenital heart defects, syndactyly, typical facial features and neurodevelopmental problems. Ventricular tachyarrhythmia is the leading cause of death at early age. Classical Timothy syndrome type 1 (TS1) results from a recurrent de novo CACNA1C mutation, G406R in exon 8 A. An atypical form of Timothy syndrome type 2 (TS2) is caused by mutations in G406R and G402S in the alternatively spliced exon 8. Only one individual for each exon 8 mutations has been described. In contrast to multiorgan disease caused by the mutation in G406R either in exon 8 A or 8, the G402S carrier manifested only an isolated cardiac phenotype with LQTS and cardiac arrest. We describe a teenage patient resuscitated from ventricular fibrillation and treated with an implantable cardioverter defibrillator. She has no other organ manifestations, no syndactyly, normal neurodevelopment and her QTc has ranged between 440-480 ms. There is no family history of arrhythmias or sudden death. Targeted oligonucleotide-selective sequencing (OS-Seq) of channelopathy genes revealed a de novo substitution, G402S in exon 8 of CACNA1C. Direct sequencing of blood and saliva derived DNA showed an identical mutation peak suggesting ubiquitous expression in different tissues. The phenotype of our patient and the previously described patient show an isolated arrhythmia disease with no other organ manifestations of classical Timothy syndrome.

Deletion of TOP3ß, a component of FMRP-containing mRNPs, contributes to neurodevelopmental disorders.

Implicating particular genes in the generation of complex brain and behavior phenotypes requires multiple lines of evidence. The rarity of most high-impact genetic variants typically precludes the possibility of accruing statistical evidence that they are associated with a given trait. We found that the enrichment of a rare chromosome 22q11.22 deletion in a recently expanded Northern Finnish sub-isolate enabled the detection of association between TOP3B and both schizophrenia and cognitive impairment. Biochemical analysis of TOP3ß revealed that this topoisomerase was a component of cytosolic messenger ribonucleoproteins (mRNPs) and was catalytically active on RNA. The recruitment of TOP3ß to mRNPs was independent of RNA cis-elements and was coupled to the co-recruitment of FMRP, the disease gene product in fragile X mental retardation syndrome. Our results indicate a previously unknown role for TOP3ß in mRNA metabolism and suggest that it is involved in neurodevelopmental disorders.

Maps of open chromatin guide the functional follow-up of genome-wide association signals: application to hematological traits.

Turning genetic discoveries identified in genome-wide association (GWA) studies into biological mechanisms is an important challenge in human genetics. Many GWA signals map outside exons, suggesting that the associated variants may lie within regulatory regions. We applied the formaldehyde-assisted isolation of regulatory elements (FAIRE) method in a megakaryocytic and an erythroblastoid cell line to map active regulatory elements at known loci associated with hematological quantitative traits, coronary artery disease, and myocardial infarction. We showed that the two cell types exhibit distinct patterns of open chromatin and that cell-specific open chromatin can guide the finding of functional variants. We identified an open chromatin region at chromosome 7q22.3 in megakaryocytes but not erythroblasts, which harbors the common non-coding sequence variant rs342293 known to be associated with platelet volume and function. Resequencing of this open chromatin region in 643 individuals provided strong evidence that rs342293 is the only putative causative variant in this region. We demonstrated that the C- and G-alleles differentially bind the transcription factor EVI1 affecting PIK3CG gene expression in platelets and macrophages. A protein-protein interaction network including up- and down-regulated genes in Pik3cg knockout mice indicated that PIK3CG is associated with gene pathways with an established role in platelet membrane biogenesis and thrombus formation. Thus, rs342293 is the functional common variant at this locus; to the best of our knowledge this is the first such variant to be elucidated among the known platelet quantitative trait loci (QTLs). Our data suggested a molecular mechanism by which a non-coding GWA index SNP modulates platelet phenotype.

Mutant CHUK and severe fetal encasement malformation.

We report an autosomal recessive lethal syndrome characterized by multiple fetal malformations, the most obvious anomalies being the defective face and seemingly absent limbs, which are bound to the trunk and encased under the skin. We identified the molecular defect that causes this syndrome, using a combined strategy of gene-expression arrays, candidate-gene analysis, clinical studies, and genealogic investigations. A point mutation in two affected fetuses led to the loss of the conserved helix­loop­helix ubiquitous kinase (CHUK), also known as IκB kinase α. CHUK has an essential role in the development of skin epidermis and its derivatives, along with various other morphogenetic events. (Funded by the Academy of Finland and others.).

Mutation spectrum of Meckel syndrome genes: one group of syndromes or several distinct groups?

Meckel syndrome (MKS) is a lethal malformation syndrome that belongs to the group of disorders that are associated with primary cilia dysfunction. Total of five genes are known to be involved in the molecular background of MKS. Here we have systematically analyzed all these genes in a total of 29 MKS families. Seven of the families were Finnish and the rest originated from elsewhere in Europe. We found 12 novel mutations in 13 families. Mutations in the MKS genes are also found in other syndromes and it seems reasonable to assume that there is a correlation between the syndromes and the mutations. To obtain some supportive information, we collected all the previously published mutations in the genes to see whether the different syndromes are dictated by the nature of the mutations. Based on this study, mutations play a role in the clinical phenotype, given that the same allelic combination of mutations has never been reported in two clinically distinct syndromes.

Identification of CC2D2A as a Meckel syndrome gene adds an important piece to the ciliopathy puzzle.

Meckel syndrome (MKS) is a lethal malformation disorder characterized classically by encephalocele, polycystic kidneys, and polydactyly. MKS is also one of the major contributors to syndromic neural tube defects (NTDs). Recent findings have shown primary cilia dysfunction in the molecular background of MKS, indicating that cilia are critical for early human development. However, even though four genes behind MKS have been identified to date, they elucidate only a minor proportion of the MKS cases. In this study, instead of traditional linkage analysis, we selected 10 nonrelated affected fetuses and looked for the homozygous regions shared by them. Based on this strategy, we identified the sixth locus and the fifth gene, CC2D2A (MKS6), behind MKS. The biological function of CC2D2A is uncharacterized, but the corresponding polypeptide is predicted to be involved in ciliary functions and it has a calcium binding domain (C2). Immunofluorescence staining of patient's fibroblast cells demonstrates that the cells lack cilia, providing evidence for the critical role of CC2D2A in cilia formation. Our finding is very significant not only to understand the molecular background of MKS, but also to obtain additional information about the function of the cilia, which can help to understand their significance in normal development and also in other ciliopathies, which are an increasing group of disorders with overlapping phenotypes.

Mutations in mRNA export mediator GLE1 result in a fetal motoneuron disease.

The most severe forms of motoneuron disease manifest in utero are characterized by marked atrophy of spinal cord motoneurons and fetal immobility. Here, we report that the defective gene underlying lethal motoneuron syndrome LCCS1 is the mRNA export mediator GLE1. Our finding of mutated GLE1 exposes a common pathway connecting the genes implicated in LCCS1, LCCS2 and LCCS3 and elucidates mRNA processing as a critical molecular mechanism in motoneuron development and maturation.

MKS1, encoding a component of the flagellar apparatus basal body proteome, is mutated in Meckel syndrome.

Meckel syndrome (MKS) is a severe fetal developmental disorder reported in most populations. The clinical hallmarks are occipital meningoencephalocele, cystic kidney dysplasia, fibrotic changes of the liver and polydactyly. Here we report the identification of a gene, MKS1, mutated in MKS families linked to 17q. Mks1 expression in mouse embryos, as determined by in situ hybridization, agrees well with the tissue phenotype of MKS. Comparative genomics and proteomics data implicate MKS1 in ciliary functions.

Evaluation of RAD50 in familial breast cancer predisposition.

The genes predisposing to familial breast cancer are largely unknown, but 5 of the 6 known genes are involved in DNA damage repair. RAD50 is part of a highly conserved complex important in recognising, signalling and repairing DNA double-strand breaks. Recently, a truncating mutation in the RAD50 gene, 687delT, was identified in 2 Finnish breast cancer families. To evaluate the contribution of RAD50 to familial breast cancer, we screened the whole coding region for mutations in 435 UK and 46 Finnish familial breast cancer cases. We identified one truncating mutation, Q350X, in one UK family. We screened a further 544 Finnish familial breast cancer cases and 560 controls for the 687delT mutation, which was present in 3 cases (0.5%) and 1 control (0.2%). Neither Q350X nor 687delT segregated with cancer in the families in which they were identified. Functional analyses suggested that RAD50 687delT is a null allele as there was no detectable expression of the mutant protein. However, the wild-type allele was retained and expressed in breast tumors from mutation carriers. The abundance of the full-length RAD50 protein was reduced in carrier lymphoblastoid cells, suggesting a possible haploinsufficiency mechanism. These data indicate that RAD50 mutations are rare in familial breast cancer and either carry no, or a very small, increased risk of cancer. Altogether, these results suggest RAD50 can only be making a very minor contribution to familial breast cancer predisposition in UK and Finland.

Breast cancer patients with p53 Pro72 homozygous genotype have a poorer survival.

PURPOSE: The p53 R72P polymorphism has been suggested to play a role in many cancers, including breast cancer. Our aim was to evaluate association of R72P with breast cancer risk as well as histopathologic features of the breast tumors and survival. EXPERIMENTAL DESIGN: The germ line R72P genotype was defined among 939 Finnish familial and 888 unselected breast cancer patients and 736 healthy population controls. The clinical and biological variables were tested for association by univariate analysis and the effects of several variables on survival by Cox's proportional hazards regression model. RESULTS: The distribution of the genotypes was similar in all groups studied, suggesting no association with breast cancer risk. Unselected breast cancer patients with 72P homozygous genotype presented significantly more often with lobular carcinoma, whereas R72 allele carriers had a significantly higher frequency of ductal carcinomas (P = 0.004). No significant association with other histopathologic variables, like tumor grade, hormone receptor status (estrogen and progesterone receptors), or tumor-node-metastasis stage, was observed. Survival analysis showed that unselected breast cancer patients with 72P homozygous genotype had significantly poorer survival than patients with other genotypes (P = 0.003). This effect on survival was independent of p53 expression in the tumors and multivariate analysis showed that 72P homozygous genotype was overall an independent prognostic factor (risk ratio of death, 2.1; 95% confidence interval, 1.4-3.3; P = 0.001). CONCLUSIONS: These results suggest no effect of either R72P allele on breast cancer risk but a significantly reduced survival for 72P homozygous breast cancer patients. The finding of codon 72 genotype as an independent prognostic marker for breast cancer warrants further studies.