Whole Exome Sequencing in 16p13.11 Microdeletion Patients Reveals New Variants Through Deductive and Systems Medicine Approaches.

The 16p13.11 microdeletion, whose prevalence in the general population is about 0.04%, is known in literature as a predisposition factor to neurodevelopmental disorders, being found in about 0.13% of patients with schizophrenia, in 0.5-0.6% of patient with epilepsy, cognitive impairment, autism spectrum disorder (ASD) and aggressiveness. The goal of this study was to identify a specific gene set pattern unique for the affected patients in comparison with other familial components. Due to the incomplete penetrance of this copy number variant (CNV), we studied by whole exome sequencing (WES), with particular regard of 850 SFARI genes, three families with an affected member carrier of inherited 16p13.11 and 16p13.11p12.3 microdeletion and one family with an affected member with a de novo 16p13.11 microdeletion. By combining a deductive approach together with personalized network models, we identified gene signatures potentially capable of explaining the clinical phenotype. Candidate variants in genes of interest were identified as possibly involved in determining the neurological phenotype of the four patients, such as compound heterozygosity in CECR2, variants in MTOR and RICTOR genes, compound heterozygous single nucleotide variants in the LRRK2 gene. Moreover, genes present in the microdeletion region were partially present as central nodes, with a focus on NDE1. No additional pathogenetic or uncertain CNVs were found in all four patients. No significant variants were detected in genes included in the microdeletion in patients 1, 2 and 3, excluding the finding of unmasked recessive variants. In conclusion, WES is a fundamental tool in the genetic investigation of patients having a predisposing variant, which is not sufficient to define the clinical phenotype. Moreover, the analysis of WES data using Systems medicine tools, such as personalized network models, led to the prioritization of genes on a high throughput scale and to discover variants in genes that were not prioritized at first.

Gene Set Enrichment Analysis of Interaction Networks Weighted by Node Centrality.

Gene set enrichment analysis (GSEA) is a powerful tool to associate a disease phenotype to a group of genes/proteins. GSEA attributes a specific weight to each gene/protein in the input list that depends on a metric of choice, which is usually represented by quantitative expression data. However, expression data are not always available. Here, GSEA based on betweenness centrality of a protein-protein interaction (PPI) network is described and applied to two cases, where an expression metric is missing. First, personalized PPI networks were generated from genes displaying alterations (assessed by array comparative genomic hybridization and whole exome sequencing) in four probands bearing a 16p13.11 microdeletion in common and several other point variants. Patients showed disease phenotypes linked to neurodevelopment. All networks were assembled around a cluster of first interactors of altered genes with high betweenness centrality. All four clusters included genes known to be involved in neurodevelopmental disorders with different centrality. Moreover, the GSEA results pointed out to the evidence of "cell cycle" among enriched pathways. Second, a large interaction network obtained by merging proteomics studies on three neurodegenerative disorders was analyzed from the topological point of view. We observed that most central proteins are often linked to Parkinson's disease. The selection of these proteins improved the specificity of GSEA, with "Metabolism of amino acids and derivatives" and "Cellular response to stress or external stimuli" as top-ranked enriched pathways. In conclusion, betweenness centrality revealed to be a suitable metric for GSEA. Thus, centrality-based GSEA represents an opportunity for precision medicine and network medicine.

Targeted re-sequencing in malformations of cortical development: genotype-phenotype correlations.

