Autistic disorder and chromosome 15q11-q13: construction and analysis of a BAC/PAC contig

Autistic disorder (AD) is a neurodevelopmental disorder that affects approximately 2-10/10,000 individuals. Chromosome 15q11-q13 has been implicated in the genetic etiology of AD based on (1) cytogenetic abnormalities; (2) increased recombination frequency in this region in AD versus non-AD families; (3) suggested linkage with markers D15S156, D15S219, and D15S217; and (4) evidence for significant association with polymorphisms in the gamma-aminobutyric acid receptor subunit B3 gene (GABRB3). To isolate the putative 15q11-q13 candidate AD gene, a genomic contig and physical map of the approximately 1.2-Mb region from the GABA receptor gene cluster to the OCA2 locus was generated. Twenty-one bacterial artificial chromosome (BAC) clones, 32 P1-derived artificial chromosome (PAC) clones, and 2 P1 clones have been isolated using the markers D15S540, GABRB3, GABRA5, GABRG3, D15S822, and D15S217, as well as 34 novel markers developed from the end sequences of BAC/PAC clones. In contrast to previous findings, the markers D15S822 and D15S975 have been localized within the GABRG3 gene, which we have shown to be approximately 250 kb in size. NotI and numerous EagI restriction enzyme cut sites were identified in this region. The BAC/PAC genomic contig can be utilized for the study of genomic structure and the identification and characterization of genes and their methylation status in this autism candidate gene region on human chromosome 15q11-q13.     

Autistic Disorder and Chromosome 15q11–q13: Construction and Analysis of a BAC/PAC Contig

Autistic disorder (AD) is a neurodevelopmental disorder that affects approximately 2–10/10,000 individuals. Chromosome 15q11–q13 has been implicated in the genetic etiology of AD based on (1) cytogenetic abnormalities; (2) increased recombination frequency in this region in AD versus non-AD families; (3) suggested linkage with markers D15S156, D15S219, and D15S217; and (4) evidence for significant association with polymorphisms in the γ-aminobutyric acid receptor subunit B3 gene (GABRB3). To isolate the putative 15q11–q13 candidate AD gene, a genomic contig and physical map of the approximately 1.2-Mb region from the GABA receptor gene cluster to the OCA2 locus was generated. Twenty-one bacterial artificial chromosome (BAC) clones, 32 P1-derived artificial chromosome (PAC) clones, and 2 P1 clones have been isolated using the markers D15S540, GABRB3, GABRA5, GABRG3, D15S822, and D15S217, as well as 34 novel markers developed from the end sequences of BAC/PAC clones. In contrast to previous findings, the markers D15S822 and D15S975 have been localized within the GABRG3 gene, which we have shown to be approximately 250 kb in size. NotI and numerous EagI restriction enzyme cut sites were identified in this region. The BAC/PAC genomic contig can be utilized for the study of genomic structure and the identification and characterization of genes and their methylation status in this autism candidate gene region on human chromosome 15q11–q13.     

Drosophila Cyclodiene Resistance Gene Shows Conserved Genomic Organization with Vertebrate γ-Aminobutyric AcidA Receptors

Genomic clones from the Rdl locus of Drosophila, whose mutant phenotype is resistant to cyclodiene insecticides and picrotoxin, were characterized by restriction mapping and partial sequencing to determine intron/exon structure. The coding region of the gene comprises nine identified exons and spans greater than 25 kb of genomic DNA. The structure of the Drosophila Rdl receptor subunit was compared with those of vertebrate gamma-aminobutyric acid subtype A (GABAA) receptors and nicotinic acetylcholine receptors (nAChRs). The first six introns in Rdl show positions similar to those in vertebrate GABAA receptors, whereas the last two differ. It is interesting that the last intron appears to be in a position similar to that in nAChRs. These results are examined in relation to the proposal, based on amino acid identities, that Rdl codes for a novel class of GABAA receptor subunit more closely related to glycine receptors, and the possible place of Rdl in the lineage of the receptor superfamily is discussed.     

Mapping of four translocation breakpoints within the Duchenne muscular dystrophy gene

