Probing the functional impact of sequence variation on p53-DNA interactions using a novel microsphere assay for protein-DNA binding with human cell extracts.

The p53 tumor suppressor regulates its target genes through sequence-specific binding to DNA response elements (REs). Although numerous p53 REs are established, the thousands more identified by bioinformatics are not easily subjected to comparative functional evaluation. To examine the relationship between RE sequence variation -- including polymorphisms -- and p53 binding, we have developed a multiplex format microsphere assay of protein-DNA binding (MAPD) for p53 in nuclear extracts. Using MAPD we measured sequence-specific p53 binding of doxorubicin-activated or transiently expressed p53 to REs from established p53 target genes and p53 consensus REs. To assess the sensitivity and scalability of the assay, we tested 16 variants of the p21 target sequence and a 62-multiplex set of single nucleotide (nt) variants of the p53 consensus sequence and found many changes in p53 binding that are not captured by current computational binding models. A group of eight single nucleotide polymorphisms (SNPs) was examined and binding profiles closely matched transactivation capability tested in luciferase constructs. The in vitro binding characteristics of p53 in nuclear extracts recapitulated the cellular in vivo transactivation capabilities for eight well-established human REs measured by luciferase assay. Using a set of 26 bona fide REs, we observed distinct binding patterns characteristic of transiently expressed wild type and mutant p53s. This microsphere assay system utilizes biologically meaningful cell extracts in a multiplexed, quantitative, in vitro format that provides a powerful experimental tool for elucidating the functional impact of sequence polymorphism and protein variation on protein/DNA binding in transcriptional networks.

Noncanonical DNA motifs as transactivation targets by wild type and mutant p53.

Sequence-specific binding by the human p53 master regulator is critical to its tumor suppressor activity in response to environmental stresses. p53 binds as a tetramer to two decameric half-sites separated by 0-13 nucleotides (nt), originally defined by the consensus RRRCWWGYYY (n = 0-13) RRRCWWGYYY. To better understand the role of sequence, organization, and level of p53 on transactivation at target response elements (REs) by wild type (WT) and mutant p53, we deconstructed the functional p53 canonical consensus sequence using budding yeast and human cell systems. Contrary to early reports on binding in vitro, small increases in distance between decamer half-sites greatly reduces p53 transactivation, as demonstrated for the natural TIGER RE. This was confirmed with human cell extracts using a newly developed, semi-in vitro microsphere binding assay. These results contrast with the synergistic increase in transactivation from a pair of weak, full-site REs in the MDM2 promoter that are separated by an evolutionary conserved 17 bp spacer. Surprisingly, there can be substantial transactivation at noncanonical (1/2)-(a single decamer) and (3/4)-sites, some of which were originally classified as biologically relevant canonical consensus sequences including PIDD and Apaf-1. p53 family members p63 and p73 yielded similar results. Efficient transactivation from noncanonical elements requires tetrameric p53, and the presence of the carboxy terminal, non-specific DNA binding domain enhanced transactivation from noncanonical sequences. Our findings demonstrate that RE sequence, organization, and level of p53 can strongly impact p53-mediated transactivation, thereby changing the view of what constitutes a functional p53 target. Importantly, inclusion of (1/2)- and (3/4)-site REs greatly expands the p53 master regulatory network.

Discovery and verification of functional single nucleotide polymorphisms in regulatory genomic regions: current and developing technologies.

The most common form of genetic variation, single nucleotide polymorphisms or SNPs, can affect the way an individual responds to the environment and modify disease risk. Although most of the millions of SNPs have little or no effect on gene regulation and protein activity, there are many circumstances where base changes can have deleterious effects. Non-synonymous SNPs that result in amino acid changes in proteins have been studied because of their obvious impact on protein activity. It is well known that SNPs within regulatory regions of the genome can result in disregulation of gene transcription. However, the impact of SNPs located in putative regulatory regions, or rSNPs, is harder to predict for two primary reasons. First, the mechanistic roles of non-coding genomic sequence remain poorly defined. Second, experimental validation of the functional consequences of rSNPs is often slow and laborious. In this review, we summarize traditional and novel methodologies for candidate rSNPs selection, in particular in silico techniques that aid in candidate rSNP selection. Additionally we will discuss molecular biological techniques that assess the impact of rSNPs on binding of regulatory machinery, as well as functional consequences on transcription. Standard techniques such as EMSA and luciferase reporter constructs are still widely used to assess effects of rSNPs on binding and gene transcription; however, these protocols are often bottlenecks in the discovery process. Therefore, we highlight novel and developing high-throughput protocols that promise to aid in shortening the process of rSNP validation. Given the large amount of genomic information generated from a multitude of re-sequencing and genome-wide SNP array efforts, future focus should be to develop validation techniques that will allow greater understanding of the impact these polymorphisms have on human health and disease.

