Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63–0.75], p = 1.86×10−18), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21–1.43], p = 3.87×10−11). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50–0.67], p = 1.17×10−12) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53–0.72], p = 8.88×10−10). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41–0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54–1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Tissue-specific alternative splicing is an important mechanism for providing spatiotemporal protein diversity. Here we show that an in-frame splice mutation in BBS8, one of the genes involved in pleiotropic Bardet-Biedl syndrome (BBS), is sufficient to cause nonsyndromic retinitis pigmentosa (RP). A genome-wide scan of a consanguineous RP pedigree mapped the trait to a 5.6 Mb region; subsequent systematic sequencing of candidate transcripts identified a homozygous splice-site mutation in a previously unknown BBS8 exon. The allele segregated with the disorder, was absent from controls, was completely invariant across evolution, and was predicted to lead to the elimination of a 10 amino acid sequence from the protein. Subsequent studies showed the exon to be expressed exclusively in the retina and enriched significantly in the photoreceptor layer. Importantly, we found this exon to represent the major BBS8 mRNA species in the mammalian photoreceptor, suggesting that the encoded 10 amino acids play a pivotal role in the function of BBS8 in this organ. Understanding the role of this additional sequence might therefore inform the mechanism of retinal degeneration in patients with syndromic BBS or other related ciliopathies.

Combinatorial regulation by transcription factor complexes is an important feature of eukaryotic gene regulation. Here, we propose a new method for identification of interactions between transcription factors (TFs) that relies on the relationship of their binding sites, and we test it using Saccharomyces cerevisiae as a model system. The algorithm predicts interacting TF pairs based on the co-occurrence of their binding motifs and the distance between the motifs in promoter sequences. This allows investigation of interactions between TFs without known binding motifs or expression data. With this approach, 300 significant interactions involving 77 TFs were identified. These included more than 70% of the known protein–protein interactions. Approximately half of the detected interacting motif pairs showed strong preferences for particular distances and orientations in the promoter sequences. These one dimensional features may reflect constraints on allowable spatial arrangements for protein–protein interactions. Evidence for biological relevance of the observed characteristic distances is provided by the finding that target genes with the same characteristic distances show significantly higher co-expression than those without preferred distances. Furthermore, the observed interactions were dynamic: most of the TF pairs were not constitutively active, but rather showed variable activity depending on the physiological condition of the cells. Interestingly, some TF pairs active in multiple conditions showed preferences for different distances and orientations depending on the condition. Our prediction and characterization of TF interactions may help to understand the transcriptional regulatory networks in eukaryotic systems.