There are over 20 females with Duchenne or Becker muscular dystrophy (DMD or BMD) who have X-autosome translocations that break the X chromosome within band Xp21. Several of these translocations have been mapped with genomic probes to regions throughout the large (approximately 2000 kb) DMD gene. In this report, a cDNA clone from the 5' end of the gene was used to further map the breakpoints in four X-autosome translocations. A t(X;21) translocation in a patient with BMD and a t(X;1) translocation in a patient with DMD were found to break within a large 110-kb intron between exons 7 and 8. Two other DMD translocations, t(X;5) and t(X;11), were found to break between the first and the second exon of the gene within a presumably large intron (greater than 100 kb). These results demonstrate that all four translocations have disrupted the DMD gene and make it possible to clone and sequence the breakpoints. This will in turn determine whether these translocations occur by chance in these large introns or whether there are sequences that predispose to translocations.     

Conserved Organization of γ-Aminobutyric AcidAReceptor Genes: Cloning and Analysis of the Chicken β4-Subunit Gene

A series of genomic clones containing DNA that encodes the chicken gamma-aminobutyric acidA (GABAA) receptor beta 4 subunit have been isolated. These have been restriction mapped and partially sequenced to determine the structural organization and the size of the beta 4-subunit gene. This gene, which comprises nine exons, spans more than 65 kb. The organization of the chicken GABAA receptor beta 4-subunit gene has been compared to that of the murine GABAA receptor delta-subunit gene and to those of the genes that encode other members of the ligand-gated ion-channel superfamily, namely muscle and neuronal nicotinic acetylcholine receptors (AChRs). Although the positions of the intron/exon boundaries of GABAA receptor subunit genes are seen to be highly conserved, there are significant differences between the genes that encode GABAA receptor and AChR subunits. These results are discussed in relation to the proposal that this superfamily of ligand-gated ion-channel receptor genes arose by duplication of an ancestral receptor gene.     

Identification of significant association and gene-gene interaction of GABA receptor subunit genes in autism

Autism is a common neurodevelopmental disorder with a significant genetic component. Existing research suggests that multiple genes contribute to autism and that epigenetic effects or gene-gene interactions are likely contributors to autism risk. However, these effects have not yet been identified. Gamma-aminobutyric acid (GABA), the primary inhibitory neurotransmitter in the adult brain, has been implicated in autism etiology. Fourteen known autosomal GABA receptor subunit genes were studied to look for the genes associated with autism and their possible interactions. Single-nucleotide polymorphisms (SNPs) were screened in the following genes: GABRG1, GABRA2, GABRA4, and GABRB1 on chromosome 4p12; GABRB2, GABRA6, GABRA1, GABRG2, and GABRP on 5q34-q35.1; GABRR1 and GABRR2 on 6q15; and GABRA5, GABRB3, and GABRG3 on 15q12. Intronic and/or silent mutation SNPs within each gene were analyzed in 470 white families with autism. Initially, SNPs were used in a family-based study for allelic association analysis--with the pedigree disequilibrium test and the family-based association test--and for genotypic and haplotypic association analysis--with the genotype-pedigree disequilibrium test (geno-PDT), the association in the presence of linkage (APL) test, and the haplotype family-based association test. Next, with the use of five refined independent marker sets, extended multifactor-dimensionality reduction (EMDR) analysis was employed to identify the models with locus joint effects, and interaction was further verified by conditional logistic regression. Significant allelic association was found for markers RS1912960 (in GABRA4; P = .01) and HCV9866022 (in GABRR2; P = .04). The geno-PDT found significant genotypic association for HCV8262334 (in GABRA2), RS1912960 and RS2280073 (in GABRA4), and RS2617503 and RS12187676 (in GABRB2). Consistent with the allelic and genotypic association results, EMDR confirmed the main effect at RS1912960 (in GABRA4). EMDR also identified a significant two-locus gene-gene effect model involving RS1912960 in GABRA4 and RS2351299 in GABRB1. Further support for this two-locus model came from both the multilocus geno-PDT and the APL test, which indicated a common genotype and haplotype combination positively associated with disease. Finally, these results were also consistent with the results from the conditional logistic regression, which confirmed the interaction between GABRA4 and GABRB1 (odds ratio = 2.9 for interaction term; P = .002). Through the convergence of all analyses, we conclude that GABRA4 is involved in the etiology of autism and potentially increases autism risk through interaction with GABRB1. These results support the hypothesis that GABA receptor subunit genes are involved in autism, most likely via complex gene-gene interactions.     

Tandem linkage of human CSF-1 receptor (c-fms) and PDGF receptor genes

