- 1Unit of Medical Genetics, IRCCS Istituto “Giannina Gaslini”, Genoa, Italy
- 2Laboratory of Veterinary Physiology, School of Veterinary Medicine, Faculty of Veterinary Science, Nippon Veterinary and Life Science University, Musashinoi, Japan
- 3Department of Molecular Carcinogenesis, Medical University of Łódź, Łódź, Poland
- 4Fetal and Perinatal Pathology Unit, IRCCS Istituto “Giannina Gaslini”, Genoa, Italy
- 5Department of Precision Medicine, University of Campania “Luigi Vanvitelli”, Naples, Italy
- 6Telethon Institute of Genetics and Medicine (TIGEM), Naples, Italy
- 7Neuroradiology Unit, IRCCS Istituto “Giannina Gaslini”, Genoa, Italy
- 8Fetal Medicine and Surgery Unit, IRCCS Istituto “Giannina Gaslini”, Genoa, Italy
- 9Medical Genetics Unit, Maria Paternò Arezzo Hospital, Ragusa, Italy
- 10Pediatric Neurology and Muscular Diseases Unit, IRCCS Istituto “Giannina Gaslini”, Genoa, Italy
- 11Department of Neurosciences, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DiNOGMI), University of Genoa, Genoa, Italy
- 12Department of Surgical Sciences and Integrated Diagnostics (DISC), Pathology Division of Anatomic Pathology, University of Genoa, Genoa, Italy
- 13Department of Neuromuscular Diseases, Queen Square Institute of Neurology, University College London, London, United Kingdom
Mutations in the WWOX gene cause a broad range of ultra-rare neurodevelopmental and brain degenerative disorders, associated with a high likelihood of premature death in animal models as well as in humans. The encoded Wwox protein is a WW domain-containing oxidoreductase that participates in crucial biological processes including tumor suppression, cell growth/differentiation and regulation of steroid metabolism, while its role in neural development is less understood. We analyzed the exomes of a family affected with multiple pre- and postnatal anomalies, including cerebellar vermis hypoplasia, severe neurodevelopmental impairment and refractory epilepsy, and identified a segregating homozygous WWOX mutation leading to a premature stop codon. Abnormal cerebral cortex development due to a defective architecture of granular and molecular cell layers was found in the developing brain of a WWOX-deficient human fetus from this family. A similar disorganization of cortical layers was identified in lde/lde rats (carrying a homozygous truncating mutation which disrupts the active Wwox C-terminal domain) investigated at perinatal stages. Transcriptomic analyses of Wwox-depleted human neural progenitor cells showed an impaired expression of a number of neuronal migration-related genes encoding for tubulins, kinesins and associated proteins. These findings indicate that loss of Wwox may affect different cytoskeleton components and alter prenatal cortical development, highlighting a regulatory role of the WWOX gene in migrating neurons across different species.
Introduction
WWOX (MIM 605131) encodes a WW domain-containing oxidoreductase known to be involved in several cellular and biological processes, including tumor suppression, cell growth and differentiation, regulation of steroid metabolism, induction of apoptosis and central nervous system (CNS) development (Bednarek et al., 2001; Aldaz et al., 2014; Iatan et al., 2014; Tanna and Aqeilan, 2018). Biallelic pathogenic WWOX variants have been implicated in WWOX-related epileptic encephalopathy (WOREE syndrome; MIM 616211), a very rare and severe neurological disorder characterized by absence of language development and acquisition of walking, drug-resistant seizures, ophthalmological involvement and a high likelihood of premature death (Tabarki et al., 2015a, b; Johannsen et al., 2018; Piard et al., 2018). Consistent with the human disease phenotypes, Wwox knockout (KO) mice exhibit ataxic gait and epilepsy since birth and die prematurely (Mallaret et al., 2014). Lethal dwarfism with epilepsy (lde/lde) rats, carrying a spontaneous homozygous frameshift mutation in the gene, show a high incidence of ataxia, audiogenic seizures, dwarfism and risk of premature death (Suzuki et al., 2009), and display pathological features of abnormal myelination and a marked reduction in astrocytes and microglia (Tochigi et al., 2019). To date, however, no study has assessed the impact of pathogenic WWOX variants on the development of cerebral cortex, and the neuropathological features of WOREE syndrome have not been investigated before. As a result, disease mechanisms implicating loss of Wwox in CNS development and diseases have not been fully understood yet. In this study, we report a WOREE syndrome family with a 4-year-old male proband presenting wide congenital anomalies, postnatal microcephaly and developmental epileptic encephalopathy, and a male fetus diagnosed by prenatal imaging with cerebellar vermis hypoplasia and terminated at 21th gestational week. A homozygous segregating WWOX c.790C > T mutation, predicted to introduce a premature stop codon at protein position 264, was found by analyzing the exomes of the family. Histopathological studies of the fetal brain revealed an abnormal formation of the developing cerebral cortex and similarly, altered cortical layers were identified in the developing brains of the lde/lde rats. Analysis of transcriptome of human neural progenitor cells (hNPC) and neurons with silenced Wwox revealed a reduced expression of different cytoskeleton-related genes. Our findings expand pre- and postnatal features associated with pathogenic WWOX variants and indicate that Wwox deficiency may affect neuronal migration and impair cerebral cortication, ultimately contributing to WOREE-related clinical and molecular phenotypes.
Materials and Methods
Patient data were reviewed for the medical history, the genetic and radiological investigations. Whole exome sequencing (WES) was performed after written informed consent was obtained from parents, according to “Giannina Gaslini” Children’s Hospital Institution review board.
Genetic Studies
Genomic DNA was isolated from 1 ml of peripheral blood using QIAamp DNA Blood Midi kit (Qiagen). Genomic DNA was enriched with SureSelect Clinical research exome 54Mb (Agilent Technologies). Whole exome sequencing (WES) runs were performed on the HiSeq1000 platform (Illumina Inc., San Diego, CA, United States) as 100-base paired-end runs, using a standard Illumina pipeline as previously described (Salpietro et al., 2019). Paired-end reads were mapped to the reference human genome sequence (GRch37/hg19) using CLC Bio Genomics Workbench 7.5.1 software (CLC Bio, Aarhus, Denmark), after removal of duplicates. Single-nucleotide polymorphisms (SNPs) and short deletion or insertion (indels) variants were called by CLC Bio Workbench software using the specific variant calling plugin and dbSNP147 and gnomAD databases. The variants were filtered for in-house variants controls. In accordance with the pedigree and phenotype, our filtering strategy prioritized rare (<1% in public databases, including 1000 Genomes project and gnomAD) biallelic and X-linked hemizygous coding variants and/or variants located in genes previously implicated in microcephaly and developmental epileptic encephalopathies. In order to investigate the presence of homozygosity regions, we performed homozygosity mapping analyzing the proband variant call format (VCF) WES data on HomozygosityMapper1. Validation, parental origin of the resulting variants and family segregation studies were performed by Sanger sequencing.
Neuropathological Studies
After expulsion of the WWOX-deficient male fetus at the 21th gestational week, his whole body was fixed in 10% buffered formalin approximately 6 h after termination. Autopsy was performed 8 weeks later. During autopsy, a sample of brain tissue (parietal coronal section) was taken, after the brain hemispheres had been exposed by means of cutting the scissurae of the parietal and frontal bones. We obtained samples of brain tissue from a control male fetus at the 21th gestational week by using the same method and fixed them 4.5 h after death in 10% buffered formalin for 2 weeks. For microscopy, tissue blocks were fixed in formalin and embedded in paraffin in accordance with standard procedures (Jansen et al., 2016). Sections (7 μ) were deparaffinized and hydrated in decreasing ethanol concentrations. Haematoxylin and eosin (H&E; Ventana HE 600 system) was performed in order to obtain a histo-architectural overview. Periodic Acid Schiff (PAS; Ventana Staining Kit 860-014) staining was performed by means of Benchmark special stains (Roche), in order to document the fine structure of the glial trajectories (Mora et al., 2007). Two further immunohistochemical reactions using CD31 (mouse monoclonal antibody JC70 Cell Marque - Roche) and GFAP (Glial Fibrillary Acidic Protein EP672Y rabbit monoclonal antibody 760-4345 Roche) were performed to investigate the superficial vascular network and processes of radial glial cells (Bär, 1980; Arisio et al., 2002; Jansen et al., 2016).
lde/lde Rat Experiments
Inbred Lethal Dwarfism with Epilepsy rat strain carrying loss-of-function mutation on Wwox gene has been maintained in a clean conventional animal room under a controlled light-dark cycle (14:10 h), with a certified diet (CR-LPF; Oriental Yeast Co., Ltd., Tokyo, Japan) and water ad libitum (Suzuki et al., 2007). All animal experiments and animal care were approved by Animal Care and Use Committee of Nippon Veterinary and Life Science University (Protocol #2019K-28, 31 March 2019), and were conducted in accordance with the Guidelines of the Animal Care and Use Committee of Nippon Veterinary and Life Science University. Pregnant rats obtained from mating between +/lde rats were intraperitoneally injected with 5-bromo-2′-deoxyuridine (BrdU, 100 mg/kg of body weight) (Sigma-Aldrich, St. Louis, MO, United States) at embryonic day (E) 16.5 to evaluate the migration of late-born neuron (Hatanaka et al., 2004; Fan et al., 2008). Resulted pups (both sexes) were dissected under the isoflurane anesthesia at postnatal day 1 (P1). Before removing the brain, whole body was perfused with phosphate-buffered saline (PBS) containing 10 unit/mL heparin followed by 4% buffered paraformaldehyde. After 48 h fixation with paraformaldehyde, isolated brains were embedded in O.T.C compound according to general protocol as previously described (Tochigi et al., 2017). Twenty-μm-thick sequential sagittal cryosections were prepared from 1 to 2 mm lateral region to the midline. Almost corresponding sections were used for comparison. Nissl staining were performed for histological examination. Air-dried cryosections were microwaved with 10 mM sodium citrate (pH 6.0) for 1 min for antigen retrieval. Anti-Satb2 (1:500, clone SATBA4B10, Abcam plc, Cambridge, United Kingdom) and anti-Tbr1 (1:5000, ab31940, Abcam plc) were used to detect layer specific-neuron subtypes (Hevner et al., 2001; Britanova et al., 2005). The positive signals of BrdU were visualized by denaturing with HCl treatment and following immunohistochemical procedure previously described (Tochigi et al., 2017). All histological images were taken on a Biozero BZ-X710 fluorescence microscope using software BZ-X analyzer Ver. 1.3.1.1 (Keyence, Tokyo, Japan). To evaluate the distribution of Satb2, Tbr1, and BrdU-positive cells, cerebral wall was divided evenly into 8 areas from below the marginal zone (MZ) to the intermediate zone (IZ) or 10 areas from below the MZ to the ventricular zone (VZ). The number of positive cells and area (mm2) that contains these cells were measured using analytical software Image J (Ver. 1.46r, NIH, Bethesda, MD, United States) (Tochigi et al., 2019). Data from each 3 +/+ and lde/lde rats genotyped by PCR were used for statistical analysis.
Transcriptomic Analyses of hNPC and Neurons With Silenced Wwox
To analyze the transcriptomic changes caused by WWOX depletion in human neuronal progenitor cells (hNPC), and their putative connection to observed disease phenotypes, we employed the data (deposited in GEO Database with accession number GSE126075), obtained in the experiment of in vitro WWOX silencing in H9-derived hNPC as previously described by some of us (Kośla et al., 2019). WWOX was silenced by a shRNA lentiviral delivery system and the transcriptome was analyzed by CAGE (cap analysis of gene expression) method. Bioinformatic analysis was performed using transcriptomic data from hNPC with silenced WWOX before and after differentiation to neurons in relation to control cells with unaltered WWOX level. The raw data were processed and analyzed to investigate expression of genes associated with neural differentiation and migration, using a selection from Molecular Signature Database of Broad Institute2. From the generated genomic datasets, genes with a fold change >2 were chosen for further investigation using a signaling pathway analysis. The Gene Set Enrichment Analysis (GSEA) allowed for indication of leading-edge genes contribution to neural differentiation and migration (t-test with a weighted scoring scheme and a permutation type regarding phenotype, significance threshold set as FDR < 0.25). Global gene expression analysis was performed using NOISeq software with NOISeq-sim algorithm, as previously described (Tarazona et al., 2015) Volcano plots show the differential gene expression profiles of control neural progenitors versus WWOX-depleted neural progenitors and control neurons versus WWOX-depleted neurons. P significance was calculated according to NOISeq-sim algorithm. Neural crest specific genes discussed in the manuscript are marked in volcano plots. Immunostaining of WWOX-deficient brain was also performed (Supplementary Notes) to confirm hints from transcriptomic studies.
Results
Antenatal and Postnatal Features of WOREE Syndrome Associated to Homozygous p.Arg264Ter Variant
Patient II.1 was born at 38 weeks after an uneventful pregnancy as the first child to Sicilian non consanguineous healthy parents. There was no history of neurological and/or genetic diseases in the family. At birth growth parameters were within normal limits. Bilateral clubfoot and cryptorchidism were noticed since the first day of life. During the neonatal period, poor spontaneous movements and joint contractures were also observed. At the age of 45 days, epilepsy started with generalized tonic-clonic seizures occurring daily and lately evolving into multidrug-resistant tonic and myoclonic seizures. At the age of 3 years, a percutaneous endoscopic gastrostomy (PEG) tube was placed due to feeding difficulties and failure to thrive. At the time of this report, the boy is 4 years old and he has not acquired any developmental milestones. His occipitofrontal circumference (OFC) is 43 cm (<3.8 SD). He shows some distinctive craniofacial features, including a round hypotonic face and short neck (Figure 1D). He shows no response to visual stimuli and fundoscopy examination revealed bilateral pale optic disks, suggestive of optic atrophy. Neurological examination showed severe truncal hypotonia, arthrogryposis, spastic tetraparesis associated with some dystonic movements, brisk deep tendon reflexes and clonus. During early infancy, extensive diagnostic and neurometabolic workup yielded normal results and cerebrospinal fluid analysis was unremarkable. Electroencephalogram (EEG) studies identified a burst suppression-like pattern in repeated evaluations. Brain magnetic resonance imaging (MRI) performed at the age of 9 months and 2.4 years demonstrated reduction of brain volume, T2 signal alterations of the parietal periventricular and frontal white matter, and hypoplasia of the cerebellar vermis and corpus callosum (Figure 1E). Patient II.3 was a fetus terminated at 21th gestational week, after the third pregnancy of this Sicilian couple (Figure 1A). During pregnancy, Tesla-3D ultrasound scan revealed increased diastolic blood pressure (>90° percentile), growth parameters were compatible with gestational age. Nuchal lucency was normal. Fetal brain MRI documented mild cerebellar vermis hypoplasia (Figure 1F), that was confirmed on high resolution post-mortem MRI performed after the fetus was terminated at 21th week of gestation.
Figure 1. Family tree, genetic studies, WOREE-associated clinico-radiological features, Wwox protein and WOREE-associated mutations. (A) Pedigree from Family. (B) Electropherograms of carrier parents and index case (II.1) with the c.790C > T homozygous variant introducing a p.Arg264Ter stop codon. (C) Homozygosity mapping reveal a homozygous bloc in the WWOX gene 16q23.2 region. (D) Distinctive craniofacial features in WOREE syndrome due to homozygous p.Arg264Ter variant, including round hypotonic face and short neck. (E) Brain MRI of Patient II.1 performed at 2.4 years of age. (E’) Sagittal T1-weighted image reveals hypoplasia of the corpus callosum (empty arrow) and inferior cerebellar vermis (arrowhead). (E”) Axial T2-weighted images demonstrates mild atrophy of the frontal lobes with associated bilateral white matter hyperintensity (arrows). (F) Fetal MRI (F’) and high-resolution post-mortem MRI (F”) studies of Patient II.3 performed at 21 gestational weeks demonstrate mild hypoplasia of the cerebellar vermis (black arrows). Note the slightly increased thickness of nuchal subcutaneous tissues on fetal MRI (white arrow). The laminar organization of the cerebral hemispheres and cortical gyration are appropriate for the gestational age (not shown). (G) Wwox protein structure and domains and WOREE-associated mutations identified so far.
Identification of the Homozygous p.Arg264Ter Variant
Whole-exome sequencing (WES) of the index case (II-1, Figure 1A) generated a total of 71,369,726 unique reads, with an average on target depth over 150, and >95% of the target bases covered at least 10X. The homozygous c.790C > T; p.Arg264Ter variant in WWOX (NM_016373) was found within a single block of homozygosity consisting in a 1.58 Mb region on chromosome 16 (chr16: 78387189-79966223), as shown by homozygosity mapping analysis (Figure 1C). No other plausible rare variants were identified in the WES data and the c.790C > T WWOX mutation emerged as the explanation for the disease pathophysiology, being fully consistent with the phenotype of WOREE syndrome and segregating within the family (Figure 1B). The p.Arg264Ter variant is predicted highly deleterious (CADD score = 55) and was not observed in approximately 6,100 individuals of European and African American ancestry (NHLBI Exome Sequencing Project) and was reported as pathogenic in three individuals on ClinVar (accession number: VCV000241105). The p.Arg264Ter homozygous variant was confirmed to segregate with the disease in the family (Figures 1A,B). Prenatal testing by Sanger sequencing of individual II.3 identified the p.Arg264Ter variant in the homozygous state, leading to termination of the fetus at 21th week of gestation.
Defective Architecture of Cortical Layers in WWOX-Deficient Human Developing Brain
Apart of cerebellar vermis hypoplasia, there were no anatomical defects or malformations detected at the macroscopic evaluation of the brain of the fetus affected with WOREE syndrome (Supplementary Figures S1A,B). In the histological HE staining of WWOX-deficient developing human brain we observed anomalous migration of the external granular layer within the molecular layer and also the latter appeared as not homogeneous (Nishimura, 1983; Figure 2A; Individual II.3). A parietal section of the cerebral cortex of three fetuses of the same gestational period (sampled as above) showed correct cortical development (Figure 2A; Ctrl). Disorganization of irregularly distributed glial trajectories was also observed with GFAP staining (black arrows; Figure 2B, Individual II.3). In the cerebral cortex, a thin vascular network (CD 31-positive) consisting of abnormally thin vessels within the areas below the molecular layer (with marked reduction of these in the context of the external granular layer) was also evident (black arrows; Figure 2C, Individual II.3). Incorrect migration of the external granular layer and inhomogeneous molecular layer were also evident using the PAS staining, which highlighted the texture of the trajectories below the molecular layer of the cortex (black arrows; Figure 2D, Individual II.3).
Figure 2. Histopathological studies of WWOX-deficient human brain. (A) H&E staining shows an incorrect migration of the external granular layer within the molecular layer in patient II.3 (right panel) compared to a normal fetus (left panel) at the same gestational age. Images at 10X and 20X magnification (molecular layer: ML, external granular layer: EGL, cortical gray: CG). Scale bar: 50 μm. (B) Histological staining with GFAP 10× (upper panel) and 20× (lower panel) in which a disorganization of the irregularly distributed trajectories is observed (black arrows). Scale bars: 50 μm. (C) Histological staining with CD-31 at 10× (upper panel) and 20× (lower panel), in which thinning of the vessels is observed in the context of the external granular layer (black arrows). (D) Histological staining with Shift Reactive Periodic Acid (PAS) at 10× (upper panel) and 20× (lower panel) show that the vascular structures of the cortex are irregularly distributed and irregularly branched (black arrows). Scale bars 50 μm.
Abnormal Neural Migration and Cortical Layer Formation in lde/lde Rats
There were no obvious anatomical defects or significant difference in weight and three-dimensional lengths of brain between +/+ and lde/lde rats at P1 (data not shown). At low magnification of Nissl staining in sagittal brain section, there was no significant difference between +/+ and lde/lde in the thickness of cortical wall (from MZ to VZ/SVZ) (Figures 3A–E). In addition, the development of cerebellum was delayed in lde/lde as shown by reduced number of foliation (Figures 3A,B, arrowheads). Two developing cortical layers (II-IV and V/VI) were recognized in the cortical plate of +/+, whereas boundary between II-IV and V/VI was unclear in lde/lde cortical plate due to irregular cell distribution (Figures 3C,D). To confirm the defects in lde/lde cerebral wall, we performed immunohistochemistry with Satb2, a layer II-IV specific marker (Britanova et al., 2005), and Tbr1, a layer V/VI specific marker (Hevner et al., 2001). Two developing cortical layers (II-IV and V/VI) were clearly recognized as positive signals against Satb2 and Tbr1. There was no significant difference between +/+ and lde/lde in the distribution of Tbr1- positive cells (Figures 3F,H), whereas the density of Satb2-positive cells was significantly decreased in the surface layer of lde/lde cortical plate (Figures 3F,G, area number 1 corresponding to cortical layer II). This defect was accompanied by increase in Satb2-positive cells of the bottom layer in lde/lde cerebral wall (Figures 3F,G, area number 7 corresponding to intermediate zone (IZ). Considerable number of Satb2-positive cells and Tbr1-positive cells were mixed in unclear boundary between II-IV and V/VI as shown in Nissl staining of lde/lde cortical plate (Figure 3F). Since these results clearly indicated that migration of late-born neuron was impaired in lde/lde cortical plate, BrdU was injected fetus at E16.5 when the late-born neuron actively proliferate (Hatanaka et al., 2004; Fan et al., 2008). When BrdU-incorporated cells were traced to P1, we found altered distribution of BrdU-positive cells in lde/lde cerebral walls, revealing delayed migration of late-born neurons (Figures 3I,J). These results indicate that Wwox deficiency impairs prenatal neuronal migration which is required for normal cortical layer formation at birth in rats.
Figure 3. Cortical layer formation and neural migration in lde/lde rat brain compared to normal one. (A–D) Sagittal brain sections at P1 rats were subjected to Nissl staining. Squares as labeled C and D in panels (A,B) indicated identical areas as shown in panels (C,D), respectively. Arrowheads indicate delayed foliation in lde/lde cerebellum. (E) The thicknesses of cerebral wall (MZ to VZ/SVZ) is shown as mean ± S.D. from three independent experiments. Student’s t-test is used to determine statistical significance (∗p < 0.05). (F) Sections are co-immunostained with anti-Satb2 (green) anti-Tbr1 (red). (I) Distribution of migrating neuronal cells labeled with BrdU (green) in cerebral cortex. Nuclei are counterstained with DAPI (blue). Representative images from three independent experiments are shown. The cell densities of Satb2 (G), Tbr1 (H), and BrdU-positive cells (J) in each area are plotted as histograms. Values are shown as the mean ± SD. The asterisk indicates a significant difference between +/+ and lde/lde in the same area by Student’s t-test (∗p < 0.05, ∗∗p < 0.01). Experimental scheme to label migrating cells by BrdU is illustrated in panel (J). MZ, marginal zone; IZ, intermediate zone; VZ, ventricular zone: SVZ, subventricular zone.
Wwox-Depleted Neurons Show Reduced Expression of Different Neural Migration-Related Genes
Transcriptomic analysis of WWOX-silenced human neural progenitors was performed before and after differentiation into mature neurons and focused on genes whose expression was significantly altered after WWOX depletion (Log2FC < >1). Within this dataset, we identified 20 neural migration-related genes of microtubule architecture and function with the high fold change (>2) of expression (Supplementary Table S1 and Supplementary Figure S2). The relative rank of the genes is presented in volcano plots (Figure 4). Among these, we found tubulin TUBA1A gene expression lowered more than 10 times in Wwox-depleted neurons compared to control neurons. Immunostaining in the WWOX-deficient human brain revealed a slight reduction of tubulin-1-alpha protein, although not significant (Supplementary Figure S3). Other tubulin-related genes which expression was found to be dramatically reduced include TUBB3 (>15-fold reduced) and TUBB2B (>3-fold reduced expression). Mutations of all three above genes are implicated in a wide array of different human neuronal migration disorders associated with impaired heterodimerization of microtubules. In addition, we found three kinesin-1 heavy chain genes (KIF5A, KIF5B, and KIF5C) showing a 2–3-fold reduced gene expression in Wwox-depleted neurons compared to normal neurons. Also, the MAP2 gene (encoding a protein involved in stabilizing microtubules by crosslinking each over and with intermediate filaments) expression was found to be 3-fold reduced in Wwox-depleted neurons.
Figure 4. Wwox-depleted neurons show reduced expression of different neural migration-related genes. The differential gene expression profiles in control hNPC vs. WWOX-depleted hNPC and control neurons vs. WWOX-depleted neurons represented in Volcano Plot. P – significance was calculated according to NOISeq-sim algorithm.
Discussion
The WWOX gene is an important transcriptional regulator, known to modulate the activity of a number of transcriptional factors and molecular pathways (including Jun, AP2gamma, NFkappaB, and ErbB4) and implicated in tumor suppressor activity, steroid-hormone metabolism regulation and brain development (Aldaz et al., 2014; Iatan et al., 2014). The gene encodes Wwox, an oxidoreductase that contains two distinct domains, (a) an amino-terminus showing high sequence conservation to the WW domain family of proteins and displaying two WW domains, and (b) a central 283 amino acids domain homologous to the Short-Chain Dehydrogenase/Reductase (SDR) superfamily (Figure 1G). Wwox shows highest expression in thyroid, testis ovary, brain and cerebellum (Chen et al., 2004; Nunez et al., 2006; Chang et al., 2007). WWOX pathogenic biallelic variants are implicated in a spectrum of neurological phenotypes, ranging between autosomal recessive spinocerebellar ataxia type 12 (SCAR12; MIM 614322) to the most severe WOREE syndrome (Gribaa et al., 2007; Abdel-Salam et al., 2014; Mignot et al., 2015). The few SCAR12 individuals described so far harbor two homozygous hotspot missense variants (Pro47Thr and Gly372Arg; NM_016373.4). According to functional studies, these variants cause only partial loss-of-function due to conformational changes that alters Wwox WW1/SDR domains ability to interact with other proteins (Mallaret et al., 2014). WOREE individuals exhibiting absent developmental milestones and profound neurological impairment usually harbor more severe WWOX recessive variants, including nonsense and frameshift variants (Mignot et al., 2015), deletions often involving multiple exons of the gene (Piard et al., 2018), and non-synonymous variants sometimes resulting in full knockout of the protein via impaired WWOX translation (Johannsen et al., 2018). In our family, we identified a homozygous p.Arg264Ter variant, resulting in a loss of normal Wwox function either through protein truncation (and disruption of the active C-terminal short-chain dehydrogenase/reductase SDR domain), or nonsense-mediated mRNA decay. Clinical phenotype of the index case (II.1, Figure 1A) was consistent with other individuals reported in the literature and included postnatal microcephaly, absent developmental milestones, visual impairment, tetraparesis, dystonic movements and drug-resistant polymorphic seizures. Moreover, mild cerebral atrophy, white matter signal alterations, and hypoplasia of the cerebellar vermis and corpus callosum were also noted. Interestingly, prenatal features such as reduced fetal movements or congenital clubfoot (occurred in individual II.1) have been occasionally described as part of the WOREE syndrome spectrum (Abdel-Salam et al., 2014; Ben-Salem et al., 2015). A large study on 20 WOREE individuals suggests that antenatal signs, including increased nuchal translucency and hands or feet anomalies, may occur in up to 25% of cases (Piard et al., 2018). Our index case (II.1) was also born with cryptorchidism and, to the best of our knowledge, this feature was not described before as part of the WOREE spectrum. However, testicular atrophy, Leydig cell dysfunction and heterogeneous gonadal abnormalities are usually observed in animal models of Wwox deficiency (Ludes-Meyers et al., 2007; Suzuki et al., 2007). These observations can be explained by the important role of the protein in gonadal development, being its SRD enzymatic domain likely involved in the regulation of sex-steroid metabolism in testis and ovaries (Nunez et al., 2006; Aqeilan et al., 2007, 2009). Notably, WOREE-associated prenatal brain anomalies have been described in 2015 in one fetus carrying a 0.6 homozygous 16q23.1 deletion encompassing WWOX (Valduga et al., 2015). In this case, autopsy identified thin corpus callosum and polymicrogyria, but confirmatory brain histopathology studies were not performed. The observation of prenatal brain malformations in WOREE syndrome is consistent with animal studies showing abundant WWOX expression in the CNS of mouse embryos, suggesting an important role of the gene in CNS development (Chen et al., 2004). Despite the above studies, disease mechanisms by which loss of Wwox leads to WOREE-associated CNS alterations have not been fully understood yet. Different hypotheses have been proposed, including abnormal tau phosphorylation, mitochondrial dysfunction, altered cell differentiation, and neuronal apoptosis (Ramos and Aldaz, 2006; Aldaz et al., 2014; Liu et al., 2018). To investigate the impact of Wwox deficiency on human brain developmental processes, we performed brain histopathological analyses of the 21th week aborted fetus (II.3) with WOREE syndrome and identified abnormal architecture of the developing cerebral cortex. Development of the cerebral cortex is a very dynamic and highly organized process, involving a series of complex morphogenetic events, which can be impaired in different neurodevelopmental disorders (Fernández-Marmiesse et al., 2019; Seto and Eiraku, 2019; Magdalon et al., 2020). Following division of progenitor cells in the ventricular zone, neurons undergo several physiological changes and migrate outward toward the cortical plate, where they differentiate and integrate into functional neural circuits (Romero et al., 2018). These events are tightly regulated by complex intra- and extra- cellular molecular pathways, including Notch, Wnt, and Hedgehog signaling (Kuzniecky, 2015; Jiang and Nardelli, 2016; Parrini et al., 2016). Importantly, migration of neurons to their final destination sites also requires changes in cell polarity and motility. A pivotal role in modulating these processes is played by cytoskeleton components such as actin filaments, microtubules, and motor proteins (including dyneins and kinesins) which facilitate migrating neurons directional movements (Budday et al., 2015; Jiang and Nardelli, 2016). In human brain fetuses, migration processes begin between the 7th and 8th week of gestation, when neuroblasts migrate along radial extensions of glial cells, to detach when they reach the proper layer of the developing cerebral cortex (Mora et al., 2007). In the present study, HE staining of the WWOX-deficient brain revealed abnormal arrangement of cell layers with defective migration of the cerebral external granular layer into the molecular layer, that appeared not homogeneous. In addition, cerebral cortex PAS-stain revealed irregular migration trajectories of glial cells, known to be important in the guidance of migrating neurons. Interestingly, it has been previously showed by some of us (KK, AKB) that Wwox silencing in neural progenitor cells may alter the expression of a large number of genes (e.g., DCLK, NEFM, and NEFL) involved in neurofilament assembly, cytoskeleton organization and chromatin remodeling (Kośla et al., 2019). In the present study, we performed detailed analysis of gene expression of Wwox-depleted neural progenitors and differentiated post-mitotic neurons. This showed a significantly reduced expression of several tubulins, microtubule- associated proteins and microtubule motor kinesins (Figure 4 and Supplementary Table S1). Mutations in genes encoding tubulins may affect neural cell differentiation and neuron migration, resulting in a wide array of brain developmental defects, including lissencephaly-pachygyria and polymicrogyria (Kumar et al., 2010; Poirier et al., 2013; Belvindrah et al., 2017; Juric-Sekhar and Hevner, 2019). These results suggest that abnormal neuronal migration and defective cerebral cortex development can occur as a result of Wwox deficiency. Previously, the clinical features of the lde/lde rats, characterized by dwarfism, postnatal lethality, epilepsy and ataxic gait, were related to a spontaneous homozygous 13-bp deletion within the exon 9 of WWOX, resulting in a premature truncation of the Wwox C-terminus (Suzuki et al., 2009). Wide neuropathological alterations have been described in the postnatal cerebral cortex of lde/lde rats, including severe hypomyelination (with a reduced number of mature oligodendrocytes) and a significant reduction in cell populations of astrocytes and microglia (Tochigi et al., 2019). In the present study, histological and immunohistological analyses of developing cerebral cortex of lde/lde rats investigated at perinatal age revealed a defective migration of immature neuron and delayed formation of cortical layers, overlapping the identified human pathological features. In conclusion, our study indicates that loss of Wwox could alter a number of transcriptionally regulated genes implicated in cytoskeleton organization and microtubule assembly, leading to cellular and sub-cellular defects of cerebral cortication. These findings highlight a conserved role of WWOX in regulating neuronal migration and CNS development across different species. Experimental studies using WOREE individuals-derived forebrain organoids, which can model critical human fetal developmental processes, could help to understand the various phases of cortical development regulated by WWOX, and to assess the different molecular pathways implicated in the pathological features of this severe brain developmental disorder.
Data Availability Statement
The datasets generated for this study can be found in the GEO database (GSE126075) and the ClinVar database (VCV000241105).
Ethics Statement
The studies involving human participants were reviewed and approved by “Giannina Gaslini” Children’s Hospital Institutional Review Board. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. The animal study was reviewed and approved by Animal Care and Use Committee of Nippon Veterinary and Life Science University. Written informed consent was obtained from the individual(s), and minor(s)’ legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author Contributions
MI: acquisition and interpretation of genetic data, and drafting of the manuscript. AA and SB: immunohistochemistry, immunofluorescence, interpretation data, and drafting of the manuscript. YT: animal model studies, interpretation of data, and drafting of the manuscript. KK and AB: transcriptomic analyses. FB and EF: histopathological studies. AT and VN: NGS analysis. MSe and DP: brain magnetic resonance imaging. LM, AR, and MSc: clinical phenotyping, assessment, and follow-up. GB: sanger sequencing validation and segregation studies. PS and CM: critical revision of manuscript and obtaining funding. FZ: interpretation of genetic data, reviewing of the manuscript, critical revision of manuscript, and obtaining funding. HS: supervision of animal model studies and reviewing of the manuscript. VS: design of the study, final acquisition, interpretation of clinical data, and drafting and reviewing the manuscript. All authors contributed to the manuscript and approved the submitted version.
Funding
This work was developed within the framework of the DINOGMI Department of Excellence of MIUR 2018-2022 (law 232, year 2016) and was supported by the National Institute for Health Research University College London Hospitals Biomedical Research Centre and National Science Centre, Poland, grant no. 2015/17/D/NZ2/01989. This work has been supported by a grant of the Italian Ministry of Health (RF-2016-02361949 to FZ).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We acknowledge the support from Enea-Eroi in ricerca onlus (www.eneaonlus.org). We also gratefully acknowledge the family for their participation to this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2020.00644/full#supplementary-material
Footnotes
References
Abdel-Salam, G., Thoenes, M., Afifi, H. H., Korber, F., Swan, D., and Bolz, H. J. (2014). The supposed tumor suppressor gene WWOX is mutated in an early lethal microcephaly syndrome with epilepsy, growth retardation and retinal degeneration. Orphanet J. Rare Dis. 9:12. doi: 10.1186/1750-1172-9-12
Aldaz, C. M., Ferguson, B. W., and Abba, M. C. (2014). WWOX at the crossroads of cancer, metabolic syndrome related traits and CNS pathologies. Biochim. Biophys. Acta. 1846, 188–200. doi: 10.1016/j.bbcan.2014.06.001
Aqeilan, R. I., Hagan, J. P., de Bruin, A., Rawahneh, M., Salah, Z., Gaudio, E., et al. (2009). Targeted ablation of the WW domain-containing oxidoreductase tumor suppressor leads to impaired steroidogenesis. Endocrinology. 150, 1530–1535. doi: 10.1210/en.2008-1087
Aqeilan, R. I., Trapasso, F., Hussain, S., Costinean, S., Marshall, D., Pekarsky, Y., et al. (2007). Targeted deletion of Wwox reveals a tumor suppressor function. PNAS 104, 3949–3954. doi: 10.1073/pnas.0609783104
Arisio, R., Bonissone, M., Piccoli, E., and Panzica, G. (2002). Central nervous system microangioarchitecture in the human foetus. Adv. Clin. Path 6, 125–129.
Bär, T. (1980). The Vascular System of the Cerebral Cortex: Advances in Anatomy, Embryology and Cell Biology, Vol. 59. Berlin: Springer.
Bednarek, A. K., Keck-Waggoner, C. L., Daniel, R. L., Laflin, K. J., Bergsagel, P. L., Kiguchi, K., et al. (2001). WWOX, the FRA16D gene, behaves as a suppressor of tumor growth. Cancer Res. 61, 8068–8073.
Belvindrah, R., Natarajan, K., Shabajee, P., Bruel-Jungerman, E., Bernard, J., Goutierre, M., et al. (2017). Mutation of the α-tubulin Tuba1a leads to straighter microtubules and perturbs neuronal migration. J. Cell Biol. 216, 2443–2461. doi: 10.1083/jcb.201607074
Ben-Salem, S., Al-Shamsi, A. M., John, A., Ali, B. R., and Al-Gazali, L. (2015). A novel whole exon deletion in WWOX gene causes early epilepsy, intellectual disability and optic atrophy. J. Mol. Neurosci. 56, 17–23. doi: 10.1007/s12031-014-0463-8
Britanova, O., Akopov, S., Lukyanov, S., Gruss, P., and Tarabykin, V. (2005). Novel transcription factor Satb2 interacts with matrix attachment region DNA elements in a tissue−specific manner and demonstrates cell−type−dependent expression in the developing mouse CNS. Eur. J. Neurosci. 21, 658–668. doi: 10.1111/j.1460-9568.2005.03897.x
Budday, S., Steinmann, P., and Kuhl, E. (2015). Physical biology of human brain development. Front. Cell Neurosci. 9:257. doi: 10.3389/fncel.2015.00257
Chang, N. S., Hsu, L. J., Lin, Y. S., Lai, F. J., and Sheu, H. M. (2007). WW domain-containing oxidoreductase: a candidate tumor suppressor. Trends Mol. Med. 13, 12–22. doi: 10.1016/j.molmed.2006.11.006
Chen, S. T., Chuang, J. I., Wang, J. P., Tsai, M. S., Li, H., and Chang, N. S. (2004). Expression of WW domain-containing oxidoreductase WOX1 in the developing murine nervous system. Neuroscience 124, 831–839. doi: 10.1016/j.neuroscience.2003.12.036
Fan, X., Kim, H. J., Bouton, D., Warner, M., and Gustafsson, J. A. (2008). Expression of liver X receptor beta is essential for formation of superficial cortical layers and migration of later-born neurons. Proc. Natl. Acad. Sci. U.S.A. 105, 13445–13450. doi: 10.1073/pnas.0806974105
Fernández-Marmiesse, A., Roca, I., Díaz-Flores, F., Cantarín, V., Pérez-Poyato, M. S., Fontalba, A., et al. (2019). Rare variants in 48 genes account for 42% of cases of epilepsy with or without neurodevelopmental delay in 246 pediatric patients. Front. Neurosci. 13:1135. doi: 10.3389/fnins.2019.01135
Gribaa, M., Salih, M., Anheim, M., Lagier-Tourenne, C., H’mida, D., Drouot, N., et al. (2007). A new form of childhood onset, autosomal recessive spinocerebellar ataxia and epilepsy is localized at 16q21-q23. Brain 130, 1921–1928. doi: 10.1093/brain/awm078
Hatanaka, Y., Hisanaga, S., Heizmann, C. W., and Murakami, F. (2004). Distinct migratory behavior of early- and late-born neurons derived from the cortical ventricular zone. J. Comp. Neurol. 479, 1–14. doi: 10.1002/cne.20256
Hevner, R. F., Shi, L., Justice, N., Hsueh, Y., Sheng, M., Smiga, S., et al. (2001). Tbr1 regulates differentiation of the preplate and layer 6. Neuron 29, 353–366. doi: 10.1016/s0896-6273(01)00211-2
Iatan, I., Choi, H. Y., Ruel, I., Reddy, M. V., Kil, H., Lee, J., et al. (2014). The WWOX gene modulates high-density lipoprotein and lipid metabolism. Circ. Cardiovasc. Genet. 7, 491–504. doi: 10.1161/CIRCGENETICS.113.000248
Jansen, A. C., Robitaille, Y., Honavar, M., Mullatti, N., Leventer, R. J., Andermann, E., et al. (2016). The histopathology of polymicrogyria: a series of 71 brain autopsy studies. Dev. Med. Child. Neurol. 58, 39–48. doi: 10.1111/dmcn.12840
Jiang, X., and Nardelli, J. (2016). Cellular and molecular introduction to brain development. Neurobiol. Dis. 92, 3–17. doi: 10.1016/j.nbd.2015.07.007
Johannsen, J., Kortüm, F., Rosenberger, G., Bokelmann, K., Schirmer, M. A., Denecke, J., et al. (2018). A novel missense variant in the SDR domain of the WWOX gene leads to complete loss of WWOX protein with early-onset epileptic encephalopathy and severe developmental delay. Neurogenetics 19, 151–156. doi: 10.1007/s10048-018-0549-5
Juric-Sekhar, G., and Hevner, R. F. (2019). Malformations of cerebral cortex development: molecules and mechanisms. Annu. Rev. Pathol. 14, 293–318. doi: 10.1146/annurev-pathmechdis-012418-012927
Kośla, K., Płuciennik, E., Styczeñ-Binkowska, E., Nowakowska, M., Orzechowska, M., and Bednarek, A. K. (2019). The WWOX gene influences cellular pathways in the neuronal differentiation of human neural progenitor cells. Front. Cell Neurosci. 13:391. doi: 10.3389/fncel.2019.00391
Kumar, R. A., Pilz, D. T., Babatz, T. D., Cushion, T. D., Harvey, K., Topf, M., et al. (2010). TUBA1A mutations cause wide spectrum lissencephaly (smooth brain) and suggest that multiple neuronal migration pathways converge on alpha tubulins. Hum. Mol. Genet. 19, 2817–2827. doi: 10.1093/hmg/ddq182
Kuzniecky, R. (2015). Epilepsy and malformations of cortical development: new developments. Curr. Opin. Neurol. 28, 151–157. doi: 10.1097/WCO.0000000000000175
Liu, C. C., Ho, P. C., Lee, I. T., Chen, Y. A., Chu, C. H., Teng, C. C., et al. (2018). WWOX phosphorylation, signaling, and role in neurodegeneration. Front. Neurosci. 12:563. doi: 10.3389/fnins.2018.00563
Ludes-Meyers, J. H., Kil, H., Nuñez, M. I., Conti, C. J., Parker-Thornburg, J., Bedford, M. T., et al. (2007). WWOX hypomorphic mice display a higher incidence of B-cell lymphomas and develop testicular atrophy. Genes Chromosomes Cancer 46, 1129–1136. doi: 10.1002/gcc.20497
Magdalon, J., Mansur, F., Teles, E., Silva, A. L., de Goes, V. A., Reiner, O., et al. (2020). Complement system in brain architecture and neurodevelopmental disorders. Front. Neurosci. 14:23. doi: 10.3389/fnins.2020.00023
Mallaret, M., Synofzik, M., Lee, J., Sagum, C. A., Mahajnah, M., Sharkia, R., et al. (2014). The tumour suppressor gene WWOX is mutated in autosomal recessive cerebellar ataxia with epilepsy and mental retardation. Brain 137, 411–419. doi: 10.1093/brain/awt338
Mignot, C., Lambert, L., Pasquier, L., Bienvenu, T., Delahaye-Duriez, A., Keren, B., et al. (2015). WWOX-related encephalopathies: delineation of the phenotypical spectrum and emerging genotype phenotype correlation. J. Med. Genet. 52, 61–70. doi: 10.1136/jmedgenet-2014-102748
Mora, M., Gaggero, G., Scandale, A., and Fulcheri, E. (2007). Modified Kluver-Barrera staining for the study and diagnosis of fetal encephalopathies. Ital. J. Anat. Embryol. 112, 111–116.
Nunez, M. I., Ludes-Meyers, J., and Aldaz, C. M. (2006). WWOX protein expression in normal human tissues. J. Mol. Histol. 37, 115–125. doi: 10.1007/s10735-006-9046-5
Parrini, E., Conti, V., Dobyns, W. B., and Guerrini, R. (2016). Genetic basis of brain malformations. Mol. Syndromol. 7, 220–233. doi: 10.1159/000448639
Piard, J., Hawkes, L., Milh, M., Villard, L., Borgatti, R., Romaniello, R., et al. (2018). The phenotypic spectrum of WWOX-related disorders: 20 additional cases of WOREE syndrome and review of the literature. Genet. Med. 21, 1308–1318. doi: 10.1038/s41436-019-0460-y
Poirier, K., Saillour, Y., Fourniol, F., Francis, F., Souville, I., Valence, S., et al. (2013). Expanding the spectrum of TUBA1A-related cortical dysgenesis to Polymicrogyria. Eur. J. Hum. Genet. 21, 381–385. doi: 10.1038/ejhg.2012.195
Ramos, D., and Aldaz, C. M. (2006). WWOX, a chromosomal fragile site gene and its role in cancer. Adv. Exp. Med. Biol. 587, 149–159. doi: 10.1007/978-1-4020-5133-3_14
Romero, D. M., Bahi-Buisson, N., and Francis, F. (2018). Genetics and mechanisms leading to human cortical malformations. Semin. Cell Dev. Biol. 76, 33–75. doi: 10.1016/j.semcdb.2017.09.031
Salpietro, V., Dixon, C. L., Guo, H., Bello, O. D., Vandrovcova, J., Efthymiou, S., et al. (2019). AMPA receptor GluA2 subunit defects are a cause of neurodevelopmental disorders. Nat. Commun. 10:3094. doi: 10.1038/s41467-019-10910-w
Seto, Y., and Eiraku, M. (2019). Human brain development and its in vitro recapitulation. Neurosci. Res. 138, 33–42. doi: 10.1016/j.neures.2018.09.011
Suzuki, H., Katayama, K., Takenaka, M., Amakasu, K., Saito, K., and Suzuki, K. (2009). A spontaneous mutation of the WWOX gene and audiogenic seizures in rats with lethal dwarfism and epilepsy. Genes Brain Behav. 8, 650–660. doi: 10.1111/j.1601-183X.2009.00502.x
Suzuki, H., Takenaka, M., and Suzuki, K. (2007). Phenotypic characterization of spontaneously mutated rats showing lethal dwarfism and epilepsy. Comp. Med. 57, 360–369.
Tabarki, B., Al Mutairi, F., and Al Hashem, A. (2015a). The fragile site WWOX gene and the developing brain. Exp. Biol. Med. (Maywood) 240, 400–402. doi: 10.1177/1535370214561952
Tabarki, B., AlHashem, A., AlShahwan, S., Alkuraya, F. S., Gedela, S., and Zuccoli, G. (2015b). Severe CNS involvement in WWOX mutations: description of five new cases. Am. J. Med. Genet. A 167A, 3209–3213. doi: 10.1002/ajmg.a.37363
Tanna, M., and Aqeilan, R. I. (2018). Modeling WWOX Loss of Function in vivo: what have we learned? Front. Oncol. 8:420. doi: 10.3389/fonc.2018.00420
Tarazona, S., Furió-Tarí, P., Turrà, D., Pietro, A. D., Nueda, M. J., Ferrer, A., et al. (2015). Data quality aware analysis of differential expression in RNA-seq with NOISeq R/Bioc package. Nucleic Acids Res. 43:e140. doi: 10.1093/nar/gkv711
Tochigi, Y., Iwasaki, Y., Sano, M., Yasuda, H., Katayama, K., and Suzuki, H. (2017). Critical roles of Astrin in the mitosis of immature rat Sertoli cells. Biochem. Biophys. Res. Commun. 486, 958–964. doi: 10.1016/j.bbrc.2017.03.137
Tochigi, Y., Takamatsu, Y., Nakane, J., Nakai, R., Katayama, K., and Suzuki, H. (2019). Loss of WWOX causes defective development of cerebral cortex with hypomyelination in a rat model of lethal dwarfism with epilepsy. Int. J. Mol. Sci. 20:3596. doi: 10.3390/ijms20143596
Keywords: WWOX, WOREE syndrome, neuropathology, animal model, developing brain, cytoskeleton
Citation: Iacomino M, Baldassari S, Tochigi Y, Kośla K, Buffelli F, Torella A, Severino M, Paladini D, Mandarà L, Riva A, Scala M, Balagura G, Accogli A, Nigro V, Minetti C, Fulcheri E, Zara F, Bednarek AK, Striano P, Suzuki H and Salpietro V (2020) Loss of Wwox Perturbs Neuronal Migration and Impairs Early Cortical Development. Front. Neurosci. 14:644. doi: 10.3389/fnins.2020.00644
Received: 05 March 2020; Accepted: 25 May 2020;
Published: 11 June 2020.
Edited by:
David F. Clayton, Queen Mary University of London, United KingdomReviewed by:
Martin W. Breuss, University of California, San Diego, United StatesMichael John Gambello, Emory University, United States
Copyright © 2020 Iacomino, Baldassari, Tochigi, Kośla, Buffelli, Torella, Severino, Paladini, Mandarà, Riva, Scala, Balagura, Accogli, Nigro, Minetti, Fulcheri, Zara, Bednarek, Striano, Suzuki and Salpietro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vincenzo Salpietro, di5zYWxwaWV0cm9AdWNsLmFjLnVr
†These authors have contributed equally to this work