A second generation human haplotype map of over 3.1 million SNPs

We describe the Phase II HapMap, which characterizes over 3.1 million human single nucleotide polymorphisms (SNPs) genotyped in 270 individuals from four geographically diverse populations and includes 25-35% of common SNP variation in the populations surveyed. The map is estimated to capture untyped common variation with an average maximum r2 of between 0.9 and 0.96 depending on population. We demonstrate that the current generation of commercial genome-wide genotyping products captures common Phase II SNPs with an average maximum r2 of up to 0.8 in African and up to 0.95 in non-African populations, and that potential gains in power in association studies can be obtained through imputation. These data also reveal novel aspects of the structure of linkage disequilibrium. We show that 10-30% of pairs of individuals within a population share at least one region of extended genetic identity arising from recent ancestry and that up to 1% of all common variants are untaggable, primarily because they lie within recombination hotspots. We show that recombination rates vary systematically around genes and between genes of different function. Finally, we demonstrate increased differentiation at non-synonymous, compared to synonymous, SNPs, resulting from systematic differences in the strength or efficacy of natural selection between populations.

The DNA sequence and biological annotation of human chromosome 1

The reference sequence for each human chromosome provides the framework for understanding genome function, variation and evolution. Here we report the finished sequence and biological annotation of human chromosome 1. Chromosome 1 is gene-dense, with 3,141 genes and 991 pseudogenes, and many coding sequences overlap. Rearrangements and mutations of chromosome 1 are prevalent in cancer and many other diseases. Patterns of sequence variation reveal signals of recent selection in specific genes that may contribute to human fitness, and also in regions where no function is evident. Fine-scale recombination occurs in hotspots of varying intensity along the sequence, and is enriched near genes. These and other studies of human biology and disease encoded within chromosome 1 are made possible with the highly accurate annotated sequence, as part of the completed set of chromosome sequences that comprise the reference human genome.

The DNA sequence of the human X chromosome

The human X chromosome has a unique biology that was shaped by its evolution as the sex chromosome shared by males and females. We have determined 99.3% of the euchromatic sequence of the X chromosome. Our analysis illustrates the autosomal origin of the mammalian sex chromosomes, the stepwise process that led to the progressive loss of recombination between X and Y, and the extent of subsequent degradation of the Y chromosome. LINE1 repeat elements cover one-third of the X chromosome, with a distribution that is consistent with their proposed role as way stations in the process of X-chromosome inactivation. We found 1,098 genes in the sequence, of which 99 encode proteins expressed in testis and in various tumour types. A disproportionately high number of mendelian diseases are documented for the X chromosome. Of this number, 168 have been explained by mutations in 113 X-linked genes, which in many cases were characterized with the aid of the DNA sequence.

DNA sequence and analysis of human chromosome 9

Chromosome 9 is highly structurally polymorphic. It contains the largest autosomal block of heterochromatin, which is heteromorphic in 6-8% of humans, whereas pericentric inversions occur in more than 1% of the population. The finished euchromatic sequence of chromosome 9 comprises 109,044,351 base pairs and represents >99.6% of the region. Analysis of the sequence reveals many intra- and interchromosomal duplications, including segmental duplications adjacent to both the centromere and the large heterochromatic block. We have annotated 1,149 genes, including genes implicated in male-to-female sex reversal, cancer and neurodegenerative disease, and 426 pseudogenes. The chromosome contains the largest interferon gene cluster in the human genome. There is also a region of exceptionally high gene and G + C content including genes paralogous to those in the major histocompatibility complex. We have also detected recently duplicated genes that exhibit different rates of sequence divergence, presumably reflecting natural selection.

The DNA sequence and comparative analysis of human chromosome 10

The finished sequence of human chromosome 10 comprises a total of 131,666,441 base pairs. It represents 99.4% of the euchromatic DNA and includes one megabase of heterochromatic sequence within the pericentromeric region of the short and long arm of the chromosome. Sequence annotation revealed 1,357 genes, of which 816 are protein coding, and 430 are pseudogenes. We observed widespread occurrence of overlapping coding genes (either strand) and identified 67 antisense transcripts. Our analysis suggests that both inter- and intrachromosomal segmental duplications have impacted on the gene count on chromosome 10. Multispecies comparative analysis indicated that we can readily annotate the protein-coding genes with current resources. We estimate that over 95% of all coding exons were identified in this study. Assessment of single base changes between the human chromosome 10 and chimpanzee sequence revealed nonsense mutations in only 21 coding genes with respect to the human sequence.