A 5' untranslated exon of the human CSF-1 receptor gene (c-fms) is separated by a 26 kb intron from the 32 kb receptor coding sequences. Nucleotide sequence analysis of cloned genomic DNA revealed that the 3' end of the PDGF receptor gene is located less than 0.5 kb upstream from this exon. Similarities in chromosomal localization, organization, and encoded amino acid sequences suggest that the genes encoding the CSF-1 and PDGF receptors arose through duplication. The as yet unidentified c-fms promoter/enhancer sequences may be confined to the nucleotides separating the two genes or could potentially lie within the PDGF receptor gene itself.     

Association Analysis of GABRB3 Promoter Variants with Heroin Dependence

GABRB3 encoding the β3 subunit of GABAA receptor has been implicated in multiple neuropsychiatric disorders, including substance abuse. Previous studies reported that SNPs at the 5′ regulatory region of GABRB3 could regulate GABRB3 gene expression and associated with childhood absence epilepsy (CAE). The study aimed to investigate whether SNPs at the 5′ regulatory region of GABRB3 were associated with heroin dependence in our population. We first re-sequenced 1.5 kb of the 5′regulatory region of GABRB3 gene to examine the SNP profile in the genomic DNA of 365 control subjects. Then, we conducted a case-control association analysis between 576 subjects with heroin dependence (549 males, 27 females) and 886 controls (472 males, 414 females) by genotyping the rs4906902 as a tag SNP. We also conducted a reporter gene assay to assess the promoter activity of two major haplotypes derived from SNPs at this region. We detected 3 common SNPs (rs4906902, rs8179184 and rs20317) at this region that had strong pair-wise linkage disequilibrium. The C allele of rs4906902 was found to be associated with increased risk of heroin dependence (odds ratio:1.27, p = 0.002). Two major haplotypes (C-A-G and T-G-C) derived from these 3 SNPs accounted for 99% of this sample, and reporter gene activity assay showed that haplotype C-A-G that contained the C allele of the tag SNP rs4906902 had higher activity than haplotype T-G-C. Our data suggest that GABRB3 might be associated with heroin dependence, and increased expression of GABRB3 might contribute to the pathogenesis of heroin dependence.     

Localization of the gene encoding the GABAA receptor β3 subunit to the Angelman/Prader-Willi region of human chromosome 15

Deletions of the proximal long arm of chromosome 15 (bands 15q11q13) are found in the majority of patients with two distinct genetic disorders, Angelman syndrome (AS) and Prader-Willi syndrome (PWS). The deleted regions in the two syndromes, defined cytogenetically and by using cloned DNA probes, are similar. However, deletions in AS occur on the maternally inherited chromosome 15, and deletions in PWS occur on the paternally derived chromosome 15. This observation has led to the suggestion that one or more genes in this region show differential expression dependent on parental origin (genetic imprinting). No genes of known function have previously been mapped to this region. We show here that the gene encoding the GABAA (gamma-aminobutyric acid) receptor beta 3 subunit maps to the AS/PWS region. Deletion of this gene (GABRB3) was found in AS and PWS patients with interstitial cytogenetic deletions. Evidence of beta 3 gene deletion was also found in an AS patient with an unbalanced 13;15 translocation but not in a PWS patient with an unbalanced 9;15 translocation. The localization of this receptor gene to the AS/PWS region suggests a possible role of the inhibitory neurotransmitter GABA in the pathogenesis of one or both of these syndromes.     

Variations in GABRA2, encoding the alpha 2 subunit of the GABA(A) receptor, are associated with alcohol dependence and with brain oscillations

Alcoholism is a complex disease with both genetic and environmental risk factors. To identify genes that affect the risk for alcoholism, we systematically ascertained and carefully assessed individuals in families with multiple alcoholics. Linkage and association analyses suggested that a region of chromosome 4p contained genes affecting a quantitative endophenotype, brain oscillations in the beta frequency range (13-28 Hz), and the risk for alcoholism. To identify the individual genes that affect these phenotypes, we performed linkage disequilibrium analyses of 69 single-nucleotide polymorphism (SNPs) within a cluster of four GABA(A) receptor genes, GABRG1, GABRA2, GABRA4, and GABRB1, at the center of the linked region. GABA(A) receptors mediate important effects of alcohol and also modulate beta frequencies. Thirty-one SNPs in GABRA2, but only 1 of the 20 SNPs in the flanking genes, showed significant association with alcoholism. Twenty-five of the GABRA2 SNPs, but only one of the SNPs in the flanking genes, were associated with the brain oscillations in the beta frequency. The region of strongest association with alcohol dependence extended from intron 3 past the 3' end of GABRA2; all 43 of the consecutive three-SNP haplotypes in this region of GABRA2 were highly significant. A three-SNP haplotype was associated with alcoholism, with P=.000000022. No coding differences were found between the high-risk and low-risk haplotypes, suggesting that the effect is mediated through gene regulation. The very strong association of GABRA2 with both alcohol dependence and the beta frequency of the electroencephalogram, combined with biological evidence for a role of this gene in both phenotypes, suggest that GABRA2 might influence susceptibility to alcohol dependence by modulating the level of neural excitation.     

