p53 responsive elements in human retrotransposons.

Long interspersed nuclear elements-1 (L1s) are highly repetitive DNA elements that are capable of altering the human genome through retrotransposition. To protect against L1 retroposition, the cell downregulates the expression of L1 proteins by various mechanisms, including high-density cytosine methylation of L1 promoters and DICER-dependent destruction of L1 mRNAs. In this report, a large number of p53 responsive elements, or p53 DNA binding sites, were detected in L1 elements within the human genome. At least some of these p53 responsive elements are functional and can act to increase the levels of L1 mRNA expression. The p53 protein can directly bind to a short 15-nucleotide sequence within the L1 promoter. This p53 responsive element within L1 is a recent addition to evolution, appearing approximately 20 million years ago. This suggests an interplay between L1 elements, which have a rich history of causing changes in the genome, and the p53 protein, the function of which is to protect against genomic changes. To understand these observations, a model is proposed in which the increased expression of L1 mRNAs by p53 actually increases, rather than decreases, the genomic stability through amplification of p53-dependent processes for genomic protection.

Genotypic association analysis using discordant-relative-pairs.

In practice, family-based design has been widely used in disease-gene association analysis. The major advantage of such design is that it is not subject to spurious association due to population structure such as population stratification (PS) and admixture. A disadvantage is that parental genotypes are hard to obtain if the disease is late onset for which a discordant-relative-pair design is useful. Designs of such kind include full-sib-pair, half-sib-pair, first-cousin-pair, and so on. The closer the relatedness of the pair, the less possible that they are subject to population stratification. On the other hand, the association test using close relative-pairs may be less powerful due to over-matching. Trade-off between these two factors (population structure and over-matching) is the major concern of this study. Some tests, namely McNemar's test, matched Cochran-Armitage trend tests (CATTs), matched maximum efficient robust test (MERT), and Bhapkar's test, are used for testing disease-gene association based on relative-pair designs. These tests are shown to be valid in the presence of PS but not admixture. Numerical studies show that the McNemar's test, additive CATT, MERT, and Bhapkar's test are robust in power, but none of them is uniformly more powerful than the others. In most simulations, the power of any of the tests increases as the pair is more distant. The proposed methods are applied to two real examples.

p53 tumor suppressor protein regulates the levels of huntingtin gene expression.

The p53 protein is a transcription factor that integrates various cellular stress signals. The accumulation of the mutant huntingtin protein with an expanded polyglutamine tract plays a central role in the pathology of human Huntington's disease. We found that the huntingtin gene contains multiple putative p53-responsive elements and p53 binds to these elements both in vivo and in vitro. p53 activation in cultured human cells, either by a temperature-sensitive mutant p53 protein or by gamma-irradiation (gamma-irradiation), increases huntingtin mRNA and protein expression. Similarly, murine huntingtin also contains multiple putative p53-responsive elements and its expression is induced by p53 activation in cultured cells. Moreover, gamma-irradiation, which activates p53, increases huntingtin gene expression in the striatum and cortex of mouse brain, the major pathological sites for Huntington's disease, in p53+/+ but not the isogenic p53-/- mice. These results demonstrate that p53 protein can regulate huntingtin expression at transcriptional level, and suggest that a p53 stress response could be a modulator of the process of Huntington's disease.

Laryngeal muscle fibre types.

The internal laryngeal muscles have evolved to subserve the highly specialized functions of airways protection, respiration, and phonation. Their contractile properties, histochemistry, biochemical properties, myosin heavy chain (MyHC) expression and their regulation by nerves and hormones are reviewed and compared with limb muscle fibres. Cricothyroid, the vocal cord tensor, is limb-like in MyHC composition and fibre type properties, while the vocal fold abductor and adductors are allotypically different, with capacity for expressing an isoform of MyHC that is kinetically faster than the fastest limb MyHC. In rats and rabbits the faster isoform is the extraocular (EO) MyHC, while in carnivores, it is the IIB MyHC. These adaptations enable the abductor and adductor muscles to remain always faster than the cricothyroid as the latter changes in speed during evolution to match changing metabolic and respiratory rates in relation to scaling with body mass. Such phylogenetic plasticity is vital to the airways protection and respiratory functions of these muscles. The posterior cricoarythenoid, the abductor muscle, is tonically driven during expiration, and consequently has a slower fibre type profile than the principal adductor, the thyroarythenoid. The human thyroarythenoid appears not to express EO or IIB MyHC significantly, but is unique in expressing the slow-tonic MyHC. The concepts of allotype and phylogenetic plasticity help to explain differences in fibre type between limb and laryngeal muscles and between homologous laryngeal muscles in different species. Laryngeal muscle fibres exhibit physiological plasticity as do limb muscles, being subject to neural and hormonal modulation.

