The 15q11.2 BP1-BP2 Microdeletion (Burnside–Butler) Syndrome: In Silico Analyses of the Four Coding Genes Reveal Functional Associations with Neurodevelopmental Disorders

The 15q11.2 BP1-BP2 microdeletion (Burnside–Butler) syndrome is emerging as the most frequent pathogenic copy number variation (CNV) in humans associated with neurodevelopmental disorders with changes in brain morphology, behavior, and cognition. In this study, we explored functions and interactions of the four protein-coding genes in this region, namely NIPA1, NIPA2, CYFIP1, and TUBGCP5, and elucidate their role, in solo and in concert, in the causation of neurodevelopmental disorders. First, we investigated the STRING protein-protein interactions encompassing all four genes and ascertained their predicted Gene Ontology (GO) functions, such as biological processes involved in their interactions, pathways and molecular functions. These include magnesium ion transport molecular function, regulation of axonogenesis and axon extension, regulation and production of bone morphogenetic protein and regulation of cellular growth and development. We gathered a list of significantly associated cardinal maladies for each gene from searchable genomic disease websites, namely MalaCards.org: HGMD, OMIM, ClinVar, GTR, Orphanet, DISEASES, Novoseek, and GeneCards.org. Through tabulations of such disease data, we ascertained the cardinal disease association of each gene, as well as their expanded putative disease associations. This enabled further tabulation of disease data to ascertain the role of each gene in the top ten overlapping significant neurodevelopmental disorders among the disease association data sets: (1) Prader–Willi Syndrome (PWS); (2) Angelman Syndrome (AS); (3) 15q11.2 Deletion Syndrome with Attention Deficit Hyperactive Disorder & Learning Disability; (4) Autism Spectrum Disorder (ASD); (5) Schizophrenia; (6) Epilepsy; (7) Down Syndrome; (8) Microcephaly; (9) Developmental Disorder, and (10) Peripheral Nervous System Disease. The cardinal disease associations for each of the four contiguous 15q11.2 BP1-BP2 genes are NIPA1- Spastic Paraplegia 6; NIPA2—Angelman Syndrome and Prader–Willi Syndrome; CYFIP1—Fragile X Syndrome and Autism; TUBGCP5—Prader–Willi Syndrome. The four genes are individually associated with PWS, ASD, schizophrenia, epilepsy, and Down syndrome. Except for TUBGCP5, the other three genes are associated with AS. Unlike the other genes, TUBGCP5 is also not associated with attention deficit hyperactivity disorder and learning disability, developmental disorder, or peripheral nervous system disease. CYFIP1 was the only gene not associated with microcephaly but was the only gene associated with developmental disorders. Collectively, all four genes were associated with up to three-fourths of the ten overlapping neurodevelopmental disorders and are deleted in this most prevalent known pathogenic copy number variation now recognized among humans with these clinical findings.

. 15q11.2 BP1-BP2 microdeletion (Burnside-Butler) syndrome region found at the proximal end of Prader-Willi syndrome (PWS) / Angelman syndrome (AS) regions within the 15q11-q13 Type I deletion depicting the location and order of the four protein-coding genes therein: NIPA1, NIPA2, CYFIP1, and TUBGCP5 within the 15q11.2 region distal to the centromere and proximal to the imprinted PWS/AS genes. The bottom enlarged horizontal chromatin figure exclusively depicts NIPA1, NIPA2, CYFIP1, and TUBGCP5 genes in the BP1-BP2 region.
The location, as well as the order of genes and transcripts (e.g., snoRNAs), are also shown in Figure 1 that causes PWS/AS. Those genes that are imprinted and paternally expressed causing PWS (in blue) and those causing AS (in red) are also imprinted and maternally expressed. The location and the relative size of the 15q11.2 BP1-BP2 microdeletion region, the typical larger 15q11-q13 Type I deletion involving breakpoints BP1 and BP3, and the typical smaller 15q11-q13 Type II deletion involving breakpoints BP2 and BP3 in both PWS and AS are illustrated. IC: imprinting center controlling the activity of imprinted genes in the 15q11-q13 region.

Overview of the Four Genes in the 15q11.2 BP1-BP2 Region
All four syntenic and bi-allelically conserved expressed genes in the 15q11.2 region between breakpoints BP1 and BP2 are functionally predicted to interact with each other along with seven other genes. The predicted STRING functional interactions network encompasses these four genes as illustrated in Figure 2.
The STRING diagram and protein interactions involving all four genes of this microdeletion syndrome are depicted in Figure 2 and their predicted functions are presented in Table 1. The identified biological processes (GO) and molecular functions (GO) are presented in Table 2 (https: //version11.stringdb.org/cgi/network.pl?taskId=XJOUKPc8icHM).
Protein Network (STRING) Interactions Encompassing the Four Genes: NIPA1, NIPA2, CYFIP1 and TUBGCP5 genes. The predicted STRING functional interactions network encompasses these four genes as illustrated in Figure 2, while the individual predicted functions are shown in Table 1. The related biological processes (GO) and molecular functions (GO) are shown in Table 2.
Protein Network (STRING) Interactions Encompassing the Four Genes: NIPA1, NIPA2, CYFIP1, and TUBGCP5 The STRING diagram and protein interactions involving all four genes of this microdeletion syndrome are depicted in Figure 2 above and their predicted functions are presented in Table 1. The identified biological processes (GO) and molecular functions (GO) are presented in Table 2 (https://version11.stringdb.org/cgi/network.pl?taskId=XJOUKPc8icHM).  (Table 1) and 34 edges and predicted functional interactions, such as, biological process (GO) and molecular function (GO) are given in Table 2, as designated in STRING: 9606.ENSP00000337452 (STRING Consortium 2019). Network nodes represent proteins (Table 1) with splice isoforms or post-translational modifications collapsed into each node for all proteins produced by a single protein-coding gene. Edges represent protein-protein associations that are considered specific and meaningful, i.e., proteins jointly contribute to a shared function ( Table 2).  (Table 1) and 34 edges and predicted functional interactions, such as, biological process (GO) and molecular function (GO) ( Table 2), as designated in STRING: 9606.ENSP00000337452 (STRING Consortium 2019). Network nodes represent proteins (Table 1) with splice isoforms or post-translational modifications collapsed into each node for all proteins produced by a single protein-coding gene. Edges represent protein-protein associations that are considered specific and meaningful, i.e., proteins jointly contribute to a shared function, , such as, Biological Process (GO) and Molecular Function (GO) ( Table 2).

Co-expression Experiments Databases
Text mining [Homology]

TUBGCP5
Gamma-tubulin complex component 5; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

CYFIP1
Cytoplasmic FMR1-interacting protein 1; Component of the CYFIP1-EIF4E-FMR1 complex which binds to the mRNA cap and mediates translational repression. In the CYFIP1-EIF4E-FMR1 complex this subunit is an adapter between EIF4E and FMR1. Promotes the translation repression activity of FMR1 in brain probably by mediating its association with EIF4E and mRNA (By similarity). Regulates formation of membrane ruffles and lamellipodia. Plays a role in axon outgrowth. Binds to F-actin but not to RNA. Part of the WAVE complex that regulates actin filament reorganization.  ). NIPA1 has two alternative splice isoforms; isoform 2 differs from the canonical sequence by missing the first 75 amino acids (GeneCards & UniProtKB/Swiss-Prot). It is a multi-pass membrane protein recruited to the cell membrane in response to low extracellular magnesium (GeneCards & UniProtKB/Swiss-Prot). This gene is widely expressed with highest levels in neuronal tissues and overexpressed in the spinal cord (20.8), frontal cortex (17.0), and fetal liver (7.3) (Protein differential expression in normal tissues from Human Integrated Protein Expression Database (HIPED) for NIPA1 Gene: https://www.genecards.org/Guide/GeneCard#protein-differential-expression; https://www.genecards.org/cgi-bin/carddisp.pl?gene=NIPA1&keywords=NIPA1). NIPA1 protein plays a potential role in nervous system development and maintenance and is most ubiquitously expressed in the brain with a mean RPKM (Reads Per Kilobase per Million reads placed) of 16.56 + 3.055 (https://www.ncbi.nlm.nih.gov/gene/123606/?report=expression; https: //www.gtexportal.org/home/gene/NIPA1) [14]. The multiple transmembrane domains of this protein localize to endosomes and plasma membrane. This protein is recruited in response to low extracellular magnesium and functions in Mg2 transport (Q7RTP0-NIPA1_HUMAN) [15]. Subcellular localization with immunofluorescence shows that endogenous NIPA1 protein associates with early endosomes and the cell surface in a variety of neuronal and epithelial cells. As expected of a magnesium-responsive gene, altered magnesium concentration leads to redistribution between the endosomal compartment and the plasma membrane; high magnesium results in diminished cell surface NIPA1 whereas low magnesium leads to accumulation in early endosomes and recruitment to the plasma membrane [15].

YES
An important paralog of this gene, NIPAL1 is also a magnesium ion transmembrane transporter. In addition, there are three other protein-coding transcript variants (paralogs): NIPAL2, NIPAL3, and NIPAL4, which are also magnesium ion transmembrane transporters. Other genes also perform a similar magnesium transport function, namely, MAGT1, MMGT1, and MRS2 (GeneCards.Org). Although NIPA1 protein acts as a Mg (2+) transporter, it can also transport other divalent cations such as Fe(2+), Sr(2+), Ba(2+), Mn(2+), and Co(2+), but to a much less extent than Mg(2+) (UniProtKB/Swiss-Prot Summary). Among its other related pathways are the transport of glucose and other sugars, bile salts and organic acids, metal ions, and amine compounds.
Only mutations (not haploinsufficiency due to 15q11.2 BP1-BP2 microdeletions including the NIPA1 gene) have thus far been associated with autosomal dominant spastic paraplegia 6 / SPG6linked hereditary spastic paraplegia: HSP (https://www.genecards.org/cgi-bin/carddisp.pl?gene= NIPA1&keywords=NIPA1). Spastic paraplegia 6 is a neurodegenerative disorder characterized by slow, gradual, progressive weakness and spasticity of the lower limbs. Rate of progression and the severity of symptoms are quite variable with initial symptoms of difficulty with balance, weakness, and stiffness in the legs, muscle spasms, and dragging the toes when walking. In some forms of the disorder, bladder symptoms (such as incontinence) may appear, or weakness and stiffness may be present to other parts of the body (UniProtKB/Swiss-Prot) [16][17][18]. Two variants in the NIPA1 gene, namely, p.Thr45Arg (VAR_023440; SNP ID: rs104894496) and p.Gly106Arg (VAR_023441; SNP ID: rs104894490) are disease-causing variations for spastic paraplegia type 6 (SPG6), autosomal dominant form. A rare NIPA1 deletion was also found in a patient with pervasive developmental disorder not otherwise specified and with mild intellectual disability [19]; this deletion is also linked to autism (SFARI.org).
The mouse NIPA1 mutants, p.Thr39Arg and p.Gly100Arg, corresponding to the respective human mutants are associated with hereditary spastic paraplegia (HSP) showing a loss-of-function when expressed in oocytes and altered trafficking in transfected COS7 cells [15]. The NIPA1 gene normally encodes an Mg2+ transporter protein and the loss-of-function of NIPA1, due to abnormal trafficking of the mutated protein provides the basis for the HSP phenotype [15,16,20]. Only abnormal trafficking of the mutated protein is causative for HSP and not from its deficiency as evident by Prader-Willi syndrome (PWS) or Angelman syndrome (AS) with only one copy of the NIPA1 gene in the typical 15q11-q13 type I deletion in both PWS or AS involving breakpoints BP1 and BP3 including the four genes in the 15q11.2 BP1-BP2 region. Other genes or transcripts, both imprinted and biallelic, do not have HSP [16,21]. Key findings in PWS and AS develop from errors of imprinting dependent on the parent of origin [22,23]. Therefore, SPG6 linked HSP is more likely caused by a dominant-negative effect or a toxic gain-of-function mechanism rather than a loss-of-function [16,17,21,24,25].
A perusal of genes involved in numerous other forms of spastic paraplegia indicates causation by genetic defects in various other metabolic pathways unrelated to magnesium transporter protein deficiency (see GeneCards.org). The syntenic NIPA2 gene, which also codes a magnesium ion transmembrane transporter protein and along with NIPA1 in the same region are both deleted in the Type I deletion but not deleted in the smaller typical 15q11-q13 Type II deletion in PWS or AS. However, patients with the larger Type I deletion do not have spastic paraplegia (GeneCards.org). Mutations in the NIPA2 gene have been reported in generalized epilepsy and childhood absence epilepsy [10,11,26]. In addition to the magnesium ion transmembrane transporter activity, NIPA1 protein inhibits Bone Morphogenetic Protein signaling by regulating the endosomal trafficking and degradation of type 2 BMP receptors (BMPR2) in Drosophila and HeLa cells [27,28].
Bone Morphogenetic Proteins (BMPs) are a group of signaling molecules that belong to the Transforming Growth Factor-β (TGFβ) superfamily of proteins. Initially discovered for their ability to induce bone formation, BMPs are now known to play crucial roles in all organ systems. BMPs are important in embryogenesis and development, and also in the maintenance of adult tissue homeostasis. Another relevant role of BMP in the neurological system is neurogenesis. Neural defects are associated with loss of BMP function in mouse models [29]. For instance, BMP11 is involved in spinal cord neurogenesis and secreted from neurons themselves serving as an inhibitory signal in the generation of new neurons from progenitors in the olfactory epithelium [30].
In rodents, BMP signaling is upregulated following lesions of the corticospinal track and suppressing upregulation that promotes regrowth of axons [31]. Targeting specific BMP receptor subunits for therapeutic purposes may provide an approach for manipulating gliosis and enhancing functional outcomes after spinal cord injury [32]. Therefore, it is no surprise that BMP signaling cuts across all hereditary spastic paraplegia (HSP) categories and are widely implicated in neurodegenerative diseases [33]. Among the HSP proteins, NIPA1 is best characterized mechanistically to inhibit BMP signaling [15]. In summary, these studies support abnormal BMP signaling in many cases probably resulting from abnormal BMP receptor trafficking and could be a unifying pathogenic mechanism for some forms of hereditary spastic paraplegia [34].
Further highlighting the importance of BMP regulation in neural development is the role of BMP7 in corticogenesis; BMP7 deletions result in reduced cortical thickening and impaired neurogenesis [35]. Interestingly, BMPR type 1A (BMPR1A) is important in the establishment of neurons involved in regulating feeding behavior [32]. Similarly, NIPA1 was shown to inhibit BMP signaling by regulating the endosomal trafficking and degradation of type 2 BMP receptors (BMPR2) in Drosophila and HeLa cells [27,28]. Perhaps, haploinsufficiency of NIPA1 protein in cases of 15q11-q13 BP1 or BP2/BP3 deletions might partially affect its inhibition of BMP signaling by regulating endosomal trafficking and degradation of type 2 BMP receptors, thereby affecting neurons. This might be the causal factor for the reported developmental and language delay, neurobehavioral disturbances, and psychiatric problems, such as autism, seizures, and schizophrenia with occasional mild dysmorphic features seen to a varying degree in patients with the 15q11.2 BP1-BP2 deletion (Burnside-Butler) syndrome [1,36]. Biological processes (GO), molecular functions (GO), cellular components with KEGG, and Reactome pathways are summarized in Table 3. In Table Table 4 examples of diseases or disorders which have been recognized and associated with the NIPA1 gene when disturbed are collated and ranked by order using MalaCards (MalaCards.org) from GeneCards (GeneCards.org). STRING functional interactions with nodes and edges (see Figure 3). Further illustrations of interactions and functions of NIPA1 with relationship to other proteins can be found in Table Table 5.    The interacting proteins with NIPA1 using the STRING interactions network are illustrated in Figure 3 with descriptions of their nodes and edges. Further illustrations of interactions and functions of NIPA1 with relationship to other proteins can be found in Table 5. Network nodes represent proteins with splice isoforms or posttranslational modifications collapsed into each node for all proteins produced by a single proteincoding gene. Edges represent protein-protein associations that are considered specific and meaningful (i.e., proteins jointly contribute to a shared function (Table 3). STRING interactants and their functions related to other proteins associated with NIPA1 gene is listed in Table 5. (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins). Network nodes represent proteins with splice isoforms or post-translational modifications collapsed into each node for all proteins produced by a single protein-coding gene. Edges represent protein-protein associations that are considered specific and meaningful (i.e., proteins jointly contribute to a shared function, such as, Biological Process (GO) and Molecular Function (GO) ( Table 3). STRING interactants and their functions related to other proteins associated with NIPA1 gene are listed in Table 5. (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins).

TGFB1
Transforming growth factor beta 1; Multifunctional protein that controls proliferation, differentiation and other functions in many cell types.

BMP4
Bone morphogenetic protein-4; Induces cartilage and bone formation. Also act in mesoderm induction, tooth development, limb formation and fracture repair.

BMP8B
Bone morphogenetic protein 8B; Induces cartilage and bone formation. Plays a role in calcium regulation and bone homeostasis; Belongs to the TGF-beta family.

MUSK
Muscle, skeletal receptor tyrosine-protein kinase; Receptor tyrosine kinase which plays a central role in the formation and the maintenance of the neuromuscular junction (NMJ), the synapse between the motor neuron and the skeletal muscle.

BMPR2
Bone morphogenetic protein receptor type 2; On ligand binding, forms a receptor complex consisting of two type II and two type I transmembrane serine/threonine kinases.,

BMP2
Bone morphogenetic protein 2; Induces cartilage and bone formation. Stimulates the differentiation of myoblasts into osteoblasts via the EIF2AK3-EIF2A-ATF4 pathway

BMP7
Bone morphogenetic protein 7; Induces cartilage and bone formation. May be the osteoinductive factor responsible for the phenomenon of epithelial osteogenesis. Plays a role in calcium regulation and bone homeostasis.

BMPR1B
Bone morphogenetic protein receptor type-1B; On ligand binding, forms a receptor complex consisting of two type II and two type I transmembrane serine/threonine kinases.

GDF2
Growth/differentiation factor 2; Potent circulating inhibitor of angiogenesis. Signals through the type I activin receptor ACVRL1. Signaling through SMAD1 in endothelial cells requires TGF-beta coreceptor endoglin.
The chromosomal band location for the NIPA2 gene is 15q11.2 with genomic location (GRCh38/hg38) at 22,838,641-22,868,384 with a size of 29,744 bases in a plus-strand orientation. This gene is protein-coding specifically for NIPA magnesium transporter 2 (https://www.genecards.org/cgi-bin/ carddisp.pl?gene=NIPA2&keywords=nipa2). When disturbed (deleted), this gene is associated with both Angelman and Prader-Willi syndromes due to specific parent of origin deletions of the 15q11-q13 region. The typical 15q11-q13 type I deletion also includes the other three protein-coding genes in the 15q11.2 BP1-BP2 region, namely, NIPA1, CYFIP1, and TUBGCP5. In the case of the smaller typical 15q11-q13 type II deletion, these four genes remain intact in the region ( Figure 1 The NIPA2 protein has the following attributes: Size: 360 amino acids; Molecular mass:39185 Da. Quaternary structure: No Data Available. This protein has two isoforms; isoform 2 differs from the canonical sequence as follows: 47-66: GQGGHAYLKEWLWWAGLLSM → V. Similar to NIPA1, NIPA2 is also a multi-pass membrane protein found in early endosomes recruited to the cell membrane in response to low extracellular magnesium (UniProtKB/Swiss-Prot for NIPA2 Gene) ( Table 6).
Unlike NIPA1, this gene is overexpressed in B-lymphocytes and the placenta ((UniProtKB/Swiss-Prot). In addition to PWS and AS involvement for this gene, it is also associated with childhood absence epilepsy [26], childhood electroclinical syndrome [10,11], and possibly autosomal recessive congenital ichthyosis (MalaCards/GeneCards) and important biological processes and functions (see Table 7). Among related pathways of this gene are miscellaneous transport and binding events and transport of glucose and other sugars, bile salts and organic acids, metal ions and amine compounds. Gene Ontology (GO) annotations related to this gene also include magnesium ion transmembrane transporter activity. An important paralog of this gene is NIPAL1 (GeneCards). Chai et al. [9] determined that the NIPA2 gene contains 7 exons and spans 29 kb. The coding region extends between exons 3 and 7, and alternative splicing utilizes alternate exons 2 and 2b which results in multiple transcript variants. Pseudogenes of this gene are found on chromosomes 3, 7, and 21 (Entrez Gene: 81614).
Similar to NIPA1, this multi-pass cell membrane protein is localized in the plasma membrane and early endosomes (Gene Ontology (GO)-GO:0005769; GO:0005886) and recruited to the cell membrane in response to low extracellular magnesium (GeneCards; ECO:0000250). Paralogs for the NIPA2 gene are NIPA1, NIPAL1, NIPAL4, NIPAL2, and NIPAL3 (PMID: 28514442; [35]). Although NIPA2 is highly conserved (PMID: 14508708; [9]) and expressed in parts of the brain as well as in most other organ systems, it is significantly overexpressed in B-lymphocytes and the placenta. In addition, five variants have been reported for this gene with either loss or gain of function or both (GeneCards). Table Table 8 describes the STRING protein interactions and their functions with this gene and description while Table 6 shows biological processes, molecular functions, cellular components, and KEGG and Reactome pathways attributed to the NIPA2 gene. To examine NIPA2 protein interactions and connectivity with other genes and their interrelated proteins, see the STRING interaction network displayed in Figure 4. To examine NIPA2 protein interactions and connectivity with other genes and their interrelated proteins, see the STRING interaction network displayed in Figure 4. Table 8 describes the STRING protein interactions and their functions with NIPA2 and description  Network nodes represent proteins with splice isoforms or post-translational modifications collapsed into each node for all proteins produced by a single protein-coding gene. Edges represent protein-protein associations that are considered specific and meaningful (i.e., proteins jointly contribute to a shared function, such as, Biological Process (GO) and Molecular Function (GO) ( Table 6). NIPA2 gene protein Interactants and their functions are listed in Table 8. (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins). STRING Protein Interactants and Their Functions for NIPA2 gene are listed in Table 8.

ACVR1C
Activin receptor type-1C; Serine/threonine protein kinase which forms a receptor complex on ligand binding.

SMAD2
Mothers against decapentaplegic homolog 2; also known as SMAD family member 2; Receptor-regulated SMAD (R-SMAD) that is an intracellular signal transducer and transcriptional modulator activated by TGF-beta (transforming growth factor) and activin type 1 receptor kinases.

ACVR1
Activin receptor type 1; On ligand binding, forms a receptor complex consisting of two type II and two type I transmembrane serine/threonine kinases.

UXS1
UDP-glucuronic acid decarboxylase 1; Catalyzes the NAD-dependent decarboxylation of UDP-glucuronic acid to UDP-xylose. Necessary for the biosynthesis of the core tetrasaccharide in glycosaminoglycan biosynthesis.

SMAD3
Mothers against decapentaplegic homolog 3; also known as SMAD family member 3; Receptor-regulated SMAD (R-SMAD) that is an intracellular signal transducer and transcriptional modulator activated by TGF-beta (transforming growth factor) and activin type 1 receptor kinases.

SMAD4
Mothers against decapentaplegic homolog 4; also known as SMAD family member 4; In muscle physiology, plays a central role in the balance between atrophy and hypertrophy.

CHUK
Inhibitor of nuclear factor kappa-B kinase subunit alpha; Serine kinase that plays an essential role in the NF-kappa-B signaling pathway.

BMPR1A
Bone morphogenetic protein receptor type-1A; On ligand binding, forms a receptor complex consisting of two type II and two type I transmembrane serine/threonine kinases.

BMPR2
Bone morphogenetic protein receptor type-2; On ligand binding, forms a receptor complex consisting of two type II and two type I transmembrane serine/threonine kinases.

TGFBR1
TGF-beta receptor type-1; Transmembrane serine/threonine kinase forming with the TGF-beta type II serine/threonine kinase receptor. Regulating a plethora of physiological and pathological processes including cell cycle arrest in epithelial and hematopoietic cells, control of mesenchymal cell proliferation and differentiation, wound healing, extracellular matrix production, and immunosuppression.

SMAD9
Mothers against decapentaplegic homolog 9; also known as SMAD family member 9; Transcriptional modulator activated by BMP (bone morphogenetic proteins) type 1 receptor kinase.

BMP7
Bone morphogenetic protein 7; Induces cartilage and bone formation. May be the osteoinductive factor responsible for the phenomenon of epithelial osteogenesis. Plays a role in calcium regulation and bone homeostasis; Bone morphogenetic proteins.

BMPR1B
Bone morphogenetic protein receptor type-1B; On ligand binding, forms a receptor complex consisting of two type II and two type I transmembrane serine/threonine kinases.

SMAD1
Mothers against decapentaplegic homolog 1; also known as SMAD family member 1; Transcriptional modulator activated by BMP (bone morphogenetic proteins) type 1 receptor kinase.

ACVR1B
Activin receptor type-1B; Transmembrane serine/threonine kinase activin type-1 receptor forming an activin receptor complex with activin receptor type-2 (ACVR2A or ACVR2B). Transduces the activin signal from the cell surface to the cytoplasm and is thus regulating many physiological and pathological processes including neuronal differentiation and neuronal survival, hair follicle development and cycling.
The chromosomal band location for the CYFIP1 gene is 15q11.2 with the genomic location (GRCh38/hg38) at 22,867,052-22,981,063 with a size of 114,012 bases in a minus strand orientation (Genecards.org). The CYFIP1 protein has the following attributes: Size: 1253 amino acids; Molecular mass: 145,182 Da. Quaternary structure: Component of the WAVE1 complex composed of ABI2, CYFIP1 or CYFIP2, BRK1, NCKAP1, and WASF1/WAVE1. CYFIP1 and CYFIP2 are part of the Wiskott-Aldrich syndrome protein-family verprolin-homologous protein (WAVE) complex that regulates actin polymerization at synapses and CYFIP1 protein has 3 described isoforms produced by alternative splicing (https://www.uniprot.org/uniprot/Q7L576) (see Table 9). This protein is required for neuronal and bristle development in Drosophila [37]. Through the CYFIP1 protein, the fragile X syndrome protein represses activity-dependent translation, and thus implicated in fragile X syndrome [38].
The CYFIP1 gene is highly expressed in the perinuclear region and enriched in synaptosomes. It is also enriched in membrane ruffles and at the tips of lamellipodia with the following subcellular localizations (GeneCards & UniProtKB/Swiss-Prot). The CYFIP1 gene is widely expressed in all tissues (GenecCards). In addition to this gene's association with fragile X syndrome, it is also associated with autism (see Table 10). A large chromosomal deletion including this gene is associated with increased risk of schizophrenia and epilepsy in human patients and reduced expression of this gene has also been observed in various human cancers, the encoded protein may inhibit tumor invasion (Entrez Gene summary: RefSeq, May 2017). Among this gene's related pathways are the regulation of actin dynamics for phagocytic cup formation and signaling by Rho GTPases (Table 9). Gene Ontology (GO) annotations related to this gene include Rac GTPase binding. An important paralog of this gene is CYFIP2 (GeneCards). Table Table 10 lists putative associated diseases for CYFIP1 Gene. STRING protein interactions with their functions and descriptions related to CYFIP1 are presented in Table 11. An illustration of the identified interactions can also be visualized in Figure 5.  Protein-Protein Inter-Relationships Involving CYFIP1 Gene. Network nodes represent proteins with splice isoforms or posttranslational modifications collapsed into each node for all proteins produced by a single proteincoding gene. Edges represent protein-protein associations that are considered specific and meaningful (i.e., proteins jointly contribute to a shared function) ( Table 9). (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins). CYFIP1 gene protein-protein interactants and their functions for CYFIP1 are listed in Table  11.

CYFIP1
Cytoplasmic FMR1-interacting protein 1; Component of the CYFIP1-EIF4E-FMR1 complex which binds to the mRNA cap and mediates translational repression. In the CYFIP1-EIF4E-FMR1 complex this subunit is an adapter between EIF4E and FMR1. Promotes the translation repression activity of FMR1 in brain probably by mediating its association with EIF4E and mRNA (By similarity). Regulates formation of membrane ruffles and lamellipodia. Plays a role in axon outgrowth. Binds to F-actin but not to RNA. Part of the WAVE complex that regulates actin filament reorganization.

VCL
Vinculin; Actin filament (F-actin)-binding protein involved in cell-matrix adhesion and cell-cell adhesion.
Regulates cell-surface E-cadherin expression and potentiates mechanosensing by the E-cadherin complex. May also play important roles in cell morphology and locomotion; Belongs to the vinculin/alpha-catenin family.

MAPK1
Mitogen-activated protein kinase 1; Serine/threonine kinase which acts as an essential component of the MAP kinase signal transduction pathway. MAPK1/ERK2 and MAPK3/ERK1 are the 2 MAPKs which play an important role in the MAPK/ERK cascade. Depending on the cellular context, the MAPK/ERK cascade mediates diverse biological functions such as cell growth, adhesion, survival and differentiation through the regulation of transcription, translation, cytoskeletal rearrangements.

OLFM4
Olfactomedin-4; May promote proliferation of pancreatic cancer cells by favoring the transition from the S to G2/M phase. In myeloid leukemic cell lines, inhibits cell growth and induces cell differentiation and apoptosis.

CTSH
Pro-cathepsin H; Important for the overall degradation of proteins in lysosomes; Belongs to the peptidase C1 family.

RETN
Resistin; Hormone that seems to suppress insulin ability to stimulate glucose uptake into adipose cells (By similarity). Potentially links obesity to diabetes (By similarity). Promotes chemotaxis in myeloid cells. Network nodes represent proteins with splice isoforms or post-translational modifications collapsed into each node for all proteins produced by a single protein-coding gene. Edges represent protein-protein associations that are considered specific and meaningful, i.e., proteins jointly contribute to a shared function, such as, Biological Process (GO) and Molecular Function (GO) ( Table 9). CYFIP1 gene protein-protein interactants and their functions are listed in Table 11. (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins).

CYFIP1
Cytoplasmic FMR1-interacting protein 1; Component of the CYFIP1-EIF4E-FMR1 complex which binds to the mRNA cap and mediates translational repression. In the CYFIP1-EIF4E-FMR1 complex this subunit is an adapter between EIF4E and FMR1. Promotes the translation repression activity of FMR1 in brain probably by mediating its association with EIF4E and mRNA (By similarity). Regulates formation of membrane ruffles and lamellipodia. Plays a role in axon outgrowth. Binds to F-actin but not to RNA. Part of the WAVE complex that regulates actin filament reorganization.

OLFM4
Olfactomedin-4; May promote proliferation of pancreatic cancer cells by favoring the transition from the S to G2/M phase. In myeloid leukemic cell lines, inhibits cell growth and induces cell differentiation and apoptosis.

CTSH
Pro-cathepsin H; Important for the overall degradation of proteins in lysosomes; Belongs to the peptidase C1 family.

RETN
Resistin; Hormone that seems to suppress insulin ability to stimulate glucose uptake into adipose cells (By similarity). Potentially links obesity to diabetes (By similarity). Promotes chemotaxis in myeloid cells.

ARPC3
Actin-related protein 2/3 complex subunit 3; Functions as component of the Arp2/3 complex which is involved in regulation of actin polymerization and together with an activating nucleation-promoting factor (NPF) mediates the formation of branched actin networks.

LTF
Lactotransferrin; Lactoferroxins A, B and C have opioid antagonist activity.

CTSD
Cathepsin D; Acid protease active in intracellular protein breakdown. Involved in the pathogenesis of several diseases such as breast cancer and possibly Alzheimer disease; Cathepsins.

MMP8
Neutrophil collagenase; Can degrade fibrillar type I, II, and III collagens; Belongs to the peptidase M10A family.

WIPF3
WAS/WASL-interacting protein family member 3; May be a regulator of cytoskeletal organization. May have a role in spermatogenesis (By similarity); Belongs to the verprolin family.

TNFAIP6
Tumor necrosis factor-inducible gene 6 protein; Possibly involved in cell-cell and cell-matrix interactions during inflammation and tumorigenesis.

CFP
Properdin; A positive regulator of the alternate pathway of complement. It binds to and stabilizes the C3-and C5-convertase enzyme complexes.

CHI3L1
Chitinase-3-like protein 1; Carbohydrate-binding lectin with a preference for chitin. Has no chitinase activity. May play a role in tissue remodeling and in the capacity of cells to respond to and cope with changes in their environment.

DOCK2
Dedicator of cytokinesis protein 2; Involved in cytoskeletal rearrangements required for lymphocyte migration in response of chemokines.

TCN1
Transcobalamin-1; Binds vitamin B12 with femtomolar affinity and protects it from the acidic environment of the stomach; Belongs to the eukaryotic cobalamin transport proteins family.

ORM1
Alpha-1-acid glycoprotein 1; Functions as transport protein in the blood stream. Binds various ligands in the interior of its beta-barrel domain. Appears to function in modulating the activity of the immune system during the acute-phase reaction.

NCKAP1L
Nck-associated protein 1-like; Essential hematopoietic-specific regulator of the actin cytoskeleton (Probable). Controls lymphocyte development, activation, proliferation and homeostasis, erythrocyte membrane stability, as well as phagocytosis and migration by neutrophils and macrophages. Component of the WAVE2 complex which signals downstream of RAC to stimulate F-actin polymerization.

ABI2
Abl interactor 2; May act in regulation of cell growth and transformation by interacting with nonreceptor tyrosine kinases ABL1 and/or ABL2. Part of the WAVE complex that regulates lamellipodia formation. The WAVE complex regulates actin filament reorganization via its interaction with the Arp2/3 complex.

NCKAP1
Nck-associated protein 1; Part of the WAVE complex that regulates lamellipodia formation. The WAVE complex regulates actin filament reorganization via its interaction with the Arp2/3 complex. As component of the WAVE1 complex, required for BDNF-NTRK2 endocytic trafficking and signaling from early endosomes.

WASF1
Wiskott-Aldrich syndrome protein family member 1; Downstream effector molecule involved in the transmission of signals from tyrosine kinase receptors and small GTPases to the actin cytoskeleton. Promotes formation of actin filaments. Part of the WAVE complex that regulates lamellipodia formation.

BRK1
Protein BRICK1; Involved in regulation of actin and microtubule organization. Part of a WAVE complex that activates the Arp2/3 complex. As component of the WAVE1 complex, required for BDNF-NTRK2 endocytic trafficking and signaling from early endosomes.

WASF2
Wiskott-Aldrich syndrome protein family member 2; Downstream effector molecule involved in the transmission of signals from tyrosine kinase receptors and small GTPases to the actin cytoskeleton. Promotes formation of actin filaments. Part of the WAVE complex that regulates lamellipodia formation. The WAVE complex regulates actin filament reorganization via its interaction with the Arp2/3 complex; Wiskott-Aldrich Syndrome protein family.

TUBGCP5 (Tubulin Gamma Complex Associated Protein 5) Gene
Attributes, Location, Description, Function, and Associated Disorders for TUBGCP5 Gene.
The chromosomal band location for TUBGCP5 gene is 15q11.2 with the genomic location (GRCh38/hg38) at 22,983,192-23,039,673 with a size of 56,482 bases in a minus strand orientation (Genecards.org). TUBGCP5 protein has the following attributes: Size: 1024 amino acids; Molecular mass:118321 Da; Quaternary structure: Gamma-tubulin complex which is composed of gamma-tubulin, TUBGCP2, TUBGCP3, TUBGCP4, TUBGCP5 and TUBGCP6 proteins (GeneCards & UniProtKB/Swiss-Prot). This protein has 2 isoforms produced by alternative splicing; the sequence of this isoform differs from the canonical sequence as follows: 1010-1024: LESLALSLMAGMEQS → CEYIMLKYFYLCISL. TUBGCP5 gene is widely expressed with the highest levels in the heart and skeletal muscle with moderate levels in the brain (GeneCards & UniProtKB/Swiss-Prot). An important paralog of this gene is TUBGCP6. Table 12 shows biological processes, molecular functions, cellular components with KEGG and Reactome pathways related to TUBGCP5. Among its related pathways are Nanog in Mammalian ESC Pluripotency and G-Beta Gamma Signaling. Gene Ontology (GO) annotations related to this gene include microtubule binding, and components of the cytoskeleton (Table 12). Table 13 shows diseases associated with TUBGCP5 gene disturbances based on MalaCards and arranged in descending order including Prader-Willi syndrome as the top disorder followed by schizophrenia and then autism. Table 14 shows proteins and their functions that interact with TUBGCP5, as depicted in Figure 6, illustrating the interactions of the TUBGCP5 gene with other related genes.

Protein-Protein
Inter-Relationships Involving the TUBGCP5 Gene Figure 6 illustrates the protein-protein inter-relationships or interactions of the TUBGCP5 gene with other genes. Table 14 shows biological processes, molecular functions, cellular components with KEGG and Reactome pathways related to TUBGCP5. Figure 6. STRING Protein-Protein Interaction Network for TUBGCP5 gene with functional interactions with 26 nodes (Source: STRING.org). Network nodes represent proteins with splice isoforms or post-translational modifications collapsed into each node for all proteins produced by a single protein-coding gene. Edges represent protein-protein associations that are considered specific and meaningful (i.e., proteins jointly contribute to a shared function (Table 12). TUBGCP5 gene interactants and their functions are listed in Table 14. (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins).

TUBGCP5
Gamma-tubulin complex component 5; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

HAUS4
HAUS augmin-like complex subunit 4; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

HAUS3
HAUS augmin-like complex subunit 3; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

TUBGCP6
Gamma-tubulin complex component 6; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

HAUS8
HAUS augmin-like complex subunit 8; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

CEP135
Centrosomal protein of 135 kDa; Centrosomal protein involved in centriole biogenesis. Acts as a scaffolding protein during early centriole biogenesis. Required for the targeting of centriole satellite proteins to centrosomes such as of PCM1, SSX2IP and CEP290 and recruitment of WRAP73 to centrioles. Also required for centriole-centriole cohesion during interphase.

DCTN3
Dynactin subunit 3; Together with dynein may be involved in spindle assembly and cytokinesis.

HAUS2
HAUS augmin-like complex subunit 2; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

TUBGCP4
Gamma-tubulin complex component 4; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

TUBGCP3
Gamma-tubulin complex component 3; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.  Table 12). TUBGCP5 gene interactants and their functions are listed in Table 14. (See Figure 2 for legend with a description of symbols and notes on interpreting the interaction of genes and their encoded proteins).

TUBGCP5
Gamma-tubulin complex component 5; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

HAUS4
HAUS augmin-like complex subunit 4; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

HAUS3
HAUS augmin-like complex subunit 3; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

TUBGCP6
Gamma-tubulin complex component 6; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

HAUS8
HAUS augmin-like complex subunit 8; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

CEP135
Centrosomal protein of 135 kDa; Centrosomal protein involved in centriole biogenesis. Acts as a scaffolding protein during early centriole biogenesis. Required for the targeting of centriole satellite proteins to centrosomes such as of PCM1, SSX2IP and CEP290 and recruitment of WRAP73 to centrioles. Also required for centriole-centriole cohesion during interphase.

DCTN3
Dynactin subunit 3; Together with dynein may be involved in spindle assembly and cytokinesis.

HAUS2
HAUS augmin-like complex subunit 2; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

TUBGCP4
Gamma-tubulin complex component 4; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

TUBGCP3
Gamma-tubulin complex component 3; Gamma-tubulin complex is necessary for microtubule nucleation at the centrosome.

CEP76
Centrosomal protein of 76 kDa; Centrosomal protein involved in regulation of centriole duplication. Required to limit centriole duplication to once per cell cycle by preventing centriole reduplication; Table 14. Cont.

CLASP1
CLIP-associating protein 1; Microtubule plus-end tracking protein that promotes the stabilization of dynamic microtubules. Involved in the nucleation of noncentrosomal microtubules originating from the trans-Golgi network (TGN). Required for the polarization of the cytoplasmic microtubule arrays in migrating cells towards the leading edge of the cell.

YWHAE
14-3-3 protein epsilon; Adapter protein implicated in the regulation of a large spectrum of both general and specialized signaling pathways.

CEP72
Centrosomal protein of 72 kDa; Involved in the recruitment of key centrosomal proteins to the centrosome. Provides centrosomal microtubule-nucleation activity on the gamma-tubulin ring complexes (gamma-TuRCs) and has critical roles in forming a focused bipolar spindle, which is needed for proper tension generation between sister chromatids. Involved in centriole duplication.

CEP70
Centrosomal protein of 70 kDa; Plays a role in the organization of both preexisting and nascent microtubules in interphase cells. During mitosis, required for the organization and orientation of the mitotic spindle.

PRKAR2B
cAMP-dependent protein kinase type II-beta regulatory subunit; Regulatory subunit of the cAMP-dependent protein kinases involved in cAMP signaling in cells.

PLK4
Serine/threonine-protein kinase PLK4; Serine/threonine-protein kinase that plays a central role in centriole duplication. Able to trigger procentriole formation on the surface of the parental centriole cylinder.

NINL
Ninein-like protein; Involved in the microtubule organization in interphase cells. Overexpression induces the fragmentation of the Golgi and causes lysosomes to disperse toward the cell periphery; it also interferes with mitotic spindle assembly.

CEP164
Centrosomal protein of 164 kDa; Plays a role in microtubule organization and/or maintenance for the formation of primary cilia (PC), a microtubule-based structure that protrudes from the surface of epithelial cells. Plays a critical role in G2/M checkpoint and nuclear divisions. A key player in the DNA damage-activated ATR/ATM signaling cascade.

HAUS1
HAUS augmin-like complex subunit 1; Contributes to mitotic spindle assembly, maintenance of centrosome integrity and completion of cytokinesis.

PLK1
Serine/threonine-protein kinase PLK1; Serine/threonine-protein kinase that performs several important functions throughout M phase of the cell cycle, including the regulation of centrosome maturation and spindle assembly, the removal of cohesins from chromosome arms, the inactivation of anaphase-promoting complex/cyclosome (APC/C) inhibitors, and the regulation of mitotic exit and cytokinesis. PRKACA cAMP-dependent protein kinase catalytic subunit alpha; Phosphorylates a large number of substrates in the cytoplasm and the nucleus. Regulates the abundance of compartmentalized pools of its regulatory subunits through phosphorylation of PJA2 which binds and ubiquitinates these subunits, leading to their subsequent proteolysis. Required for glucose-mediated adipogenic differentiation increase and osteogenic differentiation inhibition from osteoblasts.

MZT1
Mitotic-spindle organizing protein 1; Required for gamma-tubulin complex recruitment to the centrosome.

Discussion
Features seen in patients with the 15q11.2 BP1-BP2 deletion (Burnside-Butler) syndrome can be quite variable. Different clinical phenotypes are seen in children with this disorder and can depend on the source of the parental deletion [36]. Recent evidence of partial expression or bias of these four genes depending on the parent of origin leads to possible clinical differences in the offspring. For example, when the deletion is paternal, there is a greater risk of having congenital heart defects [36]. If the mother transmits the deletion, then there is a greater risk of intellectual disability and autism in their affected children [36]. A recent review further showed a phenotype with global and regional measures of surface area and cortical thickness in the brain with subcortical volumes and cognition differences compared to controls. This investigation consisted of 203 individuals compared with 4500 controls without the deletion collected from years 2015-2019 and the average age was 56 years [5]. This international study analyzed the largest CNV and neuroimaging results to date in those with the 15q11.2 BP1-BP2 deletions and non-deletion controls. They reported reduced brain surface area, a thicker cortex, and a smaller nucleus accumbens in the 15q11.2 BP1-BP2 deletion subjects. The significant difference in cortical thickness was more evident in the frontal, cingulate and parietal lobes. Furthermore, cognitive ability was lower for those with the deletion compared with non-deletion individuals suggesting involvement in neural plasticity and development leading to functional brain differences identified clinically [5].
As noted in Table 15, all four genes are protein-coding genes and neither were maternally nor paternally imprinted. Their encoded proteins interact with each other in crucial biological processes and molecular pathways (Figure 2, Tables 1 and 2). The predicted functional interactions encompassing all four genes included 11 nodes with 34 edges that pertain to biological processes and molecular functions utilizing searchable genomic databases involving the four genes in this 15q11.2 BP1-BP2 microdeletion syndrome (see Table 1). Among them, in addition to magnesium ion transport molecular function, significant other biological processes, such as regulation of axonogenesis and axon extension, regulation and production of bone morphogenetic protein (BMP), regulation of cellular growth and development were observed. These biological processes are relevant to variable clinical phenotypes, such as, autism, seizures, and schizophrenia which are also seen in this microdeletion syndrome. These clinical phenotypes often affect neurological development and are occasionally accompanied by mild dysmorphic features [1,[9][10][11][12][13].
Although NIPA1 and NIPA2 share related KEGG and Reactome pathways (Tables 3 and 6), their cardinal disease associations are different (Table 15): Spastic Paraplegia 6, Autosomal Dominant and Spastic Paraplegia 6 for NIPA1 and Angelman Syndrome and Prader-Willi Syndrome for NIPA2. However, as shown in Table 16, they both are also associated with PWS and AS. Unlike NIPA1 and NIPA2, CYFIP1, and TUBGCP5 genes do not share any KEGG or Reactome pathways (Tables 9 and 12). As shown in Table 15, cardinal diseases that are associated with the CYFIP1 gene are Fragile X Syndrome and Autism. The cardinal disease that is associated with the TUBGCP5 gene is Prader-Willi syndrome. However, as depicted in Table 16, they both are additionally associated with PWS, while the CYFIP1 gene alone is additionally associated with AS. Most notably, all four syntenic genes in this region are associated with Autism Spectrum Disorder; Schizophrenia; Epilepsy and Down Syndrome (Table 16). Except for TUBGC5, all three genes are associated with the 15q11.2 BP1-BP2 Deletion Syndrome with Attention Deficit Hyperactive Disorder & Learning Disability (Table 16). It should be noted that interacting proteins for CYFIP1 and TUBGCP5 proteins are quite different ( Figure 5, Table 11; Figure 6, Table 14, respectively), which is unlike the case of NIPA1 and NIPA2 (Figure 3, Table 5; Figure 4, Table 8, respectively).
It is interesting to note that none of the protein-protein interacting genes for the four genes, which total 70, are the same (Table 5, Table 8, Table 11, and Table 14). However, NIPA1 and NIPA2 proteins, given their involvement with common biological processes (GO) and molecular functions (GO) along with BMP6, BMP7, BMPR1A, BMPR1B, and BMPR2-proteins interact with both. In contrast, none of the protein-protein interactants are similar between CYFIP1 and TUBGCP5 proteins (Tables 11 and 14), thus affirming their functional dissimilarities.
It should be noted, as illustrated in STRING Figure 2 and Table 1, NIPA1, NIPA2, CYFIP1, and TUBGCP5 proteins show protein-protein interactions among themselves with the 11 nodes and 34 edges addressed earlier. As shown in Table 1, NIPA1 protein interacts with TUBGCP5, CYFIP1, and NIPA2 proteins (scores: 0.995, 0.967, and 0.941, respectively) along with seven other protein interactants, and thus predicted to cause significant biological processes (Table 2), and a significant molecular function is magnesium ion transmembrane transporter activity with a 0.0042 false discovery rate ( Table 2).
The varied biological processes and predicted functions of these four genes, as noted above, through their protein-protein interactions with seven other proteins (see Figure 2) collectively could play a role in the neurodevelopmental disorders. They complement regulation associated with lipoproteins and lipid metabolism and adipogenesis, encompassing CFHR1, CFHR3, and BMPR2 protein-protein interactions, given marked obesity is a common finding in PWS.
Examining the interaction of NIPA1 protein with its 11 other interacting proteins indicates that three-fourths of the interactions are important for developmental bone morphogenesis or multifunctional proteins that control proliferation, differentiation, and other functions in many cell types (see Figure 3; Table 5). The interactions of NIPA2 protein with 19 other proteins indicate their collective role as intracellular signal transducers and transcriptional modulators that are activated by TGFbeta, and thereby impacting bone morphogenesis (see Figure 4; Table 8). CYFIP1 protein interacts with 25 other proteins. These proteins have a wide range of activity with functions including cytoskeleton organization and actin filament binding with cell-matrix adhesion, MAP kinase signal transduction impacts cell growth, survival and differentiation, stimulation of glucose uptake in cells, intracellular protein breakdown and tissue remodeling (see Figure 5; Table 11). TUBGCP5 protein interacts with 25 other proteins. These proteins are focused on mitotic spindle formation and assembly, microtubule organization, and production of centrosomal proteins involved in the regulation of centriole duplication during cell division (see Figure 6; Table 14).

Overlapping Association of PWS and AS Specific Genes with Both Syndromes at Varying Degrees
Although UBE3A is a classical AS gene (9.31 score), it is also associated with PWS to a lesser degree (4.106 score) (https://www.malacards.org/search/results/UBE3A). Further, it should be noted that the 15q11-q13 region includes the imprinted gene cluster with two maternally expressed protein-coding genes UBE3A and ATP10A . Therefore, it may not be appropriate to assume that PWS/AS are exclusively caused only by the designated PWS and AS regional genes per se, as shown in Figure 1, since most of these genes are also overlappingly associated with both PWS and AS, although to varying degrees. Moreover, PWS/AS phenotype is also overlappingly associated with three proximally located contiguous genes within the 15q11.2 BP1-BP2 region: NIPA1, NIPA2, and CYFIP1, as depicted in Table 16. PWS is caused by loss of expression of the paternally inherited genes on chromosome 15q11.2-q13, and the cardinal features of PWS are attributable to the critical interval within the 15q11.2-q13 imprinted gene cluster, containing the small nucleolar RNA (snoRNA) SNORD116 and non-coding RNA IPW (Imprinted in Prader-Willi) exons [39][40][41]. Similarly, the cardinal features of AS are attributable to the loss of expression of maternally inherited genes, and UBE3A in particular, which causes a distinct AS neurodevelopmental disorder [42].
Although the PWS core features are attributable to the critical interval within the 15q11.2-q13 imprinted gene cluster and the AS core features are attributable to the loss of expression of maternally inherited gene UBE3A, these core features might be influenced by several other genes that lie within the broader BP1-BP3 region. The functions of all of these genes are in-turn likely modulated by hundreds of their interacting genes that lie throughout the genome, as shown in the STRING protein interactions for NAPI1, NAPI2, CYFIP1, and TUBGCP5: Figure 2, Table 1; Figure 3, Table 5; Figure 4, Table 8; Figure 5, Table 11; Figure 6, Table 14; Tables 15 and 16. One could speculate that PWS/AS, similar to that of ASD, depends on either the involvement of the PWS critical region in the imprinted segment within the 15q11.2-q13 BP2-BP3 region or the involvement of a critical gene, namely, UBE3A in the maternally imprinted segment [43] and in concert with other genes that lie within chromosome 15, not excluding the 15q11.2 BP1-BP2 region genes.
Every one of these genes along with their hundreds of interacting genes within the genome as a whole, as depicted and tabulated in this study pertaining to NIPA1, NIPA2, CYFIP1, and TUBGCP5 genes all in concert could impact on the PWS/AS signature phenotype in any given proband. Thus, each PWS/AS proband's clinical presentation is likely to be a unique genomic signature, within the broader spectrum of PWS/AS clinical presentations, which is akin to Autism Spectrum Disorder.  Table 15). Based on reports of their mutations or disturbed expression patterns, as reported in authentic web resources, such as, GeneCards.org, SAFARI.org, Gene Ontology (GO), OMIM.org, Entrez and UniProtKB/Swiss-Prot, their "collective loss or dosage duplication" is due to the BP1-BP2 microdeletion or microduplication that appears to impact only behavioral and neurological functions. These disorders include speech and motor delays, behavioral problems, seizures, and autism in affected individuals [44]. A recent review of the 15q11.2 BP1-BP2 microdeletion syndrome [1] found common phenotypic features, which included autism, developmental delay, motor and language delays, and behavioral problems.
Parental studies among these subjects demonstrated phenotypically normal carriers in several instances, and mildly affected carriers in others, complicating phenotypic association and/or causality. This could be due to either reduced penetrance or altered gene dosage on a particular genetic background [1,4]. The four non-imprinted and biallelically expressed genes, NIPA1, NIPA2, CFYIP1 and TUBGCP5, in this microdeletion were initially noted to impact the severity of clinical presentation and neurological impairment in these two classical genomic imprinting disorders, Prader-Willi and Angelman syndromes with typical 15q11-q13 deletions depending on the absence or presence of the genomic area between breakpoints BP1 and BP2 containing these four genes, which led to the recognition of this microdeletion (Burnside-Butler) syndrome [4].
Therefore, it is intriguing to note that the review [4] of over 200 15q11.2 BP1-BP2 microdeletion syndrome individuals, as meta-analysis, did not find either the NIPA1 associated cardinal disease: Spastic Paraplegia, the NIPA2 associated cardinal diseases: Angelman Syndrome and Prader-Willi Syndrome, the CYFIP1 associated cardinal disease: Fragile X Syndrome, or the TUBGCP5 associated cardinal disease: Prader-Willi Syndrome, except for the CYFIP1 associated cardinal disease: Autism, as presented in Table 15.
The review [4] also found that not all individuals with this deletion were clinically affected, but neuropsychiatric and behavior disturbances and mild dysmorphic features were associated with genomic imbalances of the 15q11.2 BP1-BP2 region, including microdeletions, but with an apparent incomplete penetrance and variable expressivity.

All Four BP1-BP2 Region Genes Are Significantly Associated with Autism Spectrum Disorder
From reported patient cohorts presenting for genetic services and microarray analysis, this microdeletion syndrome can now be recognized as the most common cytogenetic abnormality found in ASD instead [6], which is not only a cardinal CYFIP1 associated disease (see Table 15), but ASD is also consistently associated with TUBGCP5, NIPA1, and NIPA2 genes (see Table 4, Table 7, Table 10,  Table 13, and Table 16).
Therefore, it seems fair to conclude that except for PWS specific genes or transcripts, such as, snoRNAs, most other genes that flank on either side in broadly designated PWS and AS regions are all co-associated with autism. Most of the genes within the broader BP1-BP3 region are also co-associated with not only PWS and AS, but also with autism (Malacards.org; SFARI.org) [45]. Even the GABRB3 gene that is distal to the ATP10A gene is a SFARI (SFARI.org) recognized autism gene (Malacards.org; SFARI.org) [45].

15q11.2 B-P1-BP2 Microdeletion (Burnside-Butler) Syndrome: Frequency
The 15q11.2 BP1-BP2 microdeletion syndrome has a reported de novo frequency between 5%-22%, with 51% having inherited the microdeletion from an apparently unaffected parent and 35% having inherited the microdeletion from an affected parent [4,6]. These low penetrance estimates may relate to subclinical manifestations of neuropsychiatric/behavioral problems, incomplete information, or lack of detailed clinical or psychiatric studies in the parents of individuals with 15q11.2 BP1-BP2 microdeletion or members of control cohorts.
In addition to the above consistent correlation of each of the four genes that lie within the narrow proximal BP1-BP2 region with autism, in a more recent study [36] parent-of-origin effects (POE) of the 15q11.2 BP1-BP2 microdeletion were found to be associated with differences in clinical features in individuals inheriting the deletion. Among all probands studied, maternal deletions were found to be associated with epilepsy, autism spectrum disorder (p = 0.02) and macrocephaly (p = 0.016), while paternal deletions were associated with congenital heart disease (CHD) (p = 0.004) and abnormal muscular phenotypes (p < 0.05), while CHD and abnormal muscular phenotypes were seen in paternal deletions. This study not only supported POEs of this deletion, but notably included ASD, macrocephaly, and epilepsy.
The above key processes and summary of the four genes within the narrower proximal 15q11.2 BP1-BP2 region are collectively critical for normal neuronal development, plasticity, and function [1,5]. Their loss or gain (CNV) are associated with neurodevelopmental disorders, spastic paraplegia, seizures, learning, and gait disturbances with motor delay and neuro-behavioral psychiatric problems including autism, dyslexia, and schizophrenia/paranoid psychosis [1].
The most common maladies found to be individually associated with each of the four genes are: spastic paraplegia found in the top 12 of 15 disorders associated with NIPA1; PWS/AS, epilepsy and psychiatric/behavioral problems (autism or schizophrenia) in the top 10 disorders for NIPA2; fragile X syndrome, autism, schizophrenia, PWS and pervasive developmental disorder in the top 10 disorders for CYFIP1; and PWS, schizophrenia, autism, microcephaly, essential hypertension, body mass index (quantitative trait), epilepsy, and Down syndrome in the top 10 disorders associated with TUBGCP5 (see Table 16). As noted in Table 15, all four genes are neither maternally nor paternally imprinted, but their proteins interact with each other as described and illustrated (see Figure 2, Tables 1 and 2) in crucial biological processes (GO) and molecular pathways (GO). Collectively, either due to the microdeletion or microduplication, these genes predominantly affect brain morphology [1,4].
Our detailed exploration of the various aspects of each of these four genes, in solo and in concert is meant to enable an in-depth understanding of their individual and collective contribution in the causation of different neurodevelopmental phenotypes that have been variably reported in the literature. All four genes are individually associated with Prader-Willi syndrome, autism spectrum disorder, schizophrenia, epilepsy, and Down syndrome (Maladies.org; GeneCards.org). Except for the TUBGCP5 gene, all three remaining genes are associated with Angelman syndrome, although the TUBGCP5 gene is located most proximal to the Angelman syndrome gene, UBE3A. TUBGCP5 gene is also not associated with attention deficit hyperactivity disorder and learning disability, developmental disorder, and peripheral nervous system disease. CYFIP1 is the only gene that is not associated with microcephaly, but it is the only gene that is associated with a developmental disorder (see Table 16).
Thus, collectively, all four genes are associated with ten overlapping neurodevelopmental maladies up to 77.5% of the time.

Materials and Methods
We characterized the four genes (NIPA1, NIPA2, CYFIP1, and TUBGCP5) found within the 15q11.2 BP1-BP2 chromosome region of about 500 bp in size. When this region is deleted, an emerging disorder (i.e., Burnside-Butler syndrome) is being characterized leading to the importance of our description of these four genes. We undertook a literature review of hundreds of patients reported and summarized clinically with this microdeletion reported as the most common genetic defect found in consecutive patients presenting for genetic testing with high-resolution microarray analysis and autism spectrum disorder (e.g., [6]). Genome databases, search platforms for gene interaction, gene ontology (GO) biological process, molecular function, cellular components, along with KEGG and Reactome pathways using STRING network data along with sources, GeneCards.org, SAFARI.org, KEGG, PubMed, Gene Ontology (GO), OMIM.org, Entrez and UniProtKB/Swiss-Prot were summarized in tabular form and STRING Network results illustrated in figures and discussed in the text involving all four genes. An overview of each gene was also described throughout the summary process utilizing an approach reported previously [46].
The Gene Ontology Resource (GO; http://geneontology.org) provides structured, computable knowledge regarding the functions of genes and gene products. Founded in 1998, GO has become widely adopted in the life sciences, and its contents are under continual improvement, both in quantity and in quality. The Gene Ontology resource (GO; http://geneontology.org) is the most comprehensive and widely used knowledge base concerning the functions of genes. In GO, all functional knowledge is structured and represented in a form amenable to computational analysis, which is essential to support modern biological research. The GO knowledge base is structured using a formal ontology, by defining classes of gene functions (GO terms) developed over the past 20 years that have specified relations to each other [47]. GO terms include definitions, or equivalence axioms, that define the term relative to other terms in the GO or other ontologies. Their relationships can be computationally inferred using logical reasoning. The GO structure is constantly evolving in response to new scientific discoveries and are continuously refined using current biological information. The GO knowledgebase consists of ontology and annotations. As of the 5 September 2018 release (doi:10.5281/zenodo.1410625), there were ∼45000 terms in GO: 29698 biological processes, 11147 molecular functions, and 4201 cellular components, linked by almost 134000 relationships [47].
MalaCards Relevance Scores provide information used for the construction of our tables of Putative Associated Diseases (see Tables 4, 7, 10 and 13) for the readership. MalaCards is an integrated database of human maladies and their annotations modeled from the architecture and information from the popular GeneCards database of human genes. The search platform used by MalaCards to obtain these data [48] is SOLR, an open-source full-text search platform widely used for analytics based on Apache's Lucene text search API. When a term is searched by Lucene then it returns a set of scored hits. A "hit" represents a document (a MalaCard), whose fields (actual annotations) were previously indexed by Lucene.
The scoring is calculated by a Lucene defined algorithm: score(q,d) = coord(q,d) · queryNorm(q) · tf(t in d) · idf(t) · t.getBoost() · norm(t,d) t in q The factors in this formula are: • tf stands for term frequency-the more times a search term appears in a document, the higher the score • idf stands for inverse document frequency-matches on rarer terms count more than matches on common terms • coord is the coordination factor-if there are multiple terms in a query, the more terms that match, the higher the score • lengthNorm-matches on a smaller field score higher than matches on a larger field • index-time boost-if a boost was specified for a document at index time, scores for searches that match that document will be boosted. • query clause boost-a user may explicitly boost the contribution of one part of a query over another.
To ensure showing the best precision, MalaCards displays the score as (base 2 log of the score) + 10 For more about Lucene's scoring mechanism see Apache Lucene-Scoring [48]. GeneCards.org-the human gene database [49] is from the GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analysis.
STRING protein network was utilized to generate Figures 2-6 and their accompanying data on the protein-protein interactions as well as their Biological Processes (GO), Molecular Functions (GO), Cellular Component (GO) and KEGG and Reactome Pathways. The data are presented in Tables 1-3,  Table 5, Table 6, Table 8, Table 9, Table 11, Table 12, and Table 14 and are adopted from STRING CONSORTIUM 2019. Data in Table 15 were derived from the GeneCards website (GeneCards.org) for respective genes, while the data for Table 16 were derived from MalaCards diseases (Malacards.org) as seen in Table 4, Table 7, Table 10, and Table 13 in this study.