Significant differences in gene expression of GABA receptors in peripheral blood leukocytes of migraineurs

Migraine is a debilitating neurovascular disorder, with a substantial genetic component. The exact cause of a migraine attack is unknown; however cortical hyperexcitability is thought to play a role. As Gamma-aminobutyric Acid (GABA) is the major inhibitory neurotransmitter in the brain, malfunctioning of this system may be a cause of the hyperexcitability. To date, there has been limited research examining the gene expression or genetics of GABA receptors in relation to migraine. The aim of our study was to determine if GABA receptors play a role in migraine by investigating their gene expression using profile in migraine affected individuals and non-affected controls by Q-PCR. Gene expression of GABA(A) receptor subunit isoforms (GABRA3, GABRB3, GABRQ) and GABA(B) receptor 2 (GABBR2) was quantified in mRNA obtained from peripheral blood leukocytes from 28 migraine subjects and 22 healthy control subjects. Analysis of results showed that two of the tested genes, GABRA3 and GABBR2, were significantly down regulated in migraineurs (P=0.018; P=0.017), compared to controls. Results from the other tested genes did not show significant gene expression variation. The results indicate that there may be specific GABA receptor gene expression variation in migraine, particularly involving the GABRA3 and GABBR2 genes. This study also identifies GABRA3 and GABBR2 as potential biomarkers to select migraineurs that may be more responsive to GABA agonists with future investigations in this area warranted.     

Isolation, characterization, and localization of human genomic DNA encoding the beta 1 subunit of the GABAA receptor (GABRB1).

Genomic DNA that encodes the beta 1 subunit of the human gamma-aminobutyric acidA (GABAA) receptor was cloned and mapped. Exons and flanking introns (greater than 14 kb) were sequenced to determine the structural organization of the gene. The gene was localized on human chromosome 4, in bands p12-13. The beta 1 subunit is encoded by a relatively large gene (greater than 65 kb) on nine exons. In contrast to other conserved regions of the subunit polypeptide, the proposed channel-forming domain (M2) is derived from more than one exon. The organization of exons was compared with that of the genes that code for subunits of nicotinic acetylcholine receptors. There is no evidence for conservation of gene structure between these two members of the proposed gene superfamily. However, intron-exon junctions were found to be conserved precisely between subtypes of GABAA receptor subunits.     

Using PAC nested deletions to order contigs and microsatellite markers at the high repetitive sequence containing Npr3 gene locus

Highly polymorphic di- and tetranucleotide repeats in and around Npr3, a potential candidate gene for hypertension, have been identified using a novel approach. Because this chromosomal site is rich in repetitive DNA and difficult to sequence, P1 artificial chromosomes were retrofitted with a loxP transposon to map the gene sequence within a clone using a series of nested deletions. Sequences from ends of deletions 1-3 kb apart identified a (CA)(20) and a (TA)(18)-(CA)(8) repeat 8 kb upstream and within an intron of Npr3, respectively. DNA from 17 individuals was analyzed for length polymorphisms in these and eight additional repeats identified in 200 kb of working draft sequence from this region in GenBank. The sequence contigs and microsatellite repeats from GenBank were ordered using the P1-derived artificial chromosome deletion series. Several of these repeats were found to vary considerably in length in the set of genomic DNA tested. Since this site in chromosome 5p has recently been implicated in disease in studies with genetically hypertensive rats, the microsatellite markers reported here will be useful for genetic analysis and may even be implicated in the disease process in humans. We discuss how these types of data are useful for interpreting draft DNA sequence coming out of the genome projects, and the utility of deletion clones as a resource for ordering contigs and gap filling.