Identification of rare DNA sequence variants in high-risk autism families and their prevalence in a large case/control population
- Nori Matsunami†1,
- Charles H Hensel†2Email author,
- Lisa Baird1,
- Jeff Stevens1,
- Brith Otterud1,
- Tami Leppert1,
- Tena Varvil1,
- Dexter Hadley3,
- Joseph T Glessner3,
- Renata Pellegrino3,
- Cecilia Kim3,
- Kelly Thomas3,
- Fengxiang Wang3,
- Frederick G Otieno3,
- Karen Ho2,
- Gerald B Christensen4,
- Dongying Li5,
- Rytis Prekeris5,
- Christophe G Lambert4,
- Hakon Hakonarson3, 6 and
- Mark F Leppert1
© Matsunami et al.; licensee BioMed Central Ltd. 2014
Received: 11 September 2013
Accepted: 24 December 2013
Published: 27 January 2014
Genetics clearly plays a major role in the etiology of autism spectrum disorders (ASDs), but studies to date are only beginning to characterize the causal genetic variants responsible. Until recently, studies using multiple extended multi-generation families to identify ASD risk genes had not been undertaken.
We identified haplotypes shared among individuals with ASDs in large multiplex families, followed by targeted DNA capture and sequencing to identify potential causal variants. We also assayed the prevalence of the identified variants in a large ASD case/control population.
We identified 584 non-conservative missense, nonsense, frameshift and splice site variants that might predispose to autism in our high-risk families. Eleven of these variants were observed to have odds ratios greater than 1.5 in a set of 1,541 unrelated children with autism and 5,785 controls. Three variants, in the RAB11FIP5, ABP1, and JMJD7-PLA2G4B genes, each were observed in a single case and not in any controls. These variants also were not seen in public sequence databases, suggesting that they may be rare causal ASD variants. Twenty-eight additional rare variants were observed only in high-risk ASD families. Collectively, these 39 variants identify 36 genes as ASD risk genes. Segregation of sequence variants and of copy number variants previously detected in these families reveals a complex pattern, with only a RAB11FIP5 variant segregating to all affected individuals in one two-generation pedigree. Some affected individuals were found to have multiple potential risk alleles, including sequence variants and copy number variants (CNVs), suggesting that the high incidence of autism in these families could be best explained by variants at multiple loci.
Our study is the first to use haplotype sharing to identify familial ASD risk loci. In total, we identified 39 variants in 36 genes that may confer a genetic risk of developing autism. The observation of 11 of these variants in unrelated ASD cases further supports their role as ASD risk variants.
KeywordsFamilial autism Haplotype sharing DNA sequence variants Case/control study
Twin and family studies clearly demonstrate that genetics plays a major role in the etiology of ASDs [1–7]. However, family-based approaches using sib-ships or nuclear families have not resulted in the identification of genes or variants that explain a significant portion of the affected population. Similarly, genetic linkage studies have identified a number of chromosomal regions that are thought to contain genes that predispose to ASDs, but identification of the relevant gene(s) in these regions has proven difficult.
More recent analyses have revealed that many of the point mutations thought to predispose children to ASDs were not present in either parent, and thus occurred spontaneously in the affected individuals [8–12]. Such de novo variants likely account for as much as 20 to 30% of the genetic variation which results in ASDs. Additional cases are likely due to recessive inheritance of non-functional or hypo-functional alleles, but autosomal recessive inheritance is thought to explain only about 1% of autism cases [13, 14].
In addition to single nucleotide variants and small insertions/deletions that can be identified by DNA sequencing, larger deletions or duplications (copy number variants, CNVs) have been shown to play a role in the etiology of ASDs [15–27]. Despite the observed inheritance of many ASD predisposition CNVs from an unaffected parent, the lack of extended, multi-generation pedigrees has precluded a comprehensive analysis of segregation of putative ASD predisposition CNVs and single nucleotide polymorphisms (SNPs) and the characterization of other genetic factors necessary for their expression.
The large families available in Utah coupled with the willingness of family members to participate in genetic studies have resulted in the identification of a large number of disease predisposition genes for both Mendelian and complex diseases. The pedigrees used in this study were part of a 70-family linkage study published previously  and two smaller studies that evaluated a single extended pedigree in this collection of families [29, 30]. In this work, we analyzed members of 26 extended multi-generational ASD families and four two-generation multiplex ASD families by performing haplotype sharing analysis to identify chromosomal regions that might harbor ASD predisposition genes. We then used DNA capture and sequencing of all genes in shared regions and of additional candidate autism risk genes to identify SNPs that might predispose to ASD in these families. These SNPs were analyzed in a large case/control study and for segregation in these families. We also evaluated the segregation of CNVs, reported previously , in these families. Consistent with earlier studies, no single locus could account for more than a subset of the affected individuals in any extended pedigree. In particular, multiple potential risk alleles, including in some cases CNVs and SNPs, were identified in an extended pedigree, suggesting that no single variant is the genetic predisposition locus for all affected family members. The data presented here identify several genes that may harbor ASD predisposition mutations, and add to the growing list of genes that are targets for clinical DNA sequencing to aid in the understanding of an ASD genetic diagnosis. These data further suggest that in some individuals multiple genetic variants may be necessary to elicit the observed clinical characteristics and that for a complete understanding of ASD genetics, both sequence variants and CNVs must be analyzed.
A total of 386 DNA samples from 26 extended multi-generation and four two-generation Utah multiplex ASD pedigrees were used in this study. Families were ascertained and recruited using the Utah Population Database (UPDB) as described . Affection status was determined using the Autism Diagnostic Interview-Revised (ADI-R) and the Autism Diagnostic Observation Schedule (ADOS), for both the familial ASD cases and the unrelated ASD cases, as described previously . The average number of affected individuals in each pedigree is 7.9. The pedigrees described here are a subset of those described previously . Pedigree details are shown in Additional file 1: Table S1. The 55 samples used for our previous CNV discovery were included in these families. A total of 9,000 DNA samples previously described in a case/control study , including 3,000 individuals with ASD and 6,000 controls, were used to evaluate these variants in a broader population. All samples collected for the work described here were collected under methods approved by the University of Utah Institutional Review Board (IRB) (University of Utah IRB#: 6042–96) or The Children’s Hospital of Philadelphia IRB (CHOP IRB#: IRB 06–004886). Patients and their families were recruited through the University of Utah Department of Psychiatry or The Children’s Hospital of Philadelphia clinic or CHOP outreach clinics. Written informed consent was obtained from the participants or their parents using IRB approved consent forms prior to enrollment in the project. There was no discrimination against individuals or families who chose not to participate in the study. All data were analyzed anonymously and all clinical investigations were conducted according to the principles expressed in the Declaration of Helsinki.
SNP microarray genotyping
Affymetrix 250 K NspI SNP chip genotyping was carried out on all 386 DNA samples using the manufacturer’s recommended procedure. Genotypes were called by Affymetrix Genotyping Console software using the BRLMM  genotype calling algorithm. Only SNPs with call rates greater than or equal to 99% were used for further analyses. SNPs demonstrating Mendelian errors also were identified using PedCheck  and were excluded.
Shared haplotype analysis
We performed shared haplotype analysis on each pedigree to identify genomic regions that have significant sharing among the affected individuals in that pedigree. The HapShare algorithm  was used to perform haplotype phasing based on Mendelian inheritance and to identify shared genomic segments. The comparisons included N out of N affected individuals, (N-1) out of N, (N-2) out of N, (N-3) out of N, and so on (See Figure one in reference ). In two-generation pedigrees, in some cases co-segregation of haplotypes was observed in all affected individuals analyzed, but the shared regions were large, including up to half of a chromosome. Consequently, shared regions from nuclear families were not selected for sequencing unless they overlapped regions observed in additional families.
Custom targeted exome DNA sequencing
Sequence alignment and variant detection methods
Alignment and assembly
Sequence variant detection
CLC Bio Genomics Workbench (CLC Bio Inc.)
CLC Bio Genomics Workbench
In silico functional analysis was carried out initially using cSNP classifier, a preliminary program later incorporated into VAAST , to classify variants as synonymous, conservative missense, non-conservative missense, nonsense, frameshift, or splice site changes. Later, variants were re-annotated using the ANNOVAR program . The KnownGene and RefSeq gene tracks from the UCSC genome browser were used to annotate functional variants, and the LiftOver tool was used to convert human genome build 36 (GRCh36/hg18) coordinates to human genome build 37 (GRCh37/hg19) coordinates [39, 40].
Custom microarray design and array processing
Design of the custom iSelect InfiniumTM II BeadChip array (Illumina Inc., San Diego, CA, USA) including probes for 2,799 putative functional candidate SNPs and 7,134 CNV probes was described previously . The custom iSelect array was processed on 3,000 case and 6,000 control samples at the Center for Applied Genomics at The Children’s Hospital of Philadelphia (CHOP) . The same array was also used to analyze DNA from 196 Utah discovery cohort family members at the University of Utah Genomics Core facility for variant validation and analysis of SNP segregation in families.
Array data quality control
Subjects were withheld from SNP analysis if any of the following were true: 1) subsequent to genotyping, the DNA sample was of apparent poor quality, evidenced by very low call rates (N = 134); 2) the subject was identified as a trisomy-21 (N = 51); 3) the subject was outside of the central cluster of Caucasian subjects identified by principal component analysis (PCA) (N = 903) .
Relatedness estimation further indicated that some of the case subjects and controls were part of families with multiple relatives represented in the data. Re-evaluation of family structure in the sample cohorts used subsequently identified additional relationships. Subsequent association tests were therefore conducted using only one member of each known family in order to reduce the possibility of statistical confounding due to relatedness. For these tests, the subject selected from each family was the individual located nearest to the median centroid of the first two principal components. The number of subjects removed due to relatedness was 688. This resulted in a final sample set for association testing comprising 7,326 subjects, of which 1,541 were cases and 5,785 were controls.
was determined to be outside of the main cluster, and therefore non-Caucasian. Six hundred and eighty-two subjects were placed in the non-Caucasian category. A graphical representation of the results of this PCA analysis were reported previously .
Prior to association testing, SNPs were evaluated for call rate, Hardy-Weinberg equilibrium (HWE) and allele frequency. All SNPs with call rates lower than 99% were removed from further analysis. No SNPs had significant Hardy-Weinberg disequilibrium.
Laboratory confirmation of SNPs and CNVs
For molecular validation of SNPs, PCR products were first screened by LightScanner High Resolution Melt curve analysis (BioFire Diagnostics Inc., Salt Lake City, UT, USA) for the presence of sequence variants. PCR primer sequences are shown in Additional file 1: Table S3. Any samples that gave abnormal melt profiles were sequenced using the Sanger method at the University Utah Sequencing Core facility to confirm the presence of a sequence variant. For CNVs, pre- or custom-designed TaqMan copy number assays (Applied Biosystems Inc., Foster City, CA, USA) were used as described previously .
Protein binding assays
All glutathione S-transferase (GST)-tagged proteins were expressed and purified as described previously . To test Rab11FIP5 binding to various Rab GTPases, purified recombinant FIP5(490–653) or FIP5(490–653)-P652L were incubated with glutathione beads coated with GST, GST-Rab11a, GST-Rab4a or GST-Rab3a in the presence of 1 μm GMP-PNP. Beads were then washed with PBS and eluted with 1% sodium dodecyl sulphate (SDS). Eluates were then analyzed for the presence of FIP5 (490–653) by immunoblotting with anti-Rab11FIP5 antibodies. A similar assay also was used to test the ability of Rab11FIP5 (wild-type or P652L mutant) to dimerize.
Flow cytometry analysis of transferrin recycling
To test the effect of the Rab11FIP5-P652L mutant on endocytic recycling, the transferrin recycling assay was used as described previously . Briefly, HeLa cells expressing either wild-type FIP5-GFP or FIP5-GFP-P652L were incubated with transferrin conjugated to Alexa488. Cells were then washed and incubated with serum-supplemented media for varying amounts of time. The cell-associated (not recycled) Tf-Alexa488 was analyzed by flow cytometry.
Affymetrix 250 K SNP genotyping and haplotype sharing
SNP genotyping was carried out on 386 DNA samples from 26 extended multi-generation and four two-generation Utah multiplex ASD pedigrees. SNPs with no map location were not included in the analysis. The average call rate was 99.1% for the entire dataset.
Chromosomal regions selected for sequencing based on haplotype sharing
Eighteen shared haplotype regions
Affecteds sharing haplotype
6 of 6
6 of 6
5 of 7
4 of 7
4 of 7, 4 of 4
6 of 6
4 of 4, 4 of 6
4 of 7
5 of 8a
4 of 6
5 of 7
5 of 8
3 of 4, 6 of 6
4 of 4
3 of 4, 4 of 6, 5 of 8
4 of 7, 5 of 6, 3 of 4
4 of 7
5 of 6, 5 of 6
Sequence capture, sequence analysis and variant identification
We performed capture and DNA sequencing using DNA from 26 affected individuals from 11 families that showed the best sharing of genomic segments. These samples included individuals from two-generation pedigrees that had shared haplotypes overlapping regions identified in the extended pedigrees. Eight to nine million 36 base short reads were obtained from each sample. The short reads alignment against the NCBI reference human genome build 36 revealed coverage of 86 to 97% of the designed capture area, with the average read depth over the designed capture area of 30 to 47X. Because 1) the capture library was constructed in a directional manner, 2) all capture probes represented the same DNA strand, and 3) the library was sequenced only from one direction, some portions of exons lacked the sequence depth of coverage that we desired. Consequently, there could be additional variants (false negative results) that we did not detect in some of the genes. For example no variants were identified on haplotypes that segregate to all affected individuals in pedigree 10 on chromosomes 2 and 14 (Additional file 2: Figure S1, A and B, Additional file 3: Figure S2). Nonetheless, variant calling using the three methods shown in Table 1 identified over one million sequence variants called by at least two of the three methods. Analysis using cSNP classifier resulted in the detection of 2,825 putative functional candidate SNPs, including 210 nonsense variants, 1,614 non-conservative missense variants, 35 frameshift variants and 966 splice site variants.
We chose to design a custom microarray to evaluate the variants that we identified by sequencing in order to 1) interrogate the entire set of candidate functional SNPs in the discovery families for validation, and 2) to perform a large scale case/control study to determine if any of the variants identified predisposition genes important to the broad population of children with ASD (Figure 1). Following array design and manufacture, probes for 2,413 variants were created successfully. Custom microarray experiments on Utah discovery and CHOP case/control samples revealed 584 out of 2,413 putative variants to be polymorphic. The complete list of polymorphic variants is shown in Additional file 1: Table S4. The remaining array probes did not detect a non-reference sequence allele. These 1,829 variants thus were interpreted to be false positives due to the variant calling and alignment process of single end sequence data.
Sequence variants identified in families and observed in the case/control study
Fisher’s exact P
Odds ratio (Minor allele)
Odds ratio 95% lower confidence bound
Odds ratio 95% upper confidence bound
R288W, splice site
Sequence variants observed only in high-risk ASD families
Tested affecteds in pedigree
Affecteds with variant
Segregation of variants in high-risk pedigrees
Pedigree 1 (Figure 2) shows a two-generation family co-segregating a missense variant in RAB11FIP5 (Table 4). This variant is present in the mother and segregates to all three male affected children in the family, and not to the unaffected female child. RAB11FIP5 has previously been implicated as an ASD risk gene based on its disruption by a translocation observed in a 10 year-old male child with a diagnosis of pervasive developmental disorder not otherwise specified (PDD-NOS) . The variant detected in pedigree 1 results in a P652L substitution. Proline is conserved at this residue in all of the mammalian RAB11FIP5 genes sequenced to date, suggesting that it is important for protein function. A second individual, with a P652H variant, was detected in our case/control study (Table 3) using our custom microarray. Neither the P652L substitution nor the P652H substitution was observed in the ESP6500, 1,000 genomes project or dbSNP137 databases (Additional file 1: Table S5). Each of these variants was confirmed by Sanger sequencing (See Additional file 5: Figure S4, Additional file 6: Figure S5 for chromatograms). An additional affected individual of non-European descent, and thus not included in the case/control study, also carried the P652H variant (data not shown). The presence of the P652H variant in an additional individual with autism and not in any controls further supports the likelihood of variants in RAB11FIP5 contributing to autism risk.
Pedigree 2 (Figure 3) is a two-generation family with six affected individuals from two fathers. In this pedigree, five of the six affected individuals inherit a variant resulting in an I26T substitution in C14orf2. Two additional sequence variants, one each in the PDK4 and SDR39U1 genes, segregate to three and two affected individuals respectively. In addition, a CNV gain (OR = 3.37) described previously  is present in one affected individual. The C14orf2 and PDK4 variants were maternally inherited, while the C7orf10 and the CNV were either of paternal origin or occurred as de novo variants. Of the variants detected in this family, only the C7orf10 variant was observed in our case/control study. However, this variant had an odds ratio of 1.62 (95% confidence interval 1.04 to 2.53), suggesting the possibility for a role in autism predisposition in the general population.
Pedigree 3 (Figure 4) also is a two-generation family, with five male children affected with autism. In this pedigree, four of the five affected individuals exhibit maternal inheritance of an F154L variant in the KLHL6 gene. This A/G nucleotide variant also is found at the first nucleotide of an exon and thus also may affect splicing of the KLHL6 primary transcript. In addition to this variant, three of the five offspring have a paternally inherited D303H missense variant in the SPATA5L1 gene while two of five also have a maternally inherited P238L change in the ITPK1 gene. One affected child does not inherit any of these variants. Of interest, none of the variants observed in this small family were observed in any cases or controls in our population study, demonstrating that they are not common autism predisposition loci.
Pedigree 4 (Figure 5) is a six generation family with an ancestor common to all seven male children that are affected with autism. These children all are in the fifth or sixth generations of the pedigree. Linkage analysis was performed previously on this family using Affymetrix 10 K SNP genotype data [29, 30], and three regions of significant linkage were identified. These include 3q13.2-q13.31, 3q26.31-q27.3, and 20q11.21-q13.12. These three regions also were identified by haplotype sharing in this study (Figure 5, see Additional file 2: Figure S1C for chromosome 20 haplotype sharing). Four of the seven affected individuals in this family share a P49L variant that is the result of an A/G transition in the DEFB124 gene on chromosome 20q11.21, consistent with the haplotype sharing that we observed (Additional file 2: Figure S1c) and with the published linkage result. This variant was not observed in cases or controls in our population study. One affected individual in this pedigree does not share the DEFB124 variant, but instead has a chromosome 3q gain CNV, inherited from his father, that had an odds ratio of 3.74 in our previous study . The elevated odds ratio suggests that this CNV is an autism risk locus.
Two additional affected individuals in pedigree 4 do not carry any variant that we detected in our families. However, as indicated in Figure 5, each of these two individuals is descended from a marry-in spouse with a strong family history of autism, suggesting the possibility of additional undetected variants.
Finally, one affected individual who carries the DEFB124 variant carries variants in the HEPACAM2 gene (odds ratio 1.83 in our population study, Table 3), the AP1G2 gene (odds ratio 1.67, Table 3), the PYGO1 gene and the RELN gene. Neither the RELN variant nor the PYGO1 variant was observed in the case/control study (Table 4). Interestingly, a recent unpublished clinical case identified a heterozygous deletion involving 15 exons of the RELN gene in a patient diagnosed with autism and behavior/conduct disorder (Rena Vanzo, personal communication). Homozygous or compound heterozygous mutations in RELN are associated with lissencephaly [44, 45], but this RELN deletion is the first description of an individual with a developmental phenotype that may be due to haploinsufficiency at this locus.
Pedigree 5 (Figure 6) is a four-generation family with nine individuals affected with autism (seven male, two female). Two variants are of particular interest in this family. The first is a CNV including the 5′-flanking region of the NRXN1 gene. This CNV is inherited from a father who marries into the family in the second generation. This CNV segregates to three of the four descendants of this individual who are diagnosed with autism. An overlapping NRXN1 CNV was shown in our previous work to have an odds ratio of 14.96 , consistent with other publications suggesting a role for NRXN1 variants in autism, as well as other neurological disorders [46–48]. However, that CNV was shown to extend into the coding region of NRXN1, while TaqMan CNV analysis demonstrates that the CNV in pedigree 5 did not (data not shown). Thus, the significance of the NRXN1 CNV observed in this family is uncertain.
A second variant identified in this family, found on a haplotype shared by all five affected individuals in two branches of the family (Additional file 2: Figure S1c), is a C/T transition in the AKAP9 gene that results in an R3233C missense substitution. Note that none of the individuals in these two branches of the family carry the NRXN1 CNV. The AKAP9 variant was observed in 4/1541 cases and 4/5785 controls in our population study (odds ratio of 3.76, 95% confidence interval 0.94 to 15.03) (Table 3). A second missense variant in the AKAP9 gene was observed in a single affected individual in a nuclear family (pedigree 6, Additional file 7: Figure S6). This second AKAP9 variant was not observed in the case/control study (Table 4). The AKAP family of proteins has been suggested to connect different biological pathways that are involved in nervous system development .
Pedigree 5 also segregates other variants that are inherited by multiple children affected with autism. One branch of the pedigree segregates a G/C transversion in the CLMN gene that results in a P158A missense substitution. This variant yielded an odds ratio of 1.67 (95% confidence interval 0.73 to 3.84) in our case/control study, suggesting that it is an ASD risk allele. A variant in the ABP1 gene, also the result of a G/C transversion and resulting in an R345P missense substitution, was observed in two affected individuals in a single branch of the family. This variant was maternally inherited and not seen elsewhere in the pedigree. However, this variant was observed in 1/1,541 cases and 0/5,785 controls in our population study (Table 3) and was not observed in the ESP6500, 1,000 Genomes, or dbSNP137 databases (Additional file 1: Table S5), indicating that it may be a very rare ASD risk variant. Finally, a G/T transversion in the ALX1 gene that results in an R64L missense substitution was paternally inherited by a single individual. This variant also was seen in pedigree 7 (Additional file 8: Figure S7) and was observed multiple times in our population study (27/1,541 cases and 58/5,785 controls) yielding an odds ratio of 1.75 (95% confidence interval 1.11 to 2.77) (Table 3). Expression of this gene also may be increased by a downstream balanced translocation in a family with mental retardation, language delay and microcephaly that segregate with the translocation .
Pedigrees 8 to 10 are shown in Additional file 9: Figure S8, Additional file 10: Figure S9, Additional file 3: Figure S2. One of these pedigrees, pedigree 10, carried two haplotypes (chromosomes 2 and 14) segregating to all six affected individuals (Additional file 2: Figure S1a-b). Sequencing of the genes encompassed by these regions did not identify potential causal variants. This could be due to poor sequence coverage of some portions of the genes. However, sequencing of affected individuals in these families did result in the identification of variants that could be autism risk alleles. One of these variants, a G/A transition that result in a Q22X change in the MOK gene observed in a single affected individual and inherited from her father, is interesting, as it was observed in our population study and yielded an odds ratio of 3.76 (95% confidence interval 0.53 to 26.67) (Table 3). Other variants in pedigrees 8 to 10 (Additional file 9: Figure S8, Additional file 10: Figure S9, Additional file 3: Figure S2), including some only seen in Utah families and others seen in both families and in our population study also were identified. These variants are included in Table 3 and Table 4.
Functional analysis of RAB11FIP5
To uncover the functional consequences of the Rab11FIP5-P652L variant we evaluated binding of Rab11FIP5 to Rab11. Rab11 is a small monomeric GTPase that mediates Rab11FIP5 recruitment to endocytic membranes and is required for Rab11FIP5 function . As shown in Additional file 11: Figure S10A, the P652L substitution did not affect Rab11FIP5 binding to Rab11, nor did it affect its specificity toward the Rab11 GTPase. It was previously shown that Rab11FIP5 forms homodimers and that its ability to dimerize is also required for Rab11FIP5 cellular functions . Thus, we tested the effect of P652L substitution on Rab11FIP5 ability to dimerize. As shown in Additional file 11: Figure S10B, the Rab11FIP5-P652L mutant was still able to form dimers. Consistent with in vitro binding data, FIP5-GFP-P652L endocytic localization in HeLa cells was also not affected (Additional file 11: Figure S10).
It is now well established that Rab11FIP5 functions by regulating endocytic recycling . To that end, we next tested whether Rab11FIP5-P652L has an effect on recycling of transferrin receptors in HeLa cells and found that the P652L substitution did not alter recycling (Additional file 11: Figure S10H). Thus, so far we were unable to find any functional consequences of Rab11FIP5-P652L substitution, suggesting that core Rab11FIP5 properties are not affected. Further studies will be needed to uncover the kinetic or binding defects of Rab11FIP5 variants, especially within the context of neuronal cells.
Numerous studies have implicated genetics in the etiology of ASDs, and recent results have implicated hundreds of genes and chromosomal regions as candidate ASD predisposition loci. These include chromosomal rearrangements as well as sequence variants. None of these candidates accounts for more than a small portion of ASD cases, and it is clear, based on the apparent complexity of autism genetics, that many genes are still undiscovered. In multiplex ASD pedigrees, it has been suggested that inherited genetic variants are more likely to be responsible, while in simplex families de novo variants may be more important [8, 16].
We used a discovery/validation strategy based on identifying inherited genetic variants in two- to nine-generation ASD families, followed by a case/control analysis of those variants in DNA samples from unrelated children with autism and children with normal development, to identify putative familial ASD predisposition genes. Using haplotype analysis we identified shared genomic segments within the families, and used DNA sequencing and CNV analysis to identify potential causal mutations on those haplotypes. We also followed this with a large case/control study to determine if any of the variants we identified might play a role in the general population of individuals with ASD.
We showed previously that identification of CNVs in a family-based discovery cohort could identify copy number variants relevant to the general ASD population . In this paper, we applied the same strategy to the family-based identification of DNA sequence variants, and also follow up on our CNV analysis by providing information regarding segregation of some of these CNVs in our high-risk ASD families.
We identified 39 SNPs that are likely to affect protein function that have segregation patterns and ASD case allele frequencies suggestive of a role in ASD predisposition. Thirty-one of these variants result in non-conservative amino acid substitutions, five are predicted to affect splicing (three of these are predicted to affect both splicing and protein coding), and three introduce premature termination codons. Two variants were identified in the AKAP9 gene and the JMJD7 gene (or the JMJD7-PLA2G4B fusion gene), and two different variants were identified that affect the same amino acid residue in the RAB11FIP5 gene; so collectively, these SNPs identify 36 potential ASD risk genes.
With the exception of two-generation families, and consistent with our haplotype sharing results, no sequence variants or CNVs implicated as ASD predisposition loci segregate to all affected individuals in a pedigree. This is consistent with previous genetic studies, which to date have been unable to demonstrate segregation of a single ASD risk locus in an extended family (for example see ). In our pedigree 5 (Figure 6), two independent risk variants, a single nucleotide variant in AKAP9 and a deletion CNV in or near NRXN1, segregate to different branches of the family. Other risk variants also are found in individuals with ASD in this family, including two sequence variants with odds ratios greater than 1.5 in our population study. These results suggest that even in extended families that might be predicted to be segregating a single risk allele with reduced penetrance, multiple risk alleles in different ASD predisposition loci may be necessary. The results further suggest that use of specific inheritance models when evaluating autism genetics in large families should be approached with caution. Further, given our sequencing results and CNVs, the data suggest that a two-pronged approach, involving both DNA sequencing and microarray-based CNV analysis, may be necessary for a complete genetic diagnosis of children with ASDs.
Eleven of the putative autism risk variants that we identified in our high-risk families are further supported by data from our case/control study. Three of these variants each were seen in a single ASD case (out of 1,541 total cases) and in none of 5,785 controls. Familial variants that we detected in eight additional genes are more common in ASD cases than in controls, and each has an odds ratio greater than 1.5. Although these variants are rare (all have frequencies of < 0.01 in our case/control study), their identification in affected individuals in our ASD families and their increased prevalence in unrelated affected individuals support their role as potential ASD risk loci.
Biological functions/pathways of genes with variants found in children with ASDs
Previously associated with autism
TRPM1, RAB11FIP5, AKAP9, SCN3A
Previously associated with neurological disorder (other than autism)
RELN (autosomal recessive lissencephaly), ALX1 (facial clefting, micropthalmia), CCDC85C (seizures), EPB41L1 (intellectual disability)
ITPK1, CLMN, PRTG
PLA2G4B, c7orf10, PDK4, C14orf2
Inflammatory response/Immune function
DEFB124, BPI, RNF31, IRF9
The incidence/prevalence of autism has been reported to be as high as 1 in 38 children , while the estimated prevalence in the US is reported to be approximately 1 in 88 . Multiple factors, including genetic, epigenetic and environmental, are thought to be necessary for autism development. As with previous work involving familial autism, we found that no single gene variant described here can account for every case of autism in our pedigrees. This may reasonably be expected of a multigenic disorder with high prevalence in the population. It also is consistent with the high likelihood that autism represents a number of different conditions, at least at the genetic level, all under the same phenotypic umbrella called autism. Several variants that we identified, however, represent excellent candidates for significant risk factors for autism based on their known functions and the currently understood molecular pathways involved in autism. We discuss the most compelling cases below.
RAB11FIP5 is a member of a family of scaffolding proteins for the RAS GTPase, Rab11. Specifically, RAB11FIP5 has been characterized as a key player in apical endosome recycling, plasma membrane recycling and transcytosis [70, 71]. We identified a P652L variant in three affected siblings in a family of six members, in which the mother is an unaffected P652L carrier. An additional variant resulting in a P652H substitution also was detected in 1/1,541 Caucasian ASD cases and 0/5,785 Caucasian children with normal development (Table 3). These variants modify a conserved proline within the C-terminus of RAB11FIP5.
Heterozygous disruption of RAB11FIP5 was observed previously in a ten year old boy with a balanced translocation (46, XY, t(2;9)(p13;p24)) that disrupts only the RAB11FIP5 gene . This individual has a clinical diagnosis of PDD-NOS, an autism spectrum disorder. This translocation led the authors to suggest that haploinsufficiency of RAB11FIP5 contributes to the subject’s ASD. RAB11FIP5 works closely in conjunction with RAB11, and its presence has been detected in both presynaptic and post-synaptic densities where Rab11 plays a key role in determining synaptic strength in long-term depression , regulates norepinephrine transporter trafficking , carries out synaptic glutamate receptor recycling , and regulates dendritic branching in response to BDNF [75, 76]. All of these functions have been suggested to be significant contributors to the etiology of ASDs [77, 78] and further support the role of mutations in RAB11FIP5 as ASD risk alleles.
AKAP9 is a member of a family of over 50 proteins that serve as scaffolding partners for PKA, its effectors, and phosphorylation targets. AKAP9, also known as Yotiao, is chiefly expressed in the heart and brain, where the encoded protein serves as a scaffold for PKA, protein phosphatase I, NMDA receptors, the heart potassium channel subunit KCNQ1, IP3R1, and specific isoforms of adenylyl cyclase [79–83]. The subcellular localization and assembly of these multimeric protein scaffolds, mediated by AKAPs, are thought to be essential for function, since disruption of the interaction between the AKAP and its effectors leads to a loss of activity. In the case of KCNQ1, loss of interaction between AKAP9 and KCNQ1 leads to a potentially fatal heart condition, long-QT syndrome, which also arises in cases with loss of function mutations in KCNQ1 itself .
We identified two variants in the AKAP9 gene. These variants result in R3233C and R3832C substitutions in the encoded protein. These two variants were coincident with autism and were found in two unrelated extended ASD pedigrees (Figure 6, Additional file 7: Figure S6). The R3233C variant was additionally found in our case/control study. A recent meta-study of the genes identified from the five major autism GWAS studies and autism candidate genes arising from alternative methodologies, such as large scale CNV studies, placed AKAPS as a central, integral gene family linking many of the pathways identified by bioinformatics . Given its role in localizing PKA, adenylyl cyclase isoforms and NMDAR in the postsynaptic scaffold, AKAP9 represents a protein that, like its better-characterized counterpart AKAP5, could function in synaptic transmission and plasticity, glutamatergic receptor function regulation and recycling, and dendritic spine morphology .
It is notable that two of the genes (MOK, TRPM1) containing potential ASD risk alleles were partially or completely encompassed by risk CNVs observed in our previous study . This suggests that the same genes may be affected by different genetic mechanisms with the same or similar phenotypic result. The CNVs containing these genes were both copy number losses. The MOK sequence variant described here was a nonsense change, while the TRPM1 variant was a missense change. These results are consistent with the MOK and TRPM1 effects being due to haploinsufficiency at these two loci.
Our data demonstrate the complex nature of autism genetics. Although the heritability for autism is quite high, our data show that numerous genetic variants may confer risk to ASD even in a single family. This finding is consistent with the results of a whole genome sequencing study that used both a recessive model and model independent analyses to identify several potential ASD risk variants in an ASD family with two affected individuals . Consistent with the large number of potential ASD risk genes identified to date, none of the genes identified in this single multiplex ASD family overlapped with the genes identified in our study. Our study adds to this complexity by identifying sequence variants in regions of haplotype sharing in 30 high-risk ASD families of two to six generations. Our data further demonstrate that in very large multi-generation families, the likelihood of additional risk variants entering the family from individuals who marry into the pedigree is high. These results are consistent with the prevalence of autism in the population and the high heritability for autism seen in numerous studies. Future studies will be necessary to characterize the biological roles that the variants we identified play in the complex set of disorders that comprise ASD.
Our study is the first to use an empirical approach to identify shared genomic segments, followed by sequence variant detection, to identify potential ASD risk variants in a large set of autism families. We identified 39 DNA sequence variants in 36 genes that may represent ASD risk genes. Eleven of these variants may also be risk variants in unrelated ASD cases, as they were shown to have odds ratios greater than 1.5 in a large ASD case/control population. Twenty-eight of these variants were observed only in high-risk ASD families in our study. These very rare sequence changes also may identify additional ASD risk genes. Segregation patterns demonstrate that in multi-generation ASD families, multiple genetic variants, including CNV and SNP risk variants, may contribute to ASD risk.
Autism spectrum disorder
Copy number variant
Single nucleotide polymorphism
National Center for Biotechnology Information
Utah Population Database
Autism Diagnostic Interview-Revised
Autism Diagnostic Observation Schedule
The Children’s Hospital of Philadelphia
Principal component analysis
Minor allele frequency
Sodium dodecyl sulfate
Pervasive developmental disorder not otherwise specified.
All family subjects were ascertained and DNA collected with support from R01 MH 06359 from the National Institute of Mental Health and U19HD035476 from the National Institute of Child Health and Human Development. Utah DNA samples were processed with support from GCRC M01-RR025764 from the National Center for Research Resources. The authors gratefully acknowledge the resources provided by the Autism Genetic Resource Exchange (AGRE) Consortium and the participating AGRE families.
- Rosenberg RE, Law JK, Yenokyan G, McGready J, Kaufmann WE, Law PA: Characteristics and concordance of autism spectrum disorders among 277 twin pairs. Arch Pediatr Adolesc Med. 2009, 163: 907-914. 10.1001/archpediatrics.2009.98.View ArticlePubMedGoogle Scholar
- Hallmayer J, Cleveland S, Torres A, Phillips J, Cohen B, Torigoe T, Miller J, Fedele A, Collins J, Smith K, Lotspeich L, Croen LA, Ozonoff S, Lajonchere C, Grether JK, Risch N: Genetic heritability and shared environmental factors among twin pairs with autism. Arch Gen Psychiatry. 2011, 68: 1095-1102. 10.1001/archgenpsychiatry.2011.76.PubMed CentralView ArticlePubMedGoogle Scholar
- Lichtenstein P, Carlström E, Råstam M, Gillberg C, Anckarsäter H: The genetics of autism spectrum disorders and related neuropsychiatric disorders in childhood. Am J Psychiatry. 2010, 167: 1357-1363. 10.1176/appi.ajp.2010.10020223.View ArticlePubMedGoogle Scholar
- Ronald A, Hoekstra RA: Autism spectrum disorders and autistic traits: a decade of new twin studies. Am J Med Genet B Neuropsychiatr Genet. 2011, 156B: 255-274.View ArticlePubMedGoogle Scholar
- International Molecular Genetic Study of Autism Consortium (IMGSAC): A full genome screen for autism with evidence for linkage to a region on chromosome 7q. Hum Mol Genet. 1998, 7: 571-578.View ArticleGoogle Scholar
- International Molecular Genetic Study of Autism Consortium (IMGSAC): A genomewide screen for autism: strong evidence for linkage to chromosomes 2q, 7q, and 16p. Am J Hum Genet. 2001, 69: 570-581.View ArticleGoogle Scholar
- Buxbaum JD, Silverman J, Keddache M, Smith CJ, Hollander E, Ramoz N, Reichert JG: Linkage analysis for autism in a subset families with obsessive-compulsive behaviors: evidence for an autism susceptibility gene on chromosome 1 and further support for susceptibility genes on chromosome 6 and 19. Mol Psychiatry. 2004, 9: 144-150. 10.1038/sj.mp.4001465.View ArticlePubMedGoogle Scholar
- Iosifov I, Ronemus M, Levy D, Wang Z, Hakker I, Rosenbaum J, Yamrom B: Lee Y-h, Narzisi G, Leotta A, Kendall J, Grabowska E, Ma B, Marks S, Rodgers L, Stepansky A, Troge J, Andrews P, Bekritsky M, Pradhan K, Ghiban E, Kramer M, Parla J, Demeter R, Fulton LL, Fulton RS, Magrini VJ, Ye K, Darnell J, Darnell RB, Mardis ER, Wilson RK, Schatz MC, McCombie WR, Wigler M: De novo gene disruptions in children on the autistic spectrum. Neuron. 2012, 74 (2): 285-299. 10.1016/j.neuron.2012.04.009.View ArticleGoogle Scholar
- Sanders SJ, Murtha MT, Gupta AR, Murdoch JD, Raubeson MJ, Willsey AJ, Ercan-Sencicek AG, DiLullo NM, Parikshak NN, Stein JL, Walker MF, Ober GT, Teran NA, Song Y, El-Fishawy P, Murtha RC, Choi M, Overton JD, Bjornson RD, Carriero NJ, Meyer KA, Bilguvar K, Mane SM, Sestan N, Lifton RP, Günel M, Roeder K, Geschwind DH, Devlin B, State MW: Disruptive de novo point mutations, revealed by whole-exome sequencing, are strongly associated with autism spectrum disorders. Nature. 2012, 485 (7397): 237-241. 10.1038/nature10945.PubMed CentralView ArticlePubMedGoogle Scholar
- Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, Lin CF, Stevens C, Wang LS, Makarov V, Polak P, Yoon S, Maguire J, Crawford EL, Campbell NG, Geller ET, Valladares O, Schafer C, Liu H, Zhao T, Cai G, Lihm J, Dannenfelser R, Jabado O, Peralta Z, Nagaswamy U, Muzny D, Reid JG, Newsham I, Wu Y: Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012, 485 (7397): 242-245. 10.1038/nature11011.PubMed CentralView ArticlePubMedGoogle Scholar
- O’Roak BJ, Deriziotis P, Lee C, Vives L, Schwartz JJ, Girirajan S, Karakoc E, Mackenzie AP, Ng SB, Baker C, Rieder MJ, Nickerson DA, Bernier R, Fisher SE, Shendure J, Eichler EE: Exome sequencing in sporadic autism reveals a highly interconnected protein network and extreme locus heterogeneity. Nature. 2012, 485 (7397): 246-250. 10.1038/nature10989.PubMed CentralView ArticlePubMedGoogle Scholar
- O’Roak BJ, Vives L, Fu W, Egertson JD, Stanaway IB, Phelps IG, Carvill G, Kumar A, Lee C, Ankenman K, Munson J, Hiatt JB, Turner EH, Levy R, O’Day DR, Krumm N, Coe BP, Martin BK, Borenstein E, Nickerson DA, Mefford HC, Doherty D, Akey JM, Bernier R, Eichler EE, Shendure J: Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science. 2012, 338 (6114): 1619-1622. 10.1126/science.1227764.PubMed CentralView ArticlePubMedGoogle Scholar
- Lim ET, Raychaudhuri S, Sanders SJ, Stevens C, Sabo A, MacArthur DG, Neale BM, Kirby A, Ruderfer DM, Fromer M, Lek M, Liu L, Flannick J, Ripke S, Nagaswamy U, Muzny D, Reid JG, Hawes A, Newsham I, Wu Y, Lewis L, Dinh H, Gross S, Wang LS, Lin CF, Valladares O, Gabriel SB, De Pristo M, Altshuler DM, Purcell SM: Rare complete knockouts in humans: population distribution and significant role in autism spectrum disorders. Neuron. 2013, 77 (2): 235-242. 10.1016/j.neuron.2012.12.029.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu TW, Chahrour MH, Coulter ME, Jiralerspong S, Okamura-Ikeda K, Ataman B, Schmitz-Abe K, Harmin DA, Adli M, Malik AN, D’Gama AM, Lim ET, Sanders SJ, Mochida GH, Partlow JN, Sunu CM, Felie JM, Rodriguez J, Nasir RH, Ware J, Joseph RM, Hill RS, Kwan BY, Al-Saffar M, Mukaddes NM, Hashmi A, Balkhy S, Gascon GG, Hisama FM, LeClair E: Using whole-exome sequencing to identify inherited causes of autism. Neuron. 2013, 77 (2): 259-273. 10.1016/j.neuron.2012.11.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Girirajan S, Brkanac Z, Coe BP, Baker C, Vives L, Vu TH, Shafer N, Bernier R, Ferrero GB, Silengo M, Warren ST, Moreno CS, Fichera M, Romano C, Raskind WH, Eichler EE: Relative burden of large CNVs on a range of neurodevelopmental phenotypes. PLoS Genet. 2011, 7: e1002334-10.1371/journal.pgen.1002334.PubMed CentralView ArticlePubMedGoogle Scholar
- Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, Walsh T, Yamrom B, Yoon S, Krasnitz A, Kendall J, Leotta A, Pai D, Zhang R, Lee YH, Hicks J, Spence SJ, Lee AT, Puura K, Lehtimäki T, Ledbetter D, Gregersen PK, Bregman J, Sutcliffe JS, Jobanputra V, Chung W, Warburton D, King MC, Skuse D, Geschwind DH, Gilliam TC: Strong association of de novo copy number mutations with autism. Science. 2007, 316: 445-449. 10.1126/science.1138659.PubMed CentralView ArticlePubMedGoogle Scholar
- Marshall CR, Noor A, Vincent JB, Lionel AC, Feuk L, Skaug J, Shago M, Moessner R, Pinto D, Ren Y, Thiruvahindrapduram B, Fiebig A, Schreiber S, Friedman J, Ketelaars CE, Vos YJ, Ficicioglu C, Kirkpatrick S, Nicolson R, Sloman L, Summers A, Gibbons CA, Teebi A, Chitayat D, Weksberg R, Thompson A, Vardy C, Crosbie V, Luscombe S, Baatjes R: Structural variation of chromosomes in autism spectrum disorder. Am J Hum Genet. 2008, 82: 477-488. 10.1016/j.ajhg.2007.12.009.PubMed CentralView ArticlePubMedGoogle Scholar
- Christian SL, Brune CW, Sudi J, Kumar RA, Liu S, Karamohamed S, Badner JA, Matsui S, Conroy J, McQuaid D, Gergel J, Hatchwell E, Gilliam TC, Gershon ES, Nowak NJ, Dobyns WB, Cook EH: Novel submicroscopic chromosomal abnormalities detected in autism spectrum disorder. Biol Psychiatry. 2008, 63: 1111-1117. 10.1016/j.biopsych.2008.01.009.PubMed CentralView ArticlePubMedGoogle Scholar
- Glessner JT, Wang K, Cai G, Korvatska O, Kim CE, Wood S, Zhang H, Estes A, Brune CW, Bradfield JP, Imielinski M, Frackelton EC, Reichert J, Crawford EL, Munson J, Sleiman PM, Chiavacci R, Annaiah K, Thomas K, Hou C, Glaberson W, Flory J, Otieno F, Garris M, Soorya L, Klei L, Piven J, Meyer KJ, Anagnostou E, Sakurai T: Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature. 2009, 459: 569-573. 10.1038/nature07953.PubMed CentralView ArticlePubMedGoogle Scholar
- Bucan M, Abrahams BS, Wang K, Glessner JT, Herman EI, Sonnenblick LI, Alvarez Retuerto AI, Imielinski M, Hadley D, Bradfield JP, Kim C, Gidaya NB, Lindquist I, Hutman T, Sigman M, Kustanovich V, Lajonchere CM, Singleton A, Kim J, Wassink TH, McMahon WM, Owley T, Sweeney JA, Coon H, Nurnberger JI, Li M, Cantor RM, Minshew NJ, Sutcliffe JS, Cook EH: Genome-wide analyses of exonic copy number variants in a family-based study point to novel autism susceptibility genes. PLoS Genet. 2009, 5: e1000536-10.1371/journal.pgen.1000536.PubMed CentralView ArticlePubMedGoogle Scholar
- Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, Regan R, Conroy J, Magalhaes TR, Correia C, Abrahams BS, Almeida J, Bacchelli E, Bader GD, Bailey AJ, Baird G, Battaglia A, Berney T, Bolshakova N, Bölte S, Bolton PF, Bourgeron T, Brennan S, Brian J, Bryson SE, Carson AR, Casallo G, Casey J, Chung BH, Cochrane L, Corsello C: Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010, 466: 368-372. 10.1038/nature09146.PubMed CentralView ArticlePubMedGoogle Scholar
- Szatmari P, Paterson AD, Zwaigenbaum L, Roberts W, Brian J: Mapping autism risk loci using genetic linkage and chromosomal rearrangements. Nat Genet. 2007, 39: 319-328. 10.1038/ng1985.View ArticlePubMedGoogle Scholar
- Sanders SJ, Ercan-Sencicek AG, Hus V, Luo R, Murtha MT, Moreno-De-Luca D, Chu SH, Moreau MP, Gupta AR, Thomson SA, Mason CE, Bilguvar K, Celestino-Soper PB, Choi M, Crawford EL, Davis L, Wright NR, Dhodapkar RM, DiCola M, DiLullo NM, Fernandez TV, Fielding-Singh V, Fishman DO, Frahm S, Garagaloyan R, Goh GS, Kammela S, Klei L, Lowe JK, Lund SC: Multiple recurrent de novo CNVs, including duplications of the 7q11.23 Williams syndrome region, are strongly associated with autism. Neuron. 2011, 70: 863-885. 10.1016/j.neuron.2011.05.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Weiss LA, Shen Y, Korn JM, Arking DE, Miller DT, Fossdal R, Saemundsen E, Stefansson H, Ferreira MA, Green T, Platt OS, Ruderfer DM, Walsh CA, Altshuler D, Chakravarti A, Tanzi RE, Stefansson K, Santangelo SL, Gusella JF, Sklar P, Wu BL, Daly MJ: Autism consortium: association between microdeletion and microduplication at 16p11.2 and autism. N Engl J Med. 2008, 358: 667-675. 10.1056/NEJMoa075974.View ArticlePubMedGoogle Scholar
- Morrow EM, Yoo SY, Flavell SW, Kim TK, Lin Y, Hill RS, Mukaddes NM, Balkhy S, Gascon G, Hashmi A, Al-Saad S, Ware J, Joseph RM, Greenblatt R, Gleason D, Ertelt JA, Apse KA, Bodell A, Partlow JN, Barry B, Yao H, Markianos K, Ferland RJ, Greenberg ME, Walsh CA: Identifying autism loci and genes by tracing recent shared ancestry. Science. 2008, 321: 218-223. 10.1126/science.1157657.PubMed CentralView ArticlePubMedGoogle Scholar
- Jacquemont ML, Sanlaville D, Redon R, Raoul O, Cormier-Daire V, Lyonnet S, Amiel J, Le Merrer M, Heron D, de Blois MC, Prieur M, Vekemans M, Carter NP, Munnich A, Colleaux L, Philippe A: Array-based comparative genomic hybridization identifies high frequency of cryptic chromosomal rearrangements in patients with syndromic autism spectrum disorders. J Med Genet. 2006, 43: 843-849. 10.1136/jmg.2006.043166.PubMed CentralView ArticlePubMedGoogle Scholar
- Matsunami N, Hadley D, Hensel CH, Christensen GB, Kim C, Frackelton E, Thomas K, da Silva RP, Stevens J, Baird L, Otterud B, Ho K, Varvil T, Leppert T, Lambert CG, Leppert M, Hakonarson H: Identification of rare recurrent copy number variants in high-risk autism families and their prevalence in a large ASD population. PLoS One. 2013, 8 (1): e52239-10.1371/journal.pone.0052239.PubMed CentralView ArticlePubMedGoogle Scholar
- Allen-Brady K, Robison R, Cannon D, Varvil T, Villalobos M, Pingree C, Leppert MF, Miller J, McMahon WM, Coon H: Genome-wide linkage in Utah autism pedigrees. Mol Psychiatry. 2010, 15 (10): 1006-1015. 10.1038/mp.2009.42.PubMed CentralView ArticlePubMedGoogle Scholar
- Coon H, Matsunami N, Stevens J, Miller J, Pingree C, Camp NJ, Thomas A, Krasny L, Lainhart J, Leppert MF, McMahon W: Evidence for linkage on chromosome 3q25-27 in a large autism extended pedigree. Hum Hered. 2005, 60 (4): 220-226. 10.1159/000090546.View ArticlePubMedGoogle Scholar
- Allen-Brady K, Miller J, Matsunami N, Stevens J, Block H, Farley M, Krasny L, Pingree C, Lainhart J, Leppert M, McMahon WM, Coon H: A high-density SNP genome-wide linkage scan in a large autism extended pedigree. Mol Psychiatry. 2009, 14 (6): 590-600. 10.1038/mp.2008.14.View ArticlePubMedGoogle Scholar
- BRLMM: An improved genotype calling method for the GeneChip® human mapping 500K array Set.http://media.affymetrix.com/support/technical/whitepapers/brlmm_whitepaper.pdf,
- O’Connell JR, Weeks DE: PedCheck: a program for identification of genotype incompatibilities in linkage analysis. Am J Hum Genet. 1998, 63 (1): 259-266. 10.1086/301904.PubMed CentralView ArticlePubMedGoogle Scholar
- Arrington CB, Bleyl SB, Matsunami N, Bowles NE, Leppert TI, Demarest BL, Osborne K, Yoder BA, Byrne JL, Schiffman JD, Null DM, DiGeronimo R, Rollins M, Faix R, Comstock J, Camp NJ, Leppert MF, Yost HJ, Brunelli L: A family-based paradigm to identify candidate chromosomal regions for isolated congenital diaphragmatic hernia. Am J Med Genet A. 2012, 158A (12): 3137-3147. 10.1002/ajmg.a.35664.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18 (11): 1851-1858. 10.1101/gr.078212.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Hillier LW, Marth GT, Quinlan AR, Dooling D, Fewell G, Barnett D, Fox P, Glasscock JI, Hickenbotham M, Huang W, Magrini VJ, Richt RJ, Sander SN, Stewart DA, Stromberg M, Tsung EF, Wylie T, Schedl T, Wilson RK, Mardis ER: Whole-genome sequencing and variant discovery in C. elegans. Nat Methods. 2008, 5 (2): 183-188. 10.1038/nmeth.1179.View ArticlePubMedGoogle Scholar
- Yandell M, Huff C, Hu H, Singleton M, Moore B, Xing J, Jorde LB, Reese MG: A probabilistic disease-gene finder for personal genomes. Genome Res. 2011, 21 (9): 1529-1542. 10.1101/gr.123158.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang K, Li M, Hakonarson H: ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38 (16): e164-10.1093/nar/gkq603.PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12 (6): 996-1006. 10.1101/gr.229102. Article published online before print in May 2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer LR, Zweig AS, Hinrichs AS, Karolchik D, Kuhn RM, Wong M, Sloan CA, Rosenbloom KR, Roe G, Rhead B, Raney BJ, Pohl A, Malladi VS, Li CH, Lee BT, Learned K, Kirkup V, Hsu F, Heitner S, Harte RA, Haeussler M, Guruvadoo L, Goldman M, Giardine BM, Fujita PA, Dreszer TR, Diekhans M, Cline MS, Clawson H, Barber GP: The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013, 41 (Database issue): D64-D69.PubMed CentralView ArticlePubMedGoogle Scholar
- Junutula JR, Schonteich E, Wilson GM, Peden AA, Scheller RH, Prekeris R: Molecular characterization of Rab11 interactions with members of the family of Rab11-interacting proteins. J Biol Chem. 2004, 279 (32): 33430-33437. 10.1074/jbc.M404633200.View ArticlePubMedGoogle Scholar
- Peden AA, Schonteich E, Chun J, Junutula JR, Scheller RH, Prekeris R: The RCP-Rab11 complex regulates endocytic protein sorting. Mol Biol Cell. 2004, 15 (8): 3530-3541. 10.1091/mbc.E03-12-0918.PubMed CentralView ArticlePubMedGoogle Scholar
- Roohi J, Tegay DH, Pomeroy JC, Burkett S, Stone G, Stanyon R, Hatchwell E: A de novo apparently balanced translocation [46, XY, t(2;9)(p13;p24)] interrupting RAB11FIP5 identifies a potential candidate gene for autism spectrum disorder. Am J Med Genet B Neuropsychiatr Genet. 2008, 147B (4): 411-417. 10.1002/ajmg.b.30755.View ArticlePubMedGoogle Scholar
- Hong SE, Shugart YY, Huang DT, Shahwan SA, Grant PE, Hourihane JO, Martin ND, Walsh CA: Autosomal recessive lissencephaly with cerebellar hypoplasia is associated with human RELN mutations. Nat Genet. 2000, 26 (1): 93-96. 10.1038/79246.View ArticlePubMedGoogle Scholar
- Zaki M, Shehab M, El-Aleem AA, Abdel-Salam G, Koeller HB, Ilkin Y, Ross ME, Dobyns WB, Gleeson JG: Identification of a novel recessive RELN mutation using a homozygous balanced reciprocal translocation. Am J Med Genet A. 2007, 143A (9): 939-944. 10.1002/ajmg.a.31667.View ArticlePubMedGoogle Scholar
- Béna F, Bruno DL, Eriksson M, van Ravenswaaij-Arts C, Stark Z, Dijkhuizen T, Gerkes E, Gimelli S, Ganesamoorthy D, Thuresson AC, Labalme A, Till M, Bilan F, Pasquier L, Kitzis A, Dubourgm C, Rossi M, Bottani A, Gagnebin M, Sanlaville D, Gilbert-Dussardier B, Guipponi M, van Haeringen A, Kriek M, Ruivenkamp C, Antonarakis SE, Anderlid BM, Slater HR, Schoumans J: Molecular and clinical characterization of 25 individuals with exonic deletions of NRXN1 and comprehensive review of the literature. Am J Med Genet B Neuropsychiatr Genet. 2013, 162B (4): 388-403.View ArticlePubMedGoogle Scholar
- Nag A, Bochukova EG, Kremeyer B, Campbell DD, Muller H, Valencia-Duarte AV, Cardona J, Rivas IC, Mesa SC, Cuartas M, Garcia J, Bedoya G, Cornejo W, Herrera LD, Romero R, Fournier E, Reus VI, Lowe TL, Farooqi IS, Mathews CA, McGrath LM, Yu D, Cook E, Wang K, Scharf JM, Pauls DL, Freimer NB, Plagnol V, Ruiz-Linares A: Tourette syndrome association international consortium for genetics: CNV analysis in Tourette syndrome implicates large genomic rearrangements in COL8A1 and NRXN1. PLoS One. 2013, 8 (3): e59061-10.1371/journal.pone.0059061.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaaf CP, Boone PM, Sampath S, Williams C, Bader PI, Mueller JM, Shchelochkov OA, Brown CW, Crawford HP, Phalen JA, Tartaglia NR, Evans P, Campbell WM, Tsai AC, Parsley L, Grayson SW, Scheuerle A, Luzzi CD, Thomas SK, Eng PA, Kang SH, Patel A, Stankiewicz P, Cheung SW: Phenotypic spectrum and genotype-phenotype correlations of NRXN1 exon deletions. Eur J Hum Genet. 2012, 20 (12): 1240-1247. 10.1038/ejhg.2012.95.PubMed CentralView ArticlePubMedGoogle Scholar
- Poelmans G, Franke B, Pauls DL, Glennon JC, Buitelaar JK: AKAPs integrate genetic findings for autism spectrum disorders. Transl Psychiatry. 2013, 3: e270-10.1038/tp.2013.48.PubMed CentralView ArticlePubMedGoogle Scholar
- Liao HM, Fang JS, Chen YJ, Wu KL, Lee KF, Chen CH: Clinical and molecular characterization of a transmitted reciprocal translocation t(1;12)(p32.1;q21.3) in a family co-segregating with mental retardation, language delay, and microcephaly. BMC Med Genet. 2011, 12: 70-PubMed CentralView ArticlePubMedGoogle Scholar
- Schonteich E, Wilson GM, Burden J, Hopkins CR, Anderson K, Goldenring JR, Prekeris R: The Rip11/Rab11-FIP5 and kinesin II complex regulates endocytic protein recycling. J Cell Sci. 2008, 121 (Pt 22): 3824-3833.PubMed CentralView ArticlePubMedGoogle Scholar
- Kilpinen H, Ylisaukko-oja T, Rehnström K, Gaál E, Turunen JA, Kempas E, von Wendt L, Varilo T, Peltonen L: Linkage and linkage disequilibrium scan for autism loci in an extended pedigree from Finland. Hum Mol Genet. 2009, 18 (15): 2912-2921. 10.1093/hmg/ddp229.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen CP, Lin SP, Chern SR, Chen YJ, Tsai FJ, Wu PC, Wang W: Array-CGH detection of a de novo 2.8 Mb deletion in 2q24.2→q24.3 in a girl with autistic features and developmental delay. Eur J Med Genet. 2010, 53 (4): 217-220. 10.1016/j.ejmg.2010.03.006.View ArticlePubMedGoogle Scholar
- Uz E, Alanay Y, Aktas D, Vargel I, Gucer S, Tuncbilek G, von Eggeling F, Yilmaz E, Deren O, Posorski N, Ozdag H, Liehr T, Balci S, Alikasifoglu M, Wollnik B, Akarsu NA: Disruption of ALX1 causes extreme microphthalmia and severe facial clefting: expanding the spectrum of autosomal-recessive ALX-related frontonasal dysplasia. Am J Hum Genet. 2010, 86 (5): 789-796. 10.1016/j.ajhg.2010.04.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Mori N, Kuwamura M, Tanaka N, Hirano R, Nabe M, Ibuki M, Yamate J: Ccdc85c encoding a protein at apical junctions of radial glia is disrupted in hemorrhagic hydrocephalus (hhy) mice. Am J Pathol. 2012, 180 (1): 314-327. 10.1016/j.ajpath.2011.09.014.View ArticlePubMedGoogle Scholar
- Hamdan FF, Gauthier J, Araki Y, Lin DT, Yoshizawa Y, Higashi K, Park AR, Spiegelman D, Dobrzeniecka S, Piton A, Tomitori H, Daoud H, Massicotte C, Henrion E, Diallo O, Shekarabi M, Marineau C, Shevell M, Maranda B, Mitchell G, Nadeau A, D’Anjou G, Vanasse M, Srour M, Lafrenière RG, Drapeau P, Lacaille JC, Kim E, S2D Group:Excess of de novo deleterious mutations in genes associated with glutamatergic systems in nonsyndromic intellectual disability. Am J Hum Genet. 2011, 88 (3): 306-316. 10.1016/j.ajhg.2011.02.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Majerus PW, Wilson DB, Zhang C, Nicholas PJ, Wilson MP: Expression of inositol 1,3,4-trisphosphate 5/6-kinase (ITPK1) and its role in neural tube defects. Adv Enzyme Regul. 2010, 50 (1): 365-372. 10.1016/j.advenzreg.2009.10.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Marzinke MA, Clagett-Dame M: The all-trans retinoic acid (atRA)-regulated gene Calmin (Clmn) regulates cell cycle exit and neurite outgrowth in murine neuroblastoma (Neuro2a) cells. Exp Cell Res. 2012, 318 (1): 85-93. 10.1016/j.yexcr.2011.10.002.View ArticlePubMedGoogle Scholar
- Wong YH, Lu AC, Wang YC, Cheng HC, Chang C, Chen PH, Yu JY, Fann MJ: Protogenin defines a transition stage during embryonic neurogenesis and prevents precocious neuronal differentiation. J Neurosci. 2010, 30 (12): 4428-4439. 10.1523/JNEUROSCI.0473-10.2010.View ArticlePubMedGoogle Scholar
- Ghosh M, Loper R, Gelb MH, Leslie CC: Identification of the expressed form of human cytosolic phospholipase A2beta (cPLA2beta): cPLA2beta3 is a novel variant localized to mitochondria and early endosomes. J Biol Chem. 2006, 281 (24): 16615-16624. 10.1074/jbc.M601770200.View ArticlePubMedGoogle Scholar
- Sherman EA, Strauss KA, Tortorelli S, Bennett MJ, Knerr I, Morton DH, Puffenberger EG: Genetic mapping of glutaric aciduria, type 3, to chromosome 7 and identification of mutations in c7orf10. Am J Hum Genet. 2008, 83 (5): 604-609. 10.1016/j.ajhg.2008.09.018.PubMed CentralView ArticlePubMedGoogle Scholar
- Korotchkina LG, Patel MS: Site specificity of four pyruvate dehydrogenase kinase isoenzymes toward the three phosphorylation sites of human pyruvate dehydrogenase. J Biol Chem. 2001, 276 (40): 37223-37229. 10.1074/jbc.M103069200.View ArticlePubMedGoogle Scholar
- Meyer B, Wittig I, Trifilieff E, Karas M, Schägger H: Identification of two proteins associated with mammalian ATP synthase. Mol Cell Proteomics. 2007, 6 (10): 1690-1699. 10.1074/mcp.M700097-MCP200.View ArticlePubMedGoogle Scholar
- Jarczak J, Kościuczuk EM, Lisowski P, Strzałkowska N, Jóźwik A, Horbańczuk J, Krzyżewski J, Zwierzchowski L, Bagnicka E: Defensins: natural component of human innate immunity. Hum Immunol. 2013, 74 (9): 1069-1079. 10.1016/j.humimm.2013.05.008.View ArticlePubMedGoogle Scholar
- Holweg A, Schnare M, Gessner A: The bactericidal/permeability-increasing protein (BPI) in the innate defense of the lower airways. Biochem Soc Trans. 2011, 39 (4): 1045-1050. 10.1042/BST0391045.View ArticlePubMedGoogle Scholar
- Tokunaga F, Iwai K: Linear ubiquitination: a novel NF-κB regulatory mechanism for inflammatory and immune responses by the LUBAC ubiquitin ligase complex. Endocr J. 2012, 59 (8): 641-652. 10.1507/endocrj.EJ12-0148.View ArticlePubMedGoogle Scholar
- Nguyen H, Hiscott J, Pitha PM: The growing family of interferon regulatory factors. Cytokine Growth Factor Rev. 1997, 8 (4): 293-312. 10.1016/S1359-6101(97)00019-1.View ArticlePubMedGoogle Scholar
- Kim YS, Leventhal BL, Koh YJ, Fombonne E, Laska E, Lim EC, Cheon KA, Kim SJ, Kim YK, Lee H, Song DH, Grinker RR: Prevalence of autism spectrum disorders in a total population sample. Am J Psychiatry. 2011, 168 (9): 904-912.View ArticlePubMedGoogle Scholar
- Center for Disease Control and Prevention.http://www.cdc.gov/ncbddd/autism/data.html,
- Prekeris R, Klumperman J, Scheller RH: A Rab11/Rip11 protein complex regulates apical membrane trafficking via recycling endosomes. Mol Cell. 2000, 6 (6): 1437-1448. 10.1016/S1097-2765(00)00140-4.View ArticlePubMedGoogle Scholar
- Hales CM, Griner R, Hobdy-Henderson KC, Dorn MC, Hardy D, Kumar R, Navarre J, Chan EK, Lapierre LA, Goldenring JR: Identification and characterization of a family of Rab11-interacting proteins. J Biol Chem. 2001, 276 (42): 39067-39075. 10.1074/jbc.M104831200.View ArticlePubMedGoogle Scholar
- Fernández-Monreal M, Brown TC, Royo M, Esteban JA: The balance between receptor recycling and trafficking toward lysosomes determines synaptic strength during long-term depression. J Neurosci. 2012, 32 (38): 13200-13205. 10.1523/JNEUROSCI.0061-12.2012.View ArticlePubMedGoogle Scholar
- Matthies HJ, Moore JL, Saunders C, Matthies DS, Lapierre LA, Goldenring JR, Blakely RD, Galli A: Rab11 supports amphetamine-stimulated norepinephrine transporter trafficking. J Neurosci. 2010, 30 (23): 7863-7877. 10.1523/JNEUROSCI.4574-09.2010.PubMed CentralView ArticlePubMedGoogle Scholar
- van der Sluijs P, Hoogenraad CC: New insights in endosomal dynamics and AMPA receptor trafficking. Semin Cell Dev Biol. 2011, 22 (5): 499-505. 10.1016/j.semcdb.2011.06.008.View ArticlePubMedGoogle Scholar
- Park M, Salgado JM, Ostroff L, Helton TD, Robinson CG, Harris KM, Ehlers MD: Plasticity-induced growth of dendritic spines by exocytic trafficking from recycling endosomes. Neuron. 2006, 52 (5): 817-830. 10.1016/j.neuron.2006.09.040.PubMed CentralView ArticlePubMedGoogle Scholar
- Lazo OM, Gonzalez A, Ascaño M, Kuruvilla R, Couve A, Bronfman FC: BDNF regulates Rab11-mediated recycling endosome dynamics to induce dendritic branching. J Neurosci. 2013, 33 (14): 6112-6122. 10.1523/JNEUROSCI.4630-12.2013.PubMed CentralView ArticlePubMedGoogle Scholar
- Penzes P, Cahill ME, Jones KA, VanLeeuwen JE, Woolfrey KM: Dendritic spine pathology in neuropsychiatric disorders. Nat Neurosci. 2011, 14 (3): 285-293. 10.1038/nn.2741.PubMed CentralView ArticlePubMedGoogle Scholar
- Ebert DH, Greenberg ME: Activity-dependent neuronal signaling and autism spectrum disorder. Nature. 2013, 493 (7432): 327-337. 10.1038/nature11860.PubMed CentralView ArticlePubMedGoogle Scholar
- Piggot J, Shirinyan D, Shemmassian S, Vazirian S, Alarcón M: Neural systems approaches to the neurogenetics of autism spectrum disorders. Neuroscience. 2009, 164 (1): 247-256. 10.1016/j.neuroscience.2009.05.054.View ArticlePubMedGoogle Scholar
- Lin L, Sun W, Kung F, Dell’Acqua ML, Hoffman DA: AKAP79/150 impacts intrinsic excitability of hippocampal neurons through phospho-regulation of A-type K + channel trafficking. J Neurosci. 2011, 31 (4): 1323-1332. 10.1523/JNEUROSCI.5383-10.2011.PubMed CentralView ArticlePubMedGoogle Scholar
- Westphal RS, Tavalin SJ, Lin JW, Alto NM, Fraser ID, Langeberg LK, Sheng M, Scott JD: Regulation of NMDA receptors by an associated phosphatase-kinase signaling complex. Science. 1999, 285 (5424): 93-96. 10.1126/science.285.5424.93.View ArticlePubMedGoogle Scholar
- Marx SO, Kurokawa J, Reiken S, Motoike H, D’Armiento J, Marks AR, Kass RS: Requirement of a macromolecular signaling complex for beta adrenergic receptor modulation of the KCNQ1-KCNE1 potassium channel. Science. 2002, 295 (5554): 496-499. 10.1126/science.1066843.View ArticlePubMedGoogle Scholar
- Tu H, Tang TS, Wang Z, Bezprozvanny I: Association of type 1 inositol 1,4,5-trisphosphate receptor with AKAP9 (Yotiao) and protein kinase A. J Biol Chem. 2004, 279 (18): 19375-19382. 10.1074/jbc.M313476200.View ArticlePubMedGoogle Scholar
- Chen L, Marquardt ML, Tester DJ, Sampson KJ, Ackerman MJ, Kass RS: Mutation of an A-kinase-anchoring protein causes long-QT syndrome. Proc Natl Acad Sci USA. 2007, 104 (52): 20990-20995. 10.1073/pnas.0710527105.PubMed CentralView ArticlePubMedGoogle Scholar
- Keith DJ, Sanderson JL, Gibson ES, Woolfrey KM, Robertson HR, Olszewski K, Kang R, El-Husseini A, Dell’acqua ML: Palmitoylation of A-kinase anchoring protein 79/150 regulates dendritic endosomal targeting and synaptic plasticity mechanisms. J Neurosci. 2012, 32 (21): 7119-7136. 10.1523/JNEUROSCI.0784-12.2012.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi L, Zhang X, Golhar R, Otieno FG, He M, Hou C, Kim C, Keating B, Lyon GJ, Wang K, Hakonarson H: Whole-genome sequencing in an autism multiplex family. Mol Autism. 2013, 4 (1): 8-10.1186/2040-2392-4-8.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.