PURPOSE: Malformations of cortical development (MCD) are a phenotypically and genetically heterogeneous group of disorders, for which the diagnostic rate of genetic testing in a clinical setting remains to be clarified. In this study we aimed to assess the diagnostic rate of germline and pathogenic variants using a custom panel in a heterogeneous group of subjects with MCD and explore genotype-phenotype correlations. METHODS: A total of 84 subjects with different MCD were enrolled. Genomic DNA was isolated from peripheral blood. Fifty-nine tartget genes were assessed using a custom next-generation sequencing (NGS) panel. RESULTS: Genetic causes were identified in one-fourth of our cohort (21.4 %). Overall, we identified 19 pathogenic or likely pathogenic single-nucleotide variants in 11 genes among 18 subjects, including PAFAH1B1 (LIS1) (n = 3), TUBA1A (n = 3), DYNC1H1 (n = 3), ACTG1 (n = 2), TUBB2B (n = 1), TUBB3 (n = 1), DCX (n = 1), FLNA (n = 1), LAMA2 (n = 1), POMGNT2 (n = 1) and VLDLR (n = 1). The diagnostic yield was higher in patients with lissencephaly/pachygyria (60 %) (p = 0.001), cobblestone malformation (50 %), and subcortical band heterotopia (SBH) (40 %). Furthermore, five out of six subjects with suspect tubulinopathies on imaging harboured pathogenic variants in tubulin genes. Overall, germline pathogenic variants were more likely to be identified if MCD were diffuse (p = 0.002) and associated with other central nervous system malformations (p = 0.029). Moderate to severe intellectual disability was also more commonly associated with pathogenic variants (p = 0.044). CONCLUSION: Customized gene panels may support the diagnostic work-up for some specific MCD, especially when these are diffuse, bilateral and associated with other brain malformations.

Testing single/combined clinical categories on 5110 Italian patients with developmental phenotypes to improve array-based detection rate.

BACKGROUND: Chromosomal microarray analysis (CMA) is nowadays widely used in the diagnostic path of patients with clinical phenotypes. However, there is no ascertained evidence to date on how to assemble single/combined clinical categories of developmental phenotypic findings to improve the array-based detection rate. METHODS: The Italian Society of Human Genetics coordinated a retrospective study which included CMA results of 5,110 Italian patients referred to 17 genetics laboratories for variable combined clinical phenotypes. RESULTS: Non-polymorphic copy number variants (CNVs) were identified in 1512 patients (30%) and 615 (32%) present in 552 patients (11%) were classified as pathogenic. CNVs were analysed according to type, size, inheritance pattern, distribution among chromosomes, and association to known syndromes. In addition, the evaluation of the detection rate of clinical subgroups of patients allowed to associate dysmorphisms and/or congenital malformations combined with any other single clinical sign to an increased detection rate, whereas non-syndromic neurodevelopmental signs and non-syndromic congenital malformations to a decreased detection rate. CONCLUSIONS: Our retrospective study resulted in confirming the high detection rate of CMA and indicated new clinical markers useful to optimize their inclusion in the diagnostic and rehabilitative path of patients with developmental phenotypes.

The p.Gly130Val mutation in the GJB2 gene: A familiar case of autosomal dominant non-syndromic hearing loss.

Several forms of sensorineural hearing loss (SNHL) have been imputated to connexins mutations and prevalently to connexin 26 (Cx26), codified by the GJB2 gene (gap junction protein, beta 2). Here, we report the first familiar case (heterozygous p. G130V mutation) of non-syndromic (without any dermatological manifestation) dominant profound SNHL. Proband was a 6-years-old male with post-lingual bilateral profound SNHL, clinically identified at the age of 3 with diagnosis of severe SNHL. We confirm that the p. G130V variant of the GJB2 gene is causative of autosomal dominant form of SNHL, although it is not always associated with the presence of skin diseases.

gDNA qPCR is statistically more reliable than mRNA analysis in detecting leukemic cells to monitor CML.

