Introduction Phelan-McDermid syndrome (PMS) is caused by SHANK3 haploinsufficiency. Its wide phenotypic variation is attributed partly to the type and size of 22q13 genomic lesion (deletion, unbalanced translocation, ring chromosome), partly to additional undefined factors. We investigated a child with severe global neurodevelopmental delay (NDD) compatible with her distal 22q13 deletion, complicated by bilateral perisylvian polymicrogyria (BPP) and urticarial rashes, unreported in PMS.
Methods Following the cytogenetic and array-comparative genomic hybridization (CGH) detection of a r(22) with SHANK3 deletion and two upstream duplications, whole-genome sequencing (WGS) in blood and whole-exome sequencing (WES) in blood and saliva were performed to highlight potential chromothripsis/chromoanagenesis events and any possible BPP-associated variants, even in low-level mosaicism.
Results WGS confirmed the deletion and highlighted inversion and displaced order of eight fragments, three of them duplicated. The microhomology-mediated insertion of partial Alu-elements at one breakpoint junction disrupted the topological associating domain joining NFAM1 to the transcriptional coregulator TCF20. WES failed to detect BPP-associated variants.
Conclusions Although we were unable to highlight the molecular basis of BPP, our data suggest that SHANK3 haploinsufficiency and TCF20 misregulation, both associated with intellectual disability, contributed to the patient’s NDD, while NFAM1 interruption likely caused her skin rashes, as previously reported. We provide the first example of chromoanasynthesis in a constitutional ring chromosome and reinforce the growing evidence that chromosomal rearrangements may be more complex than estimated by conventional diagnostic approaches and affect the phenotype by global alteration of the topological chromatin organisation rather than simply by deletion or duplication of dosage-sensitive genes.
- topologically associating domains
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
The Phelan-McDermid syndrome (PMS; MIM#606232) is a neurodevelopmental disorder caused by the deletion or interruption of SHANK3 gene at chromosome 22q13,1–5 leading to haploinsufficiency of SHANK3 protein, which is a key structural component of the glutamatergic postsynaptic density.6 Patients with PMS mainly exhibit neurological and behavioural symptoms including hypotonia, intellectual disability (ID), impaired language development (delayed or absent speech), behavioural features consistent with disorders of the autistic spectrum and seizures.
The penetrance and expressivity of these symptoms may be variable, with different degrees of severity at different ages3 7–9 as well as congenital malformations,4 5 even among patients with the uniformly-sized recurrent 22q13.3 deletion2 or interstitial deletions affecting only SHANK3. 3 10
It is still uncertain whether the phenotypic variability is caused by variants regarding the genes and expression controlling elements present on the corresponding region of the intact chromosome 22, by a genome rearrangement more complex than estimated or even by a cryptic mosaic condition of the 22q rearrangement. The latter mechanism may particularly apply to the PMS individuals in whom the distal 22q deletion lies in a ring chromosome. Indeed, rings are intrinsically unstable, resulting in the continuous production of cells having duplicated rings, broken portions of the original ring or even missing the entire chromosome.
In this study, we investigated a patient with a de novo ring chromosome 22, present in all the 50 metaphases analysed, who showed a complex cortical brain abnormality and unusual dermatological features, in addition to the core phenotype of PMS. Since array-CGH had revealed a distal 22q deletion and two non-contiguous duplications, we performed whole-genome sequencing (WGS) to investigate a possible event of chromothripsis/chromoanagenesis. WGS confirmed this hypothesis, highlighting eight breakpoints along the chromosome 22 long arm and disordered assembly of the resulting fragments, with disruption of topological associating domains (TAD) and breakage of four genes.
A 3-year-old female child came to our attention after a ring chromosome 22 was detected by cytogenetics when she was 8 months old. The patient was the first child of healthy unrelated parents with uneventful family history. She was born by caesarean section to a 37-year-old mother and 32-year-old father at 37 weeks of gestational age because of umbilical cord wrapped around her neck. Birth weight was 2000 g (10th–25th centile), length 44.5 cm (25th centile), cranial circumference 31 cm (25th centile); APGAR scores were 8/9 at 1’/5’, respectively. At birth, she was initially transferred to the intensive care unit and then to the neonatal pathology unit for 28 days. Hypoglycaemia and thrombocytopaenia were solved in 4 and 12 days, respectively, while her feeding difficulties, characterised by poor sucking, required gavage. In the first months of life, generalised hypotonia and delayed psychomotor development were observed. She acquired head control at 12 months, sat at the age of 13 months, was standing at 20 months and walked with support at the age of 24 months. She uttered her first words at 12 months but her parents refer no subsequent progression of language. Hypersomnia was present since the age of 2 years.
At the age of 22/12 years (26 months), she began to manifest atopic dermatitis characterised by recurrent pruritic urticarial rushes on the face and the trunk in addition to persistent limb pilaris keratosis. Food-specific IgE antibody tests resulted negative and her family history of atopy was unremarkable. At last evaluation, at the age of 3 years, her weight was 12.25 kg (25th centile), height 95.5 cm (75th centile), cranial circumference 50.50 cm (75th centile). She showed ptosis of her left eye and long face but no major dysmorphic features. Griffiths scale scores revealed severe ID (IQ: 24, mental age equivalent 8 months) and language was absent. Poor eye contact and stereotypic hand movements characterised her neurological picture. Diffuse mild muscular hypotrophy and weakness were present, with trunk hypotonia, unsteadiness of orthostatic position and ataxic gait with support. Bilateral coxa valga, hyporeflexia in the lower limbs, together with bilateral Babinski sign, were present. Orofacial hypotonia, with drooling and bruxism, and severe swallowing disorder were observed. She suffered episodes of gastro-oesophageal reflux. Kidney ultrasound, urinalysis and function blood tests were normal. Atopic dermatitis was still present (figure 1A). The brain 3T MRI showed complex brain anomalies including bilateral perisylvian polymicrogyria (BPP) (figure 1B). The white matter volume was diffusely reduced. The third and lateral ventricles were markedly enlarged, probably because of a partial stenosis of the aqueduct caused by a large cyst compressing the lamina quadrigemina. The corpus callosum was thinned. An arachnoid cyst with mass effect compressed the right cerebellar hemisphere. EEG recording showed slowing of background activity during wakefulness, while the physiologic features (K-complex and spindles) during sleep were rare with irregular morphologies. Slow waves, of greater amplitude than the background rhythm, were present in the posterior areas, bilaterally, without side predominance. On electromyography (EMG) derivation, isolated myoclonic jerks were recorded without EEG correlates (figure 1C). She never presented episodes of epileptic febrile or non-febrile seizures and no antiepileptic drugs were administered so far.
Genomic DNA was extracted from proband’s, parents’ and maternal grandparents’ blood and saliva with standard protocols.
Array-CGH analysis was performed on proband’s DNA from blood using an Agilent Human Genome CGH Microarray Kit 4×180 k with an overall median probe space of 13 Kb. All nucleotide positions refer to the Human Genome, February 2009 Assembly (hg19). Data analysis was performed using Agilent Cytogenomics V.184.108.40.206. High-resolution G-banding metaphase chromosomes and interphase nuclei were obtained from PHA-stimulated blood lymphocytes cultured for 72 hours, according to standard protocols. Dual and three-colour FISH experiments with prokaryotic artificial chromosome (PAC) probes, selected according to the UCSC Human Genome Browser (GRCh37/hg19) and obtained from the human library RPCI-4 and 3, were performed as previously described.11 The pan-telomeric peptide nucleic acid (PNA) probe (PNA FISH kit/Cy3, Dako Denmark A/S) which recognises the consensus sequence (TTAGGG)n of human telomeres was hybridised according to the manufacturer’s instructions.
Whole-genome sequencing (WGS) and Sanger confirmation of breakpoint junctions
Sequencing libraries to perform WGS were prepared with Illumina’s TruSeq DNA PCR Free kit (San Diego, California, USA) according to manufacturer’s protocol. We obtained an average coverage of 42X in the father’s sample and 37X in the proband’s sample. Reads from the fastq files were mapped to the human reference genome GRCh37/hg19 using Isaac Genome Alignment Software (V.iSAAC-03.16.06.06).12 Coverage graphs were performed plotting the average coverage of 1000 bp windows sliding over the whole chromosome length, in order to identify large deletions and duplications in each chromosome. Structural variants were called using Lumpy (V.0.2.12) and Manta (V.0.29.6).13 14 Then these structural variants were visualised and manually checked in Integrative Genomics Viewer (IGV) genome browser (V.2.3.72) in order to reconstruct chromosomes showing chromothripsis events. Sanger sequencing was performed in order to verify the structural variations predicted by WGS findings. Primer3 software was used to design primers for sequencing de novo breakpoint junctions of structural variations. PCR was performed by using Taq Polymerase (Invitrogen) with elongation times of 1 min. Sequencing reads were analysed manually to characterise the exact breakpoints.
Whole-exome sequencing (WES)
To test for any possible constitutional or mosaic mutation that could explain the polymicrogyria, proband’s DNA from both blood and saliva was analysed by OneSeq protocol, at an average coverage of 194.3x and 134.7x, respectively, on an Illumina HiSeq2500, as previously described.15 After reads alignment, variant calling, annotation and filtering for possible artefacts was done as previously reported.16
The sequences of polymicrogyria candidate genes (online supplementary table 1 reported in Parrini et al 17) were carefully examined together with PI4KA.18 Sanger sequencing was performed to validate the only variant possibly associated with polymicrogyria, SRPX2 at Xq22.1, by using standard methods, extending the analysis to the three-generation family members and performing X-inactivation investigation.
Topologically associating domains (TADs) analysis
Search for TADs was performed by using the web-based 3D Genome Browser.19
The patient’s karyotype was 46,XX,r(22)(::p11.2>q13::)dn with no evidence of ring instability in the 50 metaphases analysed. Array-CGH showed several copy number changes involving the ring chromosome 22, potentially representing a complex chromosome rearrangement. The 22q13 band showed a multiple rearrangement pattern: a duplication (dup) followed by a normal copy region (nlm), a second duplication and a terminal deletion (dup-del) (figure 2). This rearranged region (dup-nlm-dup-del) of 21 Mb spanned at least four breakpoints, with a copy-number loss of approximately 8.4 Mb and copy-number gain of approximately 2.5 Mb. Parental origin analysis showed that both duplication and deletion involved the paternal chromosome 22 (online supplementary table 1).
Whole-genome sequencing (WGS)
WGS analysis, performed on blood samples from proband and father (average coverage 37X and 42X, respectively), revealed eight breakpoints in the proband, producing nine fragments from A to I at the q arm of chromosome 22 (figure 3A). According to array-CGH analysis (figure 2) and visual inspection of chromosome 22 by IGV (Integrative Genome Viewer), fragment B (22q11.23–q12.1; size 257 Kb), fragment E (22q13.1–q13.2; size 2,4 Mb) and fragment G (22q13.2; size 160 Kb) were found to be duplicated, while fragment I (22q13.2q13.33; size 8.4 Mb) was deleted. In addition to the copy number variations detected by array-CGH (fragments E, G, I), inversions and displaced order of fragments were highlighted, indicating a chromoanagenesis event. In more detail, fragment H (1 kb) appeared to be inverted and fused with fragment G (termed as GHinv fusion joint), showing a 1 kb length inversion event between the distal duplicated and deleted segments (figure 3B, online supplementary figure 1). Additionally, a copy of fragment G (160 kb) was inverted and fused with the start site of fragment D (termed as GinvD fusion joint), showing an inverted insertion event within chromosome 22 (figure 3B). Sanger sequencing validated both GHinv and GinvD rearrangements (online supplementary figure 2 and table 2) and revealed an unexpected 183 bp inverted insertion between fragments G and H of a chromosome 14 sequence (chr14:105 518 942–105 519 123) (figure 3B), consisting of two portions of an AluSq sequence, as denoted by the UCSC RepeatMasker track (online supplementary figure 3). The breakpoint characteristics of the GHinv junction involved a microhomology of 2 bp (5’-GA) in fragment G-chr14 fusion junction and a microhomology of 2 bp (5’-TC) with the insertion of a single bp (5’-G) in chr14-Fragment H fusion junction (online supplementary figure 4). The breakpoint junction of H-I (chr22:42814464) falls within another subfamily of Alu element (AluSx3 sequence, as denoted by the UCSC RepeatMasker track, online supplementary figure 5). The GinvD rearrangement had a blunt junction without microhomology (figure 3B; online supplementary figure 4). Sanger investigations could not confirm the other predicted breakpoint-fusions due to the presence of repeated sequences. Dual-colour and three-colour FISH experiments with PACs probes (figure 4C), either in metaphases or interphase nuclei, proved that one copy of fragment E was an insertion-duplication, located close to the centromere of the r(22) (figure 4D). This might explain why we failed to define the fusion point between the end of fragment E (chr22: 41 503 612) and the repeated pericentromeric sequences (online supplementary figure 7). The FISH data fit well with the proposed model predicting a fusion between fragment E and the centromeric region of chromosome 22, thus forming a closed ring structure (figure 4D). We also verified the absence of interstitial pan-telomeric sequences with a PNA probe (online supplementary figure 8). Since the 22q deletion and duplications was of paternal origin, the father’s DNA was investigated by WGS for the presence of any genomic rearrangement predisposing to the chromosome 22 abnormality in the proband. We detected two tandem duplications, located at Chr22: 23 015 179–23 080 761 and Chr22: 25 684 091–25 980 630, reported in DGV as well in our in-house WGS database (5 out of 21 unrelated cases and 1 out of 11, respectively).
Analysis of genes and TADs at breakpoint regions
Four of eight breakpoints detected by WGS disrupted four genes (figure 3B): EP300 (MIM 602700) has a role in the neurodevelopmental disorder Rubinstein-Taybi syndrome 2 (RSTS2, MIM#613684); NFAM1 (MIM 608740) is involved in immune system functioning; MYO18B (MIM 607295) is associated with Klippel-Feil syndrome 4 (KSF4, MIM#616549); GTPBP1 (MIM 602245) encodes a protein belonging to the family of GTP-binding proteins. TAD analysis showed that the GHinv breakpoint within NFAM1 (figure 2B) lies in a TAD that joins NFAM1 to TCF20 (figure 4).
WES analysis was conducted on blood and saliva samples from the proband and her parents in order to exclude single nucleotide variants that could explain the BPP, a sign never before reported in PMS. In fact, the coverage obtained by WGS (30–40×) could have been insufficient to highlight any possible mosaic variant known to be responsible for polymicrogyria.20 WES analysis revealed a heterozygous variant in the exon 8 of SRPX2 at Xq22.1, NM_014467.2:c.893G>A (NP_055282.1:p.Arg298His), in both proband and mother as well as in the maternal grandmother (online supplementary figure 9). This substitution is reported in several databases (dbSNP, ExAC, ESP, COSMIC) and has a frequency in the population of about 0.01%. Moreover, it affects a highly conserved amino-acid residue of the Sushi/SCR/CCP domain (InterPro#IPR000436) but with contrasting predictions of pathogenicity, according to various in silico tools (eg, SIFT, Mutation Taster, Align GVGD). Since the BPP could also be caused by somatic mosaic variants in PIK3CA, we tried to exclude this possibility by analysing proband’s blood and saliva DNA by OneSeq, under experimental conditions allowing an average coverage higher than 130x. The analysis gave normal results. Similarly, careful analysis, including manual IGV inspection, of the exonic sequence of PI4KA 18 and SNAP29, 17 21 both on chromosome 22 and reported in association with cortex disorders, did not reveal any variants.
We describe a complex genomic rearrangement associated with a constitutional ring chromosome 22 detected in a 3-year-old patient with severe neurodevelopmental disorder. Haploinsufficiency for SHANK3, fitting with the presence of the r(22), led to her initial framing within the PMS syndrome. However, the later MRI detection of BPP and the occurrence of recurrent episodes of scattered atopic dermatitis indicated a more complex molecular basis for her condition. Among the brain-imaging studies in patients with PMS, the most frequent abnormalities described are thinning of the corpus callosum (36%), white matter changes (39%), ventricular dilatation (32%) and arachnoid cysts (14%).22 Other findings showed that some patients with PMS have severe posterior cerebellar vermis hypoplasia and enlarged cisterna magna.23 Indeed, our patient’s complex brain malformations included many common brain abnormalities frequently observed in PMS individuals (thin corpus callosum, reduced volume of white matter, enlarged ventricles and arachnoid cysts). However, with the limitation that not all subjects with PMS have undergone MRI analysis, polymicrogyric cortex has never been reported in patients with 22q distal deletions,5 23 24 suggesting that mechanisms other than simple haploinsufficiency for genes located at distal 22q are responsible for it. By WGS and confirmatory WES data (the latter in both blood and saliva, at a coverage >100 x), we excluded pathogenic variants in genes associated with polymicrogyria and cortex abnormalities having a frequency <1%, and in theory, also mosaic condition, as most commonly reported for PIK3R2.17 20 Indeed, the c.893G>A (p.Arg298His) variant in the in SRPX2 gene on chromosome X, coding for a protein involved in synaptogenesis of the cerebral cortex, was the only one taken into consideration. A variant in exon 4 of SRPX2 was previously described in a male patient with BPP, his unaffected mother and two maternal aunts with mild ID and no MRI abnormality.25 Similarly, our variant was present in the unaffected mother and grandmother and X-inactivation studies revealed that the active X carried the variant in the proband but not in her mother (online supplementary figure 9). Amino-acid residue Arg298 is targeted by other missense or synonymous variations (namely p.Arg298Ser, p.Arg298Cys and p.Arg298Arg), but p.Arg298His is the most frequently observed, with 11 hemizygous subjects reported in the gnomAD browser. Therefore, based on currently available information, this variant does not seem to have a clear pathogenic role.
None of the genes disrupted by the breakpoints of the ring, namely EP300, MYO18B, NFAM1 and GTBP1, was associated with cortex abnormalities, although some of them may contribute to severity of the patient’s phenotype.
Regarding NFAM1, encoding a protein regulating B-cell development, the 183 bp insertion at GHinv breakpoint within intron 1, might theoretically impair its two transcript variants. The TAD region, disrupted by the GHinv breakpoint, joins NFAM1 to TCF20 (figure 4). De novo variants in TCF20 (OMIM 603107), encoding a widely expressed transcriptional coregulator structurally and functionally related to RAI1 (OMIM 607642), have been implicated in the aetiology of autism spectrum disorders, ID/developmental delay, schizophrenia and postnatal overgrowth.26–29 Thus, misregulation of TCF20 as a consequence of TAD perturbation may likely contribute to the neurodevelopmental features observed in our patient. Two patients with haploinsufficiency of NFAM1 resulting from interstitial 22q13.2 deletion have been reported with atopic dermatitis and urticarial rash,30 31 as observed in our patient (figure 1A), but on closer examination one of them30 has a terminal deletion that does not include NFAM1. With the limitation that RNA was not available for NFAM1 expression analysis in our patient, we cannot exclude that other genes may contribute to the urticarial rashes.
EP300 frameshift mutations were reported in Rubinstein-Taybi syndrome 2, with some phenotypic overlap with our patient (speech, gross motor delay, feeding/swallowing issues). The interruption of MYO18B does not seem to interfere with the phenotype of the patient because it is associated with an autosomal recessive condition, namely the Klippel-Feil syndrome 4, while the function of GTBP1 is still unknown. Altogether, these data explain only part of the patient’s clinical picture, including the dermatological manifestation. However, it seems likely that disruption of the TAD, either as a consequence of the duplication/deletion regions or the ring conformation itself,32 played an equally important role in her clinical manifestations
Chromoanagesis and ring chromosomes
We describe the first case of a constitutional ring chromosome derived by a chromoanagenesis event, previously reported only in cancer.33 In line with the predominant occurrence of chromoanagenesis constitutional events on paternally derived chromosomes,34 our r(22) was of paternal origin and we excluded through WGS that any bona fide predisposing rearrangement was present on paternal chromosomes 22. According to the mechanisms associated with chromothripsis and chromoanasynthesis, both belonging to the novel category of chromoanagenesis events,34 35 we have to assume that the rearrangement occurred during either paternal gametogenesis or early embryogenesis by isolation of a chromosome 22 within a micronucleus or, alternatively, as a consequence of the formation of a chromatin bridge interconnecting two daughter cells.36 In the first case, chromosome mis-segregation might have led to the trapping of chromosome 22 within the micronucleus, as well described in chromothripsis, while in the second the occurrence of a dicentric chromosome 22, pter>q13::q13-pter, possibly secondary to telomeric crisis, would have promoted chromothripsis events, as previously proposed in cancer genomes.37 38 The sequences at the breakpoints of our chromosome 22 are consistent with repair-based mechanisms, although further complexity was added at the GHinv fusion point by a 183 bp insertion from chromosome 14, containing two partial Alu elements, in contrast to the full length 300 bp Alu sequence.39 The sequence at the junction of AluSq insertion, did not resemble L1 endonuclease consensus site (5′-YYYY/RR-3′ where Y=pyrimidine, R=purine) nor its additional variations,40 41 making unlikely any L1 endonuclease role in the rearrangement. According to all these features at GHinv junction, we might predict an endonuclease-independent non-classical Alu insertion42 with a potential role in DSB repair. An Alu-Alu mediated insertion may also be an option since a second Alu was detected at H-I breakpoint (online supplementary figure 5) as well as in other chromothripsis/chromoanasynthesis rearrangements.43 44
Ring chromosomes are unstable pieces of circular chromatin that, as a consequence of sister chromatid exchanges, may lead, at least in vitro, to a mosaic condition with different cells carrying the entire duplicated or deleted ring chromosome or part of it.45 This in vivo instability should result in massive cell death, in turn causing the so-called ring syndrome that, independently from the chromosome involved, would cause extreme growth failure without major malformations, with only a few or no minor anomalies and mild to moderate ID.46 Although the existence of this syndrome remains highly hypothetical, a compensatory uniparental disomy mechanism has been demonstrated in induced pluripotent stem cells (iPSCs) obtained from patient fibroblasts containing deleted ring chromosomes, with loss of cells containing the abnormal chromosome and duplication of the wild type homologue.47 Whatever the meaning of this cell-autonomous correction of ring chromosomes in vivo, it seems clear that circular chromatin structures are potentially related to decompensation of cellular homeostasis via massive gene deregulation. In our study, where no ring chromosome instability was highlighted at least in blood, we were unable to detect the molecular basis for the patient’s BPP, although we can hypothesise that it may be linked to insufficient production during embryogenesis of proteins encoded by some chromosome 22 genes with a role in the formation of the brain cortex.
Although rearrangements more complex than the classical model, which considers the circular conformation a consequence of the deletion of either chromosome end, have been reported,48 it is likely that other ring chromosomes may be generated by a mechanism such as the one we describe here. This hypothesis might also explain the genotype–phenotype discrepancies reported for other ring chromosomes, particularly those of chromosomes 14 and 20, in which the severe phenotype cannot be explained based on the deleted region.49 50
We demonstrate for the first time the occurrence of a constitutional chromoanagenesis event resulting in the formation a ring chromosome, with both copy number losses and gains, leading to a PMS phenotype associated with BPP and atopic dermatitis with urticarial rashes. The outcome of concurrent deletions, duplications, gene interruptions and TAD is a dramatically complex phenotype. This study reinforces the growing evidence that chromosomal rearrangements (1) may be not as simple as they appear by conventional diagnostic approaches and (2) affect the phenotype by global alteration of the topological chromatin organisation rather than simply by the deletion or duplication of dosage-sensitive genes.
3D Genome Browser (http://3dgenome.org)
gnomAD browser (http://gnomad.broadinstitute.org/)
We would like to gratefully acknowledge the family participating in this study.
Contributors NK, MD, LX, MR performed WGS analysis. EE performed TAD and WES analysis. FA interpreted Brain MRI analysis. CZ assessed EEG finding. CM and MGD conducted clinical neurological evaluations and managed the patient. SiB and RG performed segregation analysis of the SRPX2 variant and X-inactivation analysis. MCB interpreted cytogenetics, array-CGH and FISH investigations. SB provided technical assistance to cytogenetic analysis. MCB and OZ designed the study. MCB, OZ, NK and EE contributed significantly to this study, wrote the article and critically discussed the final editing of the manuscript.
Funding MCB was supported by grants of the Italian Ministry of Health (Ricerca Corrente 2017). OZ was supported by Telethon Italy Grant GGP13060.
Competing interests None declared.
Patient consent Obtained.
Ethics approval The study has E. Medea Scientific Institute Research Ethics Committee approval (Prot. no. 007/15-CE).
Provenance and peer review Not commissioned; externally peer reviewed.