The DNA sequence and analysis of human chromosome 13

Chromosome 13 is the largest acrocentric human chromosome. It carries genes involved in cancer including the breast cancer type 2 (BRCA2) and retinoblastoma (RB1) genes, is frequently rearranged in B-cell chronic lymphocytic leukaemia, and contains the DAOA locus associated with bipolar disorder and schizophrenia. We describe completion and analysis of 95.5 megabases (Mb) of sequence from chromosome 13, which contains 633 genes and 296 pseudogenes. We estimate that more than 95.4% of the protein-coding genes of this chromosome have been identified, on the basis of comparison with other vertebrate genome sequences. Additionally, 105 putative non-coding RNA genes were found. Chromosome 13 has one of the lowest gene densities (6.5 genes per Mb) among human chromosomes, and contains a central region of 38 Mb where the gene density drops to only 3.1 genes per Mb.

Beta-ketoacyl-acyl carrier protein synthase III from pea (Pisum sativum L.): properties, inhibition by a novel thiolactomycin analogue and isolation of a cDNA clone encoding the enzyme

A beta-ketoacyl-acyl carrier protein (ACP) synthase III (KAS III; short-chain condensing enzyme) has been partly purified from pea leaves. The enzyme, which had acetyl-CoA:ACP acyltransferase (ACAT) activity, was resolved from a second, specific, ACAT protein. The KAS III enzyme had a derived molecular mass of 42 kDa (from its cDNA sequence) and operated as a dimer. Its enzymological characteristics were similar to those of two other plant KAS III enzymes except for its inhibition by thiolactomycin. A derivative of thiolactomycin containing a longer (C8 saturated) hydrophobic side-chain (compound 332) was a more effective inhibitor of pea KAS III and showed competitive inhibition towards malonyl-ACP whereas thiolactomycin showed uncompetitive characteristics at high concentrations. This difference may be due to the better fit of compound 332 into a hydrophobic pocket at the active site. A full-length cDNA for the pea KAS III was isolated. This was expressed in Escherichia coli as a fusion protein with glutathione S-transferase in order to facilitate subsequent purification. Demonstrated activity in preparations from E. coli confirmed that the cDNA encoded a KAS III enzyme. Furthermore, the expressed KAS III had ACAT activity, showing that the latter was inherent. The derived amino acid sequence of the pea cDNA showed 81-87% similarity to that for other plant dicotyledon KAS IIIs, somewhat less for Allium porrum (leek, 71%) and for Porphyra spp. (62%), Synechocystis spp. (65%) and various bacteria (42-65%). The pea KAS III exhibited four areas of homology, three of which were around the active-site Cys(123), His(323) and Asn(353). In addition, a stretch of 23 amino acids (residues 207-229 in the pea KAS III) was almost completely conserved in the plant KAS IIIs. Modelling this stretch showed they belonged to a peptide fragment that fitted over the active site and contained segments suggested to be involved in substrate binding and in conformational changes during catalysis, as well as an arginine suggested to participate in the acid-base catalytic mechanism.

The DNA sequence and comparative analysis of human chromosome 20

The finished sequence of human chromosome 20 comprises 59,187,298 base pairs (bp) and represents 99.4% of the euchromatic DNA. A single contig of 26 megabases (Mb) spans the entire short arm, and five contigs separated by gaps totalling 320 kb span the long arm of this metacentric chromosome. An additional 234,339 bp of sequence has been determined within the pericentromeric region of the long arm. We annotated 727 genes and 168 pseudogenes in the sequence. About 64% of these genes have a 5' and a 3' untranslated region and a complete open reading frame. Comparative analysis of the sequence of chromosome 20 to whole-genome shotgun-sequence data of two other vertebrates, the mouse Mus musculus and the puffer fish Tetraodon nigroviridis, provides an independent measure of the efficiency of gene annotation, and indicates that this analysis may account for more than 95% of all coding exons and almost all genes.

A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms

We describe a map of 1.42 million single nucleotide polymorphisms (SNPs) distributed throughout the human genome, providing an average density on available sequence of one SNP every 1.9 kilobases. These SNPs were primarily discovered by two projects: The SNP Consortium and the analysis of clone overlaps by the International Human Genome Sequencing Consortium. The map integrates all publicly available SNPs with described genes and other genomic features. We estimate that 60,000 SNPs fall within exon (coding and untranslated regions), and 85% of exons are within 5 kb of the nearest SNP. Nucleotide diversity varies greatly across the genome, in a manner broadly consistent with a standard population genetic model of human history. This high-density SNP map provides a public resource for defining haplotype variation across the genome, and should help to identify biomedically important genes for diagnosis and therapy.

An SNP map of human chromosome 22

The human genome sequence will provide a reference for measuring DNA sequence variation in human populations. Sequence variants are responsible for the genetic component of individuality, including complex characteristics such as disease susceptibility and drug response. Most sequence variants are single nucleotide polymorphisms (SNPs), where two alternate bases occur at one position. Comparison of any two genomes reveals around 1 SNP per kilobase. A sufficiently dense map of SNPs would allow the detection of sequence variants responsible for particular characteristics on the basis that they are associated with a specific SNP allele. Here we have evaluated large-scale sequencing approaches to obtaining SNPs, and have constructed a map of 2,730 SNPs on human chromosome 22. Most of the SNPs are within 25 kilobases of a transcribed exon, and are valuable for association studies. We have scaled up the process, detecting over 65,000 SNPs in the genome as part of The SNP Consortium programme, which is on target to build a map of 1 SNP every 5 kilobases that is integrated with the human genome sequence and that is freely available in the public domain.

Isolation, characterisation and expression of a cDNA for pea cholinephosphate cytidylyltransferase

In plants, phosphatidylcholine is the major phospholipid in extra-plastid membranes and is synthesised mainly by the CDP-choline pathway. Evidence from studies in animals, as well as in plants, suggests that the intermediate step catalysed by cholinephosphate cytidylyltransferase (CPCT) has a major control in carbon flux to this lipid. We have isolated a full-length CPCT cDNA (designated PCT2) from Pisum sativum cv. Feltham First using an Arabidopsis probe and the polymerase chain reaction (PCR). The deduced amino acid of PCT2 is 48%, 43% and 76% identical to the rat, yeast and Brassica napus amino acid sequences, respectively. Expression of the CPCT protein in Escherichia coli confirmed the activity of the enzyme. Expression of the PCT2 mRNA in pea roots and stems was increased by treatment with 0.1 microM indole-3-acetic acid.

Biological production of semisynthetic opiates using genetically engineered bacteria

Semisynthetic derivatives of morphine and related alkaloids are in widespread clinical use. Due to the complexity of these molecules, however, chemical transformations are difficult to achieve in high yields. We recently identified the powerful analgesic hydromorphone as an intermediate in the metabolism of morphine by Pseudomonas putida M10. Here we describe the construction of recombinant strains of Escherichia coli that express morphine dehydrogenase and morphinone reductase. These strains are capable of efficiently transforming the naturally occurring alkaloids morphine and codeine to hydromorphone and the antitussive hydrocodone, respectively. Our results demonstrate the potential for recombinant DNA technology to provide biological routes for the synthesis of known and novel semisynthetic opiate drugs.

Bacterial morphine dehydrogenase further defines a distinct superfamily of oxidoreductases with diverse functional activities

Pseudomonas putida morphine dehydrogenase is shown to be closely homologous to 18 proteins, defining a superfamily within which morphine dehydrogenase particularly resembles two bacterial, 2,5-dioxo-D-gluconic acid reductases, and two eukaryotic proteins of unknown functions. Relationships within the superfamily are extensive and complex. Residue identities between protein pairs range from 29-90%. Three subgroups are proposed. Nevertheless, on the basis of residue conservations/exchanges it is suggested that the nicotinamide coenzyme binding and substrate reduction occur in all the enzymes by broadly analogous mechanisms, among which some probable differences are identified.

Nucleotide sequence and over-expression of morphine dehydrogenase, a plasmid-encoded gene from Pseudomonas putida M10