Chronic Myeloid Leukemia (CML) is a stem cell cancer that arises when t(9;22) translocation occurs in a hematopoietic stem cells. This event results in the expression of the BCR-ABL1 fusion gene, which codes for a constitutively active tyrosine kinase that is responsible for the transformation of a HSC into a CML stem cell, which then gives rise to a clonal myeloproliferative disease. The introduction of Tyrosine Kinase Inhibitors (TKIs) has revolutionized the management of the disease. However, these drugs do not seem to be able to eradicate the malignancy. Indeed, discontinuation trials (STIM; TWISER; DADI) for those patients who achieved a profound molecular response showed 50% relapsing within 12 months. We performed a comparative analysis on 15 CML patients and one B-ALL patient, between the standard quantitative reverse-transcriptase PCR (qRT-PCR) and our genomic DNA patient-specific quantitative PCR assay (gDNA qPCR). Here we demonstrate that gDNA qPCR is better than standard qRT-PCR in disease monitoring after an average follow-up period of 200 days. Specifically, we statistically demonstrated that DNA negativity is more reliable than RNA negativity in indicating when TKIs therapy can be safely stopped.

Optimizing the molecular diagnosis of GALNS: novel methods to define and characterize Morquio-A syndrome-associated mutations.

Morquio A syndrome (MPS IVA) is a systemic lysosomal storage disorder caused by the deficiency of N-acetylgalactosamine-6-sulfatase (GALNS), encoded by the GALNS gene. We studied 37 MPS IV A patients and defined genotype-phenotype correlations based on clinical data, biochemical assays, molecular analyses, and in silico structural analyses of associated mutations. We found that standard sequencing procedures, albeit identifying 14 novel small GALNS genetic lesions, failed to characterize the second disease-causing mutation in the 16% of the patients' cohort. To address this drawback and uncover potential gross GALNS rearrangements, we developed molecular procedures (CNV [copy-number variation] assays, QF-PCRs [quantitative fluorescent-PCRs]), endorsed by CGH-arrays. Using this approach, we characterized two new large deletions and their corresponding breakpoints. Both deletions were heterozygous and included the first exon of the PIEZO1 gene, which is associated with dehydrated hereditary stomatocitosis, an autosomal-dominant syndrome. In addition, we characterized the new GALNS intronic lesion c.245-11C>G causing m-RNA defects, although identified outside the GT/AG splice pair. We estimated the occurrence of the disease in the Italian population to be approximately 1:300,000 live births and defined a molecular testing algorithm designed to help diagnosing MPS IVA and foreseeing disease progression.

Mutations and thrombosis in essential thrombocythemia: prognostic interaction with age and thrombosis history.

BACKGROUND: Vascular events in essential thrombocythemia (ET) are associated with advanced age and thrombosis history. Recent information suggests additional effect from the presence of specific mutations. OBJECTIVES: To examine the influence of age and thrombosis history on the reported association between mutational status and thrombosis-free survival in ET. PATIENTS AND METHODS: Analysis was performed using a Mayo Clinic cohort of 300 ET patients, and key findings were reanalyzed by including additional 102 Italian patients. RESULTS: Among 300 Mayo patients with ET (median age 55 yr, 60% females), mutational frequencies were 53% JAK2, 32% CALR, 3% MPL, and 12% JAK2, CALR and MPL wild type. One hundred and six (35%) patients experienced arterial (n = 75) or venous (n = 43) events, before (n = 55) or after (n = 71) diagnosis. In univariate analysis, compared to JAK2-mutated cases, JAK2, CALR and MPL wild type (HR 0.31, 95% CI 0.11-0.86), and CALR-mutated (0.53, 95% CI 0.30-0.92) patients displayed better thrombosis-free survival. JAK2, CALR, and MPL wild type remained significant (P = 0.03; HR 0.32, 95% CI 0.11-0.9) during multivariable analysis that included age (P = 0.01) and thrombosis history (P = 0.0006); a favorable impact from CALR mutations was of borderline significance (P = 0.1; HR 0.62, 95% CI 0.35-1.1), but became significant (P = 0.02) when multivariable analysis including thrombosis history (P = 0.02) was performed on patients younger than 60 yr of age. CONCLUSIONS: The favorable impact of mutational status on thrombosis-free survival in ET might be most evident for JAK2, CALR, and MPL wild type patients, whereas the favorable effect from CALR mutations might be confined to young patients.

