# Unravelling the genetic causes of multiple malformation syndromes: A whole exome sequencing study of the Cypriot population

Evie Kritioti, Athina Theodosiou, Thibaud Parpaite, Angelos Alexandrou, Nayia Nicolaou, Ioannis Papaevripidou, Nina Séjourné, Bertrand Coste, Violetta Christophidou-Anastasiadou, George Tanteles, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Evie Kritioti, Athina Theodosiou, Thibaud Parpaite, Angelos Alexandrou, Nayia Nicolaou, et al. Unravelling the genetic causes of multiple malformation syndromes: A whole exome sequencing study of the Cypriot population. PLoS ONE, 2021, 16 (7), pp.e0253562. 10.1371/journal.pone.0253562. hal-03396027

HAL Id: hal-03396027

https://amu.hal.science/hal-03396027

Submitted on 22 Oct 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Unravelling the genetic causes of multiple malformation syndromes: A whole exome sequencing study of the Cypriot population

Evie Kritioti¹,²,³, Athina Theodosiou¹,³, Thibaud Parpaite⁴, Angelos Alexandrou¹, Nayia Nicolaou², Ioannis Papaevripidou¹, Nina Séjourné⁴, Bertrand Coste⁴, Violetta Christophidou-Anastasiadou⁵, George A. Tanteles²,³*, Carolina Sismani¹,³*

¹ Department of Cytogenetics and Genomics, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus, ² Clinical Genetics Clinic, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus, ³ The Cyprus School of Molecular Medicine, The Cyprus Institute of Neurology and Genetics, Nicosia, Cyprus, ⁴ Aix Marseille Université, CNRS, LNC-UMR 7291, Marseille, France, ⁵ Clinical Genetics Clinic, The Cyprus Institute of Neurology and Genetics and Archbishop Makarios III Medical Centre, Nicosia, Cyprus

© These authors contributed equally to this work.

* csismani@cing.ac.cy

<figure><img src="image_3.png" /></figure>

## Abstract

Multiple malformation syndromes (MMS) belong to a group of genetic disorders characterised by neurodevelopmental anomalies and congenital malformations. Here we explore for the first time the genetic aetiology of MMS using whole-exome sequencing (WES) in undiagnosed patients from the Greek-Cypriot population after prior extensive diagnostics workup including karyotype and array-CGH. A total of 100 individuals (37 affected), from 32 families were recruited and family-based WES was applied to detect causative single-nucleotide variants (SNVs) and indels. A genetic diagnosis was reported for 16 MMS patients (43.2%), with 10/17 (58.8%) of the findings being novel. All autosomal dominant findings occurred de novo. Functional studies were also performed to elucidate the molecular mechanism relevant to the abnormal phenotypes, in cases where the clinical significance of the findings was unclear. The 17 variants identified in our cohort were located in 14 genes (*PCNT*, *UBE3A*, *KAT6A*, *SPR*, *POMGNT1*, *PIEZ02*, *PXDN*, *KDM6A*, *PHIP*, *HECW2*, *TFAP2A*, *CNOT3*, *AGTPBP1* and *GAMT*). This study has highlighted the efficacy of WES through the high detection rate (43.2%) achieved for a challenging category of undiagnosed patients with MMS compared to other conventional diagnostic testing methods (10–20% for array-CGH and ~3% for G-banding karyotype analysis). As a result, family-based WES could potentially be considered as a first-tier cost effective diagnostic test for patients with MMS that facilitates better patient management, prognosis and offer accurate recurrence risks to the families.

<figure><img src="image_4.png" /></figure>

## OPEN ACCESS

**Citation:** Kritioti E, Theodosiou A, Parpaite T, Alexandrou A, Nicolaou N, Papaevripidou I, et al. (2021) Unravelling the genetic causes of multiple malformation syndromes: A whole exome sequencing study of the Cypriot population. PLoS ONE 16(7): e0253562. https://doi.org/10.1371/journal.pone.0253562

**Editor:** Dylan A Mordaunt, University of Adelaide, AUSTRALIA

**Received:** March 17, 2021

**Accepted:** June 8, 2021

**Published:** July 29, 2021

**Copyright:** © 2021 Kritioti et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All sequence variants are available from the ClinVar database (accession numbers: SCV001245411 and SCV001443796 - SCV001443812.).

**Funding:** This project was funded by 2017 Telethon Grant (Cyprus), RCB Bank Ltd (Cyprus) and a departmental donation. This work was also supported by the Cy-Tera Project (ΝΑΥΠΔΜΗ/ΣΤATH/0308/31), which was co-funded by the European Regional Development Fund and the

## Introduction

Multiple malformation syndromes (MMS) are a heterogeneous group of genetic disorders that cause developmental anomalies in two or more systems and comprise multiple congenital
Republic of Cyprus through the Cyprus Research Promotion Foundation (now renamed to Cyprus Research and Innovation Foundation). Thibaud Parpaite, Nina Sejourne and Bertrand Coste received funding from the European Research Council (Grant agreement No. 678610).

**Competing interests:** This does not alter our adherence to PLOS ONE policies on sharing data and materials.

malformations [1]. Some examples of MMS are Rett syndrome, Pitt-Hopkins syndrome and Sjögren-Larsson syndrome, which are primarily characterised by intellectual disability (ID), developmental delay (DD) and congenital anomalies involving different systems such as skeletal anomalies, cardiac malformations and other neurological defects. According to the European surveillance of congenital anomalies (EUROCAT) the 2000–2009 birth prevalence of congenital malformations associated with chromosomal or single gene disorders is 4 in 1000 [2].

Despite the low prevalence of single gene disorders, such rare diseases pose a high public-health concern due to their high morbidity and mortality, healthcare costs and other medical challenges. An accurate diagnosis is therefore vital for better management, prognosis and implementation of potential therapies. A molecular diagnosis can also provide information on reproductive options, recurrence risks, offer genetic counselling, prenatal diagnosis and genetic testing to the family members of the affected patients [3].

The first-tier investigations for MMS include detailed clinical evaluation, biochemical analyses, imaging and genetic conventional diagnostic procedures such as chromosome G-banding, fluorescent *in-situ* hybridisation (FISH), array-comparative genomic hybridisation (CGH) and targeted gene testing. Nevertheless, despite the extensive diagnostic workup that most MMS patients undergo, almost half of them remain undiagnosed [4], a fact that may be accounted to a number of factors, including the low diagnostic potential of conventional methods, applied to this group of patients. In fact, several studies regarding patients with multiple congenital anomalies and/or ID have set the diagnostic yield of array-CGH to be at 10–20% [5,6], while for G-banding chromosomal analysis the diagnostic rate of was ~3%.

The past few decades have seen major progress in the development of diagnostic techniques used to characterise rare genetic diseases, especially with the introduction of next generation sequencing (NGS). Through its unbiased sequencing approach, whole exome sequencing (WES) has unveiled a detection rate of approximately 30% in the molecular diagnosis of rare Mendelian disorders [7,8] and has led to the high discovery rate of novel disease-causing genes [9]. Such discoveries facilitate the elucidation of the mode of inheritance that a Mendelian disease follows and enrich our understanding of the biological mechanisms and functions of these genes, which ultimately opens up the avenue for translational research into new diagnostic and therapeutic approaches.

Family-based WES, where patients with a suspected Mendelian disorder are analysed alongside their healthy parents, has emerged as a rapid and cost-effective approach in the identification of the causal variant [8]. Trio-WES previously demonstrated to successfully reduce the number of candidate variants by 10 times [7] and increased the diagnostic yield compared to sequencing the individual proband alone by excluding a large number of non-pathogenic heterozygous rare variants that would otherwise be flagged as potentially disease-causing [8].

Large cohort studies within Europe, such as the Deciphering Developmental Disorders (DDD) study, performing WES on thousands of patients with developmental anomalies have successfully provided a molecular diagnosis to their patients, including single nucleotide variants (SNVs) and copy number variants (CNVs), and have even discovered novel disease-causing genes [7]. To maximize the chance of finding highly penetrant monogenic causes, this study focuses on patients presenting with more than one congenital manifestation and at least one neurodevelopmental anomaly that may be presented as an early onset severe phenotype.

To our knowledge, the present study is the first evaluation of WES as a tool for investigating undiagnosed MMS patients of the Greek-Cypriot population. Novel findings of family-based WES are described and the introduction of WES to the diagnostic practice in Cyprus is outlined as a tier to accelerate the elucidation of the molecular basis of MMS for each family and ending the odyssey of undiagnosed MMS patients.
# Materials and methods

## Patients

A total of 32 unrelated families with 37 affected individuals were recruited in the study (102 total individual tests) between August 2017 and June 2019. Patients were referred by Clinical Geneticists following careful evaluation to fulfil the study's criteria.