Multi-locus interactions predict risk for post-PTCA restenosis: an approach to the genetic analysis of common complex disease.

The complexity of recognizing the potential contribution of a number of possible predictors of complex disorders is increasingly challenging with the application of large-scale single nucleotide polymorphism (SNP) typing. In the search for putative genetic factors predisposing to coronary artery restenosis following balloon angioplasty, we determined genotypes for 94 SNPs representing 62 candidate genes, in a prospectively assembled cohort of 342 cases and 437 controls. Using a customized coupled-logistic regression procedure accounting for both additive and interactive effects, we identified seven SNPs in seven genes that, together, showed a statistically significant association with restenosis incidence (P <0.0001), accounting for 11.6% of overall variance observed. Among them are candidate genes for cardiovascular pathophysiology (apolipoprotein-species and NOS), inflammatory response (TNF receptor and CD14), and cell-cycle control (p53 and p53-associated protein). Our results emphasize the need to account for complex multi-gene influences and interactions when assessing the molecular pathology of multifactorial medical entities.

The p53MH algorithm and its application in detecting p53-responsive genes.

A computer algorithm, p53MH, was developed, which identifies putative p53 transcription factor DNA-binding sites on a genomewide scale with high power and versatility. With the sequences from the human and mouse genomes, putative p53 DNA-binding elements were identified in a scan of 2,583 human genes and 1,713 mouse orthologs based on the experimental data of el-Deiry et al. [el-Deiry, W. S., Kern, S. E., Pietenpol, J. A., Kinzler, K. W. & Vogelstein, B. (1992) Nat. Genet. 1, 45-49] and Funk et al. [Funk, W. D., Pak, D. T., Karas, R. H., Wright, W. E. & Shay, J. W. (1992) Mol. Cell. Biol. 12, 2866-2871] (http://linkage.rockefeller.edu/p53). The p53 DNA-binding motif consists of a 10-bp palindrome and most commonly a second related palindrome linked by a spacer region. By scanning from the 5' to 3' end of each gene with an additional 10-kb nucleotide sequence appended at each end (most regulatory DNA elements characterized in the literature are in these regions), p53MH computes the binding likelihood for each site under a discrete discriminant model and then outputs ordered scores, corresponding site positions, sequences, and related information. About 300 genes receiving scores greater than a theoretical cut-off value were identified as potential p53 targets. Semiquantitative reverse transcription-PCR experiments were performed in 2 cell lines on 16 genes that were previously unknown regarding their functional relationship to p53 and were found to have high scores in either proximal promoter or possible distal enhancer regions. Ten (approximately 63%) of these genes responded to the presence of p53.

Probing the machinery of intracellular trafficking with the atomic force microscope.

Atomic force microscopy has emerged as a powerful tool for characterizing single biological macromolecules, macromolecular assemblies, and whole cells in aqueous buffer, in real time, and at molecular-scale spatial and force resolution. Many of the central elements of intracellular transport are tens to hundreds of nanometers in size and highly dynamic. Thus, atomic force microscopy provides a valuable means of addressing questions of structure and mechanism in intracellular transport. We begin this review of recent efforts to apply atomic force microscopy to problems in intracellular transport by discussing the technical principles behind atomic force microscopy. We then turn to three specific areas in which atomic force microscopy has been applied to problems with direct implications for intracellular trafficking: cytoskeletal structure and dynamics, vesicular transport, and receptor-ligand interactions. In each case, we discuss studies which use both intact cellular elements and reconstituted models. While many technical challenges remain, these studies point to several areas where atomic force microscopy can be used to provide valuable insight into intracellular transport at exquisite spatial and energetic resolution.

Trimming, weighting, and grouping SNPs in human case-control association studies.

The search for genes underlying complex traits has been difficult and often disappointing. The main reason for these difficulties is that several genes, each with rather small effect, might be interacting to produce the trait. Therefore, we must search the whole genome for a good chance to find these genes. Doing this with tens of thousands of SNP markers, however, greatly increases the overall probability of false-positive results, and current methods limiting such error probabilities to acceptable levels tend to reduce the power of detecting weak genes. Investigating large numbers of SNPs inevitably introduces errors (e.g., in genotyping), which will distort analysis results. Here we propose a simple strategy that circumvents many of these problems. We develop a set-association method to blend relevant sources of information such as allelic association and Hardy-Weinberg disequilibrium. Information is combined over multiple markers and genes in the genome, quality control is improved by trimming, and an appropriate testing strategy limits the overall false-positive rate. In contrast to other available methods, our method to detect association to sets of SNP markers in different genes in a real data application has shown remarkable success.

Predicting properties of intrinsically unstructured proteins.

There is increasing evidence that intrinsically unstructured proteins or protein domains have important biological functions. These types of proteins may be productively analyzed using polymer theory developed to predict global physical properties of polymers. In these theories molecular detail is "coarse grained" out of the models, and replaced with a small number of parameters that characterize the polymer. This reduction in complexity allows extremely large systems to be studied. In the case of simulations, the time scales accessible also increase significantly. Here we discuss the application of polymer theory to unstructured proteins, and consider how to classify proteins within a polymer framework. We then review polymer theory that is relevant to predicting functionally important properties, such as radius of gyration, height of a polymer brush and force required to compress a polymer brush.

AFM force measurements on microtubule-associated proteins: the projection domain exerts a long-range repulsive force.

Microtubule-associated proteins (MAPs) are thought to control spacing between microtubules. We propose that the projection domain is largely unstructured and exerts a long-range repulsive force that is predominantly entropic in origin, providing a physical mechanism for maintaining spacing. To test this hypothesis, we developed an experimental system where MAPs are electrostatically end-attached to a flat surface, such that the projection domains extend away from the surface. Atomic force microscopy force measurements on this system show that projection domains exert a long-range (>100 nm) repulsive force. This force depends on the ionic strength of the solution in a way that is consistent with a polyelectrolyte polymer brush.

A genomewide screen for autism susceptibility loci.

We report the analysis of 335 microsatellite markers genotyped in 110 multiplex families with autism. All families include at least two "affected" siblings, at least one of whom has autism; the remaining affected sibs carry diagnoses of either Asperger syndrome or pervasive developmental disorder. Affected sib-pair analysis yielded multipoint maximum LOD scores (MLS) that reach the accepted threshold for suggestive linkage on chromosomes 5, X, and 19. Nominal evidence for linkage (point-wise P<.05) was obtained on chromosomes 2, 3, 4, 8, 10, 11, 12, 15, 16, 18, and 20, and secondary loci were found on chromosomes 5 and 19. Analysis of families sharing alleles at the putative X chromosomal linked locus and one or more other putative linked loci produced an MLS of 3.56 for the DXS470-D19S174 marker combination. In an effort to increase power to detect linkage, scan statistics were used to evaluate the significance of peak LOD scores based on statistical evidence at adjacent marker loci. This analysis yielded impressive evidence for linkage to autism and autism-spectrum disorders with significant genomewide P values <.05 for markers on chromosomes 5 and 8 and with suggestive linkage evidence for a marker on chromosome 19.

Electrophoretic and immunochemical evidence showing that marsupial limb muscles express the same fast and slow myosin heavy chains as eutherians.

Limb muscles of eutherian (placental) mammals express a slow and three fast isoforms of myosin heavy chain (MyHC), but little is known about marsupial MyHCs. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) of limb MyHCs from seven marsupial species, spanning two orders, revealed four components, each of which specifically cross-reacted in Western blots with a monoclonal antibody (mAb) against a corresponding eutherian MyHC. For all seven species, the relative mobility of the band identified by each mAb matched that in the rat, suggesting that the four are homologous to eutherian slow, 2B, 2X and 2A MyHCs, respectively, in the order of decreasing mobility. Immunohistochemical analysis of fast marsupial limb muscles identitied four different fiber populations whose relative fiber size spectra (IIA

Statistical multilocus methods for disequilibrium analysis in complex traits.

Hundreds of thousands of SNP markers are being generated with the purpose of carrying out case-control association studies for complex traits, which are thought to be due to multiple underlying susceptibility genes. The number of markers is typically much larger than the number of observations so that joint analysis of marker genotypes and their interactions is not feasible. We discuss a two-stage approach to first select a small subset of markers and then model the effects of the selected markers on disease. Examples of two procedures for marker selection are given with subsequent modeling of main and interaction effects. The approaches are applied to a data set with 89 SNPs in lieu of a genome screen with many more markers.

Ascertainment and anticipation in family studies.

Many human diseases show anticipation; that is, disease occurs earlier (or with greater severity) in successive generations. In a computer simulation, we assessed the degree of anticipation that one would expect to see in two-generation breast cancer families. Under reasonable assumed distributions for age at cancer onset, number of children, and mortality, we find a consistent earlier mean age at diagnosis in daughters than in mothers, but the same mean age at diagnosis in affected aunts and nieces. We compare these results with published pedigree data for familial breast cancer that show substantial anticipation in affected daughters compared to their mothers. We find that at least some anticipation is expected in human disease families even when the disease is stable and families are ascertained without obvious sampling bias. We further demonstrate that such anticipation is reduced when comparing affected children to the parents' affected siblings.

Two approaches for consolidating results from genome scans of complex traits: selection methods and scan statistics.

This work has two purposes: (i) empirically selecting levels of significance that maximize the fraction of markers close to a gene (hit rate) when performing linkage analyses of simulated data and (ii) evaluating the utility of a previously reported scan statistic on the same data. Genotype data were simulated from a trait model of seven susceptibility genes. For purpose (i), five statistics were evaluated on all marker loci in fifty replicates; two-point lod and heterogeneity lod scores maximized over dominance (mlod, mhlod), a multi-allelic TDT test, an affected sib-pair test (ASP), and a model-free test on all sib-pairs (ALL_SIBS). Within each replicate the fraction of markers (hit rate) significant at specified levels of significance and also (a) within fifty markers of, or (b) on the same chromosome as a major gene was calculated. For purpose (ii), scan statistics of length 15 were calculated for each chromosome and their empirical significance levels estimated on the basis of 500 replicates generated under no linkage. The scan statistic was applied to the mhlod scores from one replicate (Replicate 5). Empirical p-values for the scan statistic were determined by computing mhlod scores on 500 replicates of simulated null data. For purpose (i), significance levels between 0.001 and 0.01 had the greatest hit rate for all five methods and both criteria. For criterion (a) at the 0.001 level of significance, both mlod and mhlod displayed the highest hit rates, approximately 0.4 for each. For criterion (b), all methods but ALL_SIBS and ASP had hit rates ranging between 0.4 and 0.5. For purpose (ii), the scan statistic proved equally or more powerful than the single-locus statistic for two of the seven susceptibility genes while the remaining five genes were not detected.

A train of thoughts on gene mapping.

Complex traits, by definition, are the pheonotypic outcome from multiple interacting genes. The traditional analysis of association studies on complex traits is to test one locus at a time, but a better approach is to analyze all markers simultaneously. We previously proposed a two-stage approach, first selecting the influential markers and then modeling main and interaction effects of these markers. Here we introduce alternative approaches to marker selection and discuss issues regarding analytical tools for disease gene mapping, marker selection, and statistical modeling.

Characterization of DNA condensates by atomic force microscopy.

In recent years, interest in in vitro DNA condensation has been revived by efforts to develop vectors for nonviral gene therapy (1,2). One of the critical elements for successful and versatile delivery of specific genes into targeted cells is that DNA vectors of several kilobase pairs must be compressed and packaged into small particles. Ideally, the size of these particles should be similar to that of viruses, typically less than 100 nm in diameter (3,4). It is well known that several classes of multivalent cations can be used to condense DNA in vitro to form well-defined structures, particularly toroids and rods (5,6). Because of its simplicity, this type of in vitro condensation of DNA has become a widely studied system for investigating mechanisms of DNA condensation.

Cross-bridge kinetics of rabbit single extraocular and limb muscle fibers.

PURPOSE: To gain insights into the functional significance of myosin heavy-chain (MyHC) heterogeneity by comparing the mechanical kinetic properties of single rabbit extraocular muscle (EOM) fibers with those of limb fibers. EOMs are known to contain developmental and EOM-specific MyHCs in addition to those present in limb muscles, and MyHCs profoundly influence muscle mechanics. METHODS: Isometric cross-bridge kinetics were analyzed in Ca(2+)-activated single glycerinated fibers from rabbit EOM and limb fast and slow muscles at 15 degrees C by means of mechanical perturbation analysis. The plots of stiffness and phase against frequency display a characteristic frequency, f(min), at which stiffness is minimum, and phase shift is zero. The value of f(min) is independent of Ca(2+) or force level but reflects the kinetics of cross-bridge cycling. RESULTS: Analysis of 121 limb fast fibers gave f(min) values ranging from 10 to 26 Hz. f(min) for the 10 slow soleus fibers was 0.5 Hz. Analysis of 170 EOM fibers gave f(min) values in the range for fast limb fibers, but in addition yielded f(min) values below (4-9 Hz) and above (27-33 Hz) this range. CONCLUSIONS: The wider range of mechanical kinetic characteristics in EOM fibers compared with limb fibers is likely due to the expression of developmental (low f(min)) and EOM-specific (high f(min)) MyHCs in addition to isoforms present in adult limb muscles. The considerable diversity of functional characteristics in EOM fibers is likely to be important for rotating the eyeball at various speeds during tracking and for executing saccades over a wide range of angles.

The distribution of myosin heavy chain isoforms among rat extraocular muscle fiber types.

PURPOSE: To determine the distribution of myosin heavy chain isoforms in each extraocular muscle (EOM) fiber type. METHODS: Serial sections of adult rat EOMs were stained with isoform-specific monoclonal antibodies against an array of myosin heavy chains. Immunofluorescent antibody staining of whole adult rat EOMs, examined by confocal microscopy, demonstrated the longitudinal variations of isoforms along individual fibers. RESULTS: Each global fiber type reacted predominantly with a single isoform-specific antibody and showed no longitudinal variation. Two major orbital fibers were defined, and both contained multiple myosin heavy chains. Both orbital singly and multiply innervated fibers stained proximal and distal to the neuromuscular junction with antibody to embryonic myosin heavy chain, but this isoform was sharply and completely excluded from the domain of the neuromuscular junction. Orbital singly innervated fibers also contained the EOM-specific isoform at the neuromuscular junction. Orbital multiply innervated fibers did not contain the EOM-specific isoform, but additionally contained a slow isoform along their entire length. CONCLUSIONS: Adult rat EOMs show unique fiber types with arrangements of myosin heavy chain isoforms not seen in other skeletal muscles. Moreover, unique cellular mechanisms must exist to target each isoform to its proper domain along individual orbital fibers.

Scan statistics to scan markers for susceptibility genes.

Scan statistics are applied to combine information on multiple contiguous genetic markers used in a genome screen for susceptibility loci. This information may be, for example, allele sharing proportions for sib pairs or logarithm of odds (lod) scores in general small families. We focus on a dichotomous outcome variable, for example, case and control individuals or affected-affected versus affected-unaffected siblings, and suitable single-marker statistics. A significant scan statistic based on the single-marker statistics represents evidence of the presence of a susceptibility gene. For a given length of the scan statistic, we assess its significance by Monte Carlo permutation tests. Comparing P values for varying lengths of scan statistics, we treat the smallest observed P value as our statistic of interest and determine its overall significance level. We applied this method to a genome screen with autism families. The result was informative and surprising: A susceptibility region was found (genome-wide significance level, P = 0.038), which is missed with conventional approaches.