Protein domains provide a new layer of information for classifying human variations in rare diseases.

Introduction: Using the ACMG-AMP guidelines for the interpretation of sequence variants, it remains difficult to meet the criterion associated with the protein domain, PM1, which is assigned in only about 10% of cases, whereas the criteria related to variant frequency, PM2/BA1/BS1, is reported in 50% of cases. To improve the classification of human missense variants using protein domains information, we developed the DOLPHIN system (https://dolphin.mmg-gbit.eu). Methods: We used Pfam alignments of eukaryotes to define DOLPHIN scores to identify protein domain residues and variants that have a significant impact. In parallel, we enriched gnomAD variants frequencies for each domains' residue. These were validated using ClinVar data. Results: We applied this method to all potential human transcripts' variants, resulting in 30.0% being assigned a PM1 label, whereas 33.2% were eligible for a new benign support criterion, BP8. We also showed that DOLPHIN provides an extrapolated frequency for 31.8% of the variants, compared to the original frequency available in gnomAD for 7.6% of them. Discussion: Overall, DOLPHIN allows a simplified use of the PM1 criterion, an expanded application of the PM2/BS1 criteria and the creation of a new BP8 criterion. DOLPHIN could facilitate the classification of amino acid substitutions in protein domains that cover nearly 40% of proteins and represent the sites of most pathogenic variants.

Clinical and Molecular Update on the Fourth Reported Family with Hamamy Syndrome.

We report on 2 cousins, a girl and a boy, born to first-cousin Lebanese parents with Hamamy syndrome, exhibiting developmental delay, intellectual disability, severe telecanthus, abnormal ears, dentinogenesis imperfecta, and bone fragility. Whole-exome sequencing studies performed on the 2 affected individuals and one obligate carrier revealed the presence of a homozygous c.503G>A (p.Arg168His) missense mutation in IRX5 in both sibs, not reported in any other family. Review of the literature and differential diagnoses are discussed.

The lncRNA 44s2 Study Applicability to the Design of 45-55 Exon Skipping Therapeutic Strategy for DMD.

In skeletal muscle, long noncoding RNAs (lncRNAs) are involved in dystrophin protein stabilization but also in the regulation of myocytes proliferation and differentiation. Hence, they could represent promising therapeutic targets and/or biomarkers for Duchenne and Becker muscular dystrophy (DMD/BMD). DMD and BMD are X-linked myopathies characterized by a progressive muscular dystrophy with or without dilatative cardiomyopathy. Two-thirds of DMD gene mutations are represented by deletions, and 63% of patients carrying DMD deletions are eligible for 45 to 55 multi-exons skipping (MES), becoming BMD patients (BMDΔ45-55). We analyzed the genomic lncRNA presence in 38 BMDΔ45-55 patients and characterized the lncRNA localized in introns 44 and 55 of the DMD gene. We highlighted that all four lncRNA are differentially expressed during myogenesis in immortalized and primary human myoblasts. In addition, the lncRNA44s2 was pointed out as a possible accelerator of differentiation. Interestingly, lncRNA44s expression was associated with a favorable clinical phenotype. These findings suggest that lncRNA44s2 could be involved in muscle differentiation process and become a potential disease progression biomarker. Based on these results, we support MES45-55 therapy and propose that the design of the CRISPR/Cas9 MES45-55 assay consider the lncRNA sequences bordering the exonic 45 to 55 deletion.

Identification of a peripheral blood gene signature predicting aortic valve calcification.

Calcific aortic valve disease (CAVD) is a significant cause of illness and death worldwide. Identification of early predictive markers could help optimize patient management. RNA-sequencing was carried out on human fetal aortic valves at gestational weeks 9, 13, and 22 and on a case-control study with adult noncalcified and calcified bicuspid and tricuspid aortic valves. In dimension reduction and clustering analyses, diseased valves tended to cluster with fetal valves at week 9 rather than normal adult valves, suggesting that part of the disease program might be due to reiterated developmental processes. The analysis of groups of coregulated genes revealed predominant immune-metabolic signatures, including innate and adaptive immune responses involving lymphocyte T-cell metabolic adaptation. Cytokine and chemokine signaling, cell migration, and proliferation were all increased in CAVD, whereas oxidative phosphorylation and protein translation were decreased. Discrete immune-metabolic gene signatures were present at fetal stages and increased in adult controls, suggesting that these processes intensify throughout life and heighten in disease. Cellular stress response and neurodegeneration gene signatures were aberrantly expressed in CAVD, pointing to a mechanistic link between chronic inflammation and biological aging. Comparison of the valve RNA-sequencing data set with a case-control study of whole blood transcriptomes from asymptomatic individuals with early aortic valve calcification identified a highly predictive gene signature of CAVD and of moderate aortic valve calcification in overtly healthy individuals. These data deepen and broaden our understanding of the molecular basis of CAVD and identify a peripheral blood gene signature for the early detection of aortic valve calcification.

Hox-dependent coordination of mouse cardiac progenitor cell patterning and differentiation.

Perturbation of addition of second heart field (SHF) cardiac progenitor cells to the poles of the heart tube results in congenital heart defects (CHD). The transcriptional programs and upstream regulatory events operating in different subpopulations of the SHF remain unclear. Here, we profile the transcriptome and chromatin accessibility of anterior and posterior SHF sub-populations at genome-wide levels and demonstrate that Hoxb1 negatively regulates differentiation in the posterior SHF. Spatial mis-expression of Hoxb1 in the anterior SHF results in hypoplastic right ventricle. Activation of Hoxb1 in embryonic stem cells arrests cardiac differentiation, whereas Hoxb1-deficient mouse embryos display premature cardiac differentiation. Moreover, ectopic differentiation in the posterior SHF of embryos lacking both Hoxb1 and its paralog Hoxa1 results in atrioventricular septal defects. Our results show that Hoxb1 plays a key role in patterning cardiac progenitor cells that contribute to both cardiac poles and provide new insights into the pathogenesis of CHD.

Piezo1 is required for outflow tract and aortic valve development.

AIMS: During embryogenesis, the onset of circulatory blood flow generates a variety of hemodynamic forces which reciprocally induce changes in cardiovascular development and performance. It has been known for some time that these forces can be detected by as yet unknown mechanosensory systems which in turn promote cardiogenic events such as outflow tract and aortic valve development. PIEZO1 is a mechanosensitive ion channel present in endothelial cells where it serves to detect hemodynamic forces making it an ideal candidate to play a role during cardiac development. We sought to determine whether PIEZO1 is required for outflow tract and aortic valve development. METHODS AND RESULTS: By analysing heart development in zebrafish we have determined that piezo1 is expressed in the developing outflow tract where it serves to detect hemodynamic forces. Consequently, disrupting Piezo1 signalling leads to defective outflow tract and aortic valve development and indicates this gene may be involved in the etiology of congenital heart diseases. Based on these findings, we analysed genomic data generated from patients who suffer from left ventricular outflow tract obstructions (LVOTO) and identified 3 probands who each harboured potentially pathogenic variants in PIEZO1. Subsequent in vitro and in vivo assays indicates that these variants behave as dominant negatives leading to an inhibition of normal PIEZO1 mechanosensory activity. Expressing these dominant negative PIEZO1 variants in zebrafish endothelium leads to defective aortic valve development. CONCLUSION: These data indicate that the mechanosensitive ion channel piezo1 is required for outflow tract and aortic valve development.

Loss of Cajal bodies in motor neurons from patients with novel mutations in VRK1.

Distal hereditary motor neuropathies (dHMNs) are a heterogeneous group of diseases, resembling Charcot-Marie-Tooth syndromes, but characterized by an exclusive involvement of the motor part of the peripheral nervous system. Here, we describe two new compound heterozygous mutations in VRK1, the vaccinia-related kinase 1 gene, in two siblings from a Lebanese family, affected with dHMN associated with upper motor neurons (MNs) signs. The mutations lead to severely reduced levels of VRK1 by impairing its stability, and to a shift of nuclear VRK1 to cytoplasm. Depletion of VRK1 from the nucleus alters the dynamics of coilin, a phosphorylation target of VRK1, by reducing its stability through increased proteasomal degradation. In human-induced pluripotent stem cell-derived MNs from patients, we demonstrate that this drop in VRK1 levels leads to Cajal bodies (CBs) disassembly and to defects in neurite outgrowth and branching. Mutations in VRK1 have been previously reported in several neurological diseases affecting lower or both upper and lower MNs. Here, we describe a new phenotype linked to VRK1 mutations, presenting as a classical slowly progressive motor neuropathy, beginning in the second decade of life, with associated upper MN signs. We provide, for the first time, evidence for a role of VRK1 in regulating CB assembly in MNs. The observed MN defects are consistent with a length dependent axonopathy affecting lower and upper MNs, and we propose that diseases due to mutations in VRK1 should be grouped under a unique entity named `VRK1-related motor neuron disease'.

Human pre-valvular endocardial cells derived from pluripotent stem cells recapitulate cardiac pathophysiological valvulogenesis.

Genetically modified mice have advanced our understanding of valve development and disease. Yet, human pathophysiological valvulogenesis remains poorly understood. Here we report that, by combining single cell sequencing and in vivo approaches, a population of human pre-valvular endocardial cells (HPVCs) can be derived from pluripotent stem cells. HPVCs express gene patterns conforming to the E9.0 mouse atrio-ventricular canal (AVC) endocardium signature. HPVCs treated with BMP2, cultured on mouse AVC cushions, or transplanted into the AVC of embryonic mouse hearts, undergo endothelial-to-mesenchymal transition and express markers of valve interstitial cells of different valvular layers, demonstrating cell specificity. Extending this model to patient-specific induced pluripotent stem cells recapitulates features of mitral valve prolapse and identified dysregulation of the SHH pathway. Concurrently increased ECM secretion can be rescued by SHH inhibition, thus providing a putative therapeutic target. In summary, we report a human cell model of valvulogenesis that faithfully recapitulates valve disease in a dish.

The added value of WES reanalysis in the field of genetic diagnosis: lessons learned from 200 exomes in the Lebanese population.

BACKGROUND: The past few decades have witnessed a tremendous development in the field of genetics. The implementation of next generation sequencing (NGS) technologies revolutionized the field of molecular biology and made the genetic information accessible at a large scale. However, connecting a rare genetic variation to a complex phenotype remains challenging. Indeed, identifying the cause of a genetic disease requires a multidisciplinary approach, starting with the establishment of a clear phenotype with a detailed family history and ending, in some cases, with functional assays that are crucial for the validation of the pathogenicity of a mutation. METHODS: Two hundred Lebanese patients, presenting a wide spectrum of genetic disorders (neurodevelopmental, neuromuscular or metabolic disorders, etc.), sporadic or inherited, dominant or recessive, were referred, over the last three and a half years, to the Medical Genetics Unit (UGM) of Saint Joseph University (USJ). In order to identify the genetic basis of these diseases, Whole Exome Sequencing (WES), followed by a targeted analysis, was performed for each case. In order to improve the genetic diagnostic yield, WES data, generated during the first 2 years of this study, were reanalyzed for all patients who were left undiagnosed at the genetic level. Reanalysis was based on updated bioinformatics tools and novel gene discoveries. RESULTS: Our initial analysis allowed us to identify the specific genetic mutation causing the disease in 49.5% of the cases, in line with other international studies. Repeated WES analysis enabled us to increase the diagnostics yield to 56%. CONCLUSION: The present article reports the detailed results of both analysis and pinpoints the contribution of WES data reanalysis to an efficient genetic diagnosis. Lessons learned from WES reanalysis and interpretation are also shared.

Molecular diagnosis of inherited peripheral neuropathies by targeted next-generation sequencing: molecular spectrum delineation.

PURPOSE: Inherited peripheral neuropathies (IPN) represent a large heterogenous group of hereditary diseases with more than 100 causative genes reported to date. In this context, targeted next-generation sequencing (NGS) offers the opportunity to screen all these genes with high efficiency in order to unravel the genetic basis of the disease. Here, we compare the diagnostic yield of targeted NGS with our previous gene by gene Sanger sequencing strategy. We also describe several novel likely pathogenic variants. DESIGN AND PARTICIPANTS: We have completed the targeted NGS of 81 IPN genes in a cohort of 123 unrelated patients affected with diverse forms of IPNs, mostly Charcot-Marie-Tooth disease (CMT): 23% CMT1, 52% CMT2, 9% distal hereditary motor neuropathy, 7% hereditary sensory and autonomic neuropathy and 6.5% intermediate CMT. RESULTS: We have solved the molecular diagnosis in 49 of 123 patients (~40%). Among the identified variants, 26 variants were already reported in the literature. In our cohort, the most frequently mutated genes are respectively: MFN2, SH3TC2, GDAP1, NEFL, GAN, KIF5A and AARS. Panel-based NGS was more efficient in familial cases than in sporadic cases (diagnostic yield 49%vs19%, respectively). NGS-based search for copy number variations, allowed the identification of three duplications in three patients and raised the diagnostic yield to 41%. This yield is two times higher than the one obtained previously by gene Sanger sequencing screening. The impact of panel-based NGS screening is particularly important for demyelinating CMT (CMT1) subtypes, for which the success rate reached 87% (36% only for axonal CMT2). CONCLUSION: NGS allowed to identify causal mutations in a shorter and cost-effective time. Actually, targeted NGS is a well-suited strategy for efficient molecular diagnosis of IPNs. However, NGS leads to the identification of numerous variants of unknown significance, which interpretation requires interdisciplinary collaborations between molecular geneticists, clinicians and (neuro)pathologists.

Clinical profile of comorbidity of rare diseases in a Tunisian patient: a case report associating incontinentia pigmenti and Noonan syndrome.

BACKGROUND: Noonan syndrome (NS) is an autosomal dominant multisystem disorder caused by the dysregulation of several genes belonging to the RAS Mitogen Activated Protein Kinase (MAPK) signaling pathway. Incontinentia Pigmenti (IP) is an X-linked, dominantly inherited multisystem disorder. CASE PRESENTATION: This study is the first report of the coexistence of Noonan (NS) and Incontinentia Pigmenti (IP) syndromes in the same patient. We report on the clinical phenotype and molecular characterization of this patient. The patient was examined by a pluridisciplinary staff of clinicians and geneticist. The clinical diagnosis of NS and IP was confirmed by molecular investigations. The newborn girl came to our clinics due to flagrant dysmorphia and dermatological manifestations. The clinical observations led to characterize the Incontinentia Pigmenti traits and a suspicion of a Noonan syndrome association. Molecular diagnosis was performed by Haloplex resequencing of 29 genes associated with RASopathies and confirmed the NS diagnosis. The common recurrent intragenic deletion mutation in IKBKG gene causing the IP was detected with an improved PCR protocol. CONCLUSION: This is the first report in the literature of comorbidity of NS and IP, two rare multisystem syndromes.

Loss of Calmodulin- and Radial-Spoke-Associated Complex Protein CFAP251 Leads to Immotile Spermatozoa Lacking Mitochondria and Infertility in Men.

Flagella and motile cilia share a 9 + 2 microtubule-doublet axoneme structure, and asthenozoospermia (reduced spermatozoa motility) is found in 76% of men with primary ciliary dyskinesia (PCD). Nevertheless, causal genetic variants in a conserved axonemal component have been found in cases of isolated asthenozoospermia: 30% of men with multiple morphological anomalies of sperm flagella (MMAF) carry bi-allelic mutations in DNAH1, encoding one of the seven inner-arm dynein heavy chains of the 9 + 2 axoneme. To further understand the basis for isolated asthenozoospermia, we used whole-exome and Sanger sequencing to study two brothers and two independent men with MMAF. In three men, we found bi-allelic loss-of-function mutations in WDR66, encoding cilia- and flagella-associated protein 251 (CFAP251): the two brothers were homozygous for the frameshift chr12: g.122359334delA (p.Asp42Metfs∗4), and the third individual was compound heterozygous for chr12: g.122359542G>T (p.Glu111∗) and chr12: g.122395032_122395033delCT (p.Leu530Valfs∗4). We show that CFAP251 is normally located along the flagellum but is absent in men carrying WDR66 mutations and reveal a spermatozoa-specific isoform probably generated during spermatozoon maturation. CFAP251 is a component of the calmodulin- and radial-spoke- associated complex, located adjacent to DNAH1, on the inner surface of the peripheral microtubule doublets of the axoneme. In Tetrahymena, the CFAP251 ortholog is necessary for efficient coordinated ciliary beating. Using immunofluorescent and transmission electron microscopy, we provide evidence that loss of CFAP251 affects the formation of the mitochondrial sheath. We propose that CFAP251 plays a structural role during biogenesis of the spermatozoon flagellum in vertebrates.

VarAFT: a variant annotation and filtration system for human next generation sequencing data.

With the rapidly developing high-throughput sequencing technologies known as next generation sequencing or NGS, our approach to gene hunting and diagnosis has drastically changed. In <10 years, these technologies have moved from gene panel to whole genome sequencing and from an exclusively research context to clinical practice. Today, the limit is not the sequencing of one, many or all genes but rather the data analysis. Consequently, the challenge is to rapidly and efficiently identify disease-causing mutations within millions of variants. To do so, we developed the VarAFT software to annotate and pinpoint human disease-causing mutations through access to multiple layers of information. VarAFT was designed both for research and clinical contexts and is accessible to all scientists, regardless of bioinformatics training. Data from multiple samples may be combined to address all Mendelian inheritance modes, cancers or population genetics. Optimized filtration parameters can be stored and re-applied to large datasets. In addition to classical annotations from dbNSFP, VarAFT contains unique features at the disease (OMIM), phenotypic (HPO), gene (Gene Ontology, pathways) and variation levels (predictions from UMD-Predictor and Human Splicing Finder) that can be combined to optimally select candidate pathogenic mutations. VarAFT is freely available at: http://varaft.eu.

Exome sequencing identifies recurrent BCOR alterations and the absence of KLF2, TNFAIP3 and MYD88 mutations in splenic diffuse red pulp small B-cell lymphoma.

Splenic diffuse red pulp lymphoma is an indolent small B-cell lymphoma recognized as a provisional entity in the World Health Organization 2008 classification. Its precise relationship to other related splenic B-cell lymphomas with frequent leukemic involvement or other lymphoproliferative disorders remains undetermined. We performed whole-exome sequencing to explore the genetic landscape of ten cases of splenic diffuse red pulp lymphoma using paired tumor and normal samples. A selection of 109 somatic mutations was then evaluated in a cohort including 42 samples of splenic diffuse red pulp lymphoma and compared to those identified in 46 samples of splenic marginal zone lymphoma and eight samples of hairy-cell leukemia. Recurrent mutations or losses in BCOR (the gene encoding the BCL6 corepressor) - frameshift (n=3), nonsense (n=2), splicing site (n=1), and copy number loss (n=4) - were identified in 10/42 samples of splenic diffuse red pulp lymphoma (24%), whereas only one frameshift mutation was identified in 46 cases of splenic marginal zone lymphoma (2%). Inversely, KLF2, TNFAIP3 and MYD88, common mutations in splenic marginal zone lymphoma, were rare (one KLF2 mutant in 42 samples; 2%) or absent (TNFAIP3 and MYD88) in splenic diffuse red pulp lymphoma. These findings define an original genetic profile of splenic diffuse red pulp lymphoma and suggest that the mechanisms of pathogenesis of this lymphoma are distinct from those of splenic marginal zone lymphoma and hairy-cell leukemia.

Exome Sequencing Identifies Two Variants of the Alkylglycerol Monooxygenase Gene as a Cause of Relapses in Visceral Leishmaniasis in Children, in Sudan.

Background: Visceral leishmaniasis (kala-azar, KA) is the most severe form of leishmaniasis, characterized by fever, weight loss, hepatosplenomegaly, and lymphadenopathy. During an outbreak of KA in Babar El Fugara (Sudan), 5.7% of cured patients displayed relapses, with familial clustering in half the cases. Methods: We performed whole-exome sequencing on 10 relapsing individuals and 11 controls from 5 nuclear families. Results: Rare homozygous and compound-heterozygous nonsense (c.1213C > T, rs139309795, p.Arg405*) and missense (c.701A > G, rs143439626, p.Lys234Arg) mutations of the alkylglycerol monooxygenase (AGMO) gene were associated with KA relapse in 3 families. Sequencing in additional family members confirmed the segregation of these mutations with relapse and revealed an autosomal dominant mode of transmission. These mutations were detected heterozygous in 2 subjects among 100 unrelated individuals with KA who never relapsed after cure, suggesting incomplete penetrance of AGMO deficiency. AGMO is expressed in hematopoietic cells, and is strongly expressed in the liver. AGMO modulates PAF production by mouse macrophages, suggesting that it may act through the PAF/PAF receptor pathway previously shown to have anti-Leishmania activity. Conclusions: This is the first demonstration that relapses after a first episode of KA are due to differences in human genetic susceptibility and not to modifications of parasite pathogenicity.

Shared genetic predisposition in rheumatoid arthritis-interstitial lung disease and familial pulmonary fibrosis.

Despite its high prevalence and mortality, little is known about the pathogenesis of rheumatoid arthritis-associated interstitial lung disease (RA-ILD). Given that familial pulmonary fibrosis (FPF) and RA-ILD frequently share the usual pattern of interstitial pneumonia and common environmental risk factors, we hypothesised that the two diseases might share additional risk factors, including FPF-linked genes. Our aim was to identify coding mutations of FPF-risk genes associated with RA-ILD.We used whole exome sequencing (WES), followed by restricted analysis of a discrete number of FPF-linked genes and performed a burden test to assess the excess number of mutations in RA-ILD patients compared to controls.Among the 101 RA-ILD patients included, 12 (11.9%) had 13 WES-identified heterozygous mutations in the TERT, RTEL1, PARN or SFTPC coding regions. The burden test, based on 81 RA-ILD patients and 1010 controls of European ancestry, revealed an excess of TERT, RTEL1, PARN or SFTPC mutations in RA-ILD patients (OR 3.17, 95% CI 1.53-6.12; p=9.45×10-4). Telomeres were shorter in RA-ILD patients with a TERT, RTEL1 or PARN mutation than in controls (p=2.87×10-2).Our results support the contribution of FPF-linked genes to RA-ILD susceptibility.

New advances in DPYD genotype and risk of severe toxicity under capecitabine.

BACKGROUND: Deficiency in dihydropyrimidine dehydrogenase (DPD) enzyme is the main cause of severe and lethal fluoropyrimidine-related toxicity. Various approaches have been developed for DPD-deficiency screening, including DPYD genotyping and phenotyping. The goal of this prospective observational study was to perform exhaustive exome DPYD sequencing and to examine relationships between DPYD variants and toxicity in advanced breast cancer patients receiving capecitabine. METHODS: Two-hundred forty-three patients were analysed (88.5% capecitabine monotherapy). Grade 3 and grade 4 capecitabine-related digestive and/or neurologic and/or hemato-toxicities were observed in 10.3% and 2.1% of patients, respectively. DPYD exome, along with flanking intronic regions 3'UTR and 5'UTR, were sequenced on MiSeq Illumina. DPD phenotype was assessed by pre-treatment plasma uracil (U) and dihydrouracil (UH2) measurement. RESULTS: Among the 48 SNPs identified, 19 were located in coding regions, including 3 novel variations, each observed in a single patient (among which, F100L and A26T, both pathogenic in silico). Combined analysis of deleterious variants *2A, I560S (*13) and D949V showed significant association with grade 3-4 toxicity (sensitivity 16.7%, positive predictive value (PPV) 71.4%, relative risk (RR) 6.7, p<0.001) but not with grade 4 toxicity. Considering additional deleterious coding variants D342G, S492L, R592W and F100L increased the sensitivity to 26.7% for grade 3-4 toxicity (PPV 72.7%, RR 7.6, p<0.001), and was significantly associated with grade 4 toxicity (sensitivity 60%, PPV 27.3%, RR 31.4, p = 0.001), suggesting the clinical relevance of extended targeted DPYD genotyping. As compared to extended genotype, combining genotyping (7 variants) and phenotyping (U>16 ng/ml) did not substantially increase the sensitivity, while impairing PPV and RR. CONCLUSIONS: Exploring an extended set of deleterious DPYD variants improves the performance of DPYD genotyping for predicting both grade 3-4 and grade 4 toxicities (digestive and/or neurologic and/or hematotoxicities) related to capecitabine, as compared to conventional genotyping restricted to consensual variants *2A, *13 and D949V.

Exome sequencing reveals a de novo POLD1 mutation causing phenotypic variability in mandibular hypoplasia, deafness, progeroid features, and lipodystrophy syndrome (MDPL).

BACKGROUND: Mandibular hypoplasia, deafness, progeroid features, and lipodystrophy syndrome (MDPL) is an autosomal dominant systemic disorder characterized by prominent loss of subcutaneous fat, a characteristic facial appearance and metabolic abnormalities. This syndrome is caused by heterozygous de novo mutations in the POLD1 gene. To date, 19 patients with MDPL have been reported in the literature and among them 14 patients have been characterized at the molecular level. Twelve unrelated patients carried a recurrent in-frame deletion of a single codon (p.Ser605del) and two other patients carried a novel heterozygous mutation in exon 13 (p.Arg507Cys). Additionally and interestingly, germline mutations of the same gene have been involved in familial polyposis and colorectal cancer (CRC) predisposition. PATIENTS AND METHODS: We describe a male and a female patient with MDPL respectively affected with mild and severe phenotypes. Both of them showed mandibular hypoplasia, a beaked nose with bird-like facies, prominent eyes, a small mouth, growth retardation, muscle and skin atrophy, but the female patient showed such a severe and early phenotype that a first working diagnosis of Hutchinson-Gilford Progeria was made. The exploration was performed by direct sequencing of POLD1 gene exon 15 in the male patient with a classical MDPL phenotype and by whole exome sequencing in the female patient and her unaffected parents. RESULTS: Exome sequencing identified in the latter patient a de novo heterozygous undescribed mutation in the POLD1 gene (NM_002691.3: c.3209T>A), predicted to cause the missense change p.Ile1070Asn in the ZnF2 (Zinc Finger 2) domain of the protein. This mutation was not reported in the 1000 Genome Project, dbSNP and Exome sequencing databases. Furthermore, the Isoleucine1070 residue of POLD1 is highly conserved among various species, suggesting that this substitution may cause a major impairment of POLD1 activity. For the second patient, affected with a typical MDPL phenotype, direct sequencing of POLD1 exon 15 revealed the recurrent in-frame deletion (c.1812_1814del, p.S605del). CONCLUSION: Our work highlights that mutations in different POLD1 domains can lead to phenotypic variability, ranging from dominantly inherited cancer predisposition syndromes, to mild MDPL phenotypes without lifespan reduction, to very severe MDPL syndromes with major premature aging features. These results also suggest that POLD1 gene testing should be considered in patients presenting with severe progeroid features.

How to Identify Pathogenic Mutations among All Those Variations: Variant Annotation and Filtration in the Genome Sequencing Era.

High-throughput sequencing technologies have become fundamental for the identification of disease-causing mutations in human genetic diseases both in research and clinical testing contexts. The cumulative number of genes linked to rare diseases is now close to 3,500 with more than 1,000 genes identified between 2010 and 2014 because of the early adoption of Exome Sequencing technologies. However, despite these encouraging figures, the success rate of clinical exome diagnosis remains low due to several factors including wrong variant annotation and nonoptimal filtration practices, which may lead to misinterpretation of disease-causing mutations. In this review, we describe the critical steps of variant annotation and filtration processes to highlight a handful of potential disease-causing mutations for downstream analysis. We report the key annotation elements to gather at multiple levels for each mutation, and which systems are designed to help in collecting this mandatory information. We describe the filtration options, their efficiency, and limits and provide a generic filtration workflow and highlight potential pitfalls through a use case.