It is time to change thrombosis risk assessment for PV and ET?

Polycythemia vera (PV) and essential thrombocythemia (ET) are myeloproliferative neoplasms to be diagnosed according to the WHO classification. Molecular profiling must include the analysis of JAK2 (looking for the V617F point-mutation in PV and ET, screening exon 12 for mutations only in V617F-negative PV), CALR and MPL mutations (both in V617F-negative ET). The current risk stratification to predict thrombosis requires two parameters: age over 60 years and prior history of thrombosis. On the basis of these two risk factors patients can be stratified in low-risk and high-risk and receive a proper treatment. However, a modern stratification of thrombotic risk might consider "new" low-risk patients: conventional low-risk plus absence of leukocytosis from diagnosis onwards and a hematocrit level below 45% during the course of disease for PV; conventional low-risk plus absence of leukocytosis from diagnosis onwards, JAK2 negativity, CALR positivity, and absence of cardiovascular risk factors for ET.

Type 1 versus Type 2 calreticulin mutations in essential thrombocythemia: a collaborative study of 1027 patients.

CALR (calreticulin) trails JAK2 as the second most mutated gene in essential thrombocythemia (ET). Mutant CALR in ET is a result of frameshift mutations, caused by exon 9 deletions or insertions; type-1, 52-bp deletion (p.L367fs*46), and type-2, 5-bp TTGTC insertion (p.K385fs*47) variants constitute more than 80% of these mutations. The current study includes a total of 1027 patients divided into test (n = 402) and validation (n = 625) cohorts. Among the 402 ET patients in the test cohort, 227 (57%) harbored JAK2, 11 (3%) Myeloproliferative leukemia virus oncogene (MPL), and 114 (28%) CALR mutations; 12% were wild-type for all three mutations (i.e., triple-negative). Among the 114 patients with CALR mutations, 51 (45%) displayed type-1 and 44 (39%) type-2 variants; compared to mutant JAK2, both variants were associated with higher platelet and lower hemoglobin and leukocyte counts. However, male sex was associated with only type-1 (P = 0.005) and younger age with type-2 (P = 0.001) variants. Notably, platelet count was significantly higher in type-2 vs. type-1 CALR-mutated patients (P = 0.03) and the particular observation was validated in the validation cohort that included 111 CALR-mutated ET patients (P = 0.002). These findings, coupled with the recent demonstration of preferential expression of mutant and wild-type CALR in megakaryocytes, suggest differential effects of CALR variants on thrombopoiesis.

Isolation and characterization of multipotent cells from human fetal dermis.

We report that cells from human fetal dermis, termed here multipotent fetal dermal cells, can be isolated with high efficiency by using a nonenzymatic, cell outgrowth method. The resulting cell population was consistent with the definition of mesenchymal stromal cells by the International Society for Cellular Therapy. As multipotent fetal dermal cells proliferate extensively, with no loss of multilineage differentiation potential up to passage 25, they may be an ideal source for cell therapy to repair damaged tissues and organs. Multipotent fetal dermal cells were not recognized as targets by T lymphocytes in vitro, thus supporting their feasibility for allogenic transplantation. Moreover, the expansion protocol did not affect the normal phenotype and karyotype of cells. When compared with adult dermal cells, fetal cells displayed several advantages, including a greater cellular yield after isolation, the ability to proliferate longer, and the retention of differentiation potential. Interestingly, multipotent fetal dermal cells expressed the pluripotency marker SSEA4 (90.56 ± 3.15% fetal vs. 10.5 ± 8.5% adult) and coexpressed mesenchymal and epithelial markers (>80% CD90(+)/CK18(+) cells), coexpression lacking in the adult counterparts isolated under the same conditions. Multipotent fetal dermal cells were able to form capillary structures, as well as differentiate into a simple epithelium in vitro, indicating skin regeneration capabilities.