Pseudomonas putida M10 was originally isolated from factory waste liquors by selection for growth on morphine. The NADP(+)-dependent morphine dehydrogenase that initiates morphine catabolism is encoded by a large plasmid of 165 kb. Treatment of P. putida M10 with ethidium bromide led to the isolation of a putative plasmid-free strain that was incapable of growth on morphine. The structural gene for morphine dehydrogenase, morA, has been located on the plasmid by oligonucleotide hybridization, by coupled transcription-translation of cloned restriction fragments and by nucleotide sequence analysis and is contained within a 1.7 kb SphI fragment that has been cloned into Escherichia coli. The cloned dehydrogenase enzyme is expressed at high levels in E. coli resulting in a 65-fold increase in morphine dehydrogenase activity in cell-free extracts compared with P. putida M10. Morphine dehydrogenase was rapidly purified to homogeneity, as judged by SDS/PAGE, by a one-step affinity chromatography procedure on Mimetic Orange 3 A6XL. The properties of the purified enzyme were identical with those previously reported for P. putida M10 morphine dehydrogenase. The morA gene was sequenced and the deduced amino acid sequence confirmed by N-terminal amino acid sequencing of the over-expressed protein. The predicted amino acid sequence of morA, deduced from the nucleotide sequence, indicated that morphine dehydrogenase did not belong to the non-metal-requiring short-chain class of dehydrogenases, but was more closely related to the aldo-ketoreductases.

Pea chloroplast genes encoding a 4 kDa polypeptide of photosystem I and a putative enzyme of C1 metabolism

The nucleotide sequence of 3.2 kbp of pea chloroplast DNA located upstream from the petA gene for cytochrome f, and previously reported to contain the gene for a photosystem I polypeptide, has been determined. Three open reading frames of 587, 40 and 157 codons have been identified. Orf40 encodes a highly conserved, hydrophobic, membrane-spanning polypeptide, and is identified as the gene psaI for the 4 kDa subunit of photosystem I. Orf587 is an extended version of the gene zfpA previously identified as encoding a conserved putative zinc-finger protein. The product of orf587 shows extensive homology to an unidentified open reading frame cotranscribed with a gene for folate metabolism in Escherichia coli and local homology to a region of the beta subunit of rat mitochondrial propionyl-CoA carboxylase. It is suggested that the product of orf587 is an enzyme of C1 metabolism and is unlikely to be a regulatory DNA-binding protein. Orf157 potentially encodes an unidentified basic protein, but the protein sequence is not conserved in other plants.

Molecular cloning and structural analysis of the phosphate translocator from pea chloroplasts and its comparison to the spinach phosphate translocator

Using an 5'-AvaII fragment of the spinach (Spinacia oleracea L.) phosphate translocator cDNA as a probe for a hybridization screening of a pea (Pisum sativum L.) cDNA library we have cloned and sequenced a cDNA clone coding for the phosphate translocator precursor protein from pea chloroplasts. The full-length cDNA clone comprises 42 base pairs (bp) at the 5'-non-coding region, a 1206-bp coding region corresponding to a polypeptide of 402 amino-acid residues (relative molecular mass 43 671) and 244 bp at the non-coding 3'-region. Determination of the N-terminal sequence of the phosphate translocator from both pea and spinach chloroplasts revealed that the transit peptides consist of 72 and 80 amino-acid residues, respectively. These transit peptides are different from those of other chloroplastic transit peptides in that they both contain an amphiphilic α-helix which is located either in close proximity to the processing site in pea or at the N-terminus in spinach. The mature proteins from pea and spinach both contain about 87% identical amino-acid residues and about seven putative membrane-spanning α-helices. Some of these α-helices have an amphiphilic character and might serve to form a hydrophilic translocation channel through the membrane. The in-vitro synthesized pea precursor protein is directed to the chloroplast and inserted into the chloroplast envelope membrane.

An open reading frame encoding a putative haem-binding polypeptide is cotranscribed with the pea chloroplast gene for apocytochrome f

