Genetic investigation of chromosome 5q GABAA receptor subunit genes in schizophrenia

We previously performed a genome-wide linkage scan in Portuguese schizophrenia families that identified a risk locus on chromosome 5q31-q35. This finding was supported by meta-analysis of 20 other schizophrenia genome-wide scans that identified 5q23.2-q34 as the second most compelling susceptibility locus in the genome. In the present report, we took a two-stage candidate gene association approach to investigate a group of gamma-aminobutyric acid (GABA) A receptor subunit genes (GABRA1, GABRA6, GABRB2, GABRG2, and GABRP) within our linkage peak. These genes are plausible candidates based on prior evidence for GABA system involvement in schizophrenia. In the first stage, associations were detected in a Portuguese patient sample with single nucleotide polymorphisms (SNPs) and haplotypes in GABRA1 (P=0.00062-0.048), GABRP (P=0.0024-0.042), and GABRA6 (P=0.0065-0.0088). The GABRA1 and GABRP findings were replicated in the second stage in an independent German family-based sample (P=0.0015-0.043). Supportive evidence for association was also obtained for a previously reported GABRB2 risk haplotype. Exploratory analyses of the effects of associated GABRA1 haplotypes on transcript levels found altered expression of GABRA6 and coexpressed genes of GABRA1 and GABRB2. Comparison of transcript levels in schizophrenia patients and unaffected siblings found lower patient expression of GABRA6 and coexpressed genes of GABRA1. Interestingly, the GABRA1 coexpressed genes include synaptic and vesicle-associated genes previously found altered in schizophrenia prefrontal cortex. Taken together, these results support the involvement of the chromosome 5q GABAA receptor gene cluster in schizophrenia, and suggest that schizophrenia-associated haplotypes may alter expression of GABA-related genes.

Support for involvement of neuregulin 1 in schizophrenia pathophysiology

Schizophrenia is a common, multigenic psychiatric disorder. Linkage studies, including a recent meta-analysis of genome scans, have repeatedly implicated chromosome 8p12-p23.1 in schizophrenia susceptibility. More recently, significant association with a candidate gene on 8p12, neuregulin 1 (NRG1), has been reported in several European and Chinese samples. We investigated NRG1 for association in schizophrenia patients of Portuguese descent to determine whether this gene is a risk factor in this population. We tested NRG1 markers and haplotypes for association in 111 parent-proband trios, 321 unrelated cases, and 242 control individuals. Associations were found with a haplotype that overlaps the risk haplotype originally reported in the Icelandic population ("Hap(ICE)"), and two haplotypes located in the 3' end of NRG1 (all P<0.05). However, association was not detected with Hap(ICE) itself. Comparison of NRG1 transcript expression in peripheral leukocytes from schizophrenia patients and unaffected siblings identified 3.8-fold higher levels of the SMDF variant in patients (P=0.039). Significant positive correlations (P<0.001) were found between SMDF and HRG-beta 2 expression and between HRG-gamma and ndf43 expression, suggesting common transcriptional regulation of NRG1 variants. In summary, our results suggest that haplotypes across NRG1 and multiple NRG1 variants are involved in schizophrenia.

Genomewide linkage analysis of bipolar disorder by use of a high-density single-nucleotide-polymorphism (SNP) genotyping assay: a comparison with microsatellite marker assays and finding of significant linkage to chromosome 6q22

We performed a linkage analysis on 25 extended multiplex Portuguese families segregating for bipolar disorder, by use of a high-density single-nucleotide-polymorphism (SNP) genotyping assay, the GeneChip Human Mapping 10K Array (HMA10K). Of these families, 12 were used for a direct comparison of the HMA10K with the traditional 10-cM microsatellite marker set and the more dense 4-cM marker set. This comparative analysis indicated the presence of significant linkage peaks in the SNP assay in chromosomal regions characterized by poor coverage and low information content on the microsatellite assays. The HMA10K provided consistently high information and enhanced coverage throughout these regions. Across the entire genome, the HMA10K had an average information content of 0.842 with 0.21-Mb intermarker spacing. In the 12-family set, the HMA10K-based analysis detected two chromosomal regions with genomewide significant linkage on chromosomes 6q22 and 11p11; both regions had failed to meet this strict threshold with the microsatellite assays. The full 25-family collection further strengthened the findings on chromosome 6q22, achieving genomewide significance with a maximum nonparametric linkage (NPL) score of 4.20 and a maximum LOD score of 3.56 at position 125.8 Mb. In addition to this highly significant finding, several other regions of suggestive linkage have also been identified in the 25-family data set, including two regions on chromosome 2 (57 Mb, NPL = 2.98; 145 Mb, NPL = 3.09), as well as regions on chromosomes 4 (91 Mb, NPL = 2.97), 16 (20 Mb, NPL = 2.89), and 20 (60 Mb, NPL = 2.99). We conclude that at least some of the linkage peaks we have identified may have been largely undetected in previous whole-genome scans for bipolar disorder because of insufficient coverage or information content, particularly on chromosomes 6q22 and 11p11.

Multiple lineages of R1 retrotransposable elements can coexist in the rDNA loci of Drosophila

R1 non-long terminal repeat retrotransposable elements insert specifically into the 28S rRNA genes of arthropods. One aspect of R1 evolution that has been difficult to explain is the presence of divergent lineages of R1 in the rDNA loci of the same species. Multiple lineages should compete for a limited number of insertion sites, in addition to being subject to the concerted evolution processes homogenizing the rRNA genes. The presence of multiple lineages suggests either the ability of the elements to overcome these factors and diverge within rDNA loci, or the introduction of new lineages by horizontal transmission. To address this issue, we attempted to characterize the complete set of R1 elements in the rDNA locus from five Drosophila species groups (melanogaster, obscura, testacea, quinaria, and repleta). Two major R1 lineages, A and B, that diverged about 100 MYA were found to exist in Drosophila. Elements of the A lineage were found in all 35 Drosophila species tested, while elements of the B lineage were found in only 11 species from three species groups. Phylogenetic analysis of the R1 elements, supported by comparison of their rates of nucleotide sequence substitution, revealed that both the A and the B lineages have been maintained by vertical descent. The B lineage was less stable and has undergone numerous, independent elimination events, while the A lineage has diverged into three sublineages, which were, in turn, differentially stable. We conclude that while the differential retention of multiple lineages greatly complicates its phylogenetic history, the available R1 data continue to be consistent with the strict vertical descent of these elements.