Next Article in Journal
CRISPR-Cas9 Mediated TSPO Gene Knockout alters Respiration and Cellular Metabolism in Human Primary Microglia Cells
Next Article in Special Issue
Network-Based Integrative Analysis of Genomics, Epigenomics and Transcriptomics in Autism Spectrum Disorders
Previous Article in Journal
Determination of the Role and Active Sites of PKC-Delta-Like from Lamprey in Innate Immunity
Previous Article in Special Issue
Risk Factors for Unhealthy Weight Gain and Obesity among Children with Autism Spectrum Disorder
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Functioning Autism with Missense Mutations in Synaptotagmin-Like Protein 4 (SYTL4) and Transmembrane Protein 187 (TMEM187) Genes: SYTL4- Protein Modeling, Protein-Protein Interaction, Expression Profiling and MicroRNA Studies

1
Departments of Psychiatry & Behavioral Sciences and Pediatrics, University of Kansas Medical Center, Kansas City, KS 66160, USA
2
Department of Pediatric Neurology, Hospital Universitario Quirón, 28223 Madrid, Spain
3
Genomics and Medicine, NIM Genetics, 28108 Madrid, Spain
4
Department of Biochemistry and Molecular Biology, University of Kansas Medical Center, Kansas City, KS 66160, USA
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(13), 3358; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133358
Submission received: 25 March 2019 / Revised: 10 June 2019 / Accepted: 17 June 2019 / Published: 9 July 2019

Abstract

:
We describe a 7-year-old male with high functioning autism spectrum disorder (ASD) and maternally-inherited rare missense variant of Synaptotagmin-like protein 4 (SYTL4) gene (Xq22.1; c.835C>T; p.Arg279Cys) and an unknown missense variant of Transmembrane protein 187 (TMEM187) gene (Xq28; c.708G>T; p. Gln236His). Multiple in-silico predictions described in our study indicate a potentially damaging status for both X-linked genes. Analysis of predicted atomic threading models of the mutant and the native SYTL4 proteins suggest a potential structural change induced by the R279C variant which eliminates the stabilizing Arg279-Asp60 salt bridge in the N-terminal half of the SYTL4, affecting the functionality of the protein’s critical RAB-Binding Domain. In the European (Non-Finnish) population, the allele frequency for this variant is 0.00042. The SYTL4 gene is known to directly interact with several members of the RAB family of genes, such as, RAB27A, RAB27B, RAB8A, and RAB3A which are known autism spectrum disorder genes. The SYTL4 gene also directly interacts with three known autism genes: STX1A, SNAP25 and STXBP1. Through a literature-based analytical approach, we identified three of five (60%) autism-associated serum microRNAs (miRs) with high predictive power among the total of 298 mouse Sytl4 associated/predicted microRNA interactions. Five of 13 (38%) miRs were differentially expressed in serum from ASD individuals which were predicted to interact with the mouse equivalent Sytl4 gene. TMEM187 gene, like SYTL4, is a protein-coding gene that belongs to a group of genes which host microRNA genes in their introns or exons. The novel Q236H amino acid variant in the TMEM187 in our patient is near the terminal end region of the protein which is represented by multiple sequence alignments and hidden Markov models, preventing comparative structural analysis of the variant harboring region. Like SYTL4, the TMEM187 gene is expressed in the brain and interacts with four known ASD genes, namely, HCFC1; TMLHE; MECP2; and GPHN. TMM187 is in linkage with MECP2, which is a well-known determinant of brain structure and size and is a well-known autism gene. Other members of the TMEM gene family, TMEM132E and TMEM132D genes are associated with bipolar and panic disorders, respectively, while TMEM231 is a known syndromic autism gene. Together, TMEM187 and SYTL4 genes directly interact with recognized important ASD genes, and their mRNAs are found in extracellular vesicles in the nervous system and stimulate target cells to translate into active protein. Our evidence shows that both these genes should be considered as candidate genes for autism. Additional biological testing is warranted to further determine the pathogenicity of these gene variants in the causation of autism.

1. Introduction

Whole exome sequencing (WES) and occasionally whole-genome sequencing (WGS) are increasingly used in clinical practice for diagnosis, medical intervention and prognosis [1,2,3]. High heritability estimates and family studies have supported a definite role of genetics in autism spectrum disorder (ASD). In neurodevelopmental disorders, particularly ASD and intellectual disability have diagnostic rates using WES which fluctuates in part from differences in clinical features from one study to another with rates up to 50% [4,5,6,7]. We describe a 7-year-old male with high-functioning autism spectrum disorder presenting for genetic services and whole exome sequencing following a normal microarray analysis, and variants were found in two X-linked genes: SYTL4 and TMEM187.

1.1. Synaptotagmin-Like Protein 4 (SYTL4) Gene