Expression profiling of substantia nigra in Parkinson disease, progressive supranuclear palsy, and frontotemporal dementia with parkinsonism.

BACKGROUND: Parkinson disease (PD) is characterized by loss of dopaminergic neurons in the substantia nigra. Genes contributing to rare mendelian forms of PD have been identified, but the genes involved in the more common idiopathic PD are not well understood. OBJECTIVES: To identify genes important to PD pathogenesis using microarrays and to investigate their potential to aid in diagnosing parkinsonism. DESIGN: Microarray expression analysis of postmortem substantia nigra tissue. PATIENTS: Substantia nigra samples from 14 unrelated individuals were analyzed, including 6 with PD, 2 with progressive supranuclear palsy, 1 with frontotemporal dementia with parkinsonism, and 5 control subjects. MAIN OUTCOME MEASURES: Identification of genes significantly differentially expressed (P<.05) using Affymetrix U133A microarrays. RESULTS: There were 142 genes that were significantly differentially expressed between PD cases and controls and 96 genes that were significantly differentially expressed between the combined progressive supranuclear palsy and frontotemporal dementia with parkinsonism cases and controls. The 12 genes common to all 3 disorders may be related to secondary effects. Hierarchical cluster analysis after exclusion of these 12 genes differentiated 4 of the 6 PD cases from progressive supranuclear palsy and frontotemporal dementia with parkinsonism. CONCLUSIONS: Four main molecular pathways are altered in PD substantia nigra: chaperones, ubiquitination, vesicle trafficking, and nuclear-encoded mitochondrial genes. These results correlate well with expression analyses performed in several PD animal models. Expression analyses have promising potential to aid in postmortem diagnostic evaluation of parkinsonism.

Genomic convergence to identify candidate genes for Parkinson disease: SAGE analysis of the substantia nigra.

Genomic convergence is a multistep approach that combines gene expression with genomic linkage to identify and prioritize susceptibility genes for complex disease. As a first step, we previously performed linkage analysis on 174 multiplex Parkinson's disease (PD) families, identifying five peaks for PD risk and two for genes affecting age at onset (AAO) in PD [Hauser et al., Hum Mol Genet 2003;12:671-677]. We report here the next step: serial analysis of gene expression [SAGE; Scott et al., JAMA 2001;286:2239-2242] to analyze substantia nigra tissue from three PD patients and two age-matched controls. We find 933 differentially expressed genes (P<0.05) between PD and controls, but of these, only 50 genes represented by unique SAGE tags map within our previously described PD linkage regions. Furthermore, genes encoded by mitochondrial DNA are expressed 1.5-fold higher in PD patients versus controls, without an increase in the corresponding nuclear-encoded mitochondrial components, suggesting an increase in mtDNA genomes in PD or a disjunction with nuclear expression. The next step in the genomic convergence process will be to screen these 50 high-quality candidate genes for association with PD risk susceptibility and genetic effects on AAO.

Identification of risk and age-at-onset genes on chromosome 1p in Parkinson disease.

We previously reported a linkage region on chromosome 1p (LOD = 3.41) for genes controlling age at onset (AAO) in Parkinson disease (PD). This region overlaps with the previously reported PARK10 locus. To identify the gene(s) associated with AAO and risk of PD in this region, we first applied a genomic convergence approach that combined gene expression and linkage data. No significant results were found. Second, we performed association mapping across a 19.2-Mb region centered under the AAO linkage peak. An iterative association mapping approach was done by initially genotyping single-nucleotide polymorphisms at an average distance of 100 kb apart and then by increasing the density of markers as needed. Using the overall data set of 267 multiplex families, we identified six associated genes in the region, but further screening of a subset of 83 families linked to the chromosome 1 locus identified only two genes significantly associated with AAO in PD: the gamma subunit of the translation initiation factor EIF2B gene (EIF2B3), which was more significant in the linked subset and the ubiquitin-specific protease 24 gene (USP24). Unexpectedly, the human immunodeficiency virus enhancer-binding protein 3 gene (HIVEP3) was found to be associated with risk for susceptibility to PD. We used several criteria to define significant results in the presence of multiple testing, including criteria derived from a novel cluster approach. The known or putative functions of these genes fit well with the current suspected pathogenic mechanisms of PD and thus show great potential as candidates for the PARK10 locus.

