Martina Servetti1,2†
Livia Pisciotta1,3†
Elisa Tassano4
Maria Cerminara1
Lino Nobili1,5
Silvia Boeri1,5Giulia Rosti1Margherita Lerone2
Maria Teresa Divizia2Patrizia Ronchetto4
Aldamaria Puliti1,2*- 1Department of Neurosciences, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DiNOGMI), University of Genoa, Genoa, Italy
- 2Medical Genetics Unit, IRCCS Istituto Giannina Gaslini, Genoa, Italy
- 3Child Neuropsychiatry Unit, ASST Fatebenefratelli Sacco, Milano, Italy
- 4Human Genetics Laboratory, IRCCS Istituto Giannina Gaslini, Genoa, Italy
- 5Child Neuropsychiatry Unit, Istituto Giannina Gaslini, Genoa, Italy
Neurodevelopmental disorders (NDDs) are a heterogeneous class of brain diseases, with a complex genetic basis estimated to account for up to 50% of cases. Nevertheless, genetic diagnostic yield is about 20%. Array-comparative genomic hybridization (array-CGH) is an established first-level diagnostic test able to detect pathogenic copy number variants (CNVs), however, most identified variants remain of uncertain significance (VUS). Failure of interpretation of VUSs may depend on various factors, including complexity of clinical phenotypes and inconsistency of genotype-phenotype correlations. Indeed, although most NDD-associated CNVs are de novo, transmission from unaffected parents to affected children of CNVs with high risk for NDDs has been observed. Moreover, variability of genetic components overlapped by CNVs, such as long non-coding genes, genomic regions with long-range effects, and additive effects of multiple CNVs can make CNV interpretation challenging. We report on 12 patients with complex phenotypes possibly explained by complex genetic mechanisms, including involvement of antisense genes and boundaries of topologically associating domains. Eight among the 12 patients carried two CNVs, either de novo or inherited, respectively, by each of their healthy parents, that could additively contribute to the patients’ phenotype. CNVs overlapped either known NDD-associated or novel candidate genes (PTPRD, BUD13, GLRA3, MIR4465, ABHD4, and WSCD2). Bioinformatic enrichment analyses showed that genes overlapped by the co-occurring CNVs have synergistic roles in biological processes fundamental in neurodevelopment. Double CNVs could concur in producing deleterious effects, according to a two-hit model, thus explaining the patients’ phenotypes and the incomplete penetrance, and variable expressivity, associated with the single variants. Overall, our findings could contribute to the knowledge on clinical and genetic diagnosis of complex forms of NDD.
Introduction
Neurodevelopmental disorders (NDDs) are a heterogeneous class of conditions involving the brain, including intellectual disability and autism spectrum disorder (ASD), that affect about 1–3% of children (Miller et al., 2010). The genetics of NDDs is complex and include copy number variants (CNVs) and deleterious variants in single genes (Coe et al., 2019).
Nowadays new technologies, such as whole exome (WES) and whole genome sequencing (WGS), have become more efficient and inexpensive such that it is ongoing a debate for choosing the proper first-tier diagnostic test (Srivastava et al., 2019; Alvarez-Mora et al., 2021; Arteche-Lopez et al., 2021).
Despite of that, the detection of CNVs by the array comparative genomic hybridization (array-CGH) continues to be widely used as a first-tier test, resulting in an average diagnostic yield between 15 and 20% (Miller et al., 2010). Interpretation of CNVs is complex and, despite the use of standard classification guides (Riggs et al., 2020), most identified variants remain of uncertain significance (VUS).
Failure of interpretation of VUSs may depend on different factors, including poor annotation of protein-coding genes, lack of functional information for untranslated expressed genes or for genomic intergenic regions overlapped by the CNVs. The observation of inconsistent genotype-phenotype correlations also complicates CNV interpretation. Although most NDD causative CNVs are de novo, increasing evidence indicates that CNVs with high risk for NDDs show incomplete penetrance since transmission from unaffected parents to affected children has been observed. Possible explanation for the observed inconsistency of genotype-phenotype correlation may involve complex interactions of potentially pathogenic CNVs with additional/secondary CNVs, or single nucleotide variants, that could act together to determine a disease state (Velinov, 2019). According to this hypothesis, additive effects of CNVs were for instance reported in probands carrying 16p11.2 deletions and secondary CNVs that disrupt genes associated with autism and/or intellectual disability (Girirajan et al., 2010; Duyzend et al., 2016).
The goal of this study was to clarify the contribution of microdeletions and microduplications identified, by array-CGH diagnostic test, in 12 patients with complex NDD phenotype using multiple bioinformatics resources and the synergistic effort of a team of neuropsychiatrists and geneticists.
We found potential pathogenic variants in novel candidate genes for neurodevelopmental disorders. Furthermore, CNV-mediated double-hit mechanisms seem to play a relevant role in the complex forms of neurodevelopmental disorders affecting our patients. Indeed, these co-occurring hits involved already known NDD-associated genes or hypothetical novel NDD candidate genes enriched for pathways or biological processes known to be implicated in neurodevelopment.
Methods
Clinical Assessment
The clinical assessment of the patients comprised a thorough medical history also attentive to investigating systemic symptoms and sleep, neurological-behaviour examination, developmental/cognitive assessment, and basic metabolic screening. Prior to array-CGH analysis, all patients were first evaluated by using different genetic tests for the common genetic causes of NDDs, including karyotype analysis, and tests for Fragile X syndrome. Written informed consent was obtained from the patient’s parents or legal representative. This study was reviewed and approved by the Ethics Committee of the Italian Regione Liguria (R. P. 001/2019).
CNV Detection and Annotation
Array-CGH analysis was performed on DNA samples, extracted from peripheral blood, using a whole-genome 180 K Agilent array with ∼13 Kb overall median probe spacing (Human Genome CGH Microarray, Agilent Technologies, Santa Clara, CA, USA). Data were analyzed using Agilent CytoGenomics and genomic positions reported according to the human genome assembly (GRCh37/hg19). All detected CNVs were tested for inheritance by hybridization of the parental DNA with the same array platform. De novo CNVs were confirmed either by FISH analysis or by quantitative real-time PCR as already reported (Vaccari et al., 2014; Tassano et al., 2015). To assess the clinical significance of the detected CNVs, we used several tools. Aberration segments were reviewed using GRCh37 hg19 of UCSC Genome Browser (http://genome.ucsc.edu/index.html). We annotated all detected CNVs and CNV-encompassed genes across public databases: Genomic Variants Database (DGV) (http://dgv.tcag.ca/dgv/app/home), DECIPHER (https://decipher.sanger.ac.uk/), Clinical Genome Resource (ClinGen) (https://clinicalgenome.org/), the Online Mendelian Inheritance in Man (OMIM) (http://www.omim.org), Simons Foundation Autism Research Initiative (SFARI) database (https://gene.sfari.org/), databases of mouse models (Mouse Genome Informatics, MGI) (http://www.informatics.jax.org). Literature mining was also performed. We considered only variants present in <0.1% of control individuals. We checked CNVs for the presence of dosage-sensitive genes (haploinsufficiency/triple sensitivity score) by using Clinical Genome Resource (ClinGen) consortium resource (https://dosage.clinicalgenome.org/). We also considered probability of intolerance to loss-of-function (LoF) variants (pLI) as reported in the Genome Aggregation Database (gnomAD, (https://gnomad.broadinstitute.org/). Genes with pLI scores of 0.9 or higher are extremely intolerant to heterozygous LoF variation, and thus haploinsufficient (Lek et al., 2016).
We then investigated deleted/duplicated regions for encompassing topologically associating domains (TADs), or their boundaries (TDBs), using Hi-C (http://3dgenome.fsm.northwestern.edu/). Indeed, TADs data from genome-wide higher order chromatin interaction data in human embryonic stem cells (h-ESC), h-ESC derived neural progenitor cells (H1NPC) (Dixon et al., 2015), and brain cortex (Schmitt et al., 2016), were downloaded (http://chromosome.sdsc.edu/mouse/download.html) (see Supplementary file) and mapped to hg19 coordinates using the UCSC browser. We then investigated identified CNV-genes, not previously associated with genetic diseases, for relevance to NDDs by considering different criteria: expression in the brain; reported mutations in the mouse associated with a neurodevelopmental phenotype; interaction with genes known to be associated with a neurodevelopmental disorder. To identify CNV disrupted genes with high brain expression, we used the BrainSpan database (http://www.brainspan.org/static/download.html) (Huckins et al., 2019) and extracted the list of genes with average log2 RPKM >4.5 (the top 18%), as also reported (Addis et al., 2018; Cerminara et al., 2021). We further investigated brain expression of genes at single gene isoforms using GTEx data (https://www.gtexportal.org/home/). We searched for gene-specific mouse models presenting a neurodevelopmental phenotype using MGI database. CNV-genes coding for microRNAs (miRNAs) were further investigated to detect the presence, among their possible targets, of any genes known to be associated with NDDs (Supplementary file) or to be involved in brain development and/or function by using miRDB (http://mirdb.org/), Target Scan (http://www.targetscan.org/vert_72/) and DIANA TOOLS (http://diana.imis.athena-innovation.gr/DianaTools/index.php).
Network Analysis of Candidate Genes
We used GeneCodis4 (Tabas-Madrid et al., 2012) to unveil enrichment of annotations, as already reported (Cerminara et al., 2021). Briefly, all the genes obtained from either duplicated or deleted regions were used as input in GeneCodis4 together with known NDD-associated genes, as those related to ASD and reported in SFARI database and those reported in OMIM as associated with intellectual disability (see Supplementary file). This tool allows the classification of genes according to their putative biological function by screening the Gene Ontology (GO), OMIM, Panther, and KEGG Pathways. In the analysis, the hypergeometric test was applied followed by the false discovery rate correction (FDR) with a cut-off of 5% to determine which annotations were significantly enriched. For GO analysis, various hierarchical levels of the annotation data structure were used. All genes of the top gene ontology terms, including candidate and known NDD genes, were then projected onto the STRING network (v11) (Szklarczyk et al., 2019). Edges within the STRING network were thresholded at 0.4, according to the authors’ recommendation. A graphical representation of the GeneCodis4 and STRING results was obtained by using Cytoscape tool (Shannon et al., 2003).
Pathogenicity Evaluation and Graphical Representation of CNVs
To evaluate the pathogenicity of variants we focused on the following criteria: 1) patients with multiple CNVs, enclosing NDD genes, which could have additive effects on the patients’ phenotype on the basis of results from gene enrichment and protein-protein interactions analyses; 2) patients with CNVs encompassing regions that control the expression of known NDD genes, including those with potential long-range effects; 3) patients whose CNVs encompass genes not yet reported as NDD-associated, but whose expression profile and function suggest a possible involvement in the disease, when at least another patient is reported with a similar CNV and phenotype in publically available databases. Results of this analysis were graphically reported using UCSC and specific Custom Tracks. In particular we used RefSeq Curated for gene representation and OMIM genes to underline genes already associated with disease. Decipher track was manually modified; we selected patients with CNVs similar to those of our patients, that is either deletions or duplications according to those found in our patients and overlapping approximately the same region and genes. We selected Decipher CNVs only if their related patients had a NDD phenotype (we kept out patients for whom we had no information about their phenotype). The CNVs thus selected were further investigated for being present in the same patient with other secondary CNVs overlapping additional NDD genes, and in this case reported in figures (marked with letters). Developmental Delay case track includes manually selected CNVs similar in size and gene content to those of our cases; we proceeded analogously also for Developmental Delay Control CNVs. A track reporting SFARI genes was added to underline the presence of a known ASD associated gene in the genomic region shown in figures. Brain-expressed gene track represents genes with a prevalent expression in the brain according to the BrainSpan data processed as detailed above, useful to highlight possible new NDD candidate genes.
Results
Clinical Presentation of Patients
The 12 patients we report were all born to non-consanguineous parents. After clinical evaluation and follow up, the employment of specific tests and instruments, all subjects showed a complex phenotype, each of them characterized by different neurological features and other associated impairments and comorbidities. The main clinical features were represented by a variable degree of intellectual disability, in three cases associated with a diagnosed ASD, and epilepsy. Language disorder, self-aggressive behaviours and sleep disturbances could be also observed (see Table 1).
TABLE 1. Clinical features.
Additive Effects of CNVs Involving Known NDD Genes
Four patients presented with two CNVs, overlapping known NDD genes, and inherited, respectively, by each of their healthy parents. We thought that in these cases, CNVs could additively contribute to the patient’s phenotype (Table 2).
TABLE 2. Overview of the CNVs found in the patients.
The first case, IGGAC13, is a male patient with a maternal deletion of chromosome 7 involving CNTNAP2 and a paternal deletion of chromosome 11 involving LRRC4C (Figures 1A,B). Both genes are present in the SFARI database. Few CNVs are reported to overlap these two genes in control individuals and in NDD patients according to Developmental Cases and Decipher databases. Of note, NDD patients reported in Decipher and represented in Figure 1 carried one additional CNV involving other NDD-associated genes. Indeed, a case has been previously reported with two deleterious variants, one of which disrupting LRRC4C and one affecting another SFARI-reported gene DPP6 (Maussion et al., 2017). The two variants, as in our case, were inherited one from the mother and one from the father, both unaffected.
FIGURE 1. Additive effects between CNTNAP2 and LRRC4C in patient IGGAC13. (A) Screenshots of chromosome 11 region overlapping the deletion of patient IGGAC13 (red bar IGGAC13, chr11:40,303,114-40,694,847). Decipher CNVs (a: deletion in patient 282991; b: deletion in patient 276911). (B) Screenshots of chromosome 7 region overlapping the deletion of patient IGGAC13 (red bar IGGAC13, chr7:146,464,283-146,532,437). Decipher CNVs (a: deletion in patient 319381; b: intragenic duplication in patient 255117; c: deletion in patient 250646). Only Decipher patients with additional CNVs that overlap neurodevelopmental genes were selected and reported in the figure. (C) The network depicts the protein-protein interactions (grey lines), experimentally determined interaction (red lines) and co-expression (bold black lines) of the KEGG pathway “Cell adhesion molecule,” revealed by STRING analysis on the basis of enrichment obtained with GeneCodis4, and visualized by Cytoscape tool. The Figure shows the clustering of the candidate genes with other known NDD genes. In yellow the interacting molecules directly linking the two candidate genes to each other.
CNTNAP2 encodes a member of the neurexin family which functions in the nervous system as cell adhesion molecule. LRRC4C (NGL1) is a specific binding partner for netrin G1, which is a member of the netrin family of axon guidance molecules. In fact, both molecules, CNTNAP2 and LRRC4C, are implicated in synapse formation (Um and Ko, 2017). We hypothesized that, while half dosage of each of these two molecules is compatible with a normal development, half dosage of both genes can impair synapse correct development and function. In fact, gene enrichment and String analysis showed that both genes are members of the KEGG cell adhesion molecule pathway and interact with other known NDD genes (Figure 1C; Supplementary Table S1).
The second case is a male patient (IGGAC14) that inherited from his healthy father a deletion of chromosome 17 encompassing a SFARI gene, DNAH17, and its antisense gene DNAH17-AS1. DNAH17 encodes a heavy chain associated with axonemal dynein, a member of microtubule-associated motor protein complexes. From his healthy mother, the patient inherited a duplication involving SHOX on chromosome X (Figure 2), also reported in the SFARI database.SHOX belongs to the paired homeobox family and is located in the pseudoautosomal region 1 (PAR1) of X and Y chromosomes. Enrichment of SHOX microduplications in NDD cases, particularly in those with ASD, has been described (Tropeano et al., 2016). Gene enrichment evidenced a function of SHOX in the regulation of transcription by RNA polymerase II, and of DNAH17 in microtubule-based movement, two important biological processes in neurodevelopment disorder (Figure 2; Supplementary Table S1). As shown in Figure 2A, other patients with NDD phenotype are reported in Decipher database to have duplications of SHOX and a secondary CNV involving neurodevelopmental genes. Indeed, a chromosome 17 deletion similar to that of our patient was reported in Decipher database in a case with NDD phenotype having a secondary CNV also overlapping developmental genes (Figure 2B).
FIGURE 2. Additive effects between SHOX and DNAH17 in patient IGGAC14. (A) Screenshot of chromosome X region overlapping the duplication of patient IGGAC14 (blue bar, chrX:561,527-908,212). Duplications in Decipher database (a: patient 259967; b: patient 367079; c: patient 380394; d: patient 284815; e: patient 289670; f: patient 322776; g: patient 367150; h: patient 279033; i: patient 353723). (B) Screenshot of chromosome 17 region encompassing the deletion of patient IGGAC14 (red bar, chr17:76,419,867-76,548,866). Deletions in Decipher database (a: patient 289758).
We observed a family with two affected siblings, a male (IGGAC15) and a female (IGGAC16), showing discordant phenotypes. They inherited from their unaffected mother a deletion encompassing a syndromic region at 15q11.2 (Figure 3A). The siblings also shared a duplication involving the SFARI gene SYNGAP1 (Figure 3B). Of note, the SYNGAP1 duplication was not observed in parents, at least in their blood-extracted DNA, thus leading to suppose mosaicism in either father or mother. Patient IGGAC16 also inherited a paternal, recurrent deletion involving the SFARI gene RBFOX1 (Figure 3C). Interestingly, loss of function variants of SYNGAP1 have been identified in patients with ASD and intellectual disability with or without epilepsy. However, to the best of our knowledge, no cases with duplication or gain of function variants have been reported so far. Overexpression of SynGAP was reported to block neurite outgrowth by a mechanism that involves Ras-like GTPase cascade (Tomoda et al., 2004), to produce a depression of AMPAR-mediated excitatory postsynaptic currents in neurons (Rumbaugh et al., 2006), and to regulate decreasing and increasing miniature excitatory synaptic currents in hippocampal neurons (McMahon et al., 2012). In fact, some phenotypic differences between the two siblings were observed. The male patient showed a mild phenotype characterized by a borderline cognitive functioning associated with atypical behaviours and psychomotor delay in early infancy, his sister had a more complex phenotype with mild cognitive impairment, irregular sleep rhythms, hyperactivity, coordination, and language disorders. We could speculate that the two CNVs shared by the two siblings can explain the phenotypic features shared by both of them, while the third CNV could have worsened the phenotype of the female patient.
FIGURE 3. Additive effects in a familial case, siblings IGGAC15 and IGGAC16. (A) Screenshot of chromosome 15 region encompassing the deletion shared by the two siblings (red bars, chr15:22,765,628-23,208,901). (B) Screenshot of chromosome 6 region encompassing the duplication shared by the two siblings (blue bars, chr6:33,337,684-33,669,083). (C) Screenshot of chromosome 16 region encompassing the deletion of patient IGGAC16 (red bar, chr16:6,371,455-6,568,176). Deletions reported in Decipher patients (a: patient 283262; b: patient 327250; c: patient 403649).
Addictive Effects of CNVs Involving Potential Novel Candidate Genes and New Pathogenic Mechanisms
Four patients reported double CNVs involving genes or genomic regions not yet described as NDD-associated but that could act in concert, similarly to the four cases described above, with a potential impact on the patients’ phenotype (Table 2).
The first case (IGGAC06) carried a maternal deletion of chromosome 9 and a non-maternal deletion of chromosome 11; as genetic analysis of the patient’s father was not available, we cannot determine if the chromosome 11 deletion was de novo or paternally inherited (Figures 4A,C). The chromosome 9 deletion encompassed PTPRD, which encodes a protein-tyrosine phosphatase receptor, and PTPRD-AS2, a long non-coding RNA (lncRNA) gene known as PTPRD antisense RNA 2 (head-to-head) predicted to make expression of mRNA more stable and thus enhancing PTPRD protein production. Thus, a reduced expression of PTPRD can be hypothesized because of the deletion of the first two untranslated exons of the gene and of PTPRD-AS2. A total of 7 deletions encompassing exonic regions of PTPRD are reported in controls, while several CNVs and, in particular, 16 deletions and one intragenic duplication involving coding regions of the gene are present among cases in Developmental Delay and Decipher databases with NDD phenotypes. Five among Decipher cases have an additional CNV involving brain expressed genes that, analogously to those reported in patient IGGAC06, could have an additive effect in causing patient’s phenotype. Interestingly, PTPRD is predicted to be intolerant to loss of function variants (pLI = 1) and to haploinsufficiency mechanisms (score = 0.75; value above 0.5 are predicted to be haploinsufficient) (Shihab et al., 2017). PTPRD has a role in synaptic adhesion and synapse organization and can bidirectionally induce pre- and postsynaptic differentiation of neurons by trans-synaptically binding to interleukin-1 receptor accessory protein (IL1RAP) (Yoshida et al., 2012; Yamagata et al., 2015). In the mouse, PTPRD was shown to regulate neurogenesis (Tomita et al., 2020) and mice lacking PTPRD showed impaired learning with enhanced hippocampal long-term potentiation (Uetani et al., 2000), thus suggesting a role for PTPRD in neurodevelopmental disorders. In OMIM database, PTPRD is not associated with a disorder, and its implication in human disease is still to be completely clarified. Homozygous deletions of PTPRD were observed in a patient with intellectual disability and trigonocephaly (Choucair et al., 2015) and his family members carrying the same deletions were reported to be unaffected. Instead, a heterozygous de novo splicing variant of the gene was reported in a girl with moderate nonsyndromic developmental delay (Yan et al., 2019). Furthermore, genetic studies showed an oligogenic association of PTPRD variants with obsessive-compulsive disorder (Mattheisen et al., 2015) and with restless leg syndrome (Schormair et al., 2008). Overall, results from various genetic studies are consistent with both major/oligogenic and modest/polygenic contributions of common and rare PTPRD variations in neurological, behavioural, and neurodevelopmental disorders (Uhl and Martinez, 2019). Gene enrichment and String analysis showed PTPRD interacting with many NDD genes playing a role in synaptic membrane adhesion and neuronal development (Figure 4B; Supplementary Table S1).
FIGURE 4. Additive effects between PTPRD and BUD13 in patient IGGAC06. (A) Screenshot of chromosome 9 region overlapping the deletion of patient IGGAC06 (red bar, chr9:10,286,927-11,086,155). Decipher CNVs (a: intragenic duplication in patient 333436; b: deletion in patient 412118; c: deletion in patient 337260; d: intragenic duplication in patient 260089 and e: deletion in patient 279686). (B) The network depicts the protein-protein interactions (grey lines) revealed by STRING analysis (visualized by Cytoscape tool) on the basis of the enrichment obtained with GeneCodis4 tool. The Figure shows the enriched GO terms clustering the candidate gene PTPRD and known NDD genes. (C) Screenshot of chromosome 11 region overlapping the deletion of patient IGGAC06 (red bar, chr11:116,182,875-116,677,043). (D) The network depicts the protein-protein interactions (grey lines) and experimentally determined interaction with a SFARI gene (red line) revealed by STRING analysis (visualized by Cytoscape tool) on the basis of the GO term “mRNA splicing via spliceosome”. The network indicates a direct interaction of BUD13 with the known NDD gene SF3B1.
The chromosome 11 deletion encompassed BUD13, which encodes a component of the retention and splicing (RES) complex of the spliceosomal complex. In zebrafish, the lack of BUD13 caused defects in intron splicing specifically in genes with neurodevelopmental regulatory functions, thus resulting in a decrease of differentiated neurons and brain developmental defects (Fernandez et al., 2018). GeneCodis analysis showed an enrichment of BUD13 with various NDD genes in the GeneOntology category of mRNA splicing, via spliceosome, a biological process with an important role in neurodevelopment (Figure 4D, Supplementary Table S1).
Patient IGGAC08 has a complex phenotype mainly characterized by language disorder, hyperactivity, learning difficulties, borderline cognitive level, motor skill impairment, and short stature (Table 1). His mother also had school difficulties and showed dyslexia in childhood. Mother and son share a deletion of chromosome 4q34.1 encompassing GLRA3 which encodes the alpha-3 subunit of the neuronal glycine receptor, a ligand-gated ion channel. In addition, patient IGGAC08 had a de novo deletion involving hsa-mir-4465 and a de novo duplication including INPP5A (Table 2 and Figures 5A–C). GLRA3 receptors are expressed in spinal cord, brainstem, hippocampus, amygdala, striatum, and cortex. The glycine receptor is a glycoprotein composed of 5 subunits, three α and two β subunits. The alpha subunits bind the glycine ligand while the beta subunits bind to gephyrin, a cytosolic protein required for a regulated synaptic aggregation and clustering of these receptors.
FIGURE 5. Additive effects of INPP5A, MIR4465 and GLRA3 in patient IGGAC08. (A) Screenshot of chromosome 10 region encompassing the duplication of patient IGGAC08 (blue bar, chr10:134,289,353-134,917,702). (B) Screenshot of chromosome 6 region encompassing the deletion of patient IGGAC08 (red bar, chr6:140,767,462-141,008,265). (C) Screenshot of chromosome 4 region encompassing the deletion of patient IGGAC08 (red bar, chr4:175,734,964-175,934,265).
Gephyrin has functional links with several synaptic proteins, mutations of which have been reported in various neurodevelopmental disorders (Choii and Ko, 2015; Kim et al., 2021). Half dosage of GLRA3 may contribute to impair gephyrin-mediated aggregation and post-synaptic clustering.
Loss of PTEN, a SFARI gene, can cause postsynaptic and presynaptic changes in excitatory and inhibitory connectivity. In vitro analyses have recently demonstrated that a miRNA (microRNA-301a) can downregulate PTEN after glycine receptor activation (Chen et al., 2016). Tao et al. (2019) showed that miR-4465 significantly inhibited the expression of PTEN, upregulated phosphorylated AKT and ultimately inhibited autophagy by activating mTOR in HEK293, HeLa, and SH-SY5Y cells (Tao et al., 2019). Half dosage of hsa-mir-4465 in our patient is expected to cause upregulation of PTEN, contributing to postsynaptic impairment. MiRNAs are important regulators of brain development and neuronal function, and have been associated with a variety of nervous system diseases, including ASD (Cheng et al., 2018; Wu et al., 2020).
In mice, the deletion of Inpp5a causes perinatal lethality in 90% of the homozygous mutants, early onset ataxia and relatively small stature in surviving mutants. Heterozygotes do not exhibit obvious motor coordination impairment unless challenged in motor skill test (Yang et al., 2015). Inpp5a is a downstream effector of signalling from mGlu1 receptor, one of the two members of group I metabotropic glutamate receptors, that regulates synaptic plasticity, particularly in cerebellar Purkinje cells. The mGlu1 receptor activity needs to be fine-tuned and balanced for normal motor coordination. Alteration of expression, either enhancement or decrease, of group 1 mGlu metabotropic receptors in mice is known to alter their activity, which results in defects of motor coordination (Rossi et al., 2013; Bossi et al., 2018). The deletion of INPP5A in this patient may be in part responsible for his impairment in motor skills.
IGGAC10 patient showed a complex phenotype mainly characterized by intellectual disability and epilepsy (Table 1). He carried a maternal deletion encompassing SYNCRIP and a paternal duplication involving the upstream regulatory genomic region of ADCY5. Both CNVs overlap a TDB between two flanking TADs described in different tissues, including brain cortex, which could also contribute to dysregulating the expression of implicated genes (Figures 6A,B).
FIGURE 6. Additive effects between SYNCRIP and ADCY5 in patient IGGAC10. (A) Screenshot of chromosome 6 region encompassing the deletion of patient IGGAC10 (red bar, chr6:86,338,564-86,413,285). TADs from Cortex, H1NPC and h-ESC cells. (B) Screenshot of chromosome 3 region encompassing the duplication of patient IGGAC10 (blue bar, chr3:123,113,523-123,997,093). TADs from Cortex, H1NPC and h-ESC cells.
SYNCRIP encodes the synaptotagmin-binding cytoplasmic RNA-interacting protein, a member of the cellular heterogeneous nuclear ribonucleoprotein (hnRNP) family that plays a role in multiple aspects of mRNA maturation. Missense and truncating variants of SYNCRIP were found in ASD patients with more severe phenotypes (Guo et al., 2019) and in patients with severe non-syndromic sporadic intellectual disability (Rauch et al., 2012), respectively. ADCY5 encodes a member of the membrane-bound adenylyl cyclase enzymes which mediate G-protein-coupled receptor signalling through the synthesis of the second messenger cAMP. According to SFARI database, missense and frameshift variants of ADCY5 have been reported in ASD patients. Gain of function variants of ADCY5 have been associated with a broad range of movement disorders, most notably chorea, dystonia and myoclonus (Vijiaratnam et al., 2019).
Thus, the two CNVs could cause, by a mechanism also implicating a modification of chromatin conformation, a dysregulation of both SYNCRIP and ADCY5 contributing to the complex phenotype observed in this patient.
IGGAC07 patient carried a de novo duplication of chromosome 14, and two maternally inherited duplications of chromosome 12 and chromosome 22, respectively (Table 2; Figures 7 A–C).
FIGURE 7. Additive effects of ABHD4, GNAZ and WSCD2 in patient IGGAC07. (A) Screenshot of chromosome 14 region encompassing the duplication of patient IGGAC07 (blue bar, chr14:22,360,671-23,120,435). Brain-expressed genes (see Methods). (B) Screenshot of chromosome 22 overlapping the duplication of patient IGGAC07 (blue bar, chr22:22,998,284-23,643,223). Duplications in Decipher database (a: patient 251352; b: patient 286708; c: patient 301091; d: patient 288064). Brain-expressed genes (see Methods). (C) Screenshot of chromosome 12 overlapping the duplication of patient IGGAC07 (blue bar, chr12:108,487,109-108,896,252). Brain-expressed genes (see Methods).
Interestingly, chromosome 14 duplication encompassed ABHD4, a gene recently shown to be involved in developmental anoikis, a mechanism of cell death in the prenatal brain preventing from survival of misplaced cells (Laszlo et al., 2020). Indeed, ABHD4 is involved in N-acylethanolamine biosynthesis, including the endocannabinoid molecule anandamide which has various physiological functions such as the regulation of synaptic plasticity and apoptosis. ABHD4 expression is tightly controlled spatially and temporally in the developing brain and a misregulation of its expression due to the duplication could contribute to the patient’s phenotype. Duplications encompassing ABHD4 were annotated in only one control individual, while various cases were reported both in Decipher and in Developmental Delay databases with duplications larger than that observed in our patient and a smaller one sharing with our patient ABHD4 duplication.
The maternal duplication on chromosome 22 involved a region with triplosensitivity including GNAZ, the Guanine Nucleotide-binding protein (G protein) Alpha Z polypeptide encoding gene. In this region, many duplications are reported in both Decipher and Developmental Delay databases, most of them extending over the short region containing GNAZ. Only a total of 4 control individuals are reported in either Developmental Delay Control or in DGV databases with duplications overlapping the GNAZ gene. GNAZ is a member of the Gαi subfamily of heterotrimeric G proteins, and couples to G-protein-coupled receptors triggering important pathways in neuronal development (Hultman et al., 2014). The duplication of chromosome 12 entirely encompassed WSCD2, which encodes a WSC Domain-containing protein with a still unclear function. Interestingly, genome wide association analyses found significant results for variants in WSCD2 in patients with psychiatric disorders (Lo et al., 2017). Finally, we could hypothesise that all the de novo and the inherited duplications could have jointly contributed to the NDD phenotype of patient IGGCA07.
CNV-Mediated Complex Pathogenic Mechanisms
We found 4 patients with CNVs that could affect expression of NDD genes throughout an indirect and/or long-range effect (Table 1 and Table 2).
One case (IGGAC04) presented with very severe psychomotor delay, hypotonia, and epilepsy (Table 1), carrying a de novo heterozygous deletion of chromosome 5 (Figure 8A). This deletion encompasses 6 genes, two of them expressed in the nervous system: MEF2C-AS1, which encodes a lncRNA annotated as antisense RNA of MEF2C, and ADGRV1, which encodes a member of the G-protein-coupled receptor superfamily associated with familial febrile seizures and Usher syndrome 2 (OMIM 602851). MEF2C (myocyte enhancer factor-2, OMIM 600662) is a transcription factor playing a role during development, mostly in myogenesis and neurogenesis. MEF2C haploinsufficiency has been shown to be responsible for severe cognitive deficit with stereotyped movements, epilepsy, and/or cerebral malformations (Le Meur et al., 2010). It has been recently demonstrated that, during myogenic differentiation a lncRNA, the OIP5-AS1, through its complementarity with MEF2C 3′UTR, enhances MEF2C mRNA stability thus promoting MEF2C expression and ultimately enhancing myogenesis (Yang et al., 2020). We hypothesized that during neurogenesis, analogously to OIP5-AS1, MEF2C-AS1 could enhance MEF2C mRNA stability and promote MEF2C expression and neurogenesis. In this view, the heterozygous deletion observed in patient IGGAC04 could reduce the expression of MEF2C-AS1 causing a decrease of MEF2C, thus interfering with neuronal differentiation and contributing to the patient’s phenotype. Although as a consequence of MEF2C-AS1 deletion a gene loss-of-function is likely, functional experiments would be required to exclude some compensatory mechanisms, as up-regulation of the healthy allele. The deletion of ADGRV1 could be responsible for the epilepsy phenotype. Of note, we found one patient in the Decipher database (Patient 251716) carrying a deletion overlapping the one of our patient and having similar clinical features characterized by intellectual disability, muscular hypotonia, and seizures. No deletions overlapping the same region were present among controls (Developmental Delay and DGV Controls). The chromosome 5 region deleted in these patients overlaps a TDB described in different tissues, including those relevant for neurodevelopment (e.g., Cortex), according to their relative Hi-C maps (Figure 8A), which could further complicate the fine regulation of expression of the overlapped genes, including MEF2C-AS.
FIGURE 8. Potential candidate genes and pathogenetic mechanisms in patient IGGAC04 and IGGAC01. (A) Screenshot of chromosome 5 region overlapping the deletion of patient IGGAC4 (red bar, chr5:88,232,587-90,181,774). Decipher CNVs (a: deletion in patient 251716). TADs from Cortex, H1NPC and h-ESC cells. (B) Screenshot of chromosome 16 region overlapping the duplication of patient IGGAC01 (blue bar, chr16:8,969,984-10,078,802). Decipher CNVs (a: duplication in patient 382832).
One case (IGGAC01) showed a complex phenotype mostly characterized by psychomotor delay, autistic features, hypotonia, epilepsy, gastrointestinal and sleep disturbances (Table 1), and carried a de novo heterozygous duplication of chromosome 16 (Figure 8B). The duplication involved, in particular, two OMIM reported genes, USP7, encoding the Ubiquitin Specific Peptidase 7, and GRIN2A, encoding the NR2 subunit of the N-methyl-D-aspartate (NMDA) glutamate receptor. USP7 has been shown to regulate the ubiquitination of a variety of proteins with functions in different biological processes as oxidative stress response, histone modification and regulation of chromatin remodelling (Khoronenkova et al., 2011). Deletions, truncating and missense variants of USP7 have been associated with the Hao-Fountain syndrome (OMIM 616863), characterized by speech delay, autistic spectrum disorder, attention-deficit hyperactivity disorder, sleep disturbances and gastroesophageal reflux disease (Fountain et al., 2019). All these clinical features are present in our case who has the complete duplication of USP7. Notably, in syndromic forms of ASD, duplications in addition to deletions of USP7 have been reported (Sanders et al., 2011), thus supporting the idea of dosage sensitivity, with both low and high USP7 expression causing imbalances of neuronal homeostasis. The duplication present in patient IGGAC01 partially overlapped GRIN2A. Truncating as well as missense, activating variants of GRIN2A have been associated with an idiopathic form of focal epilepsy, EEG continuous spike-and-wave during sleep (CSWS), and speech disorder (OMIM 245570) (Lemke et al., 2013) whose clinical aspects are present in our patient. Different molecular alterations of NMDA receptor subunits seem to result in a deleterious dysregulation of NMDA receptor function (Balu and Coyle, 2011). Although the effect of partial duplications on gene function are not easily predictable, some clinical features of the IGGAC01 patient suggested the occurrence of a dysregulation of GRIN2A that could have contributed to the patient’s phenotype. Of note, a patient with phenotype resembling that of patient IGGAC01 was reported in Decipher database with a similar duplication, including both USP7 and GRIN2A (Figure 8B). We thus hypothesised that the phenotype of these two patients could result from the impairment of both USP7 and GRIN2A.
IGGAC24 had a de novo duplication of chromosome 19 that encompasses various genes, 11 of which are expressed in the brain. The duplication also overlaps a TDB described in the brain cortex and other different tissues, which could contribute to dysregulating gene expression of implicated genes (Figure 9A and Table 1). Attempting to prioritize enclosed genes and define possible candidate genes, we performed a gene enrichment analysis as described above (see Methods). Two genes emerged for their functions known to play a role in neurodevelopmental disorders, namely SHISA7 and U2AF2. Indeed, in Decipher at least one NDD case with a duplication that partially overlaps the TDB and the SHISA7 gene (Patient 338729) was reported while no similar duplications, implicating brain-expressed genes, were reported among controls. We focused on SHISA7, a transmembrane protein that interacts with GABAA receptors (Han et al., 2019) and, by acting as auxiliary protein for the trafficking of AMPA receptors, regulates synaptic plasticity (Schmitz et al., 2017). Notably, spine loss has been observed in SHISA7 knockout mice and in wild type mice with SHISA7 overexpression (Schmitz et al., 2017; von Engelhardt, 2019). Thus, we hypothesised that the duplication could cause a dysregulation of SHISA7 expression, possibly also by altering chromatin conformation, and ultimately causing a dysregulation of synaptic plasticity.
FIGURE 9. Duplications of topological domain boundaries (TDBs) potentially involved in the phenotype of patient IGGAC24 and patient IGGAC23. (A) Screenshot of chromosome 19 region encompassing the duplication of patient IGGAC24 (blue bar, chr19:55,948,706-56,471,996). TADs from Cortex, H1NPC and h-ESC cells. Brain-expressed genes (see Methods). Decipher patient with reported duplication (a. patient 338729) (B) Screenshot of chromosome 10 region encompassing the duplication of patient IGGAC23 (blue bar, chr10:125,121,038-127,119,447). Brain-expressed genes (see Methods). Decipher patients with reported duplications (a: patient 283429, b: patient 378968).
Patient IGGAC23 showed a de novo duplication of chromosome 10 overlapping 7 brain-expressed genes and a TDB described in the brain cortex and other tissues (Figure 9B and Table 1). After gene enrichment analysis of these 7 genes with known NDD genes, one emerged as the most interesting, CTBP2. This gene produces two alternative transcripts for two different proteins, one of them is a component of specialized synapses (ribbon synapses) and the other one is a negative regulator of transcription. CTBP2 knock-out mice present a neurological phenotype characterized by abnormalities of brain development. A recent study reported that the overexpression of CTBP2 inhibits neuronal differentiation (Zhang et al., 2020). Of note, while no duplications similar to that of this patient were reported among controls, in Decipher a patient with autistic behaviour and a duplication shorter than the one of patient IGGC23, but still enclosing CTBP2, was reported. Thus, the duplication present in patient IGGAC23 could have caused overexpression and/or dysregulated expression of overlapped genes, among which CTBP2, whose overexpression has already been associated with impaired neuronal defects.
Discussion and Conclusion
In the frame of identifying the most likely genes and genetic mechanisms underlying complex phenotypes in neurodevelopmental disorders, and to better interpret CNV causal effects in complex cases, we performed a detailed analysis of both clinical features and CNVs of 12 patients.
In particular, we analyzed patients’ CNVs in the light of information on gene expression, gene function, and chromatin organization features of involved genomic regions. Then we assessed the presence of potential deleterious CNVs acting in concert in the same patients, according to a two-hit model emerging from the recent literature as a possible pathogenetic mechanism with important roles in neurodevelopmental disorders.
Four patients reported with a complex phenotype, i.e. not referable to a known syndromic neurodevelopmental disorder, carried a CNV that can possibly have indirect deleterious effects on neuronal expressed genes and finally on their phenotype. In one case, a deletion encompassed MEF2C-AS1, the antisense gene of a known NDD-associated gene, MEF2C. Involvement of antisense genes in NDDs has been rarely reported so far. One example is represented by PTCHD1-AS, the antisense-gene of the NDD-associated gene PTCHD1. In fact, although the molecular function of PTCHD1-AS is still to be clearly elucidated, it has been demonstrated that its disruption diminished miniature excitatory postsynaptic current frequency in iPSC-derived neurons from subjects with ASD, thus supporting a role for this long noncoding RNA in the etiology of ASD (Ross et al., 2020). Interest on regulatory elements of gene expressions, like lncRNA and regulatory enhancer elements, and their involvement in brain disorders has been recently increased and new tools have been applicated for their detection, as those implicated in ASD and Schizofrenia (Piluso et al., 2019; Alinejad-Rokny et al., 2020).
Deletions and duplications may overlap TDBs between two flanking topological associated domains causing a dysregulation of chromatin organization and, in turn, a dysregulation of implicated genes which ultimately affects the patient’s phenotype. So far, implication of CNVs in TAD boundaries have been rarely investigated in NDDs (Di Gregorio et al., 2017; Melo et al., 2020). Chromatin conformation domains can be different in various cells and tissues, and we were aware that effects of deletions and duplications deduced on the basis of TDBs data reported in databases have to be considered hypothetical, deserving further experimental validation. To tentatively achieve a more reliable data analysis of our patients, we considered only TDBs reported in more than one tissue, assuming that they could possibly represent a chromatin feature common to many tissues, including brain areas involved by NDDs. All our four cases have a CNV encompassing a TDB between two flanking topologically associated domains described in at least three tissues which could have played a role in NDD phenotype, such as brain cortex (Schmitt et al., 2016) and human embryonic stem cell-derived neurons (Dixon et al., 2015).
One important finding emerging from this evaluation of diagnostic array-CGH results is the presence of co-occurring CNVs that, according to the two-hit model, could have contributed to the patients’ phenotype.
In eight patients, we found two CNVs overlapping either known NDD genes and/or genes with potential roles in neurodevelopment. In three patients, we found one de novo and one variant inherited from unaffected parents, the remaining five patients inherited one CNV from each of their unaffected parents. In three out of eight cases, namely IGGAC06, IGGAC07, and IGGAC08, involved genes have not yet been associated with disease and we presented data in favor of their implication in neurodevelopmental disorders. Indeed, we searched the Decipher database for cases with one CNV corresponding to that presented by our patients and a secondary CNV, involving genes with a role in neurodevelopment. We found some cases, reported in the figures, that further support the hypothesis that identified co-occurring CNVs in our patients act additively causing deleterious phenotypes according to a two-hit model.
In order to gain insight into possible interaction among genes implicated in the CNVs co-occurring in the same patient and to evaluate whether these genes share molecular functions or biological processes, we performed gene enrichment and protein-protein interaction analyses (Figure 10; Supplementary Table S1). We also tested interaction with other known NDD associated genes, which helped to better define the role of new candidate genes in important neurodevelopmental pathways and biological processes. In one case, IGGAC13, the two major implicated genes, i.e., CNTNAP2 and LRRC4C, interact with each other through NRXN1 and with other NDD-associated genes and have a common function in the KEGG cell adhesion molecule pathway (Figure 1). In all the other cases, genes affected by co-occurring CNVs do not seem to directly interact, but they have a role in important, although distinct, biological processes that together contribute to proper brain development (Figure 10). We thus hypothesized that, in all these cases, the double-hit mechanism by affecting two pathways/biological processes finally determines an impairment of neurological development.
FIGURE 10. Overview of biological processes and pathways involved in patients with CNV-genes acting according to double-hit model. The network depicts the protein-protein interactions (grey lines) revealed by STRING analysis (visualized by Cytoscape tool) on the basis of the enrichment obtained with GeneCodis4 tool. The Figure shows the enriched BP terms and KEGG pathways clustering the CNV-genes of the eight patients with known NDD genes (SFARI reported genes with score 1/2 and OMIM genes correlated to NDDs). The analysis shows that all eight patients have CNV-overlapped genes involved in two known NDD-associated biological processes, except patient IGGAC13, for whom the two implicated genes appeared in the same KEGG pathway. GO terms “neuron projection development”, “synaptic membrane adhesion”, “ion trasport”, “G-protein-coupled receptor signalling pathway” and “mRNA splicing” are shared by at least two different patients.
As observable in Figure 10, couple of patients have genes involved in the same pathway/biological process and one additional gene dealing with another biological process. As an example, both IGGAC06 and IGGAC10 have one CNV each involving a different gene but with a common function in mRNA splicing, and each of the two patients has an additional CNV involving one gene with a role, respectively, in synaptic membrane adhesion and in G-coupled receptor signaling, two fundamental functions of synapse.
Of note, co-occurring secondary CNVs overlapping NDD genes were also found in patients with syndromic or recurrent gene-specific CNVs, which could modulate penetrance and/or severity of the disease. This phenomenon is well exemplified by the case of two siblings, IGGAC15 and IGGAC16, who shared two CNVs overlapping known NDD genes, including an already described deletion of chromosome 15 with reduced penetrance (Cox and Butler, 2015). In addition to these two CNVs, patient IGGAC16 also carried a deletion of chromosome 16 implicating RBFOX1, a known NDD gene, also described with reduced penetrance (Lal et al., 2013). Patient IGGAC16, the one of the two siblings who carried a total of three CNVs, exhibited a more complex phenotype compared to her brother.
Overall, our genetic analysis allowed to unveil candidate genes and to infer some synergistic and complex mechanisms that could origin the neurodevelopmental disorder in these patients. However, novel candidate genes and genes affected by co-occurring CNVs although individually associated with NDD biological processes, require further functional analyses to understand their potential interactions and ultimately to clarify their contribution to patient phenotypes.
Variability and complexity of genetic components make the functional validation of NDD genes and pathogenic mechanisms a challenge. Validation of hypothesised mechanisms in one patient can represent a research project by itself. Concerning the present cases, we could consider different levels of validation starting from verifying the effect of a CNV on expression of the implicated genes in induced neurons (iN) that could be obtained from patients’ fibroblasts. With this approach we could evaluate the effects of deletions of antisense or miRNA genes on the expression of their specific targets (as for patient IGGAC4 and IGGAC8) or of CNV involving TDB boundaries on implicated genes (as for patient IGGAC10). A further level of investigation is represented by the analysis of morphology and connections of iNs, usefull to investigate impairment of cell adhesion molecules, as for the case IGGAC13. Patient-derived iN could be used to evaluate defects on RNA splicing as in the case of IGGAC6 patient. To evaluate the effects of more complex genetic mechanisms on brain development and on behaviour, a further level of validation can be represented by studies with animal models, as mice. Indeed, as mice with double/multiple heterozygous variants can be obtained (Bossi et al., 2018), the synergistic effects of multiple variants could be also evaluated. In this view, and as mouse models for both LRRC4C and CNTNAP2 genes are available, it could be interesting to obtain a mouse strain carrying heterozygous, loss of function mutations in both genes, thus mimicking the genotype of patient IGGAC13.
Our findings further support the hypothesis that double-hit mechanisms, involving inherited and de novo CNVs, could synergistically cause, or modulate, the neurodevelopmental phenotype. Indeed, the co-occurrence of pathogenic variants has also been observed after analysing the CNVs together with the results of whole exome sequencing, underlining the importance of the genetic background for the penetrance of a deleterious variant (Pizzo et al., 2019).
Recently reported data indicate that, in NDDs, WES analyses have diagnostic yields of about 36% while chromosomal microarray analysis of less than 20% (Srivastava et al., 2019; Savatt and Myers, 2021). As variants causing NDDs can be CNVs or single nucleotide variants (SNVs), isolated or even combined, the genetic diagnostic approach has to take into consideration methods able to detect both type of variants. Indeed, several clinical laboratories are now incorporating CNV calling into WES and are able to detect multi-exon deletions and duplications (Savatt and Myers, 2021).
WGS, by definition, detects more comprehensively all classes of genetic variants and, as cost decreases, will eventually supersede WES and chromosomal array analyses in clinical testing algorithms over time.
Finally, our work describes a possible integrated framework of investigation of CNVs based on publically available resources of genomic annotations and bioinformatic analyses to prioritize candidate genes and infer pathogenic mechanisms. This approach could be used for interpretation of those CNVs that, after a first level analysis, remains of uncertain significance.
Indeed, it would be interesting to test these patients by WES to check for single nucleotide variants which, according to multi-hit model, could further contribute to their phenotype.
Data Availability Statement
The datasets presented in this study can be found in online repositories. Clinical and genetic variant details have been deposited in the Decipher database with the following accession numbers: 433116 (IGGAC01); 433118 (IGGAC04); 433122 (IGGAC06); 433123 (IGGAC07); 433125 (IGGAC08); 433130 (IGGAC10); 433133 (IGGAC13); 433135 (IGGAC14); 433138 (IGGAC15); 433139 (IGGAC16); 433172 (IGGAC23); 433173 (IGGAC24).
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Italian Regione Liguria (R.P. 001/2019). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
MS, LP, and AP participated in the design of the study. PR and ET were responsible for array-CGH experiments. PR, MC, and AP performed some array-CGH validation experiments. ML, MD, GR, LP, SB, and LN collected all clinical phenotypes and MS managed the database of clinical and CNVs data. MS, LP, MC, and AP participated in data analysis and interpretation. MS, LP, and AP wrote the manuscript. All authors critically revised the manuscript and agreed to be accountable for the content of the work.
Funding
Funding for this research was provided by Italian Ministero della Salute (Project N. RF-2016-02361949 to AP; “Cinque per Mille” and Ricerca Corrente). This work was developed within the framework of the DINOGMI Department of Excellence (MIUR 2018-2022, legge 232 del 2016).
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors wish to thank all families and referring nurses of Child Neuropsychiatric Unit for their help and collaboration.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.732002/full#supplementary-material
References
Addis, L., Sproviero, W., Thomas, S. V., Caraballo, R. H., Newhouse, S. J., Gomez, K., et al. (2018). Identification of New Risk Factors for Rolandic Epilepsy: CNV at Xp22.31 and Alterations at Cholinergic Synapses. J. Med. Genet. 55 (9), 607–616. doi:10.1136/jmedgenet-2018-105319
Alinejad-Rokny, H., Heng, J. I. T., and Forrest, A. R. R. (2020). Brain-Enriched Coding and Long Non-coding RNA Genes Are Overrepresented in Recurrent Neurodevelopmental Disorder CNVs. Cel Rep. 33 (4), 108307. doi:10.1016/j.celrep.2020.108307
Alvarez-Mora, M. I., Corominas, J., Gilissen, C., Sanchez, A., Madrigal, I., and Rodriguez-Revenga, L. (2021). Novel Compound Heterozygous Mutation in TRAPPC9 Gene: The Relevance of Whole Genome Sequencing. Genes 12 (4), 557. doi:10.3390/genes12040557
Arteche-López, A., Gómez Rodríguez, M. J., Sánchez Calvin, M. T., Quesada-Espinosa, J. F., Lezana Rosales, J. M., Palma Milla, C., et al. (2021). Towards a Change in the Diagnostic Algorithm of Autism Spectrum Disorders: Evidence Supporting Whole Exome Sequencing as a First-Tier Test. Genes 12 (4), 560. doi:10.3390/genes12040560
Balu, D. T., and Coyle, J. T. (2011). Glutamate Receptor Composition of the post-synaptic Density Is Altered in Genetic Mouse Models of NMDA Receptor Hypo- and Hyperfunction. Brain Res. 1392, 1–7. doi:10.1016/j.brainres.2011.03.051
Bossi, S., Musante, I., Bonfiglio, T., Bonifacino, T., Emionite, L., Cerminara, M., et al. (2018). Genetic Inactivation of mGlu5 Receptor Improves Motor Coordination in the Grm1 Mouse Model of SCAR13 Ataxia. Neurobiol. Dis. 109 (Pt A), 44–53. doi:10.1016/j.nbd.2017.10.001
Cerminara, M., Spirito, G., Pisciotta, L., Squillario, M., Servetti, M., Divizia, M. T., et al. (2021). Case Report: Whole Exome Sequencing Revealed Disease-Causing Variants in Two Genes in a Patient with Autism Spectrum Disorder, Intellectual Disability, Hyperactivity, Sleep and Gastrointestinal Disturbances. Front. Genet. 12, 625564. doi:10.3389/fgene.2021.625564
Chen, J., Zhuang, Y., Zhang, Z.-F., Wang, S., Jin, P., He, C., et al. (2016). Glycine Confers Neuroprotection through microRNA-301a/PTEN Signaling. Mol. Brain 9 (1), 59. doi:10.1186/s13041-016-0241-3
Cheng, Y., Wang, Z.-M., Tan, W., Wang, X., Li, Y., Bai, B., et al. (2018). Partial Loss of Psychiatric Risk Gene Mir137 in Mice Causes Repetitive Behavior and Impairs Sociability and Learning via Increased Pde10a. Nat. Neurosci. 21 (12), 1689–1703. doi:10.1038/s41593-018-0261-7
Choii, G., and Ko, J. (2015). Gephyrin: a central GABAergic Synapse Organizer. Exp. Mol. Med. 47, e158. doi:10.1038/emm.2015.5
Choucair, N., Mignon-Ravix, C., Cacciagli, P., Abou Ghoch, J., Fawaz, A., Mégarbané, A., et al. (2015). Evidence that Homozygous PTPRD Gene Microdeletion Causes Trigonocephaly, Hearing Loss, and Intellectual Disability. Mol. Cytogenet. 8, 39. doi:10.1186/s13039-015-0149-0
Coe, B. P., Stessman, H. A. F., Sulovari, A., Geisheker, M. R., Bakken, T. E., Lake, A. M., et al. (2019). Neurodevelopmental Disease Genes Implicated by De Novo Mutation and Copy Number Variation Morbidity. Nat. Genet. 51 (1), 106–116. doi:10.1038/s41588-018-0288-4
Cox, D., and Butler, M. (2015). The 15q11.2 BP1-BP2 Microdeletion Syndrome: a Review. Ijms 16 (2), 4068–4082. doi:10.3390/ijms16024068
Di Gregorio, E., Riberi, E., Belligni, E. F., Biamino, E., Spielmann, M., Ala, U., et al. (2017). Copy Number Variants Analysis in a Cohort of Isolated and Syndromic Developmental Delay/intellectual Disability Reveals Novel Genomic Disorders, Position Effects and Candidate Disease Genes. Clin. Genet. 92 (4), 415–422. doi:10.1111/cge.13009
Dixon, J. R., Jung, I., Selvaraj, S., Shen, Y., Antosiewicz-Bourget, J. E., Lee, A. Y., et al. (2015). Chromatin Architecture Reorganization during Stem Cell Differentiation. Nature 518 (7539), 331–336. doi:10.1038/nature14222
Duyzend, M. H., Nuttle, X., Coe, B. P., Baker, C., Nickerson, D. A., Bernier, R., et al. (2016). Maternal Modifiers and Parent-Of-Origin Bias of the Autism-Associated 16p11.2 CNV. Am. J. Hum. Genet. 98 (1), 45–57. doi:10.1016/j.ajhg.2015.11.017
Fernandez, J. P., Moreno-Mateos, M. A., Gohr, A., Miao, L., Chan, S. H., Irimia, M., et al. (2018). RES Complex Is Associated with Intron Definition and Required for Zebrafish Early Embryogenesis. Plos Genet. 14 (7), e1007473. doi:10.1371/journal.pgen.1007473
Fountain, M. D., Oleson, D. S., Rech, M. E., Segebrecht, L., Hunter, J. V., McCarthy, J. M., et al. (2019). Pathogenic Variants in USP7 Cause a Neurodevelopmental Disorder with Speech Delays, Altered Behavior, and Neurologic Anomalies. Genet. Med. 21 (8), 1797–1807. doi:10.1038/s41436-019-0433-1
Girirajan, S., Rosenfeld, J. A., Cooper, G. M., Antonacci, F., Siswara, P., Itsara, A., et al. (2010). A Recurrent 16p12.1 Microdeletion Supports a Two-Hit Model for Severe Developmental Delay. Nat. Genet. 42 (3), 203–209. doi:10.1038/ng.534
Guo, H., Duyzend, M. H., Duyzend, M. H., Coe, B. P., Baker, C., Hoekzema, K., et al. (2019). Genome Sequencing Identifies Multiple Deleterious Variants in Autism Patients with More Severe Phenotypes. Genet. Med. 21 (7), 1611–1620. doi:10.1038/s41436-018-0380-2
Han, W., Li, J., Pelkey, K. A., Pandey, S., Chen, X., Wang, Y.-X., et al. (2019). Shisa7 Is a GABAA Receptor Auxiliary Subunit Controlling Benzodiazepine Actions. Science 366 (6462), 246–250. doi:10.1126/science.aax5719
Huckins, L. M., Dobbyn, A., Dobbyn, A., Ruderfer, D. M., Hoffman, G., Wang, W., et al. (2019). Gene Expression Imputation across Multiple Brain Regions Provides Insights into Schizophrenia Risk. Nat. Genet. 51 (4), 659–674. doi:10.1038/s41588-019-0364-4
Hultman, R., Kumari, U., Michel, N., and Casey, P. J. (2014). Gαz Regulates BDNF-Induction of Axon Growth in Cortical Neurons. Mol. Cell Neurosci. 58, 53–61. doi:10.1016/j.mcn.2013.12.004
Khoronenkova, S. V., Dianova, , Parsons, J. L., and Dianov, G. L. (2011). USP7/HAUSP Stimulates Repair of Oxidative DNA Lesions. Nucleic Acids Res. 39 (7), 2604–2609. doi:10.1093/nar/gkq1210
Kim, S., Kang, M., Park, D., Lee, A.-R., Betz, H., Ko, J., et al. (2021). Impaired Formation of High-Order Gephyrin Oligomers Underlies Gephyrin Dysfunction-Associated Pathologies. iScience 24 (2), 102037. doi:10.1016/j.isci.2021.102037
Lal, D., Trucks, H., Møller, R. S., Hjalgrim, H., Koeleman, B. P. C., de Kovel, C. G. F., et al. (2013). Rare Exonic Deletions of theRBFOX1gene Increase Risk of Idiopathic Generalized Epilepsy. Epilepsia 54 (2), 265–271. doi:10.1111/epi.12084
László, Z. I., Lele, Z., Zöldi, M., Miczán, V., Mógor, F., Simon, G. M., et al. (2020). ABHD4-dependent Developmental Anoikis Safeguards the Embryonic Brain. Nat. Commun. 11 (1), 4363. doi:10.1038/s41467-020-18175-4
Le Meur, N., Holder-Espinasse, M., Jaillard, S., Goldenberg, A., Joriot, S., Amati-Bonneau, P., et al. (2010). MEF2C Haploinsufficiency Caused by Either Microdeletion of the 5q14.3 Region or Mutation Is Responsible for Severe Mental Retardation with Stereotypic Movements, Epilepsy And/or Cerebral Malformations. J. Med. Genet. 47 (1), 22–29. doi:10.1136/jmg.2009.069732
Lek, M., Karczewski, K. J., Minikel, E. V., Samocha, K. E., Banks, E., Fennell, T., et al. (2016). Analysis of Protein-Coding Genetic Variation in 60,706 Humans. Nature 536 (7616), 285–291. doi:10.1038/nature19057
Lemke, J. R., Lal, D., Reinthaler, E. M., Steiner, I., Nothnagel, M., Alber, M., et al. (2013). Mutations in GRIN2A Cause Idiopathic Focal Epilepsy with Rolandic Spikes. Nat. Genet. 45 (9), 1067–1072. doi:10.1038/ng.2728
Lo, M.-T., Hinds, D. A., Tung, J. Y., Franz, C., Fan, C.-C., Wang, Y., et al. (2017). Genome-wide Analyses for Personality Traits Identify Six Genomic Loci and Show Correlations with Psychiatric Disorders. Nat. Genet. 49 (1), 152–156. doi:10.1038/ng.3736
Mattheisen, M., Samuels, J. F., Wang, Y., Greenberg, B. D., Fyer, A. J., McCracken, J. T., et al. (2015). Genome-wide Association Study in Obsessive-Compulsive Disorder: Results from the OCGAS. Mol. Psychiatry 20 (3), 337–344. doi:10.1038/mp.2014.43
Maussion, G., Cruceanu, C., Rosenfeld, J. A., Bell, S. C., Jollant, F., Szatkiewicz, J., et al. (2017). Implication of LRRC4C and DPP6 in Neurodevelopmental Disorders. Am. J. Med. Genet. 173 (2), 395–406. doi:10.1002/ajmg.a.38021
McMahon, A. C., Barnett, M. W., O'Leary, T. S., Stoney, P. N., Collins, M. O., Papadia, S., et al. (2012). SynGAP Isoforms Exert Opposing Effects on Synaptic Strength. Nat. Commun. 3, 900. doi:10.1038/ncomms1900
Melo, U. S., Schöpflin, R., Acuna-Hidalgo, R., Mensah, M. A., Fischer-Zirnsak, B., Holtgrewe, M., et al. (2020). Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases. Am. J. Hum. Genet. 106 (6), 872–884. doi:10.1016/j.ajhg.2020.04.016
Miller, D. T., Adam, M. P., Aradhya, S., Biesecker, L. G., Brothman, A. R., Carter, N. P., et al. (2010). Consensus Statement: Chromosomal Microarray Is a First-Tier Clinical Diagnostic Test for Individuals with Developmental Disabilities or Congenital Anomalies. Am. J. Hum. Genet. 86 (5), 749–764. doi:10.1016/j.ajhg.2010.04.006
Piluso, G., Monteleone, P., Galderisi, S., Giugliano, T., Bertolino, A., Rocca, P., et al. (2019). Assessment of De Novo Copy-Number Variations in Italian Patients with Schizophrenia: Detection of Putative Mutations Involving Regulatory Enhancer Elements. World J. Biol. Psychiatry 20 (2), 126–136. doi:10.1080/15622975.2017.1395072
Pizzo, L., Jensen, M., Polyak, A., Rosenfeld, J. A., Mannik, K., Krishnan, A., et al. (2019). Rare Variants in the Genetic Background Modulate Cognitive and Developmental Phenotypes in Individuals Carrying Disease-Associated Variants. Genet. Med. 21 (4), 816–825. doi:10.1038/s41436-018-0266-3
Rauch, A., Wieczorek, D., Graf, E., Wieland, T., Endele, S., Schwarzmayr, T., et al. (2012). Range of Genetic Mutations Associated with Severe Non-syndromic Sporadic Intellectual Disability: an Exome Sequencing Study. The Lancet 380 (9854), 1674–1682. doi:10.1016/S0140-6736(12)61480-9
Riggs, E. R., Andersen, E. F., Cherry, A. M., Kantarci, S., Kearney, H., Patel, A., et al. (2020). Technical Standards for the Interpretation and Reporting of Constitutional Copy-Number Variants: a Joint Consensus Recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet. Med. 22 (2), 245–257. doi:10.1038/s41436-019-0686-8
Ross, P. J., Zhang, W.-B., Mok, R. S. F., Zaslavsky, K., Deneault, E., D’Abate, L., et al. (2020). Synaptic Dysfunction in Human Neurons with Autism-Associated Deletions in PTCHD1-AS. Biol. Psychiatry 87 (2), 139–149. doi:10.1016/j.biopsych.2019.07.014
Rossi, P. I. A., Musante, I., Summa, M., Pittaluga, A., Emionite, L., Ikehata, M., et al. (2013). Compensatory Molecular and Functional Mechanisms in Nervous System of the Grm1crv4 Mouse Lacking the mGlu1 Receptor: A Model for Motor Coordination Deficits. Cereb. Cortex 23 (9), 2179–2189. doi:10.1093/cercor/bhs200
Rumbaugh, G., Adams, J. P., Kim, J. H., and Huganir, R. L. (2006). SynGAP Regulates Synaptic Strength and Mitogen-Activated Protein Kinases in Cultured Neurons. Proc. Natl. Acad. Sci. 103 (12), 4344–4351. doi:10.1073/pnas.0600084103
Sanders, S. J., Ercan-Sencicek, A. G., Hus, V., Luo, R., Murtha, M. T., Moreno-De-Luca, D., et al. (2011). Multiple Recurrent De Novo CNVs, Including Duplications of the 7q11.23 Williams Syndrome Region, Are Strongly Associated with Autism. Neuron 70 (5), 863–885. doi:10.1016/j.neuron.2011.05.002
Savatt, J. M., and Myers, S. M. (2021). Genetic Testing in Neurodevelopmental Disorders. Front. Pediatr. 9, 526779. doi:10.3389/fped.2021.526779
Schmitt, A. D., Hu, M., and Ren, B. (2016). Genome-wide Mapping and Analysis of Chromosome Architecture. Nat. Rev. Mol. Cel Biol. 17 (12), 743–755. doi:10.1038/nrm.2016.104
Schmitz, L. J. M., Klaassen, R. V., Ruiperez-Alonso, M., Zamri, A. E., Stroeder, J., Rao-Ruiz, P., et al. (2017). The AMPA Receptor-Associated Protein Shisa7 Regulates Hippocampal Synaptic Function and Contextual Memory. Elife 6. doi:10.7554/eLife.24192
Schormair, B., Kemlink, D., Roeske, D., Eckstein, G., Xiong, L., Lichtner, P., et al. (2008). PTPRD (Protein Tyrosine Phosphatase Receptor Type delta) Is Associated with Restless Legs Syndrome. Nat. Genet. 40 (8), 946–948. doi:10.1038/ng.190
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Shihab, H. A., Rogers, M. F., Campbell, C., and Gaunt, T. R. (2017). HIPred: an Integrative Approach to Predicting Haploinsufficient Genes. Bioinformatics 33 (12), btx028–1757. doi:10.1093/bioinformatics/btx028
Srivastava, S., Love-Nichols, J. A., Love-Nichols, J. A., Dies, K. A., Ledbetter, D. H., Martin, C. L., et al. (2019). Meta-analysis and Multidisciplinary Consensus Statement: Exome Sequencing Is a First-Tier Clinical Diagnostic Test for Individuals with Neurodevelopmental Disorders. Genet. Med. 21 (11), 2413–2421. doi:10.1038/s41436-019-0554-6
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131
Tabas-Madrid, D., Nogales-Cadenas, R., and Pascual-Montano, A. (2012). GeneCodis3: a Non-redundant and Modular Enrichment Analysis Tool for Functional Genomics. Nucleic Acids Res. 40, W478–W483. doi:10.1093/nar/gks402
Tao, Z., Feng, C., Mao, C., Ren, J., Tai, Y., Guo, H., et al. (2019). MiR-4465 Directly Targets PTEN to Inhibit AKT/mTOR Pathway-Mediated Autophagy. Cell Stress and Chaperones 24 (1), 105–113. doi:10.1007/s12192-018-0946-6
Tassano, E., Gimelli, S., Divizia, M. T., Lerone, M., Vaccari, C., Puliti, A., et al. (2015). Thrombocytopenia-absent Radius (TAR) Syndrome Due to Compound Inheritance for a 1q21.1 Microdeletion and a Low-Frequency Noncoding RBM8A SNP: a New Familial Case. Mol. Cytogenet. 8, 87. doi:10.1186/s13039-015-0188-6
Tomita, H., Cornejo, F., Aranda-Pino, B., Woodard, C. L., Rioseco, C. C., Neel, B. G., et al. (2020). The Protein Tyrosine Phosphatase Receptor Delta Regulates Developmental Neurogenesis. Cel Rep. 30 (1), 215–228. e215. doi:10.1016/j.celrep.2019.11.033
Tomoda, T., Kim, J. H., Zhan, C., and Hatten, M. E. (2004). Role of Unc51.1 and its Binding Partners in CNS Axon Outgrowth. Genes Development 18 (5), 541–558. doi:10.1101/gad.1151204
Tropeano, M., Howley, D., Gazzellone, M. J., Wilson, C. E., Ahn, J. W., Stavropoulos, D. J., et al. (2016). Microduplications at the pseudoautosomalSHOXlocus in Autism Spectrum Disorders and Related Neurodevelopmental Conditions. J. Med. Genet. 53 (8), 536–547. doi:10.1136/jmedgenet-2015-103621
Uetani, N., Kato, K., Ogura, H., Mizuno, K., Kawano, K., Mikoshiba, K., et al. (2000). Impaired Learning with Enhanced Hippocampal Long-Term Potentiation in PTPdelta-Deficient Mice. EMBO J. 19 (12), 2775–2785. doi:10.1093/emboj/19.12.2775
Uhl, G. R., and Martinez, M. J. (2019). PTPRD: Neurobiology, Genetics, and Initial Pharmacology of a Pleiotropic Contributor to Brain Phenotypes. Ann. N.Y. Acad. Sci. 1451 (1), 112–129. doi:10.1111/nyas.14002
Um, J. W., and Ko, J. (2017). Neural Glycosylphosphatidylinositol-Anchored Proteins in Synaptic Specification. Trends Cel Biol. 27 (12), 931–945. doi:10.1016/j.tcb.2017.06.007
Vaccari, C. M., Romanini, M. V., Musante, I., Tassano, E., Gimelli, S., Divizia, M. T., et al. (2014). De Novo deletion of Chromosome 11q12.3 in Monozygotic Twins Affected by Poland Syndrome. BMC Med. Genet. 15, 63. doi:10.1186/1471-2350-15-63
Velinov, M. (2019). Genomic Copy Number Variations in the Autism Clinic-Work in Progress. Front. Cell. Neurosci. 13, 57. doi:10.3389/fncel.2019.00057
Vijiaratnam, N., Bhatia, K. P., Lang, A. E., Raskind, W. H., and Espay, A. J. (2019). ADCY5-Related Dyskinesia: Improving Clinical Detection of an Evolving Disorder. Mov Disord. Clin. Pract. 6 (7), 512–520. doi:10.1002/mdc3.12816
von Engelhardt, J. (2019). AMPA Receptor Auxiliary Proteins of the CKAMP Family. Ijms 20 (6), 1460. doi:10.3390/ijms20061460
Wu, X., Li, W., and Zheng, Y. (2020). Recent Progress on Relevant microRNAs in Autism Spectrum Disorders. Ijms 21 (16), 5904. doi:10.3390/ijms21165904
Yamagata, A., Yoshida, T., Sato, Y., Goto-Ito, S., Uemura, T., Maeda, A., et al. (2015). Mechanisms of Splicing-dependent Trans-synaptic Adhesion by PTPδ-IL1RAPL1/IL-1RAcP for Synaptic Differentiation. Nat. Commun. 6, 6926. doi:10.1038/ncomms7926
Yan, H., Shi, Z., Wu, Y., Xiao, J., Gu, Q., Yang, Y., et al. (2019). Targeted Next Generation Sequencing in 112 Chinese Patients with Intellectual Disability/developmental Delay: Novel Mutations and Candidate Gene. BMC Med. Genet. 20 (1), 80. doi:10.1186/s12881-019-0794-y
Yang, A. W., Sachs, A. J., and Nystuen, A. M. (2015). Deletion of Inpp5a Causes Ataxia and Cerebellar Degeneration in Mice. Neurogenetics 16 (4), 277–285. doi:10.1007/s10048-015-0450-4
Yang, J.-H., Chang, M.-W., Pandey, P. R., Tsitsipatis, D., Yang, X., Martindale, J. L., et al. (2020). Interaction of OIP5-AS1 with MEF2C mRNA Promotes Myogenic Gene Expression. Nucleic Acids Res. 48 (22), 12943–12956. doi:10.1093/nar/gkaa1151
Yoshida, T., Shiroshima, T., Lee, S.-J., Yasumura, M., Uemura, T., Chen, X., et al. (2012). Interleukin-1 Receptor Accessory Protein Organizes Neuronal Synaptogenesis as a Cell Adhesion Molecule. J. Neurosci. 32 (8), 2588–2600. doi:10.1523/JNEUROSCI.4637-11.2012
Keywords: variants of uncertain significance, long non-coding genes, antisense gene, oligogenic disease, array-CGH, additive effect, microdeletion, microduplication
Citation: Servetti M, Pisciotta L, Tassano E, Cerminara M, Nobili L, Boeri S, Rosti G, Lerone M, Divizia MT, Ronchetto P and Puliti A (2021) Neurodevelopmental Disorders in Patients With Complex Phenotypes and Potential Complex Genetic Basis Involving Non-Coding Genes, and Double CNVs. Front. Genet. 12:732002. doi: 10.3389/fgene.2021.732002
Received: 28 June 2021; Accepted: 03 September 2021;
Published: 21 September 2021.
Edited by:
Anjana Munshi, Central University of Punjab, IndiaReviewed by:
Patrick May, University of Luxembourg, LuxembourgGavin R. Oliver, Mayo Clinic, United States
Copyright © 2021 Servetti, Pisciotta, Tassano, Cerminara, Nobili, Boeri, Rosti, Lerone, Divizia, Ronchetto and Puliti. 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: Aldamaria Puliti, YXB1bGl0aUB1bmlnZS5pdA==
†These authors have contributed equally to this work