|
|
||||||||||||||
|
|
|||||||||||||||
ELECTRONIC LETTER |
1 Department of Pediatrics, Wayne State University, Detroit, MI 48201, USA
2 Department of Neurology, Wayne State University, Detroit, MI 48201, USA
Correspondence to:
Dr A H M Mahbubul Huq, Division of Neurology, Childrens Hospital of Michigan, Detroit, MI 48201, USA;
ahuq{at}med.wayne.edu]
Keywords: 7q36; EN2; autism
A utistic spectrum disorders (MIM 209850), which include autism, Asperger syndrome, and pervasive developmental disorder (PDD) not otherwise specified, are characterised by impairment in communications and social interactions and the presence of stereotyped behaviours. The aetiology of autistic disorder is unknown, but family and twin studies have shown a high monozygotic to dizygotic twin risk ratio and a sib relative risk between 50 and 100, suggesting that inheritance of autism is complex, but the predisposition to develop it is largely genetically determined.1,2
A number of morphological abnormalities including increased brain size and developmental abnormalities of the cerebral cortex, brainstem, and cerebellum have been reported in autism.3,4 Reported brain stem and cerebellar anomalies include hypoplasia of lobules VI and VII, decreased Purkinje cell density, olivary dysplasia, and neuronal heterotopias.35 The morphological abnormalities described in the brainstem and cerebellum of autistic subjects suggest that genes involved in cerebellar development are candidate genes in autism.
EN2, a human homologue of Drosophila engrailed gene, is a homeobox gene with an essential role in the development of the midbrain and cerebellum.6 Mice homozygous for a targeted deletion of the EN2 homeobox region were viable but showed abnormal foliation of the cerebellum.7 Petit et al8 studied two restriction fragment length polymorphisms in the EN2 gene in autism and found an association between autism and a PvuII polymorphism in the 5' region of EN2 in 100 autistic children and 100 control children.
In this study, we have attempted to replicate the association between the EN2 gene and autism using family based linkage and association studies in 196 multiplex families with autism. The EN2 gene maps to human chromosome 7q36. Information regarding a number of promoter, exon, and intron polymorphisms in the EN2 gene is available in JSNP, a database of common gene variations in the Japanese population.9 In addition, a number of unconfirmed SNPs in the EN2 gene are listed in dbSNP. We investigated a G/A SNP (JSNP ID IMS-JST089559) in exon 1 of the EN2 gene, as this SNP is located near the PvuII polymorphism studied by Petit et al,8 alters the amino acid composition of the EN2 protein, and is a common polymorphism with allele frequencies suitable for genetic association studies.
SUBJECTS AND METHODS
Family ascertainment and diagnostic criteria
DNA samples from 204 families were obtained from the Autism Genetic Resource Exchange (AGRE). AGRE, which was created by the Cure Autism Now Foundation and the Human Biological Data Interchange, is a central repository of family DNA samples for genetic studies of autism.10 Diagnoses of AGRE families with autism, Asperger syndrome, and PDD were confirmed using the Autism Diagnostic Interview-Revised protocol (ADI-R). All AGRE families included at least two affected members with a diagnosis of autism, Asperger syndrome, or PDD. Children with fragile X syndrome, tuberous sclerosis, and chromosomal anomalies were excluded. Phenotypic information on the parents was not available. However, ADI-R phenotypic data including ADI-R items and algorithm scores are available for affected sibs from the AGRE web site (www.agre.org). For quantitative trait studies, we used four ADI-R items that have shown increased phenotypic correlations between sibs as quantitative traits for linkage and association studies.
Key points
|
EN2 genotyping
We genotyped a G/A SNP (JSNP ID IMS-JST089559; JSNP: a database of common gene variations in the Japanese population) in exon 1 of the EN2 gene by a restriction enzyme based assay using AluI (New England Biolab Inc, Beverly, MA).9 The SNP alters the amino acid composition of exon 1 of the EN2 gene, changing the codon 121 from leucine (CTC) to phenylalanine (TTC). We amplified the target region with the following primers: AGCTGTCCGAGTCCGAGC (forward) and CAAGCCTGGC GAAGCAG (reverse) with a PCR programme of 95°C for five minutes, 37 PCR cycles of 95°C for 30 seconds, 58°C for 45 seconds, and 72°C for 30 seconds. For the present study, a total of 890 subjects (670 children) were available.
Statistical and genetic analysis
The data were prepared using Mega2 and ASPEX.11,12 Marker allele frequencies were obtained by counting parental genotypes. Map distances were obtained from Marshfield and Genethon. Linkage analysis was performed by GENEHUNTER 2.1, ASPEX, and MERLIN software packages.1214 Non-parametric QTL analyses, Haseman-Elston (HE) regressions, and variance component analyses were computed for the quantitative sib pair data. Children from the autistic families were classified as affected or unknown and were tested for association through the use of a transmission disequilibrium test (TDT) within the GENEHUNTER package version 2.1. Association of the quantitative traits was evaluated by the QTDT (Quantitative TDT), which uses a variance component model that partitions association into between and within family components.15
RESULTS
We had a total of 226 affected sib pairs and 207 discordant sib pairs; 181 families had two affected sibs and 15 families had three affected sibs. Eight families with identical twins were excluded from the analysis. However, these families were genotyped and served to check for genotyping errors. The program Pedcheck and MERLIN were used to find Mendelian errors.14,16 The detectable genotype errors in the sample were less than <0.1%. There was no deviation from Hardy-Weinberg disequilibrium.
We analysed sharing at the EN2 locus with the ASPEX software package. There was no overall increased identity by descent sharing in affected sib pairs at the EN2 locus, although there was increased sharing of the paternal allele in affected sib pairs compared to discordant sib pairs in a subset of families where all affected sibs have developed useful phrase speech. Sharing data for maternal and paternal alleles in affected and discordant sib pairs are shown in table 1
.
|
For QTL analysis, we performed non-parametric multipoint linkage analyses in 196 families focusing on ADI-R items for expressive language level and a composite algorithm score for stereotyped behaviour. We studied four ADI-R items that have shown increased phenotypic correlations between sibs as quantitative traits for linkage studies.18,19 ADI-R items A12 (age at first word), A13 (age at phrase speech), and A19 (useful phrase speech level) were used as language traits. ADI-R composite algorithm score D Total was used as a stereotyped behaviour trait. Age at first word refers to words repeated and consistently used for the purpose of communication with reference to a particular concept. For age at first phrase, the phrase must consist of two words, one of which must be a verb. Subjects without any word or phrase speech by interview time were coded by their current age if <8 years or 96 months if >8 years. Non-parametric QTL and Haseman Elston regression studies with GENEHUNTER and variance component analyses with MERLIN did not show any evidence of linkage with p values less than 0.05.
TDT analysis of the exon 1 variant in 196 autistic disorders families did not show any evidence of association (p=0.58). TDT using the subgroup of families with a strict diagnosis of autism or subgroups with different overall levels of language did not increase evidence of association (table 2
). We also investigated the ADI-R items A12, A13, and A19 as language traits and the ADI-R algorithm score D Total for stereotyped behaviour using the QTDT package to test association between the EN2 gene and quantitative traits. We performed 5000 permutations for each of the above traits and did not find evidence of association of EN2 with any of the quantitative traits with p values <0.1.
|
There are several possible explanations for the differences in the findings in the two studies. First, there are differences in the diagnostic definition between the two studies. In our study, diagnosis of families with autism, Asperger syndrome, and PDD were confirmed using the ADI-R, while in the original study clinical evaluation was completed using the Behavior Summarized Evaluation: BSE scale.20 More importantly, we carried out family based linkage and association studies including TDT, while Petit et al8 used a case-control approach. Case-control studies may be subject to population stratification and spurious association.21
In summary, we have failed to replicate the report by Petit et al8 suggesting that the EN2 gene on 7q36 is a genetic risk factor for autistic disorder. EN2 is a small gene within two exons contained within an 8 kb genomic region. This suggests that if a susceptibility mutation exists at EN2, our study should have detected an association between the EN2 exon 1 polymorphism and the susceptibility mutation in the EN2 gene. It should be noted, however, that the extent of linkage disequilibrium can be irregular in a number of chromosomal regions.22 The negative finding may be the result of heterogeneity among our autistic disorder families. We stratified the families according to the severity or overall level of language to identify more homogeneous subgroups, but it may be useful to subtype the autistic disorder families using other criteria. In addition, a larger sample may be needed to detect a susceptibility locus with small effect.
The EN2 gene maps to human chromosome 7q36. Although most studies have implicated the chromosome 7q2134 region in the aetiology of autism,17,2327 one study reported a novel peak with a lod score of 2.13 at marker D7S483 on 7q36.17 There is also evidence of a quantitative trait locus (QTL) for language on 7q3536.18 Although our study did not provide support for a role of the EN2 gene in autistic disorder, further studies are required to investigate the presence of an additional susceptibility locus for autistic disorders further distally on the chromosome 7q36 region where EN2 is located.
ACKNOWLEDGEMENTS
The study described here was supported in part by a grant from the Childrens Hospital of Michigan Research Foundation. We acknowledge support from Cure Autism Now and the Autism Genetic Resource Exchange (AGRE). We gratefully acknowledge the resources provided by the AGRE consortium* and the participating AGRE families.
*The AGRE Consortium. AGRE Scientific Advisory Board: Chairman: Daniel H Geschwind, UCLA, Los Angeles, CA. Members: W Ted Brown, NYS Institute for Basic Research in Developmental Disabilities, Long Island, NY; Maya Bucan, University of Pennsylvania, Philadelphia, PA; Joseph Buxbaum, Mount Sinai School of Medicine, NY; Edwin H Cook, University of Chicago, Chicago, IL; T Conrad Gilliam, Columbia Genome Center, NY; David A Greenberg, Mount Sinai Medical Center, NY; David H Ledbetter, University of Chicago, Chicago, IL; Bruce Miller, UCSF, San Francisco, CA; Stanley F Nelson, UCLA School of Medicine, Los Angeles, CA; Jonathan Pevsner, Kennedy Krieger Institute, Baltimore, MD; Jerome I Rotter, Cedars-Sinai Medical Center, Los Angeles, CA; Gerard D Schellenberg, University of Washington, Seattle, WA; Carol A Sprouse, Childrens National Medical Center, Baltimore, MD; Rudolph E Tanzi, Massachusetts General Hospital, Boston, MA; Kirk C Wilhelmsen, University of California, San Francisco, CA; and Jeremy M Silverman, Mount Sinai Medical School, New York, NY.
REFERENCES
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS | REGISTER |