The inclusion criteria were: (i) presence of congenital malformations or dysmorphic features and at least one neurodevelopmental anomaly; (ii) a prior negative diagnostic workup (including chromosomal G-banding analysis and array-CGH, at a resolution of at least 200kb and/or targeted gene testing and/or Fragile-X syndrome testing; (iii) negative history of exposure to known causative environmental factors including exposure to teratogens, viral infections and alcohol during pregnancy; (iv) availability of detailed phenotypic description including the patient's family pedigree and any family medical history if applicable; (v) availability of parental samples (vi) signed informed consent by the patient and parents/other family members (if applicable).

This study was approved by the Cyprus National Bioethics Committee under the title "Genetic and Molecular Characterisation of Greek-Cypriot Patients with Multiple Malformation Syndromes (MMS) of Unknown Aetiology" (EEBK/EII2017/19).

## Whole exome sequencing

Genomic DNA (gDNA) was extracted from peripheral blood samples using the QIAmp DNA Blood Midi Kit (Qiagen, Toronto, ON, Canada) and DNA exome library preparation was performed using the TruSeq DNA Exome library preparation kit (Illumina, San Diego, CA, USA) according to the manufacturer's protocol guidelines. The high-output flow cell and the Illumina NextSeq500 platform were used for paired-end (75x2) sequencing of twelve exomes per run.

## Bioinformatics analysis

NGS analysis was performed using an updated version of an in-house developed pipeline [10] according to the Genome Analysis Toolkit (GATK) Best Practices (version 3.3), executed on the Cy-Tera high performance computer (HPC) clusters of the Cyprus Institute (http://web.cytera.cyi.ac.cy/). Raw reads (FASTQ format) were aligned to the human reference genome (GRCh37/hg19) using the Burrows-Wheeler Aligner (BWA-MEM) software [11]. Variants were called using the GATK HaplotypeCaller and GATK Hardfiltering was performed to remove low confidence calls. Variants were annotated for their effect using the standalone script of the Variant Effect Predictor (VEP) tool (Ensembl release 81) (https://asia.ensembl.org/info/docs/tools/vep/) and further annotated and prioritised using GEMINI (Genome MINING, version 0.20.1) [12]. In-house filtering guidelines were implemented to isolate variants with a strong correlation to the phenotype. Variants with potential damaging effects were prioritised using in silico prediction tools, such as CADD (https://cadd.gs.washington.edu/) as well as PolyPhen (version 2.2.2) (http://genetics.bwh.harvard.edu/pph2/) and SIFT (version 5.2.2) (https://sift.bii.a-star.edu.sg/) for missense variants. Prioritisation was performed according to the suspected mode of inheritance of the disease and minor allele frequencies (MAF) of <1% and <5%, from publicly available databases such as gnomAD exomes r2.0.1 (https://gnomad.broadinstitute.org/) and 1000 Genomes (Phase 3) (http://phase3browser.1000genomes.org/index.html). Filtered variants were inspected with an in-house allele frequency dataset specific to the Greek-Cypriot population containing 166 exomes. Databases such as OMIM (https://omim.org/) and ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/)
aided in the identification of known pathogenic variants or established genes contributing to
known clinical phenotypes. Individual variants were further manually assessed using publicly
available databases and tools such as the public version of the Human Gene Mutation Database
(HGMD®) (version 2017.4) (http://www.hgmd.cf.ac.uk/ac/index.php) and Mutation Taster
(http://www.mutationtaster.org/) respectively. Pathogenicity verdict of variants according to
ACMG was assisted with Varsome [13].

*Peddy* tool [14] was used after each run to derive an accurate prediction of the individual's sex and to check relatedness between samples from the same family/trio. The latter feature is essential when discovering *de novo* variants.

To prevent false-positive results originating from the NGS analysis and variant filtering, the integrative genomics viewer (IGV) [15] was used to establish sufficient coverage of the area enclosing the candidate variant, GC content and presence of variants within low complexity regions.

Detection of CNVs from WES using ExomeDepth

CNV analysis was carried out using WES data of patients from our cohort who did not have a
significant finding (i.e. SNV/indel) that could explain their phenotype.

CNV detection was performed using ExomeDepth (v. 1.1.15) [16]. For this analysis an in-house R script was implemented using R version 3.6 on the Cytera cluster of the Cyprus Institute (http://web.cytera.cyi.ac.cy/).

The analysis was performed with standard parameters using bam files created previously for SNV and indel variant calling procedure as described before. For each reference set, 5-10 samples were used depending on the available unrelated samples in the batch. For potential CNV calls on X chromosome only samples from the same sex were correlated.

A typical experiment produced on average 238 CNV calls. Results were filtered using Bayes factor (BF) >50 [17] and each call was further annotated using the online VEP tool (Ensembl release 81) and BioMart (http://grch37.ensembl.org/biomart/martview/).

Confirmation and interpretation of candidate variants

All flagged potentially pathogenic variants were initially evaluated in an interdisciplinary man-
ner between the Geneticists, Bioinformaticians and Clinical Geneticists.

Polymerase chain reaction (PCR) amplification and subsequent bi-directional Sanger sequencing using the ABI 3130xl Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) were used to confirm these variants. All primers were designed using the Primer3 web tool (version 0.4.0) (http://bioinfo.ut.ee/primer3-0.4.0/). A list comprising all primers used for the amplification of all candidate variants within this study is included in S1 Table.

All candidate sequence variants were classified as pathogenic, likely pathogenic, variants of uncertain significance (VUS), likely benign and benign (Table 1), using the standards and guidelines of the American College of Medical Genetics and Genomics (ACMG) [18], assisted by the automatic classification offered by Varsome [13]. All potential causative sequence variants discovered from this study were deposited in ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar) with accession numbers SCV001443796–SCV001443812.

In the case of a confirmed molecular diagnosis the family was informed and was offered
genetic counselling accordingly.

**Additional investigations/functional studies**

**Splice site variant analysis.** For splice site variants with an unclear significance, Human Splicing Finder (http://www.umd.be/HSF/index.html) was used to predict the variants' effect
<table><caption>Table 1. Clinical and molecular findings of molecularly characterised MMS patients.</caption><thead><tr><th>Patient</th><th>Gender</th><th>Age at WES</th><th>Main clinical characteristics</th><th>Gene</th><th>Exon</th><th>Nucleotide change</th><th>Protein change</th><th>Coding impact</th><th>Inheritance</th><th>Zygosity</th><th>Novel/Known</th><th>MAF (gnomAD)</th><th>CADD score</th><th>Variant classification<sup>a</sup></th><th>OMIM Diagnosis</th></tr></thead><tbody><tr><td>1</td><td>F</td><td>8</td><td>Microcephaly, microsomy, dwarfism, mild-moderate DD, mild-moderate DD, FTT</td><td>PCNT</td><td>14</td><td>NN_006031.6: c.2407C&gt;T</td><td>p.Gln803*</td><td>Nonsense</td><td>AR</td><td>Homozygous</td><td>Novel</td><td>N/A</td><td>40</td><td>Pathogenic (PVS1, PM2, PP3, PP4)</td><td>Microcephalic osteodysplastic primordial dwarfism, type II (OMIM #210720)—AR</td></tr><tr><td>2</td><td>M</td><td>4</td><td>GDD, hypotonia, glossophysis, tapering fingers, plagiocephaly, over-folded ear helices, myopia</td><td>UBE3A</td><td>12</td><td>NN_130839.4: c.2410C&gt;T</td><td>p.Gln804*</td><td>Nonsense</td><td>De novo</td><td>Heterozygous</td><td>Novel</td><td>N/A</td><td>45</td><td>Pathogenic (PVS1, PS2, PM2, PP3, PP4)</td><td>Angelman Syndrome (OMIM #105830)—AD</td></tr><tr><td>3</td><td>F</td><td>7</td><td>GDD, FTT, static leukoencephalopathy, microcephaly, pinched nose, bifid first molar, unilateral blepharoptosis</td><td>KAT6A</td><td>16-17</td><td>NN_006766.5: c.3039 +1G&gt;A</td><td>N/A</td><td>Splice donor defect</td><td>De novo</td><td>Heterozygous</td><td>Novel</td><td>N/A</td><td>25.1</td><td>Pathogenic (PVS1, PS2, PM2, PP3, PP4)</td><td>Dystonia, dopa-responsive, due to sepiapterin reductase deficiency (OMIM #612716)—AR</td></tr><tr><td>4</td><td>F</td><td>59</td><td>Moderate-severe ID, cataplexy, obesity, tapering fingers, drowsiness, speech delay</td><td>SPR</td><td>3</td><td>NN_003124.5: c.655C&gt;T</td><td>p.Arg219*</td><td>Nonsense</td><td>AR</td><td>Homozygous</td><td>rs779204655 [19]</td><td>3.182x10<sup>5</sup></td><td>36</td><td>Pathogenic (PVS1, PM2, PP4, PP5)</td><td>Dystonia, dopa-responsive, due to sepiapterin reductase deficiency (OMIM #612716)—AR</td></tr><tr><td>5</td><td>F</td><td>48</td><td>Moderate-severe ID, cataplexy, obesity, tapering fingers, drowsiness, speech delay, rigidity, Parkinsonism</td><td>SPR</td><td>3</td><td>NN_003124.5: c.655C&gt;T</td><td>p.Arg219*</td><td>Nonsense</td><td>AR</td><td>Homozygous</td><td>rs779204655 [19]</td><td>3.182x10<sup>5</sup></td><td>36</td><td>Pathogenic (PVS1, PM2, PP4, PP5)</td><td>Dystonia, dopa-responsive, due to sepiapterin reductase deficiency (OMIM #612716)—AR</td></tr><tr><td>6</td><td>M</td><td>22</td><td>Severe ID, DD, congenital muscular dystrophy, absent speech, sensorineural hearing loss, hyperCKemia, stereotypies, severe visual impairment, cataract, hyperactivity, behavioral abnormalities, macular dystrophy</td><td>POMGNT1</td><td>4</td><td>NN_001243766.1: c.304G&gt;T</td><td>p.Glu102*</td><td>Nonsense</td><td>AR</td><td>Compound heterozygous</td><td>Novel</td><td>N/A</td><td>36</td><td>Pathogenic (PVS1, PM2, PP3, PP4)</td><td>Muscular dystrophy-dystroglycanopathy (congenital with brain and eye anomalies), type A.3 (OMIM #253280)—AR</td></tr><tr><td>7</td><td>M</td><td>18</td><td>Severe ID, DD, congenital muscular dystrophy, absent speech, sensorineural hearing loss, hyperCKemia, stereotypies, severe visual impairment, hydrocephaly, eczenematous rash, squint</td><td>POMGNT1</td><td>4</td><td>NN_001243766.1: c.304G&gt;T</td><td>p.Glu102*</td><td>Nonsense</td><td>AR</td><td>Compound heterozygous</td><td>Novel</td><td>N/A</td><td>36</td><td>Pathogenic (PVS1, PM2, PP3, PP4)</td><td>Muscular dystrophy-dystroglycanopathy (congenital with brain and eye anomalies), type A.3 (OMIM #253280)—AR</td></tr><tr><td>8</td><td>M</td><td>26</td><td>Severe ID, DD, seizures, primary microcephaly, tetraparesis, spasticity with permanent contractures, long fingers, thickened alveolar ridge, thickened gums, strabismus, coliosis, exaggerated reflexes, bilateral foot deformities</td><td>PIEZO2</td><td>24</td><td>NN_022068.3: c.3655G&gt;A</td><td>p.Vall1291Met</td><td>Missense</td><td>AR</td><td>Compound heterozygous</td><td>rs374051556</td><td>4.04x10<sup>5</sup></td><td>22</td><td>VUS (PS3, PS4, PM2, PP4, BP4, BP6)</td><td>Arthrogryposis, distal, with impaired proprioception and touch (OMIM #617146)—AR</td></tr></tbody></table>

(Continued)
<table><caption>Table 1. (Continued)</caption><thead><tr><th>Patient</th><th>Gender</th><th>Age at WES</th><th>Main clinical characteristics</th><th>Gene</th><th>Exon</th><th>Nucleotide change</th><th>Protein change</th><th>Coding impact</th><th>Inheritance</th><th>Zygosity</th><th>Novel/Known</th><th>MAF (gnomAD)</th><th>CADD score</th><th>Variant classification<sup>a</sup></th><th>OMIM Diagnosis</th></tr></thead><tbody><tr><td>9</td><td>M</td><td>26</td><td>Bilateral iris defects, extensive anterior synechiae, posterior embryotoxon, nystagmus, glaucoma, cataracts, bilateral hydronephrosis, vesicoureteric obstruction, bilateral megaureter and bipolar disorder</td><td>PXDN</td><td>16</td><td>NM_012993.3: c.2098G&gt;T</td><td>p.Gly700*</td><td>Nonsense</td><td>AR</td><td>Homozygous</td><td>Novel</td><td>N/A</td><td>41</td><td>Pathogenic (PVS1, PM2, PP3, PP4)</td><td>Anterior segment dysgenesis 7, with sclerocornea (OMIM #269400)—AR</td></tr><tr><td>10</td><td>M</td><td>27</td><td>Severe ID, FTT, microcephaly, spasticity, large prominent teeth, dysplastic ears, sparse eyebrows, high arched palate, short stature, frequent otitis media, dry skin, broad thumbs, sparse teeth, hypoplastic teeth</td><td>KDM6A</td><td>12</td><td>NM_021140.3: c.1192C&gt;T</td><td>p.Gln398*</td><td>Nonsense</td><td>De novo</td><td>Heterozygous</td><td>Novel</td><td>N/A</td><td>37</td><td>Pathogenic (PVS1, PM2, PP3, PP4)</td><td>Kabuli syndrome (OMIM #300867)—AD</td></tr><tr><td>11</td><td>M</td><td>25</td><td>DD, mild myopathy, right orchioepxy, horseshoe kidney, anger issues, shirt terminal phalanges, lactose intolerance, lipomatosis, slanting palpebral fissures, brachydactyly, tapering fingers, mild retrognathia, increased laxity, coordination issues</td><td>PHIP</td><td>31</td><td>NM_017934.7: c.3631_3634delCAAA</td><td>p.Gln1211Aspfs-13</td><td>Frameshift</td><td>De novo</td><td>Heterozygous</td><td>Novel</td><td>N/A</td><td>N/A</td><td>Pathogenic (PVS1, PM2, PP3, PP4)</td><td>Chung Jansen syndrome (OMIM #617991)—AD</td></tr><tr><td>12</td><td>F</td><td>22</td><td>Severe ID, DD, self-injurious behaviour, absent speech, hypotonia, contractures of upper limbs, hand flapping, scoliosis, microstomia, broad nasal tip, squint, irregular menses, facial dysmorphic features</td><td>HECW2</td><td>20</td><td>NM_020760.4: c.3542C&gt;G</td><td>p.Ala1181Gly</td><td>Missense</td><td>De novo</td><td>Heterozygous</td><td>Novel</td><td>N/A</td><td>24.9</td><td>Likely pathogenic (P2, PM2, PP3, PP4)</td><td>Neurodevelopmental disorder with hypotonia, seizures, and absent language (OMIM #617268)—AD</td></tr><tr><td>13</td><td>M</td><td>3</td><td>GDD, ID, camptodactyly, upplanting palpebral fissures, microstomia, coloboma, sparse hair</td><td>HECW2</td><td>20</td><td>NM_020760.4: c.3587A&gt;G</td><td>p.Lys1196Arg</td><td>Missense</td><td>De novo</td><td>Heterozygous</td><td>Novel</td><td>N/A</td><td>27.5</td><td>Pathogenic (P2, PM1, PM2, PP3, PP4)</td><td>Branchiofacial syndrome (OMIM #113620)—AD</td></tr><tr><td>14</td><td>M</td><td>9</td><td>GDD, ID, speech delay, short stature, craniocystosis, strabismus, anal artesia, hypopadias, small chin, 5th finger clinodactyly, dysplastic kidneys, scoliosis</td><td>TFAP2A</td><td>4</td><td>NM_003220.2: c.763A&gt;G</td><td>p.Arg255Gly</td><td>Missense</td><td>De novo</td><td>Heterozygous</td><td>rs121909574 [22–24]</td><td>N/A</td><td>24</td><td>Pathogenic (P2, PM1, PM2, PP2, PP3, PP4, PP5)</td><td>Branchiofacial syndrome (OMIM #113620)—AD</td></tr><tr><td>15</td><td>F</td><td>Deceased</td><td>GDD, ID, FTT, cerebellar atrophy, microcephaly, tremors, chilblains, erythroderma, feeding difficulties, hypotonia, seizures, spinal muscular atrophy, death in childhood</td><td>AGTPBP1</td><td>10</td><td>NM_015239.2: c.820_821delCA</td><td>p.Gln274GlufsTer17</td><td>Missense</td><td>De novo</td><td>Heterozygous</td><td>rs1238165628 [22–24]</td><td>3.99x10<sup>-6</sup></td><td>33</td><td>Likely pathogenic (P2, PM2, PP2, PP3)</td><td>Intellectual developmental disorder with speech delay, autism, and dysmorphic facies (OMIM #618672)—AD</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1.** (Continued)</caption><thead><tr><th>Patient</th><th>Gender</th><th>Age at WES</th><th>Main clinical characteristics</th><th>Gene</th><th>Exon</th><th>Nucleotide change</th><th>Protein change</th><th>Coding impact</th><th>Inheritance</th><th>Zygosity</th><th>Novel/Known</th><th>MAF (gnomAD)</th><th>CADD score</th><th>Variant classification<sup>a</sup></th><th>OMIM Diagnosis</th></tr></thead><tbody><tr><td>16</td><td>F</td><td>44</td><td>Severe ID, absent speech, epileptic encephalopathy, coarse facial features, asymmetric skull, macroglossia, hand biting, stereotypies</td><td>GAMT</td><td>2</td><td>NM_138924.3: c.32/G&gt;A</td><td>p.Lys109 =</td><td>Splice donor defect/ Synonymous</td><td>AR</td><td>Homozygous</td><td>rs80338735 [25--30]</td><td>1.86x10<sup>-4</sup></td><td>14.7</td><td>Pathogenic (PVS1, PM2, PP4, PP5, BP4)</td><td>Cerebral creatine deficiency syndrome 2 (OMIM #612736)—AR</td></tr></tbody></table>

Abbreviations: F—Female; M—Male; ID—Intellectual Disability; DD—Developmental Delay; FTT—Failure to Thrive; GDD—Global Developmental Delay; AR—Autosomal Recessive; AD—Autosomal Dominant.

<sup>a</sup> Variant classification according to the current ACMG guidelines.

https://doi.org/10.1371/journal.pone.0253562.t001
on splicing. In cases with a predicted positive effect, RNA was extracted from Epstein-Barr Virus (EBV) transformed lymphoblastoid cell lines using the RNeasy midi kit (Qiagen) and underwent reverse transcription using the ProtoScript® First Strand cDNA Synthesis Kit (New England Biolabs, Ipswich, MA, USA). PCR amplifications were performed using custom designed primers (S2 Table) and Sanger sequencing was performed using the ABI 3130xl Genetic Analyzer (Applied Biosystems).

**PIEZ02 constructs, site-directed mutagenesis and transformations.** Mutations

Val1219Met and Arg1575Glu were introduced in full-length human *PIEZ02* cDNA by site-
directed mutagenesis using Phusion® High-Fidelity DNA Polymerase (New England Biolabs)
and primers Piezo2-V1219M-pcDNA-F, Piezo2-V1219M-pcDNA-R, Piezo2-R1575Q-
pcDNA-F and Piezo2-R1575Q-pcDNA-R (S3 Table) according to a modified version of the
QuickChange II XL Site-Directed Mutagenesis Kit (Agilent) (conditions are available upon
request). The primer pair for the mutagenesis of Arg1575Glu was designed manually as to con-
tain the altered base mid-sequence and the primer pair for Val1219Met was designed using the
QuickChange® Primer Design Program (Agilent) [31]. The mutated constructs were trans-
formed in DH5 alpha turbo competent *E. coli* cells. Plasmid cDNA was purified using the Qia-
gen Plasmid Maxi kit according to the manufacturer's protocol. Sanger sequencing was used
to confirm the successful mutagenesis of both variants followed by sequential digestion using
*Apal* and *NotI*.

**Cell culture and transient transfection.** PIEZO1-deficient HEK293T cells (Human embryonic kidney) [32] were grown in Dulbecco's Modified Eagle Medium containing 4.5 mg.mL⁻¹ glucose, 10% (vol/vol) FBS, 100 U.mL⁻¹ (1%) penicillin/streptomycin. Cells were plated onto poly-D-lysine-coated 12-mm round glass coverslips (Corning) in 24-well plates and transfected using lipofectamine 2000 (Invitrogen) according to the manufacturer's instruction. Each well was transfected with 600 ng.mL⁻¹ of specified constructs. For double mutant experiments, cells were co-transfected with WT + Arg1575Glu constructs or Val1219-Met + Arg1575Glu constructs at a ratio of 1:2 (300 ng.mL⁻¹ + 600 ng.mL⁻¹ plasmid per condition). All transfections were supplemented with 300 ng.mL⁻¹ of GFP plasmid to visualise transfected cells and GFP positive cells were recorded 24–48 hours later.

**Electrophysiology.** Patch-clamp experiments were performed in standard whole-cell recordings using an Axopatch 200B amplifier (Axon Instruments). Patch pipettes had a resistance of 1.5–2.5 MΩ when filled with an internal solution consisting of (in mM) 140 CsCl, 10 Hepes, 5 EGTA, 1 CaCl₂, 1 MgCl₂, 4 MgATP, and 0.4 Na₂GTP (pH adjusted to 7.3 with CsOH). The extracellular solution consisted of (in mM) 133 NaCl, 3 KCl, 1 MgCl₂, 10 Hepes, 2.5 CaCl₂, 10 glucose (pH adjusted to 7.3 with NaOH). All experiments were done at room temperature. Currents were sampled at 20 kHz and filtered at 2 kHz. Leak currents before mechanical stimulations were subtracted off-line from the current traces. Clampfit 10.7 (Molecular Devices) software was used to analyse recordings and biophysical parameters. PRISM 7.0 (GraphPad) software was used for statistical analysis (non-parametric Krustal-Wallis test). All chemicals were sourced from Sigma-Aldrich.

**Mechanical stimulation.** Mechanical stimulation was achieved using a firepolished glass pipette (tip diameter 3–4 µm) positioned at an angle of 80° to the cell being recorded. Downward movement of the probe toward the cell was driven by a Clampex controlled piezo-electric crystal microstage (E625 LVPZT Controller/Amplifier, Physik Instrumente). The probe had a velocity of 0.7 µm.ms⁻¹ during the ramp segment of the command for forward motion, and the stimulus was applied for 150 ms. A series of mechanical steps in 0.5 µm increments was applied every 10 s. Inward MA currents were recorded at a holding potential of -80 mV.
# Results

## Patients

Thirty-two unrelated families with 37 patients met the inclusion criteria of this study as outlined in the Methods section. The patient group included 16 children (43.2%) and 21 adults (56.7%). The mean age of patients at the time of referral for WES was 21 (ranging from 3.5 to 59), with a nearly equal male to female ratio (18:19). The vast majority of the individuals included within this cohort are from the European population and more specifically from Greek-Cypriot origin. The observed ancestry composition of our cohort was estimated with a randomized Principal Component Analysis (PCA) as shown in Fig 1.

The most common clinical characteristics within our cohort were ID and DD. A wide spectrum of congenital anomalies, dysmorphic features and other neurological deficits were presented including short stature, microcephaly and seizures respectively.

Three (9.4%) of the recruited families were self-reported consanguineous. All parents included in the cohort were apparently healthy.

## Whole exome sequencing

**Coverage/Depth.** Based on the 100 samples that underwent WES within this cohort, the mean target coverage was 59X, with 89% of the target bases achieving at least 10X coverage.

## Molecular results

**Detection rate.** Family-based WES and subsequent analysis performed on 32 families revealed 17 potentially pathogenic SNVs or small indels in 16 MMS patients, achieving an overall detection rate of 43.2% (16/37). All molecular diagnoses and clinical information are summarised in Table 1. The clinical information of the recruited MMS patients without a molecular diagnosis is available upon request.

<figure><img src="image_5.png" /><figcaption>**Fig 1.** Principal Component Analysis (PCA) using *peddy* tool, of samples (n = 166) from in-house dataset (in purple) projected to those from 1000 genomes. The samples of this study are included in the dataset.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0253562.g001
<table><caption>Table 2. Distribution of Mendelian mode of inheritance, novelty of variants and variant type of index cases with a molecular diagnosis.</caption><thead><tr><th></th><th>Number of diagnosed patients (percentage %)</th></tr></thead><tbody><tr><td><strong>Mode of inheritance</strong></td><td></td></tr><tr><td><strong>Dominant</strong></td><td>7 (43.8%)</td></tr><tr><td>Autosomal dominant</td><td>6 (37.5%)</td></tr><tr><td>X-linked dominant</td><td>1 (6.3%)</td></tr><tr><td><strong>Recessive</strong></td><td>9 (56.3%)</td></tr><tr><td>Homozygous</td><td>6 (37.5%)</td></tr><tr><td>Compound heterozygous</td><td>3 (18.8%)</td></tr><tr><td>X-linked recessive</td><td>0</td></tr><tr><td></td><td><strong>Number of candidate pathogenic variants (percentage %)</strong></td></tr><tr><td><strong>Variant novelty</strong></td><td></td></tr><tr><td>Known</td><td>7 (41.2%)</td></tr><tr><td>Novel</td><td>10 (58.8%)</td></tr><tr><td><strong>Type of mutation</strong></td><td></td></tr><tr><td>Nonsense</td><td>6 (35.3%)</td></tr><tr><td>Missense</td><td>7 (41.2%)</td></tr><tr><td>Frameshift</td><td>2 (11.8%)</td></tr><tr><td>Splice defect</td><td>2 (11.8%)</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0253562.t002

## Candidate variants

All 17 candidate variants had a strong correlation with the patient's phenotype (Table 1).

Based on the distribution of the confirmed candidate variants, missense (7/17, 41.1%) were most abundant, followed closely by nonsense (6/17, 35.3%) variants. For the detection of potentially causative variants, *in silico* tools such as PolyPhen2 and SIFT predicted pathogenicity for all 7 flagged missense mutations and most mutations had a CADD score of more than 20, suggesting deleteriousness (Table 1). Based on the ACMG guidelines for the classification of sequence variants, the majority of our findings were classified as pathogenic (11/17, 64.7%). Based on our in-house dataset, the flagged variants detected in the 16 MMS patients were not identified in any other individuals referred to the lab. An autosomal recessive mode of inheritance was the most common between the diagnosed MMS patients (9/16, 56.3%). All candidate variants were identified in genes known to cause Mendelian monogenic disorders with the majority of the variants being described as novel (10/17, 58.8%), meaning they have not been reported before in the literature or in publicly available databases. An analytical summary of these statistics is illustrated in Table 2.

All 7 autosomal dominant conditions resulted from *de novo* mutations, whereas autosomal recessive conditions resulted from homozygous and compound heterozygous mutations. All patients with *de novo* disease-causative variants were born to apparently healthy, non-consanguineous parents.

Patient 12 had 2 heterozygous *de novo* variants in gene *HECW2* causing neurodevelopmental disorder with hypotonia, seizures, and absent language (OMIM #617268). The variants were not present in neither of the parents upon Sanger confirmation and when observing our NGS data the two variants were located on the same reads (S2 Fig).

Of the recessive diagnosed disorders, variant p.Arg219* observed in *SPR* (patients 4 and 5) has been previously identified in patients with similar phenotypes [19]. Although some of the other variants including compound heterozygous variants p.Glu102* and p.Arg129Trp found
in *POMGNT1* (patients 6 and 7) and p.Val1219Met and p.Arg1575Gln identified in *PIEZO2* (patient 8) are present in databases such as ClinVar and gnomAD (Table 1), they have not been phenotypically described in the literature.

Of the conditions that were caused by monogenic homozygous mutations, one was found in *PCNT* in patient 1, from a quartet (two affected siblings and both parents) with consanguineous parents. Another was identified in *PXDN* (patient 9) within a trio with parents who were third cousins and the other four cases were found in *SPR* (patients 4 and 5), *AGTPBP1* (patient 15) and *GAMT* (patient 16), as mentioned above, within different trios with apparently unrelated parents.

All of the above variants following a recessive mode of inheritance demonstrated clear segregation from apparently unaffected carrier parents to their affected offspring.

Collectively, 15 genes (*PCNT*, *UBE3A*, *KAT6A*, *SPR*, *POMGNT1*, *PIEZO2*, *PXDN*, *KDM6A*, *PHIP*, *HECW2*, *TFAP2A*, *CNOT3*, *AGTPBP1* and *GAMT*) with presumed pathogenic SNVs were identified within this cohort.

## CNV findings from WES data

Upon CNV analysis of all patients with a negative result from family-based WES, no significant CNVs were identified that could explain the patients' clinical presentations.

## Further experimental investigations

A predicted pathogenic splice site variant c.(3039+1G>A) in *KAT6A* was identified in patient 3. Sanger sequencing of cDNA derived from patient lymphoblastoid cell lines revealed that the functional effect of the splice variant resulted in a 19bp intron retention at the donor site of exon 16 (Fig 2), which is predicted to result in the protein's truncation.

For patient 8, two missense variants were identified in a compound heterozygous form in gene *PIEZO2*. *PIEZO2* homo-trimers form mechanically activated (MA) ion channels involved in touch and proprioception [33–36], whose mutations, either gain- or loss-of-function, have been linked to MMS with overlapping clinical phenotypes [37]. We generated human *PIEZO2* constructs with p.(Val1219Met) maternally inherited mutation and p.(Arg1575Gln) paternally inherited mutation to test functional consequences of these mutations on PIEZO2 channel activity. Heterologous expression of these mutants was done in HEK-P1 KO cells to avoid contamination by endogenous MA channels and compared to WT channels (Fig 3). Expression of Val1219Met mutant induces MA currents similar to WT PIEZO2 regarding current amplitude, inactivation kinetics and activation threshold. Interestingly, expression of Arg1575Gln mutant does not induce any detectable MA channel activity (Fig 3B and 3D).

Since Val1219Met mutation seems to have no effect on properties of recorded MA currents, we tested if its association with Arg1575Gln mutant could induce MA channels with altered functional properties. However, co-expression of Val1219Met and Arg1575Gln mutants induced MA currents with similar characteristic than co-expression of WT and Arg1575Gln mutant (Fig 3C and 3E), arguing that Val1219Met associated to Arg1575Gln mutation do not induce MA channels with altered biophysical properties. Altogether, these experiments show that heterologous expression of Val1219Met mutant induces MA currents similar to WT currents, and that Arg1575Gln mutation is deleterious and induces loss-of-function.

## Discussion

Here we present the first genomic investigation performed on MMS patients from the Greek-Cypriot population. This study highlights the high detection rate of 43.2% obtained by WES
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<table><thead><tr><th></th><th colspan="4">Exon 16</th><th colspan="12">Exon 17</th></tr></thead><tbody><tr><th>Wild type</th><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>G</td><td>A</td><td>A</td><td>A</td><td>C</td><td>C</td><td>A</td><td>T</td><td>T</td><td>T</td><td>C</td><td>T</td><td>C</td><td>C</td><td>A</td><td>C</td><td>C</td><td>G</td><td>A</td><td>A</td><td>G</td><td>G</td><td>A</td><td>G</td><td>A</td><td>T</td><td>T</td></tr><tr><th>Mutation</th><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>G</td><td>A</td><td>T</td><td>A</td><td>A</td><td>C</td><td>G</td><td>C</td><td>A</td><td>T</td><td>T</td><td>T</td><td>C</td><td>T</td><td>T</td><td>T</td><td>C</td><td>C</td><td>A</td><td>G</td><td>A</td><td>A</td><td>A</td><td>C</td><td>C</td><td>A</td><td>G</td><td>A</td></tr></tbody></table>

**Fig 2.** Analysis of the c.3039+1G>A splice variant. A—IGV image illustrating presence of de novo splice variant c.(3039+1G>A) in *KAT6A*, identified in patient 3 and absent in parents. B—Electropherograms showing presence of c.(3039+1G>A) in DNA of proband (G/A) and absence in both parents (G/G) C—RT-PCR of *KAT6A* and actin (ACTB) cDNA transcripts of patient 3 (P) and a healthy control (C). D—The functional effect of c.(3039+1G>A) splice variant in patient 3, which was identified by sequencing the cDNA product from RNA extracted from patient's lymphoblastoid cell lines. The splice donor variant resulted in the intron retention of 19bp and a frameshift change.

https://doi.org/10.1371/journal.pone.0253562.g002
<figure><img src="image_11.png" /><figcaption>**Fig 3.** A—Schematic topology of PIEZO2 protein with position of the two identified mutation. B, C—Typical example of mechanically-activated currents recorded in HEK293T-P1KO cells transiently transfected with each specified constructs (600ng expression plasmid per condition) (B) or co-transfected with WT+R1575Q or V1219M+R1575Q at a ratio of 1:2 (300ng+600ng plasmid per condition) (C). Holding potential is -80mV. Top traces represent the displacement of the stimulation probe, bottom traces represent the recorded currents. Traces corresponding to mechanical stimulations incremented by 0.5 µm are superimposed. D, E—Average maximal current amplitude (left panel), inactivation time-constant (middle panel) and threshold for mechanical current activation (right panel). Cells were recorded 24 to 48 hours post transfection transfection (n = 38, 28 and 14 for WT, V1219M and R1575Q, respectively; n = 16 and 22 for WT +R1575Q and V1219M+R1575Q, respectively; mean±s.e.m.; Mann-Whitney test).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0253562.g003

for a challenging group of undiagnosed patients, identifying a total of 17 potentially causative variants.

A range of inheritance patterns, types of mutations and a high number of novel findings were identified (Table 2). Fifty-three percent of the total pathogenic variants involved autosomal recessive disorders, which demonstrates the high potential of family-based WES in diagnosing recessive diseases in patients who have already undergone an extensive diagnostic workup. The number of detected autosomal recessive disorders is significantly higher compared to other similar studies accounting for 27% and 34.3% in patient groups of 416 and 2000 respectively [38,39]. This increase in recessive diseases detected within the Cypriot population may be accounted to the small population size of Cyprus as well as the small sample size of this study. Due to the high prevalence of founder effect mutations observed in regional studies, especially those with small population sizes, the study's group will resume the study in order to identify the true frequency of potential founder effect mutations within the Cypriot population. Specifically, a high frequency (>0.05) of a founder effect mutation related to β-
thalassemia has already been extensively demonstrated within the Cypriot population in the past [40]. Concurrently, the lower number of *de novo* variants detected may account for their novelty and for the fact that they may be present in genes that have not been linked to a disease before, stressing the importance of depositing novel molecular findings to freely-accessible databases. Doing so will subsequently increase knowledge within the scientific community and accelerate diagnosis for other MMS patients. We also observed that missense mutations were abundant in our cohort (43.2%) as it has also been observed in a recent study on patients with neurodevelopmental disorders and multiple congenital anomalies (40.8%) [41]. Even though the frequencies of autosomal recessive disorders and missense variants detected in this study are higher or similar to other larger studies, a more reliable frequency can be achieved with the continuity of this study and subsequent increase of the study's sample size.

Literature regarding rare genetic syndromes is limited, therefore even in the cases where already known disease-causing variants are identified, they offer valuable additions to limited pre-existing knowledge. In our study, the identification of already known disease-causing variants has helped expand the spectrum of phenotypes. More specifically, of the 7 known mutations identified within our cohort only three have been characterised both molecularly and phenotypically in the past and are present in OMIM. These mutations were detected in genes *SPR*, *TFAP2A* and *GAMT* [19,22,23,25-30]. The current study contributed to the phenotypic spectrum of the disorder caused by mutations found within *SPR*, especially since the affected siblings with a variant p.(Arg219*) in *SPR* are significantly older than the majority of previously reported cases. Obesity and a Prader-Willi-like manifestation are present in both of our *SPR* patients, both features being novel to patients with sepiapterin reductase deficiency (SRD). Such findings can ultimately offer a new and faster avenue for the clinical diagnosis of younger SRD patients.

Novel phenotypes were also identified in patient 15 who is a carrier of a homozygous *AGTPBP1* variant p.(Gln274GlufsTer17), which is a novel variant within a known disease-associated gene. The patient's clinical manifestations included novel phenotypes such as erythroderma and an abnormal metabolic profile, which have not been reported before in patients with recessive mutations in *AGTPBP1*.

Known variants with a strong genotype-phenotype correlation within well-established disease-causing genes were easily identified as the primary flagged variants. Such findings include two compound heterozygous variants in *POMGNT1* in two male siblings (patients 6 and 7), with typical MRI findings of patients with *POMGNT1* mutations (S1 Fig), a homozygous variant in *PCNT* in patient 1 who was subsequently diagnosed with microcephalic osteodysplastic primordial dwarfism, type II (MOPDII) and patient 2 with a variant in *UBE3A* causative for Angelman syndrome. Although the latter two are novel variants, other flagged variants such as ones identified within *SPR* and *TFAP2A* have already been characterised in the literature before, making the diagnoses of these patients easier and marking the end of the diagnostic odyssey of these families much faster. For variants found within genes that have recently been linked to a genetic disorder and for variants with uncertain significance, such as ones located within *KAT6A* and *PIEZO2*, further work was carried out to assess their pathogenicity.

Even though RNA sequencing can demonstrate the direct effect of a mutation on the RNA or even protein level, WES provided insight to investigate the effect of splice variants at an RNA-level without the need to perform RNA sequencing. This was more specifically observed in patient 3, with a splice donor variant in *KAT6A* (c.3039+1G>A) that led to an intron retention of 19bp.

*PIEZO2* mutations in humans lead to disease of variable clinical phenotypes with overlapping clinical features including arthrogryposis syndrome [37]. The two previously reported missense variants identified in *PIEZO2* from patient 8 were tested functionally. The paternally
inherited variant p.(Arg1575Gln) was functionally inactive, as its heterologous expression induced no detectable MA channel activity. Arg1575Gln is located within a sequence connecting Beam and Latch domains of the protein whose alteration could affect gating of Piezo2 channels [42,43] or could induce mis-fold of the protein leading to its degradation. No functional difference was identified in MA currents induced by the maternally inherited variant p. (Val1219Met) whose mutation lies in the TM21-TM22 extracellular loop [43]. Since PIEZO channels are formed by homo-trimeric assembly, co-expression experiments were done with deleterious p.(Arg1575Gln) in excess. Co-expression of p.(Arg1575Gln) with WT induces MA currents similar to WT currents suggesting no dominant negative effect. Co-expression of both mutants induces MA currents similar to WT or p.(Val1219Met) currents, excluding an interaction-dependent mechanism between mutants. It cannot be ruled out that the p. (Val1219Met) mutation may have effects on channel activity that are not detected due to sensitivity limitations of our functional assay. Another hypothesis is that p.(Val1219Met) mutation could alter PIEZO properties that are not related to its ion channel activity or that are dependent from the cellular context, such as interaction with partners absent in heterologous expression system. Nevertheless, patient 8 clinical phenotype is reminiscent of the various clinical phenotypes associated with biallelic loss-of-function or dominant gain-of-function PIEZO2 mutations [44–46] suggesting that biallelic p.(Val1219Met) and p.(Arg1575Gln) mutations alter PIEZO2 function leading to MMS, which could therefore be causative for the patient’s phenotype.

It has been shown that approximately 2% of all *de novo* SNVs are multinucleotide variants (MNVs) occurring 20bp-5kb from one another [47,48] and constitute an important, sometimes undetected cause of genetic disorders. Based on a recent study on 6688 trio-WES cases, *de novo* MNVs were more frequently observed within genes associated with developmental disorders and were suggested to be more damaging than SNVs, even when the MNVs were both missense mutations [49]. This appears to be the case of patient 12, where two missense variants 45bp apart were identified in *HECW2*. Both variants were classified as *likely pathogenic* and their combination appears to be deleterious. This class of SNV variants is very challenging and special attention is required for their interpretation and detection.

Since the use of WES can give different detection rates for different patient groups, this study was also able to highlight the importance of patient selection and setting strict criteria in achieving a high detection rate (43.2%). Receiving a detailed phenotypic description from the clinicians and having stricter patient inclusion criteria can ultimately lead to higher detection rates being achieved using WES.

Despite a high detection rate, the limitations of WES as a diagnostic tool are associated to its incomplete coverage of exonic regions, insufficient sequencing of areas with high GC content, low complexity regions, stretches of repeats, inability in detecting mosaicism in the parents due to the limited depth of coverage of the technique and the challenges in detecting CNVs.

CNVs not detected by array-CGH or WES are still left undiscovered, which could account for a large number of cases that are left undiagnosed. Recent advances in CNV detection using WES data have proven to be beneficial and have increased the overall diagnostic rate of various studies [50,51]. Based on the literature, the diagnostic rates of CNV detection on patients using WES data vary widely among different tools [16]. In this cohort, all patients had already undergone an extensive diagnostic workup until the time of CNV analysis, leaving little room for the detection of pathogenic CNVs.

In order to increase the rate of accurate genetic diagnosis for MMS patients, future advances in higher read depths, the use of whole genome sequencing (WGS), long-read sequencing and reanalysis of WES data could potentially make it possible to accurately identify pathogenic SNVs and CNVs in a single experiment.
**Conclusions**

This is the first genomic study on MMS patients from the Greek-Cypriot population. WES provided a potential genetic diagnosis for 43.2% of a cohort of MMS patients, with a previous negative diagnostic workup, and included the discovery of rare variants, most of which are novel. The high diagnostic yield and efficacy of WES observed in this study suggest that WES could even be considered as a first-tier cost effective diagnostic test over conventional diagnostic approaches for MMS patients when a specific disorder is not suspected. The high detection rate indicates that WES is a valuable tool for the efficient diagnosis of patients with MMS and highlights the importance of selecting patients that have more than one congenital anomaly and/or neurodevelopmental disorder or an early onset severe phenotype. Lastly, we propose that the genes in which flagged variants were identified in this study could be added in clinical exome sequencing (CES) and related panels for MMS. For the undiagnosed cases, WES data reanalysis and WGS will be carried out in the near future.

**Supporting information**

S1 Fig. Example of strong genotype-phenotype correlation illustrated through MRI findings of patients 6 and 7 showing hydrocephalus and vermian hypoplasia; typical characteristics of patients with mutations in *POMGNT1*.

(TIF)

S2 Fig. IGV image illustrating presence of *de novo* variants p.(Ala1181Gly) and p. (Lys1196Arg) in *HECW2*, in patient 13, on same reads.

(TIF)

S1 Table. List of primer pairs used in Sanger sequencing for validation of flagged variants. The synthesis scale and shipping condition of all primers were 0.02 µmol and diss. 100 µM respectively. [F: forward; R: reverse]

(DOCX)

S2 Table. Primers used for cDNA sequencing for confirmation of functional effect of splice change in patient 3.

(DOCX)

S3 Table. Primers used for site-directed mutagenesis in human *PIEZ02* pcDNA3.1.

(DOCX)

S1 Raw images.

(PDF)

**Acknowledgments**

The authors would like to thank all the families for their participation in the study. Authors thank Virginie Penalba, Panayiotis Myrianthopoulos, Dr Pavlos Fanis and Dr Anthi Demetria-dou for technical assistance and Dr Ludmila Kousoulidou for her help with this manuscript.

**Author Contributions**

**Conceptualization:** George A. Tanteles, Carolina Sismani.

**Formal analysis:** Evie Kritioti, Athina Theodosiou, Thibaud Parpaite.

**Investigation:** Evie Kritioti, Thibaud Parpaite, Angelos Alexandrou, Nayia Nicolaou, Ioannis Papaevripidou, Nina Séjourné.
**Methodology:** Evie Kritioti, Thibaud Parpaite, Angelos Alexandrou, Ioannis Papaevripidou, Bertrand Coste.

**Project administration:** George A. Tanteles, Carolina Sismani.

**Resources:** Violetta Christophidou-Anastasiadou, George A. Tanteles, Carolina Sismani.

**Software:** Athina Theodosiou.

**Supervision:** George A. Tanteles, Carolina Sismani.

**Writing – original draft:** Evie Kritioti, Thibaud Parpaite, Bertrand Coste.

**Writing – review & editing:** Evie Kritioti, Athina Theodosiou, George A. Tanteles, Carolina Sismani.

## References

1. Segen JC. Concise Dictionary of Modern Medicine [Internet]. 2nd ed. McGraw-Hill; 2006 [cited 2020 Mar 26]. https://medical-dictionary.thefreedictionary.com/multiple+malformation+syndromes.

2. Moorthie S, Blencowe H, Darlison MW, Lawn J, Morris JK, Modell B, et al. Estimating the birth prevalence and pregnancy outcomes of congenital malformations worldwide. J Community Genet. 2018 Oct 1; 9(4):387-96. https://doi.org/10.1007/s12687-018-0384-2 PMID: 30218347

3. Schieppati A, Henter JI, Daina E, Aperia A. Why rare diseases are an important medical and social issue. Vol. 371, The Lancet. 2008. p. 2039-41. https://doi.org/10.1016/S0140-6736(08)60872-7 PMID: 18555915

4. Shashi V, McConkie-Rosell A, Rosell B, Schoch K, Vellore K, McDonald M, et al. The utility of the traditional medical genetics diagnostic evaluation in the context of next-generation sequencing for undiagnosed genetic disorders. Genet Med. 2014; 16(2):176-82. https://doi.org/10.1038/gim.2013.99 PMID: 23928913

5. Sagoo GS, Butterworth AS, Sanderson S, Shaw-Smith C, Higgins JPT, Burton H. Array CGH in patients with learning disability (mental retardation) and congenital anomalies: updated systematic review and meta-analysis of 19 studies and 13,926 subjects. Genet Med. 2009 Mar; 11(3):139-46. https://doi.org/10.1097/GIM.0b013e318194ee8f PMID: 19367186

6. Vianna GS, Medeiros PFV, Alves AF, Silva TO, Jehee FS. Array-CGH analysis in patients with intellectual disability and/or congenital malformations in Brazil. 2016.

7. Wright CF, Fitzgerald TW, Jones WD, Clayton S, McRae JF, van Kogelenberg M, et al. Genetic diagnosis of developmental disorders in the DDD study: a scalable analysis of genome-wide research data. Lancet. 2015; 385(9975):1305-14. https://doi.org/10.1016/S0140-6736(14)61705-0 PMID: 25529582

8. Farwell KD, Shahmizadi L, El-Khechen D, Powis Z, Chao EC, Tippin Davis B, et al. Enhanced utility of family-centered diagnostic exome sequencing with inheritance model-based analysis: results from 500 unselected families with undiagnosed genetic conditions. Genet Med. 2015 Jul 13; 17(7):578-86. https://doi.org/10.1038/gim.2014.154 PMID: 25356970

9. Fitzgerald TW, Gerety SS, Jones WD, van Kogelenberg M, King DA, McRae J, et al. Large-scale discovery of novel genetic causes of developmental disorders. Nature. 2014 Dec 24; 519(7542):223-8. https://doi.org/10.1038/nature14135 PMID: 25533962

10. Theodosiou A. Bioinformatics variant calling analysis pipeline (version 1.0) guidelines for Next generation sequencing (exome) of human samples using Illumina NextSeq500 platform. NGS analysis pipeline guidelines. 2017 [cited 2019 May 11]; https://training.vi-seem.eu/images/trainingMaterial/LifeSciences/NGS_analysis_pipeline_Theodosiou_CYL.pdf.

11. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009 Jul 15; 25(14):1754-60. https://doi.org/10.1093/bioinformatics/btp324 PMID: 19451168

12. Paila U, Chapman BA, Kirchner R, Quinlan AR. GEMINI: Integrative Exploration of Genetic Variation and Genome Annotations. Gardner PP, editor. PLoS Comput Biol. 2013 Jul 18; 9(7). https://doi.org/10.1371/journal.pcbi.1003153 PMID: 23874191

13. Kopanos C, Tsiolkas V, Kouris A, Chapple CE, Aguilera MA, Meyer R, et al. VarSome: the human genomic variant search engine. Bioinformatics. 2019; 35(11):1978-80. https://doi.org/10.1093/bioinformatics/bty897 PMID: 30376034
14. Pedersen BS, Quinlan AR. Who's Who? Detecting and Resolving Sample Anomalies in Human DNA Sequencing Studies with Peddy. Am J Hum Genet. 2017 Mar 2; 100(3):406–13. https://doi.org/10.1016/j.ajhg.2017.01.017 PMID: 28190455

15. Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013 Mar 1; 14(2):178–92. https://doi.org/10.1093/bib/bbs017 PMID: 22517427

16. Plagnol V, Curtis J, Epstein M, Mok KY, Stebbings E, Grigoriadou S, et al. A robust model for read count data in exome sequencing experiments and implications for copy number variant calling. Bioinformatics. 2012 Nov 1; 28(21):2747–54. https://doi.org/10.1093/bioinformatics/bts526 PMID: 22942019

17. Zanardo ÉA, Monteiro FP, Chehimi SN, Oliveira YG, Dias AT, Costa LA, et al. Application of Whole-Exome Sequencing in Detecting Copy Number Variants in Patients with Developmental Delay and/or Multiple Congenital Malformations. J Mol Diagnostics. 2020 Aug 1; 22(8):1041–9. https://doi.org/10.1016/j.jmoldx.2020.05.007 PMID: 32497716

18. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: A joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015 May 8; 17(5):405–24. https://doi.org/10.1038/gim.2015.30 PMID: 25741868

19. Dill P, Wagner M, Somerville A, Thöny B, Blau N, Weber P. Child Neurology: Paroxysmal stiffening, upward gaze, and hypotonia Hallmarks of sepiapterin reductase deficiency. Neurology. 2012; 78(4):e29–32.

20. Sframeli M, Sarkozy A, Bertoli M, Astrea G, Hudson J, Scoto M, et al. Congenital muscular dystrophies in the UK population: Clinical and molecular spectrum of a large cohort diagnosed over a 12-year period. Neuromuscul Disord. 2017; 27(9):793–803. https://doi.org/10.1016/j.nmd.2017.06.008 PMID: 28688748

21. Borisovna KO, Yurievna KA, Yurievich TK, Igorevna KO, Olegovich KD, Igorevna DA, et al. Compound heterozygous POMGNT1 mutations leading to muscular dystrophy-dystroglycanopathy type A3: A case report. BMC Pediatr. 2019 Apr 8; 19(1):98. https://doi.org/10.1186/s12887-019-1470-2 PMID: 30961548

22. Milunsky JM, Maher TA, Zhao G, Roberts AE, Stalker HJ, Zori RT, et al. TFAP2A Mutations Result in Branchio-Oculo-Facial Syndrome. Am J Hum Genet. 2008 May 9; 82(5):1171–7. https://doi.org/10.1016/j.ajhg.2008.03.005 PMID: 18423521

23. Reiber J, Sznajer Y, Posteguillo EG, Müller D, Lyonnet S, Baumann C, et al. Additional clinical and molecular analyses of TFAP2A in patients with the branchio-oculo-facial syndrome. Am J Med Genet Part A. 2010 Apr; 152(4):994–9.

24. Meshcheryakova TI, Zinchenko RA, Vasilyeva TA, Marakhonov AV., Zhylina SS, Petrova NV., et al. A clinical and molecular analysis of branchio-oculo-facial syndrome patients in russia revealed new mutations in TFAP2A. Ann Hum Genet. 2015 Mar 1; 79(2):148–52. https://doi.org/10.1111/ahg.12098 PMID: 25590586

25. Stöckler S, Isbrandt D, Hanefeld F, Schmidt B, Von Figura K. Guanidinoacetate methyltransferase deficiency: The first inborn error of creatine metabolism in man. Am J Hum Genet. 1996; 58(5):914–22. PMID: 8651275

26. Mercimek-Mahmutoglu S, Stoeckler-Ipsiroglu S, Adami A, Appleton R, Araújo HC, Duran M, et al. GAMT deficiency: Features, treatment, and outcome in an inborn error of creatine synthesis. Neurology. 2006 Aug 8; 67(3):480–4. https://doi.org/10.1212/01.wnl.0000234852.43688.bf PMID: 16855203

27. Morris AAM, Appleton RE, Power B, Isherwood DM, Abernethy LJ, Taylor RW, et al. Guanidinoacetate methyltransferase deficiency masquerading as a mitochondrial encephalopathy. J Inherit Metab Dis. 2007 Feb 14; 30(1):100–100. https://doi.org/10.1007/s10545-006-0478-2 PMID: 17171576

28. Dhar SU, Scaglia F, Li FY, Smith L, Barshop BA, Eng CM, et al. Expanded clinical and molecular spectrum of guanidinoacetate methyltransferase (GAMT) deficiency. Mol Genet Metab. 2009 Jan 1; 96(1):38–43. https://doi.org/10.1016/j.ymgme.2008.10.008 PMID: 19027335

29. Comeaux MS, Wang J, Wang G, Kleppe S, Zhang VW, Schmitt ES, et al. Biochemical, molecular, and clinical diagnoses of patients with cerebral creatine deficiency syndromes. Mol Genet Metab. 2013; 109(3):260–8. https://doi.org/10.1016/j.ymgme.2013.04.006 PMID: 23660394

30. Mercimek-Mahmutoglu S, Ndika J, Kanhai W, de Villemeur TB, Cheillan D, Christensen E, et al. Thirteen New Patients with Guanidinoacetate Methyltransferase Deficiency and Functional Characterization of Nineteen Novel Missense Variants in the GAMT Gene. Hum Mutat. 2014 Apr 1; 35(4):462–9. https://doi.org/10.1002/humu.22511 PMID: 24415674

31. Novoradovsky A, Zhang V, Ghosh M, Hogrefe H, Sorge JA, Gaasterland T. Computational Principles of Primer Design for Site Directed Mutagenesis. In: Technical Proceedings of 2005 NSTI Nanotechnology Conference and Trade Show, Anaheim. 2005. p. 532–5.
32. Lukacs V, Mathur J, Mao R, Bayrak-Toydemir P, Procter M, Cahalan SM, et al. Impaired PIEZO1 function in patients with a novel autosomal recessive congenital lymphatic dysplasia. Nat Commun. 2015 Sep 21; 6(1):1-7. https://doi.org/10.1038/ncomms9329 PMID: 26387913

33. Coste B, Mathur J, Schmidt M, Earley TJ, Ranade S, Petrus MJ, et al. Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels. Science. 2010 Oct 1; 330(6000):55-60. https://doi.org/10.1126/science.1193270 PMID: 20813920

34. Ranade SS, Woo SH, Dubin AE, Moshourab RA, Wetzel C, Petrus M, et al. Piezo2 is the major transducer of mechanical forces for touch sensation in mice. Nature. 2014 Dec 4; 516(729):121-5. https://doi.org/10.1038/nature13980 PMID: 25471886

35. Florez-Paz D, Bali KK, Kuner R, Gomis A. A critical role for Piezo2 channels in the mechanotransduction of mouse proprioceptive neurons. Sci Rep. 2016 May 17; 6(1):1-9.

36. Woo S-H, Lukacs V, de Nooij JC, Zaytseva D, Criddle CR, Francisco A, et al. Piezo2 is the principal mechanotransduction channel for proprioception. Nat Neurosci. 2015 Dec 9; 18(12):1756-62. https://doi.org/10.1038/nn.4162 PMID: 26551544

37. Alper SL. Genetic Diseases of PIEZO1 and PIEZO2 Dysfunction. Curr Top Membr. 2017 Jan 1; 79:97-134. https://doi.org/10.1016/bs.ctm.2017.01.001 PMID: 28728825

38. Nambot S, Thevenon J, Kuentz P, Duffourd Y, Tisserant E, Bruel AL, et al. Clinical whole-exome sequencing for the diagnosis of rare disorders with congenital anomalies and/or intellectual disability: Substantial interest of prospective annual reanalysis. Genet Med. 2018 Jun 1; 20(6):645-54. https://doi.org/10.1038/gim.2017.162 PMID: 29095811

39. Yang Y, Muzny DM, Xia F, Niu Z, Person R, Ding Y, et al. Molecular Findings Among Patients Referred for Clinical Whole-Exome Sequencing. JAMA. 2014 Nov 12; 312(18):1870-9. https://doi.org/10.1001/jama.2014.14601 PMID: 25326635

40. Baysal E, Indrak K, Bozkurt G, Berkalp A, Aritkan E, Old JM, et al. The β-thalassaemia mutations in the population of Cyprus. Br J Haematol. 1992; 81(4):607-9.

41. Stojanovic JR, Miletic A, Peterlin B, Maver A, Mijovic M, Borlja N, et al. Diagnostic and Clinical Utility of Clinical Exome Sequencing in Children With Moderate and Severe Global Developmental Delay/Intellectual Disability. J Child Neurol. 2020 Feb 1; 35(2):116-31. https://doi.org/10.1177/0883073819879835 PMID: 31623504

42. Taberner FJ, Prato V, Schaefer I, Schrenk-Siemens K, Heppenstall PA, Lechner SG. Structure-guided examination of the mechanogating mechanism of PIEZO2. Proc Natl Acad Sci U S A. 2019 Jul 9; 116 (28):14260-9. https://doi.org/10.1073/pnas.1905985116 PMID: 31235572

43. Wang L, Zhou H, Zhang M, Liu W, Deng T, Zhao Q, et al. Structure and mechanogating of the mammalian tactile channel PIEZO2. Nature. 2019 Sep 12; 573(7773):225-9. https://doi.org/10.1038/s41586-019-1505-8 PMID: 31435011

44. Coste B, Houge G, Murray MF, Stitziel N, Bandell M, Giovanni MA, et al. Gain-of-function mutations in the mechanically activated ion channel PIEZO2 cause a subtype of Distal Arthrogryposis. Proc Natl Acad Sci U S A. 2013 Mar 19; 110(12):4667-72. https://doi.org/10.1073/pnas.1221400110 PMID: 23487782

45. McMillin MJ, Beck AE, Chong JX, Shively KM, Buckingham KJ, Gildersleeve HIS, et al. Mutations in PIEZO2 Cause Gordon Syndrome, Marden-Walker Syndrome, and Distal Arthrogryposis Type 5. Am J Hum Genet. 2014 May 1; 94(5):734-44. https://doi.org/10.1016/j.ajhg.2014.03.015 PMID: 24726473

46. Chesler AT, Szczot M, Bharucha-Goebel D, Čeko M, Donkervoort S, Laubacher C, et al. The Role of PIEZO2 in Human Mechanosensation. N Engl J Med. 2016 Oct 6; 375(14):1355-64. https://doi.org/10.1056/NEJMoa1602812 PMID: 27653382

47. Schrider DR, Hourmozdi JN, Hahn MW. Pervasive multinucleotide mutational events in eukaryotes. Curr Biol. 2011 Jun 21; 21(12):1051-4. https://doi.org/10.1016/j.cub.2011.05.013 PMID: 21636278

48. Besenbacher S, Sulem P, Helgason A, Helgason H, Kristjansson H, Jonasdottir A, et al. Multi-nucleotide de novo Mutations in Humans. Petrov DA, editor. PLOS Genet. 2016 Nov 15; 12(11):e1006315. https://doi.org/10.1371/journal.pgen.1006315 PMID: 27846220

49. Kaplanis J, Akawi N, Gallone G, McRae JF, Prigmore E, Wright CF, et al. Exome-wide assessment of the functional impact and pathogenicity of multinucleotide mutations. Genome Res. 2019 Jul 1; 29 (7):1047-56. https://doi.org/10.1101/gr.239756.118 PMID: 31227601

50. Jiao Q, Sun H, Zhang H, Wang R, Li S, Sun D, et al. The combination of whole-exome sequencing and copy number variation sequencing enables the diagnosis of rare neurological disorders. Clin Genet. 2019 Aug 1; 96(2):140-50. https://doi.org/10.1111/cge.13548 PMID: 30945278

51. Pfundt R, del Rosario M, Vissers LELM, Kwint MP, Janssen IM, de Leeuw N, et al. Detection of clinically relevant copy-number variants by exome sequencing in a large cohort of genetic disorders. Genet Med. 2017 Jun 27; 19(6):667-75. https://doi.org/10.1038/gim.2016.163 PMID: 28574513