Association between the neuron-specific RNA-binding protein ELAVL4 and Parkinson disease.

Inflammatory processes have been implicated in the cascade of events that lead to nerve cell death. In the nervous system, a number of genes involved in inflammation pathways are regulated post-transcriptionally via the interaction of their mRNAs with specific RNA-binding Hu proteins, the vertebrate homologues of the Drosophila ELAV (for embryonic lethal abnormal vision). The gene encoding ELAVL4, a member of the Hu family of proteins, is located 2 Mb from the chromosome 1p linkage region peak for age-at-onset (AAO) of Parkinson disease (PD) (LOD = 3.41). Nine single-nucleotide polymorphisms (SNPs) in ELAVL4 were genotyped for 266 multiplex families (1,223 samples). Additional genotyping in 377 singleton families was performed for a subset of five SNPs (SNPs 1-5) that were not in linkage disequilibrium. SNP 2 (located in the first intron of ELAVL4) showed a strong significant association with AAO of PD (P = 0.006), and SNP 5 (a coding SNP in ELAVL4) showed a moderately significant association (P = 0.035). Haplotype analysis revealed that the A-C haplotype at SNPs 2 and 3 has the strongest significant association with AAO (P = 0.0001) among all combinations of two or three loci. The A-C haplotype remained significant for AAO after the inclusion of the C allele at SNP 5 to this haplotype (A-C-C haplotype, P = 0.00018). Although SNP 5 was found to associate with PD risk in the early-onset subset of PD families (at least one affected with AAO <40 years, 60 families), we believe that it is a by-product of its association with AAO. Taken together, these results suggest a potential role for ELAVL4 as a modifier gene for AAO of PD.

Complement factor H variant increases the risk of age-related macular degeneration.

Age-related macular degeneration (AMD) is a leading cause of visual impairment and blindness in the elderly whose etiology remains largely unknown. Previous studies identified chromosome 1q32 as harboring a susceptibility locus for AMD. We used single-nucleotide polymorphisms to interrogate this region and identified a strongly associated haplotype in two independent data sets. DNA resequencing of the complement factor H gene within this haplotype revealed a common coding variant, Y402H, that significantly increases the risk for AMD with odds ratios between 2.45 and 5.57. This common variant likely explains approximately 43% of AMD in older adults.

Mutations in the pleckstrin homology domain of dynamin 2 cause dominant intermediate Charcot-Marie-Tooth disease.

Charcot-Marie-Tooth (CMT) disease is a clinically and genetically heterogeneous group of peripheral neuropathies. Different chromosomal loci have been linked with three autosomal dominant, 'intermediate' types of CMT: DI-CMTA, DI-CMTB and DI-CMTC. We refined the locus associated with DI-CMTB on chromosome 19p12-13.2 to 4.2 Mb in three unrelated families with CMT originating from Australia, Belgium and North America. After screening candidate genes, we identified unique mutations in dynamin 2 (DNM2) in all families. DNM2 belongs to the family of large GTPases and is part of the cellular fusion-fission apparatus. In transiently transfected cell lines, mutations of DNM2 substantially diminish binding of DNM2 to membranes by altering the conformation of the beta3/beta4 loop of the pleckstrin homology domain. Additionally, in the Australian and Belgian pedigrees, which carry two different mutations affecting the same amino acid, Lys558, CMT cosegregated with neutropenia, which has not previously been associated with CMT neuropathies.

Targeted disruption of Drosophila Roc1b reveals functional differences in the Roc subunit of Cullin-dependent E3 ubiquitin ligases.