The nucleotide sequence of a 1 kbp region of pea chloroplast DNA upstream from the gene petA encoding apocytochrome f has been determined. An open reading frame of 231 codons (ORF231) encoding a putative membrane-spanning polypeptide is separated by 205 bp from the coding region of petA. The open reading frame is homologous to open reading frames located in a similar position with respect to petA in chloroplast DNA from Marchantia polymorpha, tobacco, rice, wheat and Vicia faba. The sequence around a conserved histidine residue in a putative membrane-spanning region of the polypeptide resembles sequences present in cytochrome b from chromaffin granules and neutrophil membranes, suggesting that the open reading frame may encode a haem-binding polypeptide, possibly a b-type cytochrome. Northern hybridisation analysis indicates the presence in pea chloroplasts of a complex pattern of transcripts containing ORF231. Large transcripts of 5.5 kb, 4.3 kb, 3.4 kb and 2.7 kb encode both ORF231 and apocytochrome f, indicating that ORF231 and petA are co-transcribed.

Two small open reading frames are co-transcribed with the pea chloroplast genes for the polypeptides of cytochrome b-559

The genes encoding the 9 kDa and 4 kDa polypeptides of cytochrome b-559 have been located in pea chloroplast DNA by coupled transcription-translation of cloned restriction fragments of chloroplast DNA in a cell-free extract of Escherichia coli and by nucleotide sequence analysis. The genes (psbE and psbF) are located approximately 1.0 kbp downstream of the gene for cytochrome f and are transcribed in the opposite direction, similar to the arrangement in the chloroplast genomes of other higher plants. Nucleotide sequence analysis of this region revealed four open reading frames encoding hydrophobic proteins of 83 (psbE), 39 (psbF), 38 and 40 amino acid residues, which are co-transcribed as a single major RNA of 1.1 kb. The 5' and 3' ends of this RNA have been located by primer extension and S1 nuclease mapping. The 5' end of the RNA is located 140 bp upstream of the initiating ATG codon of psbE and is preceded by typical chloroplast promoter sequences. The 3' end of the RNA is located approximately 515 bp downstream of the TAA stop codon of psbF close to a stable stem-loop structure.

Synthesis and assembly of the cytochrome b-f complex in higher plants

The cytochrome b-f complex is composed of four polypeptide subunits, three of which, cytochrome f, cytochrome b-563 and subunit IV, are encoded in chloroplast DNA and synthesised within the chloroplast, and the fourth, the Rieske FeS protein, is encoded in nuclear DNA and synthesised in the cytoplasm. The assembly of the cytochrome b-f complex therefore requires the interaction of subunits encoded by different genomes. A key role for the nuclear-encoded Rieske FeS protein in the assembly of the complex is suggested by a study of cytochrome b-f complex mutants. The assembly of individual subunits of the complex may be regulated by the availability of prosthetic groups. The genes for the chloroplast-encoded subunits and cDNA clones for the Rieske FeS protein have been isolated and characterised. Cytochrome f and the Rieske FeS protein are synthesised initially with N-terminal presequences required for their correct assembly within the chloroplast. The deduced amino acid sequences of the four subunits have been used to suggest models for the arrangement of the polypeptides in the thylakoid membrane.

Binding of pea cytochrome f to the inner membrane of Escherichia coli requires the bacterial secA gene product

Various sequences from the 5' end of the pea chloroplast gene for cytochrome f have been fused in the correct reading frame with lacZ, and the cellular location of the hybrid polypeptides in Escherichia coli has been examined. Hybrid polypeptides containing N-terminal parts of cytochrome f are located in the cytoplasmic membrane of E. coli. Membrane localization is most efficient when the intact signal sequence of cytochrome f is present at the N-terminal end of the fusion proteins. Fusion within the signal sequence, so that the processing site is absent, reduces the efficiency of membrane binding. Membrane insertion of fusion proteins containing signal sequences is prevented in a temperature-sensitive secA strain at the nonpermissive temperature and the hybrid proteins accumulate in the cytoplasm. This indicates that specific recognition of the chloroplast signal sequence occurs in the bacterial secretory pathway.

Structure and topology of cytochrome f in pea chloroplast membranes

A transmembrane arrangement of cytochrome f in chloroplast thylakoid membranes, with the N-terminal heme-containing region in the intrathylakoid space and a 15 amino acid C-terminal sequence in the stroma, is suggested by the amino acid sequence deduced from the nucleotide sequence of the pea chloroplast gene. This topology has been confirmed by partial proteolysis of the polypeptide in intact and disrupted thylakoid membranes and in inside-out and right-side-out vesicles of chloroplast membranes.