The personal genome project.
Explore the source record for details and available documents.
Ethics Committees, Research↗
Biomedical subjects
G M Church
Publications and source records attributed to G M Church.
Exploring the DNA-binding specificities of zinc fingers with DNA microarrays.
A key step in the regulation of networks that control gene expression is the sequence-specific binding of transcription factors to their DNA recognition sites. A more complete understanding of these DNA-protein interactions will permit a more comprehensive and quantitative mapping of the regulatory pathways within cells, as well as a deeper understanding of the potential functions of individual genes regulated by newly identified DNA-binding sites. Here we describe a DNA microarray-based method to characterize sequence-specific DNA recognition by zinc-finger proteins. A phage display library, prepared by randomizing critical amino acid residues in the second of three fingers of the mouse Zif268 domain, provided a rich source of zinc-finger proteins with variant DNA-binding specificities. Microarrays containing all possible 3-bp binding sites for the variable zinc fingers permitted the quantitation of the binding site preferences of the entire library, pools of zinc fingers corresponding to different rounds of selection from this library, as well as individual Zif268 variants that were isolated from the library by using specific DNA sequences. The results demonstrate the feasibility of using DNA microarrays for genome-wide identification of putative transcription factor-binding sites.
Amino Acid Sequence↗
Computational comparison of two draft sequences of the human genome.
We are in the enviable position of having two distinct drafts of the human genome sequence. Although gaps, errors, redundancy and incomplete annotation mean that individually each falls short of the ideal, many of these problems can be assessed by comparison. Here we present some comparative analyses of these drafts. We look at a number of features of the sequences, including sequence gaps, continuity, consistency between the two sequences and patterns of DNA-binding protein motifs.
Algorithms↗
Selection analyses of insertional mutants using subgenic-resolution arrays.
We describe a method of genome-wide analysis of quantitative growth phenotypes using insertional mutagenesis and DNA microarrays. We applied the method to assess the fitness contributions of Escherichia coli gene domains under specific growth conditions. A transposon library was subjected to competitive growth selection in Luria-Bertani (LB) and in glucose minimal media. Transposon-containing genomic DNA fragments from the selected libraries were compared with the initial unselected transposon insertion library on DNA microarrays to identify insertions that affect fitness. Genes involved in the biosynthesis of nutrients not provided in the growth medium were found to be significantly enriched in the set of genes containing negatively selected insertions. The data also identify fitness contributions of several uncharacterized genes, including putative transcriptional regulators and enzymes. The applicability of this high-resolution array selection in other species is discussed.
Culture Media↗
Identifying regulatory networks by combinatorial analysis of promoter elements.
Several computational methods based on microarray data are currently used to study genome-wide transcriptional regulation. Few studies, however, address the combinatorial nature of transcription, a well-established phenomenon in eukaryotes. Here we describe a new approach using microarray data to uncover novel functional motif combinations in the promoters of Saccharomyces cerevisiae. In addition to identifying novel motif combinations that affect expression patterns during the cell cycle, sporulation and various stress responses, we observed regulatory cross-talk among several of these processes. We have also generated motif-association maps that provide a global view of transcription networks. The maps are highly connected, suggesting that a small number of transcription factors are responsible for a complex set of expression patterns in diverse conditions. This approach may be useful for modeling transcriptional regulatory networks in more complex eukaryotes.
Cell Cycle↗
Correlation between transcriptome and interactome mapping data from Saccharomyces cerevisiae.
Genomic and proteomic approaches can provide hypotheses concerning function for the large number of genes predicted from genome sequences. Because of the artificial nature of the assays, however, the information from these high-throughput approaches should be considered with caution. Although it is possible that more meaningful hypotheses could be formulated by integrating the data from various functional genomic and proteomic projects, it has yet to be seen to what extent the data can be correlated and how such integration can be achieved. We developed a 'transcriptome-interactome correlation mapping' strategy to compare the interactions between proteins encoded by genes that belong to common expression-profiling clusters with those between proteins encoded by genes that belong to different clusters. Using this strategy with currently available data sets for Saccharomyces cerevisiae, we provide the first global evidence that genes with similar expression profiles are more likely to encode interacting proteins. We show how this correlation between transcriptome and interactome data can be used to improve the quality of hypotheses based on the information from both approaches. The strategy described here may help to integrate other functional genomic and proteomic data, both in yeast and in higher organisms.
Fungal Proteins↗
A dynamic programming approach to de novo peptide sequencing via tandem mass spectrometry.
Tandem mass spectrometry fragments a large number of molecules of the same peptide sequence into charged molecules of prefix and suffix peptide subsequences and then measures mass/charge ratios of these ions. The de novo peptide sequencing problem is to reconstruct the peptide sequence from a given tandem mass spectral data of k ions. By implicitly transforming the spectral data into an NC-spectrum graph G (V, E) where /V/ = 2k + 2, we can solve this problem in O(/V//E/) time and O(/V/2) space using dynamic programming. For an ideal noise-free spectrum with only b- and y-ions, we improve the algorithm to O(/V/ + /E/) time and O(/V/) space. Our approach can be further used to discover a modified amino acid in O(/V//E/) time. The algorithms have been implemented and tested on experimental data.
Algorithms↗
Algorithms for identifying protein cross-links via tandem mass spectrometry.
Cross-linking technology combined with tandem mass spectrometry (MS-MS) is a powerful method that provides a rapid solution to the discovery of protein-protein interactions and protein structures. We studied the problem of detecting cross-linked peptides and cross-linked amino acids from tandem mass spectral data. Our method consists of two steps: the first step finds two protein subsequences whose mass sum equals a given mass measured from the mass spectrometry; and the second step finds the best cross-linked amino acids in these two peptide sequences that are optimally correlated to a given tandem mass spectrum. We designed fast and space-efficient algorithms for these two steps and implemented and tested them on experimental data of cross-linked hemoglobin proteins. An interchain cross-link between two beta subunits was found in two tandem mass spectra. The length of the cross-linker (7.7 A) is very close to the actual distance (8.18 A) obtained from the molecular structure in PDB.
Algorithms↗
Dynamic simulation of the human red blood cell metabolic network.
We have developed a Mathematica application package to perform dynamic simulations of the red blood cell (RBC) metabolic network. The package relies on, and integrates, many years of mathematical modeling and biochemical work on red blood cell metabolism. The extensive data regarding the red blood cell metabolic network and the previous kinetic analysis of all the individual components makes the human RBC an ideal 'model' system for mathematical metabolic models. The Mathematica package can be used to understand the dynamics and regulatory characteristics of the red blood cell.
Computer Simulation↗
Aligning gene expression time series with time warping algorithms.
UNLABELLED: motivation: Increasingly, biological processes are being studied through time series of RNA expression data collected for large numbers of genes. Because common processes may unfold at varying rates in different experiments or individuals, methods are needed that will allow corresponding expression states in different time series to be mapped to one another. RESULTS: We present implementations of time warping algorithms applicable to RNA and protein expression data and demonstrate their application to published yeast RNA expression time series. Programs executing two warping algorithms are described, a simple warping algorithm and an interpolative algorithm, along with programs that generate graphics that visually present alignment information. We show time warping to be superior to simple clustering at mapping corresponding time states. We document the impact of statistical measurement noise and sample size on the quality of time alignments, and present issues related to statistical assessment of alignment quality through alignment scores. We also discuss directions for algorithm improvement including development of multiple time series alignments and possible applications to causality searches and non-temporal processes ('concentration warping').
Algorithms↗
Sequence variability and candidate gene analysis in complex disease: association of mu opioid receptor gene variation with substance dependence.
To analyze candidate genes and establish complex genotype-phenotype relationships against a background of high natural genome sequence variability, we have developed approaches to (i) compare candidate gene sequence information in multiple individuals; (ii) predict haplotypes from numerous variants; and (iii) classify haplotypes and identify specific sequence variants, or combinations of variants (pattern), associated with the phenotype. Using the human mu opioid receptor gene (OPRM1) as a model system, we have combined these approaches to test a potential role of OPRM1 in substance (heroin/cocaine) dependence. All known functionally relevant regions of this prime candidate gene were analyzed by multiplex sequence comparison in 250 cases and controls; 43 variants were identified and 52 different haplotypes predicted in the subgroup of 172 African-Americans. These haplotypes were classified by similarity clustering into two functionally related categories, one of which was significantly more frequent in substance-dependent individuals. Common to this category was a characteristic pattern of sequence variants [-1793T-->A, -1699Tins, -1320A-->G, -111C-->T, +17C-->T (A6V)], which was associated with substance dependence. This study provides an example of approaches that have been successfully applied to the establishment of complex genotype-phenotype relationships in the presence of abundant DNA sequence variation.
Adult↗
Predicting regulons and their cis-regulatory motifs by comparative genomics.
We have combined and compared three techniques for predicting functional interactions based on comparative genomics (methods based on conserved operons, protein fusions and correlated evolution) and optimized these methods to predict coregulated sets of genes in 24 complete genomes, including Saccharomyces cerevisiae, Caenorhabditis elegans and 22 prokaryotes. The method based on conserved operons was the most useful for this purpose. Upstream regions of the genes comprising these predicted regulons were then used to search for regulatory motifs in 22 prokaryotic genomes using the motif-discovery program AlignACE. Many significant upstream motifs, including five known Escherichia coli regulatory motifs, were identified in this manner. The presence of a significant regulatory motif was used to refine the members of the predicted regulons to generate a final set of predicted regulons that share significant regulatory elements.
Animals↗
The Isw2 chromatin remodeling complex represses early meiotic genes upon recruitment by Ume6p.
The ISWI class of chromatin remodeling factors exhibits potent chromatin remodeling activities in vitro. However, the in vivo functions of this class of factors are unknown at a molecular level. We have found that S. cerevisiae Isw2 complex represses transcription of early meiotic genes during mitotic growth in a parallel pathway to Rpd3-Sin3 histone deacetylase complex. This repressor function of lsw2 complex is largely dependent upon Ume6p, which recruits the complex to target genes. Nuclease digestion analyses revealed that lsw2 complex establishes nuclease-inaccessible chromatin structure near the Ume6p binding site in vivo. Based on these findings, we propose a model for the mechanism of transcriptional repression by two distinct chromatin remodeling complexes.
Binding Sites↗
Predicting ligand-binding function in families of bacterial receptors.
The three-dimensional fold of a new protein sequence can often be inferred directly from sequence homology to a protein of known structure. The function of a new protein sequence is more difficult to predict, however, since homologues can have different molecular and cellular functions. To develop and automate computational methods for determining molecular function, we have analyzed ligand-binding specificity in two related families of binding proteins. One of these families includes Escherichia coli lactose repressor and ribose-binding protein, and the other includes E. coli sulfate- and phosphate-binding proteins. These proteins have similar folds but varying specificity, binding many different small molecules, including mono- and disaccharides, purines, oxyanions, ferric iron, and polyamines. Starting from template structural alignments, alignments of over 90 sequences per family were generated by iterative database searches with hidden Markov models. Phylogenetic trees were made of full-length sequences and of subsets of residues lining the binding cleft, to determine whether subbranches of the trees correlate with ligand-binding preference. Automated analyses of residues in the binding pocket were also used to predict ligand-binding function for many uncharacterized database sequences and to identify specific side chain-ligand contacts in proteins without solved structures. Our results demonstrate the utility of anchoring functional annotation within a protein family context.
Amino Acid Sequence↗
Microarray analysis of the transcriptional network controlled by the photoreceptor homeobox gene Crx.
BACKGROUND: Terminal differentiation of many cell types is controlled and maintained by tissue- or cell-specific transcription factors. Little is known, however, of the transcriptional networks controlled by such factors and how they regulate differentiation. The paired-type homeobox transcription factor, Crx, has a pivotal role in the terminal differentiation of vertebrate photoreceptors. Mutations in the human CRX gene result in either congenital blindness or photoreceptor degeneration and targeted mutation of the mouse Crx results in failure of development of the light-detecting outer segment of photoreceptors. RESULTS: We have characterized the transcriptional network controlled by Crx by microarray analysis of gene expression in developing retinal tissue from Crx(+/+) and Crx(-/-) mice. These data were combined with analyses of gene expression in developing and adult retina, as well as adult brain. The most abundant elements of this network are ten photoreceptor-specific or -enriched genes, including six phototransduction genes. All of the available 5' regulatory regions of the putative Crx targets contain a novel motif that is composed of a head-to-tail arrangement of two Crx-binding-element-like sequences. Analysis of the 5' regions of a set of mouse and human genes suggests that this motif is specific to Crx targets. CONCLUSIONS: This study demonstrates that cDNA microarrays can be successfully used to define the transcriptional networks controlled by transcription factors in vertebrate tissue in vivo.
Animals↗
Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae.
AlignACE is a Gibbs sampling algorithm for identifying motifs that are over-represented in a set of DNA sequences. When used to search upstream of apparently coregulated genes, AlignACE finds motifs that often correspond to the DNA binding preferences of transcription factors. We previously used AlignACE to analyze whole genome mRNA expression data. Here, we present a more detailed study of its effectiveness as applied to a variety of groups of genes in the Saccharomyces cerevisiae genome. Published functional catalogs of genes and sets of genes grouped by common name provided 248 groups, resulting in 3311 motifs. In conjunction with this analysis, we present measures for gauging the tendency of a motif to target a given set of genes relative to all other genes in the genome and for gauging the degree to which a motif is preferentially located in a certain distance range upstream of translational start sites. We demonstrate improved methods for comparing and clustering sequence motifs. Many previously identified cis-regulatory elements were found. We also describe previously unidentified motifs, one of which has been verified by experiments in our laboratory. An extensive set of AlignACE runs on randomly selected sets of genes and on sets of genes whose upstream regions contain known transcription factor binding sites serve as controls.
Algorithms↗
A computational analysis of whole-genome expression data reveals chromosomal domains of gene expression.
Chromosome correlation maps display correlations between the expression patterns of genes on the same chromosome. Using these maps, we show here that adjacent pairs of genes, as well as nearby non-adjacent pairs of genes, show correlated expression independent of their orientation. We present specific examples of adjacent pairs with highly correlated expression patterns, in which the promoter of only one of the two genes contains an upstream activating sequence (UAS) known to be associated with that expression pattern. Finally, we show that genes with similar functions tend to occur in adjacent positions along the chromosomes. Our results suggest that, in certain chromosomal expression domains, an UAS can affect the transcription of genes that are not immediately downstream from it.
Cell Cycle↗
RNA expression analysis using a 30 base pair resolution Escherichia coli genome array.
We have developed a high-resolution "genome array" for the study of gene expression and regulation in Escherichia coli. This array contains on average one 25-mer oligonucleotide probe per 30 base pairs over the entire genome, with one every 6 bases for the intergenic regions and every 60 bases for the 4,290 open reading frames (ORFs). Twofold concentration differences can be detected at levels as low as 0.2 messenger RNA (mRNA) copies per cell, and differences can be seen over a dynamic range of three orders of magnitude. In rich medium we detected transcripts for 97% and 87% of the ORFs in stationary and log phases, respectively. We found that 1, 529 transcripts were differentially expressed under these conditions. As expected, genes involved in translation were expressed at higher levels in log phase, whereas many genes known to be involved in the starvation response were expressed at higher levels in stationary phase. Many previously unrecognized growth phase-regulated genes were identified, such as a putative receptor (b0836) and a 30S ribosomal protein subunit (S22), both of which are highly upregulated in stationary phase. Transcription of between 3,000 and 4,000 predicted ORFs was observed from the antisense strand, indicating that most of the genome is transcribed at a detectable level. Examples are also presented for high-resolution array analysis of transcript start and stop sites and RNA secondary structure.
Escherichia coli↗