Cullin-dependent ubiquitin ligases regulate a variety of cellular and developmental processes by recruiting specific proteins for ubiquitin-mediated degradation. Cullin proteins form a scaffold for two functional modules: a catalytic module comprised of a small RING domain protein Roc1/Rbx1 and a ubiquitin-conjugating enzyme (E2), and a substrate recruitment module containing one or more proteins that bind to and bring the substrate in proximity to the catalytic module. Here, we present evidence that the three Drosophila Roc proteins are not functionally equivalent. Mutation of Roc1a causes lethality that cannot be rescued by expression of Roc1b or Roc2 by using the Roc1a promoter. Roc1a mutant cells hyperaccumulate Cubitus interruptus, a transcription factor that mediates Hedgehog signaling. This phenotype is not rescued by expression of Roc2 and only partially by expression of Roc1b. Targeted disruption of Roc1b causes male sterility that is partially rescued by expression of Roc1a by using the Roc1b promoter, but not by similar expression of Roc2. These data indicate that Roc proteins play nonredundant roles during development. Coimmunoprecipitation followed by Western or mass spectrometric analysis indicate that the three Roc proteins preferentially bind certain Cullins, providing a possible explanation for the distinct biological activities of each Drosophila Roc/Rbx.

Fibroblast growth factor 20 polymorphisms and haplotypes strongly influence risk of Parkinson disease.

The pathogenic process responsible for the loss of dopaminergic neurons within the substantia nigra of patients with Parkinson disease (PD) is poorly understood. Current research supports the involvement of fibroblast growth factor (FGF20) in the survival of dopaminergic cells. FGF20 is a neurotrophic factor that is preferentially expressed within the substantia nigra of rat brain. The human homologue has been mapped to 8p21.3-8p22, which is within an area of PD linkage revealed through our published genomic screen. To test whether FGF20 influences risk of PD, we genotyped five single-nucleotide polymorphisms (SNPs) lying within the FGF20 gene, in a large family study. We analyzed our sample (644 families) through use of the pedigree disequilibrium test (PDT), the genotype PDT, the multilocus-genotype PDT, and the family-based association test to assess association between risk of PD and alleles, genotypes, multilocus genotypes, and haplotypes. We discovered a highly significant association of PD with one intronic SNP, rs1989754 (P=.0006), and two SNPs, rs1721100 (P=.02) and ss20399075 (P=.0008), located in the 3' regulatory region in our overall sample. Furthermore, we detected a haplotype (A-G-C-C-T) that is positively associated with risk of PD (P=.0003), whereas a second haplotype (A-G-G-G-C) was found to be negatively associated with risk of PD (P=.0009). Our results strongly support FGF20 as a risk factor for PD.

Genomic convergence: identifying candidate genes for Parkinson's disease by combining serial analysis of gene expression and genetic linkage.

We present a multifactorial, multistep approach called genomic convergence that combines gene expression with genomic linkage analysis to identify and prioritize candidate susceptibility genes for Parkinson's disease (PD). To initiate this process, we used serial analysis of gene expression (SAGE) to identify genes expressed in two normal substantia nigras (SN) and adjacent midbrain tissue. This identified over 3700 transcripts, including the three most abundant SAGE tags, which did not correspond to any known genes or ESTs. We developed high-throughput bioinformatics methods to map the genes corresponding to these tags and identified 402 SN genes that lay within five large genomic linkage regions, previously identified in 174 multiplex PD families. These genes represent excellent candidates for PD susceptibility alleles and further genomic convergence and analyses.

Drosophila Roc1a encodes a RING-H2 protein with a unique function in processing the Hh signal transducer Ci by the SCF E3 ubiquitin ligase.

Substrate specificity of SCF E3 ubiquitin ligases is thought to be determined by the F box protein subunit. Another component of SCF complexes is provided by members of the Roc1/Rbx1/Hrt1 gene family, which encode RING-H2 proteins. Drosophila contains three members of this gene family. We show that Roc1a mutant cells fail to proliferate. Further, while the F box protein Slimb is required for Cubitus interruptus (Ci) and Armadillo/beta-catenin (Arm) proteolysis, Roc1a mutant cells hyperaccumulate Ci but not Arm. This suggests that Slimb and Roc1a function in the same SCF complex to target Ci but that a different RING-H2 protein acts with Slimb to target Arm. Consequently, the identity of the Roc subunit may contribute to the selection of substrates by metazoan SCF complexes.