SYTL4 gene, also known as Granuphilin/SLP4, encodes a member of the synaptotagmin-like protein family. Members of this family are characterized by an N-terminal RAB27-binding domain and C-terminal tandem C2 domains (Figure 1 and Figure 2). The first C2 domain binds phospholipids in a calcium-independent manner, whereas the second C2 domain does not (http://omim.org/entry/300723). The encoded protein binds to specific small RAB-GTPases involved in intracellular membrane trafficking. This protein binds to RAB27 and may be involved in inhibiting dense core vesicle exocytosis. Alternate splicing results in multiple transcript variants that encode the same protein.
SYTL4 gene is relevant in neuronal system development and implicated in neurological and psychological diseases [8] (Entrez Gene ID # 94121 (Human); ID # 27359 (Mouse). This gene is highly expressed in the bed nucleus of the stria terminalis (BNTS) [9], a brain region that regulates mood, motivation for social behavior and social attachment [10]. SYTL4 expression is down regulated in the dorsal raphe nucleus from patients with major depressive disorder [11]. In a mouse model of anxiety, significant change in Sytl4 was observed among the altered protein networks in the brain proteome [12].
Targeted Null/Knockout mutant Sytl4 mammalian phenotypes for Sytl4 showed abnormal behavior and neurological phenotype [13,14,15]. This gene interacts directly with other known autism genes. Confirmatively, gene expression knock-out Sytl4 mutant mouse includes abnormal behavior and neurological problems [13,14]. Diseases already associated with SYTL4 include Epileptic Encephalopathy (early Infantile) and Branchiootic syndrome (www.genecards.org/cgi-Fbin/carddisp.pl?gene= SYTL4&keywords=sytl4). A missense mutation in SYTL4 was also recently reported in a female with autism and non-skewed X-chromosome inactivation [16].

1.2. Transmembrane Protein 187 (TEM187) Gene

The second gene in our study, TMEM187, is additionally known as ITBA1/CXORF12/DXS9878E- gene and consists of two exons which encode a multi-pass membrane protein expressed in all regions of the brain (www.genecards.org/cgi-bin/carddisp.pl?gene=TMEM187&keywords=TMEM187). This is a conserved gene that codes a 261-amino-acid protein with six transmembrane helical domains. Unlike SYTL4, not much is known about this gene, but it belongs to a group of genes which host microRNA genes in their introns or exons [17]. One of this gene’s related phenotypes is schizophrenia (GWAS catalog for TMEM187 gene: Gene relation via enhancers containing phenotype SNP: Enhancer ID: GH0XJ153980). The latest STRING- network of protein interactions for this gene reveals that it directly interacts with four known ASD genes, namely, HCFC1; TMLHE; MECP2 and GPHN. (https://gene.sfari.org/database/human-gene/).

2. Results

2.1. Genomic Study

Whole exome sequencing performed in trios revealed maternally inherited X-linked missense variants: c.835C>T, p. Arg279Cys in the SYTL4 gene and c.708G>T; and p. Gln236His in the TMEM187 gene. Research findings will be discussed that relate specifically to SYTL4 and TMEM187 gene variants seen in our patient and as emerging candidate genes for autism.

2.2. Synaptotagmin-Like 4 (SYTL4) Gene

2.2.1. Deleterious and Damaging Nature of the SYTL4- Variant

The maternally inherited variant in Synaptotagmin-like protein 4 gene (SYTL4) at chromosome-X position g.99944930G>A, c.835C>T, p. Arg279Cys [NM_001174068.1] sits within exon 9 of the SYTL4 gene located at the Xq22.1 cytoband (Figure 1). This particular SYTL4 nucleotide variant CGC⇒TGC at position 835 (exon 9; alternating) has been predicted by large-scale genomic sequencing studies, such as the 1000 Genome project, to be deleterious (SIFT score 0.01) and possibly damaging (Polyphen score 0.79), which is further validated with a MAF score of 0.0005/2. The Exome Aggregation Consortium (ExAC) database including data from 1000 genomes further indicated that the SYTL4 variant is also damaging.
Our patient is from Spain and according to Exome Aggregation Consortium, the allele frequency among European (Non-Finnish) population for this variant is 0.0004222. No homozygotes were found, but eight hemizygotes (male) were observed in 47372 alleles which presumably included both isoforms. There is no reported information from these individuals regarding findings related to ASD when examining clinical information-sharing resources such as Phenome Central, Gene Matcher or ClinVar.

2.2.2. Modeling of Native and R279C Mutant for SYTL4 Gene

There is limited high-resolution structural information about SYTL4. The crystal [18] and NMR [19] structure of C2(1) (residues 354-483) and ring domain (residues 43-105) cover only 29% of the protein primary structure. A structural comparison of SYTL4 and the R279C mutant was carried out using the I-TASSER multiple threading approach [20]. The protein structure closest to both forms of SYTL4 found in the Protein Data Bank (PDB) [21] is human synaptotagmin 2 (PDB #4P42) [22,23], with TM scores (0.658 and 0.661) and RMSD values (1.15 Å and 0.96 Å), respectively, for the native and R279C form (Figure 2).

2.2.3. Hierarchical Protein Structural Modeling Study of Both Native and R (279) C SYTL4

Both atomic models preserved the canonical C2 structures of synaptotagmins; [21,22,23], however, only the native SYTL4 preserved the extended structure of the ring domain (Figure 3) [18]. An overlay (Figure 3C) of two atomic models demonstrated good alignment of the C-terminal C2 domains, whereas the structures of the region comprising both the ring domain and the 279 amino acid change sites for the two SYTL4 forms differed considerably. In the native SYTL4 model, Arg 279 (side chain in red) is part of a large extended loop conformation that appears to be stabilized by an apparent salt bridge formed between it and Asp60 (side chain in Salmon) in the ring domain. The distance (3.116 Å) calculated between the Arginine guanidinium nitrogen and Asp carboxyl oxygen (magnified in Figure 3B) is well within the threshold distance observed for salt bridges in a comprehensive survey of crystal structures [24]. In the R279C model, Cys279 is in a beta sheet, with no apparent hydrogen bonding contacts observed within distance constraints for such interactions [25] (Figure 3A,B). Our modeling suggests a potential structural change induced by the R279C variant, which eliminates the stabilizing Arg279-Asp60 salt bridge and, therefore, leads to significant structural changes in the N-terminal half of the SYTL4.

2.2.4. STRING- Protein–Protein Interaction Network Study Reveals Direct Interaction of SYTL4 with Other Known Autism Genes

SYTL4 directly interacts with three known ASD genes, STX1A, SNAP25 and STXBP1 (https://gene.sfari.org/database/human-gene/). SYTL4 also interacts with 14 other genes, namely, FGF4, STX1B, SNAP29, RAB3A, SNAP23, RAB27B, RAB8A, SNAP47, STX4, STX19, STX3, STX11, STX2, and RAB27A, whose transcriptome-wide isoform-level multiple related gene family members are known ASD genes (https://gene.sfari.org/database/human-gene/). SYTL4 also directly interacts with several members of the RAB-family of genes, such as RAB27A, RAB27B, RAB8A, and RAB3A [26] (Figure 4).

2.2.5. SYTL4- Molecular Pathways and Associated Diseases

The molecular pathways of the SYTL4 gene includes a synaptic vesicle cycle, insulin secretion, AMPK- signaling, and SNARE interactions in vesicular transport (Table 1). Diseases associated with synaptic vesicle cycle defects include early infantile epileptic encephalopathy and defects in insulin secretion including defects in degradation of gangliosides which are abundantly expressed in the nervous system. The deficiency of certain gangliosides will affect the regenerative ability of injured hypoglossal nerves [27]. Lastly, defective SNARE interactions in vesicular transport have been implicated in CEDNIK-syndrome, which includes cerebral dysgenesis and neuropathy. All the noted pathways and their associated diseases broadly impact nervous system function which are relevant to autism (Table 1) [26].

2.2.6. SYTL4- Networks of Biological Processes

As presented in Table 2 [26], extensive biological processes of the SYTL4 gene pertain to synaptic vesicle functions, including neurotransmitter secretion and regulation of signaling, RAB protein signal transduction, glutamate secretion, neuro-muscular synaptic transmission, and axonogenesis, all of which are relevant for proper neuronal function, and thus important in autism.

2.2.7. SYTL4- Molecular Functions

The molecular function of the SYTL4 gene includes syntaxin binding, which is essential for neurotransmission (Table 3) [26] and directly interacts with syntaxin-binding protein 1 (STXBP1) (Figure 4). Mutations in STXBP1 are associated with infantile-epileptic encephalopathy-4 [28]. SYTL4 molecular gene function also relates to synaptosomal-associated protein receptor (SNAPRe) activity which regulates neurotransmitter release to ensure vesicle-to-target specificity [29]. In addition, the molecular function of the SYTL4 gene pertains to trimeric-G-protein (GDP binding protein) which plays a pivotal role in signal transduction pathways for numerous hormones and neurotransmitters [30].

2.2.8. Missense Mutation Causing R (279) C Amino Acid Change Affects the Structures of Canonical SYTL4 Gene as Well as Its Shorter Isoform

The amino acid change found in our study, R(arg)279C(cys) in exon 9 (Figure 1), is integral for the full length ‘canonical’ sequence of the SYTL4 protein with 671 amino acids. It is also part of the shorter isoform containing 349 amino acids shown in Figure 5. The R>C amino acid variation at residue 279 has been classified by the Human Genome Variation Society (HGVS) as a missense variant (variation ID # rs141441277) and determined to be least common with a Minor Allele Frequency (MAF) of 0.000529801. It is further predicted as “possible damaging” with a PolyPhen score of 0.79, while SIFT predicts the effect of this amino acid change on protein function as “deleterious” (https://gnomad.broadinstitute.org/variant/X-99944930-G-A; Variant Effect Predictor (VEP) program: https://useast.ensembl.org/info/docs/tools/vep/index.html.

2.2.9. Autism Predictive Human Serum MicroRNAs with Predicted Interaction with Mouse Sytl4 Gene

Among the total of 298 mouse Sytl4 associated/predicted miRs, three of the five miRs (60%), namely, miR181b-5p, miR320a, and miR130a-3p- are with good ASD predictive power in serum, as identified in Table 4. Apart from the five-serum miRNAs with good ASD predictive power, miR106b-5p and miR328- were up and down regulated (respectively) in the serum ASD study. Both miRs have been reported as showing altered expression among schizophrenics [31,32,33,34]. Interestingly, these two miRs have been predicted to interact with mouse Sytl4 gene (Table 4). In addition, four other miRs (miR98, miR103, miR132, and miR320) have been found to be dysregulated in the superior temporal gyrus of ASD [35] and predicted to interact with mouse Sytl4 gene. Furthermore, miR106b, miR181b-5p, miR320, and miR328 are differentially expressed in the ASD cerebellar cortex [35,36] and predicted to interact with mouse Sytl4 gene (Table 4).

2.3. Transmembrane Protein 187 (TEM187) Gene

2.3.1. Structure of TMEM187 Gene, Expression and Location of the Novel Variant

TMEM187 gene consists of two exons and encodes a multi-pass membrane protein (Figure 6).
The maternally inherited TMEM187 missense variant lies beyond the last transmembrane helix region at the near terminal end of the Pfam domain 8-245aa represented by multiple sequence alignments and hidden Markov models (https://pfam.xfam.org/family/tmem187; TMEM187 (PF15100)). The TMEM187 gene has three transcripts (splice variants) and 35 orthologs.

2.3.2. Deleterious and Damaging Nature of the Novel TMEM187 Gene Variant

Unlike the SYTL4 gene variant, the TMEM187 missense gene variant (c.708G>T; p. Gln236His) is neither found in the Exome Aggregation Consortium (ExAC) database (which includes the 1000 genome data) or in the Ensemble.org data base. The nature of this novel Glutamine(Q)236 to Histidine(H) variation in TMEM187 is considered deleterious as determined by the PROVEAN (Protein Variation Effect Analyzer v1.1) score of = −4.474 with prediction cutoff of −2.5; SIFT score as deleterious (0.05); and PolyPhen-2 score as probably damaging (0.432). (http://gnomad.broadinstitute.org/variant/X-153248221-G-T). PROVEAN prediction generated score of −4.47 is termed as deleterious (scores <−2.5 considered deleterious), while the SIFT score of 0.05 approaches 0.0 which is the most damaging score.

2.3.3. TMEM187 Gene Is Expressed in the Brain

The TMEM187 gene is ubiquitously expressed in all systems including all parts of the brain (www.uniprot.org/uniprot/Q14656; www.genecards.org/).

2.3.4. Latest STRING- Gene Interaction Network Study Reveals Direct Protein–Protein Interactions of TMEM187 with Several Other Known Autism Genes

Analysis of its STRING- network of interactions (Figure 7) reveals that it directly interacts with four known ASD genes, namely HCFC1, TMLHE, MECP2, and GPHN (https://gene.sfari.org/database/human-gene/). Moreover, three other directly interacting genes, namely UBL4A, RBM25, and AKAP4, though not known ASD genes, do show other transcriptome-wide isoform-level family member genes known as ASD genes: UBL7, RBM27, RBM8A, RBMS3, and AKAP9 (https://gene.sfari.org/database/human-gene/).

3. Discussion

Our research findings will be discussed related specifically to both the SYTL4 and TMEM187 gene variants seen in our patient and evidence for the two gene variants playing a role in the causation of high-functioning autism spectrum disorder.

3.1. Synaptotagmin-Like 4 (SYTL4) Gene

As an emerging candidate gene for autism with analyzed protein modeling, protein interactome networks, expression profiling and microRNA interactions will follow.

3.1.1. Protein Structure Altering Rare Variants Have Been Observed to Be More Frequent in Individuals with Autism

Though common variants are a large driving factors in autism spectrum disorder, the effect size of individual common variants is estimated to be small [38,39]. Therefore, the search for rare variants exhibiting a much larger individual effect is ongoing. Protein structure-altering rare variants have been observed more frequently in ASD cases, and ASD risks are increased when two rare variants may deleteriously affect both copies for an autosomal protein, or a single copy of an X-chromosomal protein among ASD males [40]. Rare hemizygous mutations on the X-chromosome are found to be more enriched in male ASD patients compared to controls. Furthermore, if rare hemizygous mutations on the X-chromosome in males alter gene expression known to be present in the brain, then the overall odds-ratio for ASD will be increased [41].

3.1.2. Deleterious and Damaging Nature of the SYTL4 Gene Variant

It is important to note that this particular SYTL4 nucleotide variant CGC⇒TGC resulting in p. Arg279Cys (Figure 1) has been predicted through large-scale genomic sequencing studies, such as the 1000 Genome project, and has been judged to be deleterious (SIFT score 0.01) and possibly damaging (Polyphen score 0.79) with a MAF score of 0.0005/2.
It should also be noted that ExAC database which includes the 1000 genomes data, also calls this SYTL4 variant damaging and deleterious based on Polyphen and SIFT, respectively. Our patient being from Spain, according to ExAC, among the European (Non-Finnish) population, the allele frequency for this variant is 0.0004222 without any homozygotes and with eight hemizygotes (male) for 47,372 alleles, which presumably includes both isoforms. There is no information available as to any of these eight hemizygotes, or any other hemizygote being reported as exhibiting ASD in any of the clinical information-sharing resources such as Phenome Central, Gene Matcher and ClinVar.

3.1.3. Randomness of X Chromosome Inactivation Could Render the Mother Asymptomatic

It is important to note that our male heterozygous carrier of the R (279) C variant is expected to be fully penetrant for the mutant allele. The proband’s mother is a heterozygous carrier of the R (279) C variant and due to randomness of X chromosome inactivation may render her as asymptomatic given that the X-linked SYTL4 gene is not over expressed among females, indicating that it does not escape inactivation [42]. The proband’s mother is expected to express equally her normal SYTL4 allele along with the mutant allele.

3.1.4. Our Modeling Results Show Large Conformational Changes Proximal to the R (279) C Amino Acid Variation

It is possible that the structure and function of both the canonical and the truncated isoform of SYTL4 protein will be affected, particularly considering their interactions at membrane surfaces. Our modeling results show large conformational changes proximal to the R (279) C amino acid variation within exon 9. The flanking regions contain two known phosphorylatable serines (YTKS (@274) VIDLR (@279) P EDVVHESGS (@289) L) as shown in Figure 5.

3.1.5. Missense Mutations Change the Size or Properties of Amino Acids Preventing the Function of Proteins

A study (reported by our co-authors: SKR, MGB) in whole exome sequencing in females with autism, a non-synonymous missense mutation (X: 99941091; C>G; p.H448D) of the SYTL4 gene was observed in a female with autism and random X-chromosome inactivation (46–54%). This patient additionally harbored four other autosomal missense gene mutations [16]. Missense mutations are of importance in understanding the structure or function of a protein since they usually occur in amino acid residues of structural or functional significance by changing the size or properties of the amino acid there by preventing the function of that protein [43,44]. The effect of such a mutation is additionally dependent on the sequence and structure context of the alteration [45]. Proteins fold according to minimum free energy [46]. Only correctly folded proteins can deliver the functional properties of a protein, and even minor changes in the size or properties of an amino acid side chain can alter or prevent the function of the protein [47]. On the other hand, even large deletions or insertions may be tolerated in numerous positions within a protein [48].
Protein function and interactions require both stability and specificity. Since most disease-causing mutations produce structural effects, the importance of a specific gene location and protein production is emphasized [47]. Structural information is needed to fully understand the effects and consequences of mutations, whether disease-causing or used purposefully to modify the properties of a protein. Three-dimensional structures and computer models have been used to elucidate disease mechanisms from specific amino acid substitutions [49,50]. For example, in a study of 4236 mutations from 436 genes, mutations at arginine and glycine residues are collectively responsible for about 30% of genetic diseases [43].

3.1.6. SYTL4 Amino Acid Change R (279) C in Exon 9: RAB-Binding Domain

The amino acid change in our study patient, SYTL4 R(279)C in exon 9 (Figure 1, Figure 3 and Figure 5), has been determined to severely affect the critical functioning of this gene’s encoded RAB protein binding region at its N-terminal [51,52], which is perhaps analogous to the effect of a significant change in Sytl4 gene expression observed among the altered protein networks in the mouse brain proteome in a mouse model of anxiety [12]. The SYTL4 RAB-Binding Domain within which lies D[Asp], at 60 takes part in the apparent salt bridge formation with R[Arg] at 279 only in the native protein configuration with the presence of mutant C[Cys] at 279 leading to the formation of an extended beta-pleated sheet (see Figure 3A,B).

3.1.7. Effect of the R[Arg]⇒C[Cys] Amino Acid Change at 279 on the functionality of the RAB-Binding Domain

Our analysis of the mutant and native SYTL4 protein structure models shows that in the native protein, arginine (R279) forms an apparent salt bridge with aspartic acid (D60) (Figure 2 and Figure 3A). This arginine is part of a large extended loop conformation that appears to be stabilized by the apparent salt bridge formed between it and Asp60 within the RAB-Binding Domain [aa 4-122] and the Ring Domain [aa 43-105] (Figure 2 and Figure 3A).
In our R (279) C mutant protein structure model (Figure 3B), cysteine (C279) is located amidst a beta sheet with no apparent hydrogen bonding contacts observed within distance constraints for such interactions. Our modeling suggests a potential structural change induced by the R (279) C variation eliminating the stabilizing Arg279-Asp60 salt bridge and leads to significant structural changes in the N-terminal half of SYTL4 (Figure 3B,C). This change could very well affect the functionality of the RAB-Binding Domain, not only for the canonical-full-length SYTL4 protein (isoform-1), but also for the truncated SYTL4 protein: isoform-2 (Figure 5).

3.1.8. Potentially Deleterious R (279) C Amino Acid Change “Likely” To Affect Its Neighboring Active Phosphorylation Sites

Given our modeling and profiling results, large conformational changes were seen proximal to the R279C amino acid variation within exon 9. Flanking regions do contain two known phosphorylatable serines (YTKS (@274) VIDLR (@279) P EDVVHESGS (@289) L) (Figure 5). The structure and function of both canonical and truncated isoforms of the SYTL4 protein may be affected, particularly considering their interactions at membrane surfaces.

3.1.9. Role of Arginine (R279) in SYTL4 Protein Structure and Function

The missense mutation identified in our study causes a change in the amino acid at 279 from R [Arginine] to C [Cysteine] (Figure 1 and Figure 2B). Arginine is a large polar amino acid with sidechains that prefers to reside in an aqueous environment and is found more commonly on the surface of a protein. SYTL4 is a protein that anchors to the cell membrane and is exposed to both sides of the membrane by aqueous environments of the luminal and cytoplasmic sides. The C2 domain of SYTL4 protein facilitates binding of this protein to cell membranes and is often found in the active centers of proteins that bind phosphorylated substrates [53,54].
The arginine (R279) residue is in-fact located in between two Serine(S) residues at 274 and 289- positions (YTKS (@274) VIDLR (@279) PEDVVHESGS (@289) L) (Figure 5) and is shown to undergo post-translational phosphorylation [55,56,57]. Arginine is capable of “salt bridging” by forming non-bonded and hydrogen-bonded paired electrostatic interactions between acidic carboxyl groups and basic amino groups in single or adjacent protein chains [58]. One important role of “salt bridging” is connecting protein subunits or joining two secondary structures to form quaternary structures where they can connect as many as five secondary structure units [58].

3.1.10. Arginine Disfavors Cysteine for Substitution

Given the above contrasting basic differences between arginine and cysteine in their properties and their role in protein structure and function, arginine disfavors cysteine for substitution, particularly in extracellular and membrane proteins, such as SYTL4. Since such substitution can be devastating to protein stability and function given the loss of arginine’s ability to create stabilizing hydrogen bonds, the substituted cysteine residue’s ability to alter the native three-dimensional conformation of the protein molecule is important [44]. Thus, the potential for a structural and functional change induced by the SYTL4 R279C amino acid change is consistent with the noted physiochemical differences between arginine and cysteine [44].

3.1.11. Dysfunction of Evolutionarily Conserved RAB-Binding GTPases Play a Role in Autism and Neuronal Disorders

Among the transcriptome-wide RAB family of genes, RAB2A, RAB11FIP5, RAB19, RAB39B, and RAB43 are also known ASD genes [59,60,61] (https://gene.sfari.org/database/human-gene/) and SYTL4 protein interacts with several other members of the RAB family of proteins including RAB3A, RAB8A, RAB27A, and RAB27B (Figure 4). Therefore, the defect in the RAB protein-binding region in the N-terminal half of the SYTL4 protein is due to the R (279) C amino acid variant and may be a causal factor for the high-functioning autism in our patient.
In addition, the SYTL4 gene directly interacts with Syntax Binding Protein 1 (STXBP1). Both SYTL4 and STXBP1 are known to interact with RAB3A. STXBP1 interaction with RAB3A promotes RAB3A dissociation from the vesicle membrane (https://www.uniprot.org/uniprot/P61764). SYTL4 protein additionally interacts with several other members of the RAB family of proteins, including RAB3A, RAB8A, RAB27A, and RAB27B. Further, RAB-binding domain of SYTL4 serves as a preferred effector binding site for the GTP-bound form of the RAB27A protein that regulates the exocytosis of secretary granules [62,63,64]. RAB27A binds to the N-terminus SLP homology domains 1 and 2 of SYTL4 protein and the C-terminal domain (Figure 1 and Figure 2A). It appears to play a role in the localization of RAB27A to specific sites in a cell [51,52,65]. Upregulation of RAB27A protein in basal forebrain neurons has been associated with mild cognitive impairment and Alzheimer’s disease [52].
Further, specific evolutionarily conserved RAB-binding GTPases function as regulators of membrane trafficking and binding and act as binary molecular switches turned on by binding GTP and off by hydrolyzing GTP to GDP [66]. Their dysfunction through mutations has also been shown to play a crucial role in causing diverse patho-physiologies including X-linked mental retardation associated with autism, epilepsy and macrocephaly. This suggests a major role for specific RAB- binding effector proteins, such as SYTL4 and interacting RAB-activating GTPases (RAB GTPases) in the maintenance of normal neuronal function [51,59,60,61,64,65,66].

3.1.12. Significance of Defect in RAB- Protein Binding Region of N-Terminal Half of SYTL4 Protein due to R (279) C Amino Acid Variant

The potentially deleterious R (279) C amino acid change could critically affect the RAB-binding domain of SYTL4 proteins. Our analysis of the mutant and native SYTL4 protein structure models shows that in the native protein, arginine (R279) forms an apparent salt bridge with aspartic acid (D60). This arginine is part of a large extended loop conformation that appears to be stabilized by the apparent salt bridge formed between it and Asp60 within the RAB-Binding Domain [aa 4-122] and the Ring Domain [aa 43-105].
In our R279C mutant model, cysteine (C279) is located amidst a beta sheet with no apparent hydrogen bonding contacts observed within distance constraints for such interactions. Our modeling thus suggests a potential structural change induced by the R279C variation which eliminates the stabilizing Arg279-Asp60 salt bridge and leads to significant structural changes in the N-terminal half of SYTL4 (Figure 3B). This could very well affect the functionality of the RAB-Binding Domain, not only for the canonical-full-length SYTL4 protein (isoform-1), but also for the truncated SYTL4 protein: isoform-2.

3.1.13. SYTL4-Protein-Protein Interactions with RAB27A and with Other RAB-Family of Genes

As seen in Figure 4, the SYTL4 protein directly interacts with several members of the RAB (Ras-Associated proteins in Brain) family of proteins, RAB3A, RAB8A, RAB27A, and RAB27B. Further, RAB- binding domain of the SYTL4 protein serves as a preferred effector binding site for the GTP-bound form of RAB27A protein that regulates the exocytosis of secretary granules [62,63,64]. Among the RAB- family of genes, RAB2A, RAB11FIP5, RAB19, RAB39B, and RAB43 are known ASD genes [59,60,61] (https://gene.sfari.org/database/human-gene/).

3.1.14. Upregulation of RAB27A Protein Associated with Mild Cognitive Impairment and Alzheimer Disease

RAB27A binds to the N-terminus SLP homology domains 1 and 2 of the SYTL4 protein, and the C-terminal domain seems to play a role in the localization of RAB27A to specific sites in a cell [51,52,65]. It important to note that upregulation of the RAB27A protein in basal forebrain neurons has been associated with mild cognitive impairment and Alzheimer disease [52].

3.1.15. Dysfunction of Conserved RAB-Binding GTPases Play a Role in X-Linked Mental Retardation with Autism

As indicated earlier, specific evolutionarily conserved RAB-binding GTPases do function as regulators of membrane trafficking and binding by acting as binary molecular switches that are turned on by binding GTP and off by hydrolyzing GTP to GDP. Their dysfunction through mutations has been shown to play a crucial role in causing diverse patho-physiologies including X-linked mental retardation associated with autism, epilepsy, and macrocephaly, suggesting a major role for specific RAB-binding effector proteins, such as the SYTL4, and interacting RAB-activating GTPases (RAB- GTPases) in the maintenance of normal neuronal function [51,52,59,60,61,66,67].

3.1.16. SYTL4 Gene Is Relevant to Neuronal System Function and Disorders

The SYTL4 gene is relevant for neuronal system development, function and behavior, and is implicated in neurological and psychological diseases (Entrez Gene ID # 94121 (Human); ID # 27359 (Mouse). The SYTL4 gene expression is down regulated in the dorsal raphe nucleus of patients with major depressive disorders [11]. In addition, in a mouse model of anxiety, significant changes in Sytl4 were observed among the altered protein networks in the brain proteome [12].

3.1.17. SYTL4 Protein Is Abundantly Expressed in the Bed Nucleus of Stria Terminalis and Is Upregulated in Male Brain

The bed nucleus of stria terminalis (BNST) is a heterogeneous complex limbic forebrain structure, which plays an important role in controlling autonomic, neuroendocrine, and behavioral responses, and is thought to serve as a key relay connecting limbic forebrain structures to hypothalamic and brainstem regions associated with autonomic and neuroendocrine functions [10,14]. Its control of physiological and behavioral activity is mediated by local action of numerous neurotransmitters [68]. Therefore, one could argue that the abundantly-expressed SYTL4 mutant protein in the bed nucleus of stria terminalis could potentially affect normal behavioral responses resulting in an autistic phenotype, given that Sytl4 is upregulated in the brain of male mice, specifically in the posteromedial area of the medial BNST [9].
Furthermore, in mouse brain, the expression of Sytl4 protein is sexually dimorphic due to the sex hormone (estrogen and testosterone)-specific control of this gene and the developmental influence of sex hormone can lead to enduring effects on brain and behavior [9,69]. Mice with targeted disruptions of the Sytl4 gene exhibit specific deficits in sex-specific behavior and deficits. Sytl4 is required for patterning male sexual behavior [9]. Given the increased risk estimates of ASD among males [70,71], the overall odds-ratio for ASD increases if the rare hemizygous mutation on the X-chromosome pertains to gene expression that is known to express in the brain [41].

3.1.18. Targeted Knockout Mutant Mammalian Phenotypes for Sytl4 Includes Abnormal Behavior and Abnormal Neurological Phenotype

The amino acid changes of SYTL4 R279C in exon 9 in our patient have been determined to severely affect the RAB protein-binding region at its N-terminal [51,52], which may be analogous to the partial or full effect of targeted Null/Knockout mutant Sytl4 mammalian phenotypes. Gomi et al. [16] found that a targeted (Null/Knockout) mutation for Granuphilin (Sytl4) in XY-male mouse with 129P2/OlaHsd- genetic background replaced exons 3–5: targeted mutation 1, Tetsuro Izumi (tm1Tiz): Grn−/Y (Sytl4 tm1Tiz/Y) particularly affected its endocrine and nervous systems with abnormal corticotroph morphology and hypersecretion of adrenocorticotropin [16,72]. Importantly, the same mammalian-targeted Sytl4 knock-outs in 129P2/OlaHsd*C3H/He (Sytl4tm1Tiz/Sytl4tm1Tiz) and in 129P2/OlaHsd (Sytl4tm1Tiz/Y) have been associated with abnormal behavior and neurological phenotypes encompassing alertness, behavioral response to light, circadian rhythm, cognition, consumption behavior, emotion/affect behavior, grooming behavior, impulsive behavior, motor capabilities/coordination/movement, sensory capabilities/reflexes/nociception, sheltering behavior, sleep behavior, vocalization, and social interaction [13,14].

3.1.19. Targeted Sytl4 Knock Out Mouse Model Studies Affirm That Defective SYTL4 Protein Function is Likely to Effectuate Neurological and Phenotypic Defects

A wide range of functions of the SYTL4 protein is noted for various cellular components, such as intracellular protein transport and positive regulation of protein secretion involving nucleoplasm, cytoplasm, endosome, centrosome, and plasma membrane which govern diverse molecular functions. These include protein and phospholipid binding, zinc and metal ion binding, clathrin binding, and neurexin protein binding. One could postulate that the observed low birth weight and height, as well as the down-slanting palpebral fissures, mild hypertelorism, thin upper lip, and pointed chin seen in our patient could relate to the abnormal SYTL4 protein. These phenotypic abnormalities are perhaps not observable in the Tetsuro Izumi (tm1Tiz) and Grn−/Y (Sytl4 tm1Tiz/Y) mice models [16,72], however, decreased body weight has been reported in targeted Sytl4 knock-out (Sytl4tm1Tiz/Sytl4tm1Tiz) mice with 129P2/OlaHsd*C3H/He genetic background. Thus, targeted Sytl4 knock out mouse model studies affirm that defective SYTL4 protein function may effectuate neurological and phenotypic defects.

3.1.20. SYTL4 Protein Directly Interacts with Proteins Known to Cause Autism

Our protein–protein interactions study of the SYTL4 gene shows that the SYTL4 protein directly interacts with three other proteins which are known to cause autism, namely, STX1A, SNAP25 and STXBP1 (Figure 4; https://gene.sfari.org/database/human-gene/). In addition, the SYTL4 protein directly interacts with FGF4, STX1B, SNAP29, SNAP23, SNAP47, RAB3A, RAB27A, RAB27B, RAB8A, STX4, STX19, STX3, STX11, and STX2 whose numerous transcriptome-wide isoforms are known as autism genes (https://gene.sfari.org/database/human-gene/).

3.1.21. SYTL4 Gene Sequence Shows Similarity to a Known Autism Gene: SYT1

Moreover, the SYTL4 gene sequence alignment shows similarity to the SYT1 gene which is a known autism gene (https://gene.sfari.org/database/human-gene/). Yet another larger Synaptotagmin gene family member, SYT17, is also a known ASD gene (https://gene.sfari.org/database/human-gene/).

3.1.22. SYTL4 Gene Sequence Alignment Shows Similarity to SYT1(Synaptotagmin 1) Gene Which is a Known ASD Gene

Although the SYTL4 gene is presented here as a candidate gene, its transcriptome-wide isoform-level related gene family members, such as SYT1, SYT17, and SYT3 are known ASD genes (https://gene.sfari.org/database/human-gene/). Moreover, the SYTL4 gene sequence alignment shows similarity to the SYT1 gene [70]. Synaptotagmins are integral membrane proteins of synaptic vesicles thought to serve as Ca (2+) sensors in the process of vesicular trafficking and exocytosis. Calcium binding to synaptotagmin-1 participates in triggering neurotransmitter release at the synapse [70].

3.1.23. Direct Protein-Protein STRING Interactions of the SYTL4 Gene with Other ASD Genes

More significantly, the SYTL4 protein directly interacts with three known ASD proteins, STX1A, STXBP1 and SNAP25 (Figure 4) (https://gene.sfari.org/database/human-gene/). STX1A (Syntaxin 1A (brain)) encodes a protein involved in the regulation of serotonergic and GABAergic systems and its expression is altered in autism [70]. Rare single gene mutations in STX1A have been implicated in ASD. This gene is located at 7q11.23. Common STX1A variants are nominally associated with high-functioning autism and Asperger syndrome (https://gene.sfari.org/database/human-gene/). This protein governs the release and uptake of extracellular vesicles in the nervous system and glial cells facilitating transcellular communication [73] and serves as a key molecule in ion channel regulation and synaptic exocytosis (https://gene.sfari.org/database/human-gene/).
The SYTL4 gene also directly interacts with Syntax binding protein 1 (STXBP1). A frameshift mutation in the STXBP1 gene has been implicated in a study of quartet families with autism spectrum disorder [74,75,76] and is a known autism gene (https://gene.sfari.org/database/human-gene/). SYTL4 recruits and binds STXBP1 to promote exocytosis [26] (http://www.genecards.org/cgi-bin/carddisp.pl?gene=STXBP1&keywords=STXBP]. It is important to note that like STXBP1, SYTL4 is also associated with early onset epileptic encephalopathy and both these genes are expressed in every part of the brain [http://www.genecards.org/). SYTL4 as well as STXBP1 are known to be expressed in the frontal cortex and brain and are essential for protein–protein interactions at synapses and the neurotransmitter release cycle in human neurons (http://www.genecards.org/). STXBP1 is overexpressed in the frontal cortex and even heterozygous mutations cause early onset epileptic encephalopathy, specifically through presynaptic impairment and autism [74,75,76]. Like STXBP1, SYTL4 is also associated with early onset epileptic encephalopathy and is expressed in brain regions (http://www.genecards.org/). Polymorphisms in the SNAP25 (Synaptosomal-associated protein, 25 kDa) gene are associated with ADHD. SNAP25+/− mice also exhibit hyperactivity and cognitive and social impairment [77]. In addition to the SYTL4 protein’s direct interactions with three known ASD proteins, STX1A, STXBP1 and SNAP25, it interacts with 14 other genes, namely, FGF4, STX1B, SNAP29, RAB3A, SNAP23, RAB27B, RAB8A, SNAP47, STX4, STX19, STX3, STX11, STX2, and RAB27A, whose numerous transcriptome-wide isoform-level multiple gene family members are known ASD genes (https://gene.sfari.org/database/human-gene/).

3.1.24. SYTL4- Molecular Pathways, Biological Processes and Molecular Functions

Significant SYTL4 gene-involved molecular pathways are the synaptic vesicle cycle, insulin secretion, AMPK-signaling pathway, and SNARE interactions (Table 1). The synaptic vesicle cycle pertains to synaptic vesicles that are filled with neurotransmitters by active transport, and the diseases that are associated with defects in this cycle are early infantile epileptic encephalopathy; centronuclear myopathy; episodic ataxias; and familial or sporadic hemiplegic migraine. Insulin secretion is regulated by several hormones and neurotransmitters, and the diseases associated with this pathway are defects in the degradation of ganglioside and type II diabetes mellitus. AMPK signaling acts as a sensor of cellular energy status while SNARE interactions mediate the docking of synaptic vesicles with the presynaptic membrane in neurons. Diseases associated with defective SNARE interactions are pseudohypoparathyroidism and cerebral dysgenesis, neuropathy, ichthyosis, and palmoplantar keratoderma or CEDNIK syndrome (Table 1).
The SYTL4 gene’s extensive biological processes (Table 2) pertain to RAB-protein signal transduction along with synaptic vesicle functions, neurotransmitter secretion, regulation of signaling, glutamate secretion, neuro-muscular synaptic transmission, and axonogenesis, which are all relevant for proper neuronal function, and thus also relevant to autism. The SYTL4 gene’s molecular functions (Table 3), biological processes (Table 2) and molecular pathways (Table 1) are indicative of its significant role in neuronal function that is meaningful in the causation of high-functioning autism in our proband. Moreover, in recent studies, disturbed SYTL4 gene function has been associated with neuropsychiatric disorders, such as autism, schizophrenia and depression as well as the immune system. It is also considered a major modulator of central nervous system function [11,12].

3.1.25. Synaptic Dysfunction in Neurodevelopmental Disorders Is Associated with Autism and Intellectual Disabilities

The SYTL4 modulates exocytosis of dense-core granules and secretion of hormones in the pancreas and the pituitary. It interacts with vesicles containing negatively charged phospholipids in a Ca (2+)-independent manner (http://www.genecards.org/). The significant SYTL4 molecular pathways are the synaptic vesicle cycle, insulin secretion, AMPK-signaling pathway, and SNARE interactions (Table 1). Like SYTL4, SYT1, STX1A, STXBP1 and SNAP25 genes play an important role in the extravesicular synaptic function of neuronal systems and neurotransmission (https://gene.sfari.org/database/human-gene/; www.genecards.org/). Thus, these proteins are generally involved in the functioning of the synaptic vesicles’ (or neurotransmitter vesicles) cycle to facilitate synaptic vesicle exocytosis by which a synaptic vesicle fuses with the plasma membrane of the pre-synaptic axon terminal and releases its contents in the synaptic cleft, which is essential for propagating nerve impulses between neurons and are constantly created by the cells to facilitate transcellular communication [73].
It should be noted that the SYTL4 mRNAs are found in extracellular vesicles and stimulate target cells to translate into active protein [78]. The SNAP25 gene encodes t-SNARE which is involved in the molecular regulation of neurotransmitter release that associates with proteins in vesicle docking and membrane fusion. It may play an important role in the synaptic function of specific neuronal systems [77] (https://gene.sfari.org/database/human-gene/). Similarly, STXBP1 protein is essential for neurotransmission and binds to syntaxin, a component of the synaptic vesicle fusion machinery (https://www.uniprot.org/uniprot/P61764).
Genes that regulate presynaptic processes ultimately affect neurotransmitter release when disturbed [79,80]. Synaptic dysfunction in neurodevelopmental disorders is associated with autism and intellectual disabilities [81,82,83,84] as noted by Baker et al. [80] who reported the first case of a rare missense variant (I368T) in the Synaptotagmin1 (SYT1) gene causing a human neurodevelopmental disorder. It has a dominant negative effect involving both synaptic vesicle exocytosis and endocytosis [80].
Significant molecular pathways, functions and biological processes (Table 1, Table 2 and Table 3) of the SYTL4 gene include synaptic vesicle cycle and fusion, exocytosis and neurotransmitter secretion, akin to the molecular functions of its directly interacting known ASD genes: STXBP1 and SNAP25. Therefore, deficiency in the functioning of the SYTL4 protein due to the R (279) C amino acid change affects the canonical structures of the SYTL4 gene as well as its shorter isoform affecting synaptic vesicle cycle and fusion, exocytosis, and neurotransmitter secretion likely to cause autism and intellectual disabilities.

3.1.26. ASD-Predictive MicroRNAs among Mouse Sytl4- Interacting MicroRNAs

Our analysis of 298 validated/predicted microRNA interactions with mouse Sytl4 gene [34,85,86,87,88] has identified three of five autism-associated serum miRs (60%), namely, miR181b-5p, miR320a, and miR130a-3p, which have good predictive power in serum [34]. Among these three miRs, miR320 has the greatest ASD predictive power reported in serum as well as in the superior temporal gyrus and cerebellar cortex of ASD individuals [34,35,37]. In addition, 5 of 13 (38%) miRs that were differentially expressed in ASD serum samples were predicted to be interacting with mouse Sytl4 gene [34]. Three of these five miRs, namely, miR130a-3p, miR181b-5p, and miR328 are predicted to be associated with the mouse Sytl4 gene. They have been shown to be differentially expressed in schizophrenics [31,32].

3.1.27. Dysregulation of miR-320—Most Predictive for ASD in Serum and Brain Tissues

Recently, miR-320, along with miR-197 in human follicular fluid were found to be associated with embryonic development potential [89]. Knocking down miR-320 in mouse oocytes negatively affects embryonic developmental potentially by inhibiting the expression of the Wnt-signaling pathway and therefore miRNAs in human follicular fluid might reflect an effect on embryo quality [89]. Furthermore, autopsy tissue sections showed concordantly dysregulated miR-320a and voltage-dependent anion channel 1 levels in HIV-1 patients suffering from mild cognitive impairment [90]. One could consider dysregulation of miR-320 as the most predictive for ASD in serum and brain tissues, since it was also found to be dysregulated in the superior temporal gyrus of ASD specimens [37]. Additionally, miR-320 and miR-197 are differentially expressed in the ASD cerebellar cortex [35,36]. Given the fact that miR-320 has been predicted to interact with mouse Sytl4 gene, it gives credence to the SYTL4 gene as a plausible new gene ASD.

3.1.28. SYTL4 Interacting miR181b-1- Being Predictive of ASD

The micro RNA second in line to miR320 is miR181b-1 as the most predictive for ASD. It not only showed good predictive power for ASD in serum [34] but is also differentially expressed in the ASD cerebellar cortex [35,36]. Unlike miR-320, miR181b-1 has additionally been found to show altered expression in the cortical regions in schizophrenia [32,33,34], supporting the contention that ASD and schizophrenia share common neurobiological features [32,33,34]. In a recent study, significant down-regulation of miRNA-181b expression in schizophrenics predicted improvement of negative symptoms to treatment, and thus miRNA-181b is predicted to serve as a potential plasma-molecular marker for antipsychotic responses [91].

3.1.29. SYTL4 Interacting miR130a-Being Predictive of ASD

Another potentially important micro RNA is miR130a (miR130a-3p) for prediction of ASD but is lacking corroboration from any ASD brain tissue studies (unlike miR320 and miR181b-1, above). However, like miRNA-181b, miR130a also showed altered expression (serum/cortical) in schizophrenia, again supporting the contention that ASD and schizophrenia share common neurobiological features [32,33,34].

3.1.30. SYTL4-Interacting miR106b and miR328 Dysregulated in ASD Cerebellar Cortex and Altered among Schizophrenics

Other SYTL4-associated miRs, miR106b (106b-3p) and miR328, have shown up and down regulation, respectively, in ASD serum and are differentially expressed in the ASD cerebellar cortex [34,35]. In addition, both miRs have also been shown to be altered among schizophrenics [32,33,34].

3.1.31. SYTL Interacting miR63, miR103, 5nd miR132 Are Dysregulated in Superior Temporal Gyrus of ASD

Other predicted SYTL4-interacting miRs, such as, miR93, miR103 (miR103-1 and miR103-2), and miR132, have been found to be dysregulated in superior temporal gyrus of ASD [37]. Thus, a total of eight microRNAs which are predicted to be associated with mouse Sytl4 are found to be altered in ASD serum and/or brain, thereby, augmenting our contention that the SYTL4 gene is a plausible new ASD candidate gene. Fifty percent (4/8) of the ASD-associated miRs (miRs106, miRs130a, miRs181b, and miRs328) that are predicated to interact with the mouse Sytl4 gene are also known to be associated with schizophrenia, supporting the contention that ASD and schizophrenia share common neurobiological features [34,92,93] (Table 4).
Fifty percent (4/8) of the ASD-associated mirRs (mirRs93, mirRs103, mirRs132, mirRs320) that are predicted to interact with mouse Sytl4 have also been determined to be dysregulated in superior temporal gyrus of ASD [37]. Yet again, 50% (4/8) of the ASD-associated mirRs (mirRs106b, mirRs181b-5p, mirRs320, mirRs328) predicted to interact with the mouse Sytl4 have also been determined to be differentially expressed in the ASD cerebellar cortex [35]. Over expression of miR142-5p, miR142-3p, miR451a, miR144-3p, and miR21-5p has been reported in ASD brain tissue along with hypomethylation of the promoter region of the miR142 gene in the same samples, suggesting dysregulation of these microRNAs [94]. However, the mouse Sytl4 gene has not been predicted to interact with any of the miRs that were found to be over expressed by Mor et al. [94]. Furthermore, these five miRs are not represented among the 13 differentially expressed miRs in ASD serum [34] (Table 4). Other studies have also not reported any of the five miRs to be over expressed by Mor et al., [94] based on specific regions of ASD brains, such as the superior temporal gyrus [37] and cerebellar cortex [35].

3.2. Transmembrane Protein 187 (TMEM187) Gene

TMEM187 is an emerging candidate gene for autism with a discussion undertaken on protein interactome networks, expression profiling and microRNA interaction studies.

3.2.1. TMEM187 Gene Belongs to a Group of Genes Which Host MicroRNA Genes in Their Introns or Exons

At the outset, it should be noted that unlike the SYTL4 gene, not much is known about the TMEM187 gene, and there is limited information available regarding its biological processes, molecular pathways and functions or microRNA interactions. The TMEM187 gene, like SYTL4, is a protein-coding gene, but belongs to a group of genes which host microRNA genes in their introns or exons [17]. However, we introduce the TMEM187 gene as an emerging candidate gene for autism with our mutation analysis of its novel missense variant c.708G>T; p. Gln236His, its STRING-protein interactome network and its expression profiling.

3.2.2. Novel TMEM187 Missense Variant c.708G>T: Glutamine(Q)236 Histidine(H)

Unlike the SYTL4 gene variant, TMEM187 Glutamine(Q)236 Histidine(H) variant in our patient is not found in the Exome Aggregation Consortium (ExAC) database. It is also not found in the listed 261 previous variants currently listed at ensembl.org.

3.2.3. Deleterious and Damaging Nature of the Novel TMEM187- Variant

Our extensive analyses of this novel variant, as detailed in the Results section, was often determined to be deleterious or damaging. Glutamine, which is a polar amino acid, was changed to histidine and this alters the protein [43,44].

3.2.4. TMEM187 Protein Is Expressed in Brain

Like SYTL4, the TMEM187 gene is ubiquitously expressed in all systems including all parts of the brain (www.uniprot.org/uniprot/Q14656; www.genecards.org/).

3.2.5. STRING–Gene Interaction Network Study Reveals Direct Protein–Protein Interactions of the TMEM187 Gene with Several Other Known Autism Genes

Although the novel X-linked TMEM187 missense gene variant c.708G>T; p.Gln236His found in our high-functioning autism patient is not known as an ASD gene, but analysis of the latest STRING network interactions reveal direct interactions with four known ASD genes, namely HCFC1, TMLHE, MECP2, and GPHN (https://gene.sfari.org/database/human-gene/) (Figure 7).

3.2.6. Significance of TMEM187 Protein-Protein Interacting Autism Genes

HCFC1 is a syndromic ASD gene that interacts with TMEM187, while TMLHE, MECP2 and GPHN are rare single gene autism genes (https://gene.sfari.org/database/human-gene/). HCFC1 is involved in control of the cell cycle with mutations in this X-linked (Xq28) gene associated with intellectual disability [26]. The two genes (TMEM187 and HCFC1) lie just 2kb apart [95]. Over expression of HCFC1 due to a variant is linked to intellectual disability [26]. A rare mutation in TMLHE has been identified with autism (https://gene.sfari.org/database/human-gene/).
The TMEM187 gene is in linkage with the MECP2 gene which is a well-known determinant of brain structure, and amino acid variations in the MECP2 protein cause micro-encephalopathy and are also associated with several neurodevelopmental disorders that affect both brain morphology and cognition [96]. Mutations in this gene underlie Rett syndrome, a well-known autism disorder (https://gene.sfari.org/database/human-gene/). Rare single gene mutations in the GPHN gene are associated with ASD and this gene encodes a neuronal assembly protein that anchors inhibitors of neurotransmitter receptors to postsynaptic cytoskeleton (https://gene.sfari.org/database/human-gene/).

3.2.7. Significance of Other Protein-Protein Interactions of TMEM187

The TMEM187 protein interacts directly with UBL4A, RBM25, and AKAP4 (Figure 7). Though these genes are not known ASD genes, their other transcriptome-wide isoform-level family member genes are known ASD genes: UBL7, RBM27, RBM8A, RBMS3 and AKAP9 (https://gene.sfari.org/database/human-gene/).
Additionally, the TMEM187 protein directly interacts with the LAGE3 (L Antigen Family Member 3) protein (Figure 7), however, LAGE3 is not a known ASD gene. It is associated with Galloway-Mowat syndrome 2, an X-linked early-onset nephrotic syndrome associated with microcephaly, central nervous system abnormalities, developmental delay, and a propensity for seizures. Brain anomalies include gyration defects such as lissencephaly, pachygyria, polymicrogyria, and cerebellar hypoplasia. Most patients show facial dysmorphism characterized by a small, narrow forehead, large/floppy ears, deep-set eyes, hypertelorism, and micrognathia (www.uniprot.org/uniprot/Q14657).

3.2.8. Other TMEM Proteins Gene Family Members Are Known Autism, Bipolar and Panic Disorder Genes

Other transmembrane proteins gene family members, such as TMEM231, are known syndromic autism genes (https://gene.sfari.org/database/human-gene/). This gene is associated with two neurological syndromes: Joubert syndrome-20 [MIM:614970] and Meckel syndrome 11 (https://gene.sfari.org/database/human-gene/; MIM:615397). Additionally, two other members (TMEM132E and TMEM132D) are known to be associated with bipolar and panic disorders [97,98].

3.2.9. X-chromosome Harbors Disproportionately Higher Number of TMEM187-Interacting Autism and Nervous System Disorder Genes: Implications for Boys vs Girls Ratio

Except for GPHN gene located on chromosome 14, all the other genes (TMEM187, HCFC1, TMLHE, MECP2, LAGE3 and SYTL4) are located exclusively on the long arm of the X-chromosome. The LAGE3 gene, whose family of genes are clustered together at Xq28, is like that of TMEM187, HCFC1, TMLHE, and MECP2. SYTL4 is also located on the X-chromosome at q22.1, proximal to the centromere. This augments the assertion that X-chromosome harbors a disproportionately higher number of ASD and nervous system disorder genes [99], and consequently, disproportionately affects more boys than girls, given that the overall odds-ratio for ASD is increased if the rare hemizygous mutation is on the X-chromosome and X-linked genes expressed in the brain [41] as is the case in all six of these X-linked genes (www.genecards.org).

4. Materials and Methods

4.1. Clinical Report

The 7-year-old male proband was the only child born to healthy young non-consanguineous parents. There was no family history of genetic disorders, malformations, epilepsy, autism, or intellectual disability. Our proband was the product of a 38-week pregnancy to a primigravida mother via an uncomplicated C-section due to a transverse presentation. The Apgar scores were 9 and 9 at 1 and 5 min, respectively. The birth weight was 3600 gm (55th percentile), length was 52 cm (85th percentile), and head circumference was 35 cm (55th percentile).
The proband was evaluated in the Department of Pediatric Neurology, Hospital Universitario Quirón, Madrid, Spain at the age of 4.5 years due to longstanding impairment in social and communicative functioning. Although an early intervention program was established in the first months of life for motor, cognitive, speech development, and social behavior, the proband exhibited mild psychomotor delay during his first years of life. He walked unsupported at 11 months but had significant problems with walking, squatting or dressing at 4.5 years of age. First bi-syllabic babbling occurred at 18 months; at the age of 3 years, he only spoke words without making sentences. His social development was markedly affected; he showed atypical behaviors, refused playing with other children, had a limited amount of interests, and eye contact was minimal.
His weight was 19 kg (65th percentile) and height was 109 cm (65th percentile). He had down-slanting palpebral fissures, mild hypertelorism, thin upper lip, and a pointed chin. Conventional genetic studies (karyotype and array comparative genomic hybridization) showed no abnormalities. The neurological exam was normal. He had impaired social interaction during the examination and lacked eye contact, had peculiar language (echolalias, verbosity, and abnormal pitch), stereotyped mannerisms, and restricted patterns of interest. Brain MRI and sleep video-EEG tests displayed normal results. Cognitive assessment using the Wechsler Preschool and Primary Test of Intelligence-III (WPPSI-III) revealed a verbal and non-verbal IQ at above-average level without significant discrepancies. The Behavior Assessment System for Children (BASC) completed by his parents and preschool teachers revealed significant problems in “social skills,” “adaptability” and “atypicality” domains.
At the age of 6 years, his neurological examination remained normal, but he had an unusually high-pitched voice with stereotyped phrases and echolalia. He tended to perseverate on repetitive interests and activities (chronology of history, borders of countries). His eye contact was inconsistent and poorly integrated with other communicative efforts. He reacted aversely to sensory stimuli (e.g., loud noise, flavors). Autism Spectrum Screening Questionnaire (ASSQ), the Autism Diagnostic Interview-Revised (ADI-R), and the Autism Diagnostic Observation Scale (ADOS, Module 3) were administered. His total score on the ASSQ was 31 and 29 according to the evaluation by the teachers and parents (high-functioning autism cut-off = 22 and 19, respectively). His ADI-R algorithm scores were 15 on the social domain (autism cut-off = 10), 11 on the communication domain for verbal children (autism cut-off = 8), and 5 on the repetitive behavior’s domain (autism cut-off = 3). His total score on the ADOS communication and social algorithm items were 6 and 9 (autism cut-off = 3 and 6, respectively). The clinical and neuropsychological evaluations were consistent with high-functioning autism. The patients allowed for us to undertake research investigations. The study was approved by the local ethics committees on January 7, 2016 and was conducted in accordance with the ethical principles of the Declaration of Helsinki and Good Clinical Practice standards. Informed consent was obtained from parents, with the child giving assent.

4.2. Genomic Investigations

Exome sequencing was performed using genomic DNA isolated (MagnaPure, Roche Applied Science, Manheim, Germany) from whole blood from the proband and parents. Libraries were prepared using the Ion AmpliSeq™ Exome Kit (Life Technologies, Carlsbad, California, USA) and quantified by qPCR. The enriched libraries were prepared using Ion Chef™ and sequenced on PI™ Chip in the Ion Proton™ System (Life Technologies) to provide >90% of amplicons covered with at least 20X. Signal processing, base calling, alignment and variant calling were performed on a Proton™ Torrent Server using the Torrent Suite™ Software (v4.4 Life Technologies, Carlsbad, CA, USA). Variants were annotated using Ion Reporter™ Software with the human genome reference assemble GRCh37 (hg19) and pedigree analysis performed using the Genetic Disease Screen (GDS) trio workflow.
Candidate variants were visualized using IGV (Integrative Genomics Viewer, Cambridge, MA, USA) and evaluated based on stringent assessments at both the gene and variant levels, taking into consideration both the patient’s phenotype and the inheritance pattern. Variants in the SYTL4 and TMEM187 genes were recognized as probable pathogenic and confirmed by Sanger sequencing. However, due to the lack of cooperation from other family members, additional testing was not available.

4.3. Modeling of Native and R279C Variant for SYTL4 Gene

The theoretical atomic models of native and R279C SYTL4 were constructed using I-TASSER [20,100,101]. With human SYTL4 sequence (UniProtKB Accession # = Q96C24) as query, multiple sequence-template alignments were initially generated by the meta-threading program LOMETS [20,102,103], followed by generation of the predicted atomic structures. Native SYTL4 matched well with several moderately high-scoring templates corresponding to synaptotagmin family members, with an estimated TM score of 0.5 ± 0.15 and RMSD of 12.2 ± 4.4 Å. The R279C variant also matched to synaptotagmin family member templates, with slightly lower scores, yielding a TM score and RMSD of 0.45 ± 0.15 and 13.5 ± 4.0 Å, respectively.
Both structures were close to the threshold (TM score >0.5) for correct topology. Molecular graphics and analyses were performed with the UCSF Chimera package. Chimera is developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco (supported by NIGMS P41-GM103311) [85].

4.4. MicroRNAs

MicroRNAs (miRNAs / miRs) play a key role in the transcriptional networks of the developing human brain, as regulators of gene expression. Autism spectrum disorder (ASD), being a complex neurodevelopmental disorder, is characterized by multiple deficits in communication, social interaction and behavior [34]. Vasu et al. [34] examined the serum expression profiles of 125 neurologically relevant miRNAs expression profiles in 55 individuals with ASD. These neurologically relevant miRNAs represented pathways involved in axon guidance, TGF-beta signaling, MAPK signaling, adherents’ junction, regulation of actin cytoskeleton, oxidative phosphorylation, hedgehog signaling, focal adhesion, mTOR signaling, and Wnt signaling [34].
Vasu et al. [34] found that only 13 miRNAs (out of the selected 125 miRs; ~10%) were differentially expressed among the 55 ASD individuals compared to the controls in their serum. Of these 13, miR151a-3p, miR181b-5p, miR320a, miR328, miR433, miR489, miR572, and miR663a were downregulated (61.5%), while miR101-3p, miR106b-5p, miR130a-3p, miR195-5p, and miR19b-3p were upregulated (38.5%) [34]. Furthermore, of these 13 miRs, only five miRs (38.5%), namely, miR181b-5p, miR320a, miR572, miR130a-3p and miR19b-3p had high values for sensitivity, specificity and area under the curve, thereby showing good predictive power for distinguishing individuals with ASD [34]. Therefore, it was decided to screen for the presence of these five miRNAs with good ASD predictive power, namely, miR181b-5p, miR320a, miR572, miR130a-3p and miR19b-3p as well as the other five miRs that were also differentially expressed among the 55 ASD individuals among the 298 validated/predicted microRNA interactions of mouse Sytl4 gene [86,87,88,104]. Additionally, we screened for the presence of other significantly ASD-associated miRs from brain-specific micro RNA expression studies by others [35,36,37].

5. Conclusions

1. TMEM187 as well as SYTL4 genes are: X-linked and both located on the long arm of the X-chromosome at Xq28 and Xq22.1, respectively;
2. Both Q236H TMEM187 and R279C SYTL4 gene variants have been predicted to be damaging or deleterious by SIFT, PolyPhen2, MutationTaster, Provean, and LRT variant calling programs;
3. Both TMEM187 and SYTL4 mRNAs are found in extracellular vesicles and stimulate target cells to translate into active protein [17], and the release and uptake of extracellular vesicles in the nervous system and glial cells provides novel mechanisms of transcellular communication [73];
4. Together, TMEM187 and SYTL4 genes directly interact with seven known ASD genes: four and three ASD genes, respectively (Figure 4 and Figure 7; https://gene.sfari.org/database/human-gene/);
5. Another transmembrane protein gene family member, TMEM231, is a known syndromic autism gene: TMEM231 gene is associated with two neurological syndromes: Joubert syndrome-20 [MIM:614970] and Meckel syndrome 11 (https://gene.sfari.org/database/human-gene/; MIM:615397);
6. The SYTl4-RAB-binding protein RAB27A is specifically associated with mild cognitive impairment and Alzheimer disease [52];
7. RAB-binding GTPases play a crucial role in causing diverse patho-physiologies including X- linked mental retardation (intellectual disability) associated with autism, epilepsy, and macrocephaly, suggesting a major role for specific RAB-binding effector proteins, such as SYTL4, and interacting RAB-activating GTPases (RAB GTPases) in the maintenance of normal neuronal function [51,59,60,61,64,65,66];
8. Two other members of the transmembrane protein gene family, TMEM132E and TMEM132D, are known to be associated with bipolar and panic disorders, respectively [97,98];
9. One of the TMEM187 genes’ related phenotypes is schizophrenia (GWAS catalog for TMEM187 gene: Gene relation via enhancers containing phenotype SNP: Enhancer ID: GH0XJ153980; https://genecards.weizmann.ac.il/geneloc-bin/display_map.pl?chr_nr=0X&range_type=gh_id&gh_id=GH0XJ153941#GH0XJ153941). Similarly, our extensive ASD predictive mouse Sytl4-interacting microRNAs study reveals that 50% of the ASD-associated miRs are known to be associated with schizophrenia (Table 4);
10. A recent large study of gene expression patterns from postmortem brain tissues found transcriptome-wide isoform-level dysregulation in ASD, schizophrenia, bipolar disorder, panic disorder, and other related neurological disorders [105], supporting our above findings in SYTL4 and TMEM187 gene variants that ASD and schizophrenia share common neurobiological features [34,92,93];
11. Both SYTL4 and the TMM187 genes are ubiquitously expressed, including in the brain [9,10,14] (www.GeneCards.Org; www.Uniprot.Org). The overall odds-ratio for ASD is increased if rare hemizygous mutations on the X-chromosome in male patients are found as these genes have known expression in the brain [41].
Given the above analytical analyses, there is evidence to support that the missense mutations seen in both the TMEM187 and SYTL4 genes, either synergistically or individually, are causal mutations for the high-functioning autism seen in our patient. Consequently, both genes are proposed as novel autism candidate genes.
It is probable that both gene variants synergistically are causative of the high-functioning autism seen in our patient. Oligogenic heterozygosity or involvement of more than one gene observed in patients, suggests a new potential mechanism in the pathogenesis of autism spectrum disorders [106], further supporting the suggestion that the multifactorial model of ASD risk or monogenic may be too simplistic even for the most penetrant causes of ASD [107]. Pathogenicity of mutations individually or synergistically would require biological assays, such as in-vivo models, for study. Meanwhile, publications of similar findings by others in the study of autism involving either of these two gene variants would lend credence to our findings and assertions.

Author Contributions

The manuscript was conceived and entirely drafted by S.K.R., including the title, figures, tables, microRNA study, interpretation, discussion, conclusion, materials and methods, and references. M.G.B. reviewed, contributed and edited the entire manuscript. The clinical studies were conducted by A.F.-J., Neuro-pediatrician at the Department of Pediatric Neurology, Hospital Universitario Quirón, Madrid, Spain. The clinical case report was by A.F.-J. The genomic studies were conducted by S.Á., at the Genomics and Medicine, NIM Genetics, Madrid, Spain. The theoretical 3D structural modeling of native SYTL4 as well as the analysis of predicted atomic threading models of the mutant and native SYTL4 proteins, and the interpretation thereof was by O.W.N., at the Department of Biochemistry & Molecular Biology, University of Kansas Medical Center. Kansas City, USA. S.K.R. was granted permissions to reproduce in this manuscript the following: (1) the gene interaction STRING networks (by M.S., Head of GeneCards Development, Weizmann Institute of Science, Israel (www.GeneCards.Org & STRING Consortium)) and (2) the SYTL4 protein secondary structure graphics by E.G., SIB Swiss Institute of Bioinformatics, Geneva, Switzerland (UniProtKB/Swiss-Prot; WWW.Uniprot.Org; Q96C24).

Funding

Owen W. Nadeau acknowledges the research funding support for this project from NIH grant # DK32953. Merlin G. Butler acknowledges National Institute of Child Health and Human Development (NICHD) grant HD02528.

Acknowledgments

Waheeda Hossain in the Departments of Psychiatry & Behavioral Sciences, University of Kansas Medical Center (Kansas City, KS., 66160 USA) reviewed and edited the manuscript. Charlotte Iannachi in the Departments of Psychiatry & Behavioral Sciences, University of Kansas Medical Center rendered valuable help with the manuscript processing. Mark Grubb (Kansas City, MO. 64171 USA) and Azhar Mohammed (Falls Church, VA. 22044 USA) assisted in the production of figures. We acknowledge the seminal contribution of the proband’s family to this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shashi, V.; McConkie-Rosell, A.; Rosell, B.; Schoch, K.; Vellore, K.; McDonald, M.; Jiang, Y.H.; Xie, P.; Need, A.; Goldstein, D.B. The utility of the traditional medical genetic diagnostic evaluation in the context of next-generation sequencing for undiagnosed genetic disorders. Genet. Med. 2014, 16, 176–182. [Google Scholar] [CrossRef] [PubMed]
  2. LePichon, J.B.; Saunders, C.; Soden, S.E. The future of next-generation sequencing in neurology. JAMA Neurol. 2015, 72, 971–972. [Google Scholar] [CrossRef] [PubMed]
  3. Xue, Y.; Ankala, A.; Wilcox, W.R.; Hegde, M.R. Solving the molecular diagnostic testing conundrum for Mendelian disorders in the era of next-generation sequencing: Single-gene, gene panel, or exome/genome sequencing. Genet. Med. 2015, 17, 444–451. [Google Scholar] [CrossRef] [PubMed]
  4. O’Roak, B.J.; Deriziotis, P.; Lee, C.; Vives, L.; Schwartz, J.J.; Girirajan, S.; Karakoc, E.; Mackenzie, A.P.; Ng, S.B.; Baker, C.; et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat. Genet. 2011, 43, 585–589. [Google Scholar] [CrossRef] [PubMed]
  5. Dixon-Salazar, T.J.; Silhavy, J.L.; Udpa, N.; Schroth, J.; Bielas, S.; Schaffer, A.E.; Olvera, J.; Bafna, V.; Zaki, M.S.; Abdel-Salam, G.H.; et al. Exome sequencing can improve diagnosis and alter patient management. Sci. Transl. Med. 2012, 4, 138ra78. [Google Scholar] [CrossRef] [PubMed]
  6. Sanders, S.J.; Murtha, M.T.; Gupta, A.R.; Murdoch, J.D.; Raubeson, M.J.; Willsey, A.J.; Ercan-Sencicek, A.G.; DiLullo, N.M.; Parikshak, N.N.; Stein, J.L.; et al. De novo mutations revealed by whole-exome sequencing are strongly associated with autism. Nature 2012, 485, 237–241. [Google Scholar] [CrossRef]
  7. Yu, T.W.; Chahrour, M.H.; Coulter, M.E.; Jiralerspong, S.; Okamura-Ikeda, K.; Ataman, B.; Schmitz-Abe, K.; Harmin, D.A.; Adli, M.; Malik, A.N.; et al. Using whole-exome sequencing to identify inherited causes of autism. Neuron 2013, 77, 259–273. [Google Scholar] [CrossRef]
  8. Vasieva, O.; Cetiner, S.; Savage, A.; Schumann, G.G.; Bubb, J.V.; Quinn, J.P. Primate specific retrotransposons, SVAs, in the evolution of networks that alter brain function in Neurons and Cognition. arXiv 2016, arXiv:1602.07642v2. [Google Scholar]
  9. Xu, X.; Coats, J.K.; Yang, C.F.; Wang, A.; Ahmed, O.M.; Alvarado, M.; Izumi, T.; Shah, N.M. Modular genetic control of sexually dimorphic behaviors. Cell 2012, 48, 596–607. [Google Scholar] [CrossRef]
  10. Lebow, M.; Chen, A. Overshadowed by the amygdala: The bed nucleus of the stria terminalis emerges as key to psychiatric disorders. Mol. Psychiatry 2016, 21, 450–463. [Google Scholar] [CrossRef]
  11. Kerman, I.A.; Bernard, R.; Bunney, W.E.; Jones, E.G.; Schatzberg, A.F.; Myers, R.M.; Barchas, J.D.; Akil, H.; Watson, S.J.; Thompson, R.C. Evidence for transcriptional factor dysregulation in the dorsal raphe nucleus of patients with major depressive disorder. Front. Neurosci. 2012, 6, 135. [Google Scholar] [CrossRef] [PubMed]
  12. Szego, E.M.; Janáky, T.; Szabó, Z.; Csorba, A.; Kompagne, H.; Müller, G.; Lévay, G.; Simor, A.; Juhász, G.; Kékesi, K.A. A mouse model of anxiety molecularly characterized by altered protein networks in the brain proteome. Eur. Neuropsychopharmacol. 2010, 20, 96–111. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, H.; Ishizaki, R.; Xu, J.; Kasai, K.; Kobayashi, E.; Gomi, H.; Izumi, T. The Rab27a effector exophilin7 promotes fusion of secretory granules that have not been docked to the plasma membrane. Mol. Biol. Cell 2013, 24, 319–330. [Google Scholar] [CrossRef] [PubMed]
  14. Xu, X. Modular genetic control of innate behaviors. Bioessays 2013, 35, 421–424. [Google Scholar] [CrossRef]
  15. Butler, M.G.; Rafi, S.K.; Hossain, W.; Stephan, D.A.; Manzardo, A.M. Whole exome sequencing in females with autism implicates novel and candidate genes. Int. J. Mol. Sci. 2015, 16, 1312–1335. [Google Scholar] [CrossRef] [PubMed]
  16. Gomi, H.; Mori, K.; Itohara, S.; Izumi, T. Rab27b is expressed in a wide range of exocytic cells and involved in the delivery of secretory granules near the plasma membrane. Mol. Biol. Cell 2007, 18, 4377–4386. [Google Scholar] [CrossRef] [PubMed]
  17. Boivin, V.; Deschamps-Francoeur, G.; Scott, M.S. Protein coding genes as hosts for noncoding RNA expression. Semin. Cell Dev. Biol. 2018, 75, 3–12. [Google Scholar] [CrossRef] [PubMed]
  18. Bonanno, J.B.; Rutter, M.; Bain, K.T.; Miller, S.; Romero, R.; Wasserman, S.; Sauder, J.M.; Burley, S.K.; Almo, S.C. Protein Data Bank in Europe (PDBe) ID: 3fdw: Crystal Structure of a C2 Domain from Human Synaptotagmin-Like Protein 4; New York SGX Research Center for Structural Genomics (NYSGXRC): New York, NY, USA, 2008. [Google Scholar]
  19. Miyamoto, K.; Sato, M.; Koshiba, S.; Inoue, M.; Kigawa, T.; Yokoyama, S. Protein Data Bank (PDB) ID: 2CSZ: Solution Structure of the RING Domain of the Synaptotagmin-Like Protein 4. RIKEN Structural Genomics/Proteomics Initiative (RSGI), 2005. Available online: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/Structure/pdb/2CSZ (accessed on 8 July 2018).
  20. Zhang, Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9, 40. [Google Scholar] [CrossRef] [PubMed]
  21. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The protein data bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed]
  22. Fuson, K.L.; Montes, M.; Robert, J.J.; Sutton, R.B. Structure of human synaptotagmin 1 C2AB in the absence of Ca2+ reveals a novel domain association. Biochemistry 2007, 46, 13041–13048. [Google Scholar] [CrossRef]
  23. Schauder, C.M.; Wu, X.; Saheki, Y.; Narayanaswamy, P.; Torta, F.; Wenk, M.R.; De Camilli, P.; Reinisch, K.M. Structure of a lipid-bound extended synaptotagmin indicates a role in lipid transfer. Nature 2014, 510, 552–555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Donald, J.E.; Kulp, D.W.; DeGrado, W.F. Salt bridges: Geometrically specific, designable interactions. Proteins 2011, 79, 898–915. [Google Scholar] [CrossRef]
  25. Gregoret, L.M.; Rade, S.D.; Fletterick, R.J.; Cohen, F.E. Hydrogen bonds involving sulfur atoms in proteins. Proteins 1991, 9, 99–107. [Google Scholar] [CrossRef] [PubMed]
  26. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef] [PubMed]
  27. Okada, M.; Itoh, M.; Haraguchi, M.; Okajima, T.; Inoue, M.; Oishi, H.; Matsuda, Y.; Iwamoto, T.; Kawano, T.; Fukumoto, S.; et al. b-series Ganglioside deficiency exhibits no definite changes in the neurogenesis and the sensitivity to Fas-mediated apoptosis but impairs regeneration of the lesioned hypoglossal nerve. J. Biol. Chem. 2002, 277, 1633–1636. [Google Scholar] [CrossRef] [PubMed]
  28. Stamberger, H.; Nikanorova, M.; Willemsen, M.H.; Accorsi, P.; Angriman, M.; Baier, H.; Benkel-Herrenbrueck, I.; Benoit, V.; Budetta, M.; Caliebe, A.; et al. STXBP1 Encephalopathy: A neurodevelopmental disorder including epilepsy. Neurology 2016, 86, 954–962. [Google Scholar] [CrossRef]
  29. Söllner, T.; Bennett, M.K.; Whiteheart, S.W.; Scheller, R.H.; Rothman, J.E. A protein assembly-disassembly pathway in vitro that may correspond to sequential steps of synaptic vesicle docking, activation, and fusion. Cell 1993, 75, 409–430. [Google Scholar] [CrossRef]
  30. Neves, S.R.; Ram, P.T.; Iyengar, R. G protein pathways. Science 2002, 296, 1636–1639. [Google Scholar] [CrossRef]
  31. Rapoport, J.; Chavez, A.; Greenstein, D.; Addington, A.; Gogtay, N. Autism spectrum disorders and childhood-onset schizophrenia: Clinical and biological contributions to a relation revisited. J. Am. Acad. Child Adolesc. Psychiatry 2009, 48, 10–18. [Google Scholar] [CrossRef]
  32. Beveridge, N.J.; Cairns, M.J. MicroRNA dysregulation in schizophrenia. Neurobiol. Dis. 2012, 46, 263–271. [Google Scholar] [CrossRef] [Green Version]
  33. Shi, W.; Du, J.; Qi, Y.; Liang, G.; Wang, T.; Li, S.; Xie, S.; Zeshan, B.; Xiao, Z. Aberrant expression of serum miRNAs in schizophrenia. J. Psychiatr. Res. 2012, 46, 198–204. [Google Scholar] [CrossRef] [PubMed]
  34. Vasu, M.; Anitha, A.; Thanseem, I.; Suzuki, K.; Yamada, K.; Takahashi, T.; Wakuda, T.; Iwata, K.; Tsujii, M.; Sugiyama, T.; et al. Serum microRNA profiles in children with autism. Mol. Autism 2014, 5, 40. [Google Scholar] [CrossRef] [PubMed]
  35. Abu-Elneel, K.; Liu, T.; Gazzaniga, F.S.; Nishimura, Y.; Wall, D.P.; Geschwind, D.H.; Lao, K.; Kosik, K.S. Heterogeneous dysregulation of microRNAs across the autism spectrum. Neurogenetics 2008, 9, 153–161. [Google Scholar] [CrossRef] [PubMed]
  36. Ghahramani Seno, M.M.; Hu, P.; Gwadry, F.G.; Pinto, D.; Marshall, C.R.; Casallo, G.; Scherer, S.W. Gene and miRNA expression profiles in autism spectrum disorders. Brain Res. 2011, 1380, 85–97. [Google Scholar] [CrossRef] [PubMed]
  37. Stamova, B.; Ander, B.P.; Barger, N.; Sharp, F.R.; Schumann, C.M. Specific regional and age-related small noncoding RNA expression patterns within superior temporal gyrus of typical human brains are less distinct in autism brains. J. Child Neurol. 2015, 30, 1930–1946. [Google Scholar] [CrossRef] [PubMed]
  38. Anney, R.; Klei, L.; Pinto, D.; Almeida, J.; Bacchelli, E.; Baird, G.; Bolshakova, N.; Bölte, S.; Bolton, P.F.; Bourgeron, T.; et al. Individual common variants exert weak effects on the risk for autism spectrum disorders. Hum. Mol. Genet. 2012, 21, 4781–4792. [Google Scholar] [CrossRef] [PubMed]
  39. Klei, L.; Sanders, S.J.; Murtha, M.T.; Hus, V.; Lowe, J.K.; Willsey, A.J.; Moreno-De-Luca, D.; Yu, T.W.; Fombonne, E.; Geschwind, D.; et al. Common genetic variants, acting additively, are a major source of risk for autism. Mol. Autism 2012, 3, 9. [Google Scholar] [CrossRef]
  40. Stein, J.L.; Parikshak, N.N.; Geschwind, D.H. Rare inherited variation in autism: Beginning to see the forest. Neuron 2013, 77, 209–211. [Google Scholar] [CrossRef]
  41. Crestani, C.C.; Alves, F.H.; Gomes, F.V.; Resstel, L.B.; Correa, F.M.; Herman, J.P. Mechanisms in the bed nucleus of the stria terminalis involved in control of autonomic and neuroendocrine functions: A review. Curr. Neuropharmacol. 2013, 11, 141–159. [Google Scholar] [CrossRef]
  42. Lopes, A.M.; Burgoyne, P.S.; Ojarikre, A.; Bauer, J.; Sargent, C.A.; Amorim, A.; Affara, N.A. Transcriptional changes in response to X chromosome dosage in the mouse: Implications for X inactivation and the molecular basis of Turner Syndrome. BMC Genom. 2010, 11, 82. [Google Scholar] [CrossRef]
  43. Vitkup, D.; Sander, C.; Church, G.M. The amino-acid mutational spectrum of human genetic disease. Genome Biol. 2003, 4, R72. [Google Scholar] [CrossRef] [PubMed]
  44. Betts, M.J.; Russell, R.B. Amino acid properties and consequences of substitutions. In Bioinformatics for Geneticists; Barnes, M.R., Gray, I.C., Eds.; John Wiley & Sons: New York, NY, USA, 2003. [Google Scholar]
  45. Khan, S.; Vihinen, M. Spectrum of disease-causing mutations in protein secondary structures. BMC Struct. Biol. 2007, 7, 56. [Google Scholar] [CrossRef] [PubMed]
  46. Poussu, E.; Vihinen, M.; Paulin, L.; Savilahti, H. Probing the alpha-complementing domain of E. coli beta-galactosidase with use of an insertional pentapeptide mutagenesis strategy based on Mu in vitro DNA transposition. Proteins 2004, 54, 681–692. [Google Scholar] [CrossRef]
  47. Vihinen, M.; Vetri, D.; Maniar, H.S.; Ochs, H.D.; Zhu, Q.; Vorechovský, I.; Webster, A.D.; Notarangelo, L.D.; Nilsson, L.; Sowadski, J.M.; et al. Structural basis for chromosome X-linked agammaglobulinemia: A tyrosine kinase disease. Proc. Natl. Acad. Sci. USA 1994, 91, 12803–12807. [Google Scholar] [CrossRef] [PubMed]
  48. Yue, P.; Li, Z.; Moult, J. Loss of protein structure stability as a major causative factor in monogenic disease. J. Mol. Biol. 2005, 353, 459–473. [Google Scholar] [CrossRef] [PubMed]
  49. Rong, S.B.; Vihinen, M. Structural basis of Wiskott-Aldrich syndrome causing mutations in the WH1 domain. J. Mol. Med. (Berl.) 2000, 78, 530–537. [Google Scholar] [CrossRef] [PubMed]
  50. Lappalainen, I.; Vihinen, M. Structural basis of ICF-causing mutations in the methyltransferase domain of DNMT3B. Protein Eng. 2002, 15, 1005–1014. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Strom, M.; Hume, A.N.; Tarafder, A.K.; Barkagianni, E.; Seabra, M.C. A family of Rab27-binding proteins. Melanophilin links Rab27a and myosin Va function in melanosome transport. J. Biol. Chem. 2002, 277, 25423–25430. [Google Scholar] [CrossRef] [PubMed]
  52. Mufson, E.J.; Counts, S.E.; Ginsberg, S.D. Gene expression profiles of cholinergic nucleus basalis neurons in Alzheimer’s disease. Neurochem. Res. 2002, 27, 1035–1048. [Google Scholar] [CrossRef] [PubMed]
  53. Zhang, D.; Aravind, L. Identification of novel families and classification of the C2 domain superfamily elucidate the origin and evolution of membrane targeting activities in eukaryotes. Gene 2010, 469, 18–30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lyakhova, T.A.; Knight, J.D. The C2 domains of granuphilin are high-affinity sensors for plasma membrane lipids. Chem. Phys. Lipids 2014, 182, 29–37. [Google Scholar] [CrossRef] [PubMed]
  55. Dephoure, N.; Zhou, C.; Villén, J.; Beausoleil, S.; Bakalarski, C.E.; Elledge, S.J.; Gygi, S.P. A quantitative atlas of mitotic phosphorylation. Proc. Natl. Acad. Sci. USA 2008, 105, 10762–10767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Bian, Y.; Song, C.; Cheng, K.; Dong, M.; Wang, F.; Huang, J.; Sun, D.; Wang, L.; Ye, M.; Zou, H. An enzyme assisted RP-RPLC approach for in-depth analysis of human liver phosphoproteome. J. Proteom. 2014, 96, 253–262. [Google Scholar] [CrossRef] [PubMed]
  57. Sharma, K.; D’Souza, R.C.; Tyanova, S.; Schaab, C.; Wiśniewski, J.R.; Cox, J.; Mann, M. Ultradeep human phosphoproteome reveals a distinct regulatory nature of Tyr and Ser/Thr-based signaling. Cell Rep. 2014, 8, 1583–1594. [Google Scholar] [CrossRef] [PubMed]
  58. Musafia, B.; Buchner, V.; Arad, D. Complex salt bridges in proteins: Statistical analysis of structure and function. J. Mol. Biol. 1995, 254, 761–770. [Google Scholar] [CrossRef] [PubMed]
  59. Mitra, S.; Cheng, K.W.; Mills, G.B. Rab GTPases implicated in inherited and acquired disorders. Semin. Cell Dev. Biol. 2010, 22, 57–68. [Google Scholar] [CrossRef] [PubMed]
  60. Zhen, Y.; Stenmark, H. Cellular functions of Rab GTPases at a glance. J. Cell Sci. 2015, 128, 3171–3176. [Google Scholar] [CrossRef] [PubMed]
  61. Zerial, M.; McBride, H. Rab proteins as membrane organizers. Nat. Rev. Mol. Cell Biol. 2001, 2, 107–117. [Google Scholar] [CrossRef] [PubMed]
  62. Fukuda, M.; Kanno, E. Analysis of the role of Rab27 effector Slp4-a/Granuphilin-a in dense-core vesicle exocytosis. Methods Enzymol. 2005, 403, 445–457. [Google Scholar]
  63. Izumi, T.; Gomi, H.; Torii, S. Functional analysis of Rab27a effector granuphilin in insulin exocytosis. Methods Enzymol. 2005, 403, 216–229. [Google Scholar]
  64. Chavas, L.M.; Ihara, K.; Kawasaki, M.; Torii, S.; Uejima, T.; Kato, R.; Izumi, T.; Wakatsuki, S. Elucidation of Rab27 recruitment by its effectors: Structure of Rab27a bound to Exophilin4/Slp2-a. Structure 2008, 16, 1468–1477. [Google Scholar] [CrossRef] [PubMed]
  65. Krzewski, K.; Cullinane, A.R. Evidence for defective Rab GTPase-dependent cargo traffic in immune disorders. Exp. Cell Res. 2013, 319, 2360–2367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Giannandrea, M.; Bianchi, V.; Mignogna, M.L.; Sirri, A.; Carrabino, S.; D’Elia, E.; Vecellio, M.; Russo, S.; Cogliati, F.; Larizza, L.; et al. Mutations in the small GTPase gene RAB39B are responsible for X-linked mental retardation associated with autism, epilepsy, and macrocephaly. Am. J. Hum. Genet. 2010, 86, 185–195. [Google Scholar] [CrossRef] [PubMed]
  67. Goitre, L.; Trapani, E.; Trabalzini, L.; Retta, S.F. The Ras superfamily of small GTPases: The unlocked secrets. Methods Mol. Biol. 2014, 1120, 1–18. [Google Scholar] [PubMed]
  68. Yang, C.F.; Shah, N.M. Representing sex in the brain, one module at a time. Neuron 2014, 8, 261–278. [Google Scholar] [CrossRef] [PubMed]
  69. Palmer, N.; Beam, A.; Agniel, D.; Eran, A.; Manrai, A.; Spettell, C.; Steinberg, G.; Mandl, K.; Fox, K.; Nelson, S.F.; et al. Association of Sex with Recurrence of Autism Spectrum Disorder among Siblings. JAMA Pediatr. 2017, 171, 1107–1112. [Google Scholar] [CrossRef] [PubMed]
  70. Durdiakova, J.; Warrier, V.; Banerjee-Basu, S.; Baron-Cohen, S.; Chakrabarti, B. STX1A and Asperger syndrome: A replication study. Mol. Autism 2014, 5, 14. [Google Scholar] [CrossRef] [PubMed]
  71. Halladay, A.K.; Bishop, S.; Constantino, J.M.; Daniels, A.M.; Koenig, K.; Palmer, K.; Messinger, D.; Pelphrey, K.; Sanders, S.J.; Singer, A.T.; et al. Sex and gender differences in autism spectrum disorder: Summarizing evidence gaps and identifying emerging areas of priority. Mol. Autism 2015, 6, 36. [Google Scholar] [CrossRef] [PubMed]
  72. Gomi, H.; Mizutani, S.; Kasai, K.; Itohara, S.; Izumi, T. Granuphilin molecularly docks insulin granules to the fusion machinery. J. Cell Biol. 2005, 171, 99–109. [Google Scholar] [CrossRef] [PubMed]
  73. Budnik, V.; Ruiz-Cañada, C.; Wendler, F. Extracellular vesicles round off communication in the nervous system. Nat. Rev. 2016, 17, 160–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Van Breevoort, D.; Snijders, A.P.; Hellen, N.; Weckhuysen, S.; van Hooren, K.W.; Eikenboom, J.; Valentijn, K.; Fernandez-Borja, M.; Ceulemans, B.; De Jonghe, P.; et al. STXBP1 promotes Weibel-Palade body exocytosis through its interaction with the Rab27A effector Slp4-a. Blood 2014, 123, 3185–3194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Yuen, R.K.; Thiruvahindrapuram, B.; Merico, D.; Walker, S.; Tammimies, K.; Hoang, N.; Chrysler, C.; Nalpathamkalam, T.; Pellecchia, G.; Liu, Y.; et al. Whole-genome sequencing of quartet families with autism spectrum disorder. Nat. Med. 2015, 21, 185–191. [Google Scholar] [CrossRef] [PubMed]
  76. Franceschini, A.; Lin, J.; von Mering, C.; Jensen, L.J. SVD-phy: Improved prediction of protein functional associations through singular value decomposition of phylogenetic profiles. Bioinformatics 2016, 32, 1085–1087. [Google Scholar] [CrossRef] [PubMed]
  77. Braida, D.; Guerini, F.R.; Ponzoni, L.; Corradini, I.; De Astis, S.; Pattini, L.; Bolognesi, E.; Benfante, R.; Fornasari, D.; Chiappedi, M.; et al. Association between SNAP-25 gene polymorphisms and cognition in autism: Functional consequences and potential therapeutic strategies. Transl. Psychiatry 2015, 5, e500. [Google Scholar] [CrossRef] [PubMed]
  78. Hong, B.S.; Cho, J.H.; Kim, H.; Choi, E.J.; Rho, S.; Kim, J.; Kim, J.H.; Choi, D.S.; Kim, Y.K.; Hwang, D.; et al. Colorectal cancer cell-derived microvesicles are enriched in cell cycle-related mRNAs that promote proliferation of endothelial cells. BMC Genom. 2009, 10, 556. [Google Scholar] [CrossRef] [PubMed]
  79. Waites, C.L.; Garner, C.C. Presynaptic function in health and disease. Trends Neurosci. 2011, 34, 326–337. [Google Scholar] [CrossRef]
  80. Baker, K.; Gordon, S.L.; Grozeva, D.; van Kogelenberg, M.; Roberts, N.Y.; Pike, M.; Blair, E.; Hurles, M.E.; Chong, W.K.; Baldeweg, T.; et al. Identification of a human synaptotagmin-1 mutation that perturbs synaptic vesicle cycling. J. Clin. Investig. 2015, 125, 1670–1678. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Geppert, M.; Goda, Y.; Hammer, R.E.; Li, C.; Rosahl, T.W.; Stevens, C.F.; Südhof, T.C. Synaptotagmin I: A major Ca2+ sensor for transmitter release at a central synapse. Cell 1994, 79, 717–727. [Google Scholar] [CrossRef]
  82. Bai, J.; Wang, P.; Chapman, E.R. C2A activates a cryptic Ca (2+)-triggered membrane penetration activity within the C2B domain of synaptotagmin I. Proc. Natl. Acad. Sci. USA 2002, 99, 1665–1670. [Google Scholar] [CrossRef]
  83. Zoghbi, H.Y.; Bear, M.F. Synaptic dysfunction in neurodevelopmental disorders associated with autism and intellectual disabilities. Cold Spring Harb. Perspect. Biol. 2012, 4, a009886. [Google Scholar] [CrossRef]
  84. Bacaj, T.; Wu, D.; Yang, X.; Morishita, W.; Zhou, P.; Xu, W.; Malenka, R.C.; Südhof, T.C. Synaptotagmin-1 and synaptotagmin-7 trigger synchronous and asynchronous phases of neurotransmitter release. Neuron 2013, 80, 947–959. [Google Scholar] [CrossRef] [PubMed]
  85. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera--a visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef]
  86. Plaisance, V.; Abderrahmani, A.; Perret-Menoud, V.; Jacquemin, P.; Lemaigre, F.; Regazzi, R. MicroRNA-9 controls the expression of Granuphilin/Slp4 and the secretory response of insulin-producing cells. J. Biol. Chem. 2006, 281, 26932–26942. [Google Scholar] [CrossRef] [PubMed]
  87. Rolland, T.; Taşan, M.; Charloteaux, B.; Pevzner, S.J.; Zhong, Q.; Sahni, N.; Yi, S.; Lemmens, I.; Fontanillo, C.; Mosca, R.; et al. A proteome-scale map of the human interactome network. Cell 2014, 159, 1212–1226. [Google Scholar] [CrossRef] [PubMed]
  88. Eppig, J.T.; Blake, J.A.; Bult, C.J.; Kadin, J.A.; Richardson, J.E. Mouse Genome Database Group. The mouse genome database (MGD): Facilitating mouse as a model for human biology and disease. Nucleic Acids Res. 2015, 43, D726–D736. [Google Scholar] [CrossRef] [PubMed]
  89. Feng, R.; Sang, Q.; Zhu, Y.; Fu, W.; Liu, M.; Xu, Y.; Shi, H.; Xu, Y.; Qu, R.; Chai, R.; et al. MiRNA-320 in the human follicular fluid-development in vitro. Sci. Rep. 2015, 5, 8689. [Google Scholar] [CrossRef] [PubMed]
  90. Fatima, M.; Prajapati, B.; Saleem, K.; Kumari, R.; Mohindar Singh Singal, C.; Seth, P. Novel insights into role of miR-320a-VDAC1 axis in astrocyte-mediated neuronal damage in neuroAIDS. Glia 2017, 65, 250–263. [Google Scholar] [CrossRef] [PubMed]
  91. Song, H.T.; Sun, X.Y.; Zhang, L.; Zhao, L.; Guo, Z.M.; Fan, H.M.; Zhong, A.F.; Niu, W.; Dai, Y.H.; Zhang, L.Y.; et al. A preliminary analysis of association between the downregulation of microRNA-181b expression and symptomatology improvement in schizophrenia patients before and after antipsychotic treatment. J. Psychiatr. Res. 2014, 54, 134–140. [Google Scholar] [CrossRef] [PubMed]
  92. Khanzada, N.S.; Butler, M.G.; Manzardo, A.M. GeneAnalytics pathway analysis and genetic overlap among autism spectrum disorder, bipolar disorder and schizophrenia. Int. J. Mol. Sci. 2017, 18, 527. [Google Scholar] [CrossRef] [PubMed]
  93. Sundararajan, T.; Manzardo, A.M.; Butler, M.G. Functional analysis of schizophrenia genes using GeneAnalytics program and integrated databases. Gene 2018, 641, 25–34. [Google Scholar] [CrossRef] [PubMed]
  94. Mor, M.; Nardone, S.; Sams, D.S.; Elliott, E. Hypomethylation of miR-142 promoter and upregulation of microRNAs that target the oxytocin receptor gene in the autism prefrontal cortex. Mol. Autism 2015, 6, 46. [Google Scholar] [CrossRef] [PubMed]
  95. Frattini, S.A.; Zucchi, I.; Patrosso, C.P. Characterization sequence/EST database screening. Genomics 1996, 34, 323–327. [Google Scholar]
  96. Joyner, A.H.; Cooper Roddey, J.; Bloss, C.S.; Bakken, T.E.; Rimol, L.M.; Melle, I.; Agartz, I.; Djurovic, S.; Topol, E.J.; Schork, N.J.; et al. A common MECP2 haplotype associates with reduced cortical surface area in humans in two independent populations. Proc. Natl. Acad. Sci. USA 2009, 106, 15483–15488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Sklar, P.; Smoller, J.W.; Fan, J.; Ferreira, M.A.; Perlis, R.H.; Chambert, K.; Nimgaonkar, V.L.; McQueen, M.B.; Faraone, S.V.; Kirby, A.; et al. Whole-genome association study of bipolar disorder. Mol. Psychiatry 2008, 13, 558–569. [Google Scholar] [CrossRef] [PubMed]
  98. Erhardt, A.; Akula, N.; Schumacher, J.; Czamara, D.; Karbalai, N.; Müller-Myhsok, B.; Mors, O.; Borglum, A.; Kristensen, A.S.; Woldbye, D.P.; et al. Replication and meta-analysis of TMEM132D gene variants in panic disorder. Transl. Psychiatry 2012, 2, e156. [Google Scholar] [CrossRef] [PubMed]
  99. Butler, M.G.; Rafi, S.K.; Manzardo, A.M. High-resolution chromosome ideogram representation of currently recognized genes for autism spectrum disorders. Int. J. Mol. Sci. 2015, 16, 6464–6495. [Google Scholar] [CrossRef] [PubMed]
  100. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef]
  101. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef]
  102. Wu S, Zhang Y: LOMETS: A local meta-threading-server for protein structure prediction. Nucleic Acids Res. 2007, 35, 3375–3382. [CrossRef]
  103. Rual, J.F.; Venkatesan, K.; Hao, T.; Hirozane-Kishikawa, T.; Dricot, A.; Li, N.; Berriz, G.F.; Gibbons, F.D.; Dreze, M.; Ayivi-Guedehoussou, N.; et al. Towards a proteome-scale map of the human protein-protein interaction network. Nature 2005, 437, 1173–1178. [Google Scholar] [CrossRef]
  104. Smith, C.M.; Finger, J.H.; Hayamizu, T.F.; McCright, I.J.; Xu, J.; Berghout, J.; Campbell, J.; Corbani, L.E.; Forthofer, K.L.; Frost, P.J. The mouse gene expression database (GXD): 2014 update. Nucleic Acids Res. 2014, 42, D818–D824. [Google Scholar] [CrossRef] [PubMed]
  105. Gandal, M.J.; Haney, J.R.; Parikshak, N.N.; Leppa, V.; Ramaswami, G.; Hartl, C.; Schork, A.J.; Appadurai, V.; Buil, A.; Werge, T.M.; et al. Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap. Science 2018, 359, 693–697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Schaaf, C.P.; Sabo, A.; Sakai, Y.; Crosby, J.; Muzny, D.; Hawes, A.; Lewis, L.; Akbar, H.; Varghese, R.; Boerwinkle, E.; et al. Oligogenic heterozygosity in individuals with high-functioning autism spectrum disorders. Hum. Mol. Genet. 2011, 20, 3366–3375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Guo, H.; Wang, T.; Wu, H.; Long, M.; Coe, B.P.; Li, H.; Xia, K. Inherited and multiple de novo mutations in autism/developmental delay risk genes suggest a multifactorial model. Mol. Autism 2018, 9, 64. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graphic representation of the secondary structure of SYTL4 gene indicating the location of Rab binding domain, ring domain, C2 domains, and the location of arginine at amino acid position 279.
Figure 1. Graphic representation of the secondary structure of SYTL4 gene indicating the location of Rab binding domain, ring domain, C2 domains, and the location of arginine at amino acid position 279.
Ijms 20 03358 g001
Figure 2. The native SYTL4 protein structure contains the C2 domains, ring domain, and an apparent salt bridge between arginine (R279) and aspartic acid (D60) as noted in (A). Arginine 279 (side chain in red) is part of a large extended loop conformation that appears to be stabilized by an apparent salt bridge formed between it and Asp60 (side chain in Salmon) in the ring domain. The distance (3.116Å) calculated between the Arg guanidinium nitrogen and Asp carboxyl oxygen is well within the threshold distance observed for salt bridges in a comprehensive survey of crystal structures as noted in (B).
Figure 2. The native SYTL4 protein structure contains the C2 domains, ring domain, and an apparent salt bridge between arginine (R279) and aspartic acid (D60) as noted in (A). Arginine 279 (side chain in red) is part of a large extended loop conformation that appears to be stabilized by an apparent salt bridge formed between it and Asp60 (side chain in Salmon) in the ring domain. The distance (3.116Å) calculated between the Arg guanidinium nitrogen and Asp carboxyl oxygen is well within the threshold distance observed for salt bridges in a comprehensive survey of crystal structures as noted in (B).
Ijms 20 03358 g002
Figure 3. Effect of the R[Arg]⇒C[Cys] amino acid change at 279 on the structure of SYTL4 Proteins. Theoretical 3D structures of native and R279C SYTL4. Hierarchical protein structural modeling of both native (A) and R279C (B) SYTL4 proteins were carried out using I-TASSER. Three X-ray crystal structures of human Synaptotagmin C2 domains (Protein Data Bank ID: 2R83, 4P42 and 3HN8) among the top 10 templates used for both SYTL4 forms. Native SYTL4 (grey ribbon trace) and R279C variant (gold ribbon trace) are shown with C2(1), C2(2) and ring domains rendered as light green, dark green and cyan traces, respectively, in Panels (A,B). Panel (A): SYTL4 native protein structure indicating the C2 domains, ring domain, and an apparent salt bridge between arginine (R279) and aspartic acid (D60). Shown magnified in the native SYTL4 structure is the apparent salt bridge formed between D60 and R279. Panel (B): SYTL4 R279C variant (gold ribbon trace) shows significant displacement of both the 279 amino acid site (salmon ribbon trace) and ring domains of native (cyan ribbon trace). Panel (C): An overlay of native (A) and R279C (B) SYTL4 structures demonstrates relatively good alignment for the C terminal C2 domains and confirms the significant displacement of both the 279 amino acid site (salmon ribbon trace) and ring domains of native (cyan ribbon trace) SYTL4.
Figure 3. Effect of the R[Arg]⇒C[Cys] amino acid change at 279 on the structure of SYTL4 Proteins. Theoretical 3D structures of native and R279C SYTL4. Hierarchical protein structural modeling of both native (A) and R279C (B) SYTL4 proteins were carried out using I-TASSER. Three X-ray crystal structures of human Synaptotagmin C2 domains (Protein Data Bank ID: 2R83, 4P42 and 3HN8) among the top 10 templates used for both SYTL4 forms. Native SYTL4 (grey ribbon trace) and R279C variant (gold ribbon trace) are shown with C2(1), C2(2) and ring domains rendered as light green, dark green and cyan traces, respectively, in Panels (A,B). Panel (A): SYTL4 native protein structure indicating the C2 domains, ring domain, and an apparent salt bridge between arginine (R279) and aspartic acid (D60). Shown magnified in the native SYTL4 structure is the apparent salt bridge formed between D60 and R279. Panel (B): SYTL4 R279C variant (gold ribbon trace) shows significant displacement of both the 279 amino acid site (salmon ribbon trace) and ring domains of native (cyan ribbon trace). Panel (C): An overlay of native (A) and R279C (B) SYTL4 structures demonstrates relatively good alignment for the C terminal C2 domains and confirms the significant displacement of both the 279 amino acid site (salmon ribbon trace) and ring domains of native (cyan ribbon trace) SYTL4.
Ijms 20 03358 g003
Figure 4. STRING- Protein-Protein Interaction Network of SYTL4. [https://version11.string-db.org/cgi/network.pl?taskId=f7upTuHlbV0A]. SYTL4- Interacting ASD Genes; SYTL4: Synaptotagmin-like protein 4; modulates exocytosis of dense-core granules and secretion of hormones in the pancreas and the pituitary. Interacts with vesicles containing negatively charged phospholipids in a Ca (2+)-independent manner; synaptotagmin-like tandem C2 proteins. STXBP1: Syntaxin-binding protein 1; may participate in the regulation of synaptic vesicle docking and fusion. Essential for neurotransmission and binds syntaxin, a component of the synaptic vesicle fusion machinery (Score: 0.746). STX1A: Syntaxin-1A; plays a role in hormone and neurotransmitter exocytosis. Potentially involved in docking of synaptic vesicles at presynaptic active zones (Score: 0.805). SNAP25: Synaptosomal-associated protein 25; t-SNARE involved in the molecular regulation of neurotransmitter release. May play an important role in the synaptic function of specific neuronal systems (Score: 0.551).
Figure 4. STRING- Protein-Protein Interaction Network of SYTL4. [https://version11.string-db.org/cgi/network.pl?taskId=f7upTuHlbV0A]. SYTL4- Interacting ASD Genes; SYTL4: Synaptotagmin-like protein 4; modulates exocytosis of dense-core granules and secretion of hormones in the pancreas and the pituitary. Interacts with vesicles containing negatively charged phospholipids in a Ca (2+)-independent manner; synaptotagmin-like tandem C2 proteins. STXBP1: Syntaxin-binding protein 1; may participate in the regulation of synaptic vesicle docking and fusion. Essential for neurotransmission and binds syntaxin, a component of the synaptic vesicle fusion machinery (Score: 0.746). STX1A: Syntaxin-1A; plays a role in hormone and neurotransmitter exocytosis. Potentially involved in docking of synaptic vesicles at presynaptic active zones (Score: 0.805). SNAP25: Synaptosomal-associated protein 25; t-SNARE involved in the molecular regulation of neurotransmitter release. May play an important role in the synaptic function of specific neuronal systems (Score: 0.551).
Ijms 20 03358 g004
Figure 5. SYTL4 amino acids sequence data for the “canonical” form and its truncated isoform 2, both encompassing the alternating exon 9 with the missense mutation resulting in R[Arg]⇒C[Cys] at 279: Amino acid sequences that are highlighted in blue are the first 335 amino acids which are common to both the full length ‘canonical’ sequence of the SYTL4 protein and its truncated isoform 2, both containing the alternating exon 9 Glu(E) 270–Met(M) 304, highlighted in yellow with the R[Arg]⇒C[Cys] at 279. The location of native arginine residue (at 279) in-between two Serine(S) amino acid residues at 274 and 289 positions (which have been shown to undergo post-translational phosphorylation) is underlined. The RAB-Binding Domain (comprised of amino acids 4 through 122) is also indicated with underline, within which lies D[Asp] at position 60 and takes part in the apparent salt bridge formation with R[Arg] at 279 in the native protein configuration. The presence of mutant C[Cys] at 279 leads to the formation of an extended beta-pleated sheet therein instead (as seen in Figure 3A,B).
Figure 5. SYTL4 amino acids sequence data for the “canonical” form and its truncated isoform 2, both encompassing the alternating exon 9 with the missense mutation resulting in R[Arg]⇒C[Cys] at 279: Amino acid sequences that are highlighted in blue are the first 335 amino acids which are common to both the full length ‘canonical’ sequence of the SYTL4 protein and its truncated isoform 2, both containing the alternating exon 9 Glu(E) 270–Met(M) 304, highlighted in yellow with the R[Arg]⇒C[Cys] at 279. The location of native arginine residue (at 279) in-between two Serine(S) amino acid residues at 274 and 289 positions (which have been shown to undergo post-translational phosphorylation) is underlined. The RAB-Binding Domain (comprised of amino acids 4 through 122) is also indicated with underline, within which lies D[Asp] at position 60 and takes part in the apparent salt bridge formation with R[Arg] at 279 in the native protein configuration. The presence of mutant C[Cys] at 279 leads to the formation of an extended beta-pleated sheet therein instead (as seen in Figure 3A,B).
Ijms 20 03358 g005
Figure 6. (A). Graphic representation of TMEM187 protein indicating the location of Q236H amino acid variation beyond the last (6th) transmembrane helical domain (191-210aa), at the near terminal end of the Pfam domain (8-245aa) (https://pfam.xfam.org/family/tmem187; TMEM187 (PF15100)). (B). Theoretical partial 3D protein structure of the encoded multi-pass transmembrane protein targeting the 15-200aa region (https://modbase.compbio.ucsf.edu/), excluding the Q236H amino acid variant near the terminal end region of the protein which is represented by multiple sequence alignments and hidden Markov models preventing comparative structural analysis of the variant harboring region.
Figure 6. (A). Graphic representation of TMEM187 protein indicating the location of Q236H amino acid variation beyond the last (6th) transmembrane helical domain (191-210aa), at the near terminal end of the Pfam domain (8-245aa) (https://pfam.xfam.org/family/tmem187; TMEM187 (PF15100)). (B). Theoretical partial 3D protein structure of the encoded multi-pass transmembrane protein targeting the 15-200aa region (https://modbase.compbio.ucsf.edu/), excluding the Q236H amino acid variant near the terminal end region of the protein which is represented by multiple sequence alignments and hidden Markov models preventing comparative structural analysis of the variant harboring region.
Ijms 20 03358 g006
Figure 7. STRING-Protein–Protein Interactions of TMEM187. (https://string-db.org/network/9606.ENSP00000358999). TMEM187- Interacting ASD Genes: HCFC1: Host cell factor 1: Involved in control of the cell cycle; Coactivator for EGR2 and GABP2 (Score: 0.747); TMLHE: Trimethyllysine dioxygenase, mitochondrial: Converts trimethyllysine (TML) into Hydroxytrimethyllysine (HTML). (Score: 0.657); MECP2: Methyl-CpG-binding protein 2: Chromosomal protein that binds to methylated DNA; Mediates transcriptional repression through interaction with histone deacetylase and the corepressor SIN3A (Score: 0.555); GPHN: Gephyrin: Microtubule-associated protein involved in membrane protein-cytoskeleton interactions (Score: 0.531).
Figure 7. STRING-Protein–Protein Interactions of TMEM187. (https://string-db.org/network/9606.ENSP00000358999). TMEM187- Interacting ASD Genes: HCFC1: Host cell factor 1: Involved in control of the cell cycle; Coactivator for EGR2 and GABP2 (Score: 0.747); TMLHE: Trimethyllysine dioxygenase, mitochondrial: Converts trimethyllysine (TML) into Hydroxytrimethyllysine (HTML). (Score: 0.657); MECP2: Methyl-CpG-binding protein 2: Chromosomal protein that binds to methylated DNA; Mediates transcriptional repression through interaction with histone deacetylase and the corepressor SIN3A (Score: 0.555); GPHN: Gephyrin: Microtubule-associated protein involved in membrane protein-cytoskeleton interactions (Score: 0.531).
Ijms 20 03358 g007
Table 1. SYTL4 Molecular Pathways and Associated Diseases.
Table 1. SYTL4 Molecular Pathways and Associated Diseases.
Pathway IDPathway DescriptionCount in Gene SetFalse Discovery RateFunctional Description
4721Synaptic vesicle cycle30.0161Communication between neurons is mediated by the release of neurotransmitter from synaptic vesicles (SVs). At the nerve terminal, SVs cycle through repetitive episodes of exocytosis and endocytosis. SVs are filled with neurotransmitters by active transport. DISEASES: Early infantile epileptic encephalopathy; Centronuclear myopathy; Episodic ataxias; Familial or sporadic hemiplegic migraine
4911Insulin secretion30.0161Insulin secretion is regulated by several hormones and neurotransmitters. Peptide hormones, such as glucagon-like peptide 1 (GLP-1), increase cAMP levels and thereby potentiate insulin secretion via the combined action of PKA and Epac2.
Acetylcholine (Ach), a major parasympathetic neurotransmitter. DISEASES: Type II diabetes mellitus; Defects in the degradation of ganglioside.
4152AMPK signaling pathway30.0325AMP-activated protein kinase (AMPK) is a serine threonine kinase that is highly conserved through evolution. AMPK system acts as a sensor of cellular energy status.
4130SNARE interactions 20.0432SNARE proteins (an acronym derived from “SNAP (Soluble NSF Attachment Protein) Receptor”). The primary role of SNARE proteins is to mediate vesicle fusion, that is, the fusion of vesicles with their target membrane-bound compartments. The best studied SNAREs are those that mediate docking of synaptic vesicles with the presynaptic membrane in neurons.
DISEASES: Pseudohypoparathyroidism and Cerebral dysgenesis, neuropathy, ichthyosis, and palmoplantar keratoderma syndrome; CEDNIK syndrome.
Table 2. SYTL4- Network: Biological Processes
Table 2. SYTL4- Network: Biological Processes
Pathway IDPathway DescriptionCount in Gene SetFalse Discovery Rate
GO:0048489synaptic vesicle transport81.14 × 10−9
GO:0097479synaptic vesicle localization81.14 × 10−9
GO:0016079synaptic vesicle exocytosis71.79 × 10−9
GO:0016082synaptic vesicle priming41.14 × 10−7
GO:0031629synaptic vesicle fusion to presynaptic membrane45.45 × 10−7
GO:0007269neurotransmitter secretion63.06 × 10−6
GO:0048167regulation of synaptic plasticity50.000171
GO:0032482RAB protein signal transduction40.000701
GO:0031630regulation of synaptic vesicle fusion to presynaptic membrane20.00156
GO:0014047glutamate secretion30.00165
GO:0007268synaptic transmission60.0115
GO:0023051regulation of signaling10.0211
GO:0050803regulation of synapse structure or activity40.0249
GO:0050804modulation of synaptic transmission40.0356
GO:0007274neuromuscular synaptic transmission20.0385
GO:0065008regulation of biological quality10.0389
GO:0007409axonogenesis50.0393
Table 3. Molecular Function (GO)
Table 3. Molecular Function (GO)
Pathway IDPathway DescriptionCount in Gene SetFalse Discovery RateFunctional Description and Associated Diseases
GO:0019905; GO:0017075Syntaxin binding65.46 × 10−7Syntaxin binding is essential for neurotransmission: syntaxin is a component of the synaptic vesicle fusion machinery.
Mutations in Syntaxin binding protein 1 (STXBP1) have been associated with infantile-epileptic encephalopathy-4 [28].
GO:0005484SNAP receptor activity40.000161SNAPRE activity also regulates neurotransmitter release to ensure vesicle-to-target specificity (SNAP receptors implicated in vesicle targeting and fusion [29].
GO:0019003GDP binding30.019The trimeric-G-protein (GTP binding proteins) play a pivotal role in the signal transduction pathways for numerous hormones and neurotransmitters [30].
Table 4. Showing the Five Autism Predictive Human Serum MicroRNAs with Predicted Interaction with Mouse Sytl4 Gene.
Table 4. Showing the Five Autism Predictive Human Serum MicroRNAs with Predicted Interaction with Mouse Sytl4 Gene.
Sytl4-miR=ASD-miRMouse Sytl4- miRs: Predicted Interactions with ASD & Schizophrenia- Associated miRsValidationReference
miR93Sytl4predictedMGI:1351606c
miR93ASDDysregulated in superior temporal gyrus of ASDStomova, et al., 2015 [37]
miR103-1; miR103-2Sytl4PredictedMGI:1351606c
miR103ASDDysregulated in superior temporal gyrus of ASDStomova, et al., 2015 [37]
miR106bSytl4PredictedMGI:1351606c
miR106b-5p (miR106b) *ASDUpregulated in ASD-serum; differentially expressed in ASD cerebellar cortexVasu, et al., 2014 [34] and Abu-Elneel K, et al., 2008 [35]
miR106bSchizophreniaAltered expression (serum/cortical)Vasu, et al., 2014 [34] and Shi, et al., 2012 [33], Beveridge and Cairns, 2012 [32].
miR130aSytl4PredictedMGI:1351606c
miR130a-3p (miR130a) **ASDGood predictive power for ASD- serumVasu, et al., 2014 [34]
miR-130aSchizophreniaAltered expression (serum/cortical)Vasu, et al., 2014 [34]; Shi, et al., 2012 [33], Beveridge and Cairns, 2012 [32]
miR132Sytl4PredictedMGI:1351606c
miR132ASDDysregulated in superior temporal gyrus of ASDStomova, et al., 2015 [37]
miR181b-1 (miR181b) ***Sytl4PredictedMGI:1351606c
miR181b-2Sytl4PredictedMGI:1351606c
miR181b-5p (miR181b/b1) ***ASDGood predictive power for ASD in serum; differentially expressed in ASD cerebellar cortexVasu, et al. 2014 [34]. Abu-Elneel K, et al. 2008 [35]. Ghahramani Seno, et al., 2011 [36].
miR181bSchizophreniaAltered expression (serum/cortical)Vasu, et al., 2014 [34]; Shi, et al., 2012 [33], Beveridge and Cairns, 2012 [32].
miR320Sytl4PredictedMGI:1351606c
miR320a (miR320) ****ASDGood predictive power for ASD in serumVasu, et al., 2014 [34]
miR320ASDDysregulated in superior temporal gyrus of ASDStomova, et al., 2015 [37]
miR320ASDDifferentially expressed in ASD cerebellar cortexAbu-Elneel K, et al. 2008 [35]
miR328Sytl4PredictedMGI:1351606c
miR328ASDDown regulated in serum; differentially expressed in ASD cerebellar cortexVasu, et al., 2014 [34] and Abu-Elneel K, et al., 2008 [35]
miR328SchizophreniaAltered expression (serum/cortical)Vasu, et al., 2014 [34] and Shi, et al., 2012 [33]; Beveridge and Cairns, 2012 [32]

Share and Cite

MDPI and ACS Style

Rafi, S.K.; Fernández-Jaén, A.; Álvarez, S.; Nadeau, O.W.; Butler, M.G. High Functioning Autism with Missense Mutations in Synaptotagmin-Like Protein 4 (SYTL4) and Transmembrane Protein 187 (TMEM187) Genes: SYTL4- Protein Modeling, Protein-Protein Interaction, Expression Profiling and MicroRNA Studies. Int. J. Mol. Sci. 2019, 20, 3358. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133358

AMA Style

Rafi SK, Fernández-Jaén A, Álvarez S, Nadeau OW, Butler MG. High Functioning Autism with Missense Mutations in Synaptotagmin-Like Protein 4 (SYTL4) and Transmembrane Protein 187 (TMEM187) Genes: SYTL4- Protein Modeling, Protein-Protein Interaction, Expression Profiling and MicroRNA Studies. International Journal of Molecular Sciences. 2019; 20(13):3358. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133358

Chicago/Turabian Style

Rafi, Syed K., Alberto Fernández-Jaén, Sara Álvarez, Owen W. Nadeau, and Merlin G. Butler. 2019. "High Functioning Autism with Missense Mutations in Synaptotagmin-Like Protein 4 (SYTL4) and Transmembrane Protein 187 (TMEM187) Genes: SYTL4- Protein Modeling, Protein-Protein Interaction, Expression Profiling and MicroRNA Studies" International Journal of Molecular Sciences 20, no. 13: 3358. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20133358

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop