Evolution of termination codons of proteins and the TAG-TGA paradox

In most eukaryotes and prokaryotes TGA is used at a significantly higher frequency than TAG as termination codon of protein-coding genes. Although this phenomenon has been recognized several years ago, there is no generally accepted explanation for the TAG-TGA paradox. Our analyses of human mutation data revealed that out of the eighteen sense codons that can give rise to a nonsense codon by single base substitution, the CGA codon is exceptional: it gives rise to the TGA stop codon at an order of magnitude higher rate than the other codons. Here we propose that the TAG-TGA paradox is due to methylation and hypermutabilty of CpG dinucleotides. In harmony with this explanation, we show that the coding genomes of organisms with strong CpG methylation have a significant bias for TGA whereas those from organisms that lack CpG methylation use TGA and TAG termination codons with similar probability.

Use of signals of positive and negative selection to distinguish cancer genes and passenger genes

A major goal of cancer genomics is to identify all genes that play critical roles in carcinogenesis. Most approaches focused on genes positively selected for mutations that drive carcinogenesis and neglected the role of negative selection. Some studies have actually concluded that negative selection has no role in cancer evolution. We have re-examined the role of negative selection in tumor evolution through the analysis of the patterns of somatic mutations affecting the coding sequences of human genes. Our analyses have confirmed that tumor suppressor genes are positively selected for inactivating mutations, oncogenes, however, were found to display signals of both negative selection for inactivating mutations and positive selection for activating mutations. Significantly, we have identified numerous human genes that show signs of strong negative selection during tumor evolution, suggesting that their functional integrity is essential for the growth and survival of tumor cells.

The SDHB Arg230His mutation causing familial paraganglioma alters glycolysis in a new Caenorhabditis elegans model

The conserved B-subunit of succinate dehydrogenase (SDH) participates in the tricarboxylic acid cycle (TCA) cycle and mitochondrial electron transport. The Arg230His mutation in SDHB causes heritable pheochromocytoma/paraganglioma (PPGL). In Caenorhabditiselegans, we generated an in vivo PPGL model (SDHB-1 Arg244His; equivalent to human Arg230His), which manifests delayed development, shortened lifespan, attenuated ATP production and reduced mitochondrial number. Although succinate is elevated in both missense and null sdhb-1(gk165) mutants, transcriptomic comparison suggests very different causal mechanisms that are supported by metabolic analysis, whereby only Arg244His (not null) worms demonstrate elevated lactate/pyruvate levels, pointing to a missense-induced, Warburg-like aberrant glycolysis. In silico predictions of the SDHA-B dimer structure demonstrate that Arg230His modifies the catalytic cleft despite the latter's remoteness from the mutation site. We hypothesize that the Arg230His SDHB mutation rewires metabolism, reminiscent of metabolic reprogramming in cancer. Our tractable model provides a novel tool to investigate the metastatic propensity of this familial cancer and our approach could illuminate wider SDH pathology.This article has an associated First Person interview with the first author of the paper.

Structure, function and disease relevance of Wnt inhibitory factor 1, a secreted protein controlling the Wnt and hedgehog pathways

Wnts and Hedgehogs (Hh) are large, lipid-modified extracellular morphogens that play key roles in embryonic development and stem cell proliferation of Metazoa. Both morphogens signal through heptahelical Frizzled-type receptors of the G-Protein Coupled Receptor family and there are several other similarities that suggest a common evolutionary origin of the Hh and Wnt pathways. There is evidence that the secreted protein, Wnt inhibitory factor 1 (WIF1) modulates the activity of both Wnts and Hhs and may thus contribute to the intertwining of these pathways. In this article, we review the structure, evolution, molecular interactions and functions of WIF1 with major emphasis on its role in carcinogenesis.

Morphological Stasis and Proteome Innovation in Cephalochordates

Lancelets, extant representatives of basal chordates, are prototypic examples of evolutionary stasis; they preserved a morphology and body-plan most similar to the fossil chordates from the early Cambrian. Such a low level of morphological evolution is in harmony with a low rate of amino acid substitution; cephalochordate proteins were shown to evolve slower than those of the slowest evolving vertebrate, the elephant shark. Surprisingly, a study comparing the predicted proteomes of Chinese amphioxus, Branchiostoma belcheri and the Florida amphioxus, Branchiostoma floridae has led to the conclusion that the rate of creation of novel domain combinations is orders of magnitude greater in lancelets than in any other Metazoa, a finding that contradicts the notion that high rates of protein innovation are usually associated with major evolutionary innovations. Our earlier studies on a representative sample of proteins have provided evidence suggesting that the differences in the domain architectures of predicted proteins of these two lancelet species reflect annotation errors, rather than true innovations. In the present work, we have extended these studies to include a larger sample of genes and two additional lancelet species, Asymmetron lucayanum and Branchiostoma lanceolatum. These analyses have confirmed that the domain architecture differences of orthologous proteins of the four lancelet species are because of errors of gene prediction, the error rate in the given species being inversely related to the quality of the transcriptome dataset that was used to aid gene prediction.

Influence of WFIKKN1 on BMP1-mediated activation of latent myostatin

The NTR domain of WFIKKN1 protein has been shown to have significant affinity for the prodomain regions of promyostatin and latent myostatin but the biological significance of these interactions remained unclear. In view of its role as a myostatin antagonist, we tested the assumption that WFIKKN1 inhibits the release of myostatin from promyostatin and/or latent myostatin. WFIKKN1 was found to have no effect on processing of promyostatin by furin, the rate of cleavage of latent myostatin by BMP1, however, was significantly enhanced in the presence of WFIKKN1 and this enhancer activity was superstimulated by heparin. Unexpectedly, WFIKKN1 was also cleaved by BMP1 and our studies have shown that the KKN1 fragment generated by BMP1-cleavage of WFIKKN1 contributes most significantly to the observed enhancer activity. Analysis of a pro-TGF-β -based homology model of homodimeric latent myostatin revealed that the BMP1-cleavage sites are buried and not readily accessible to BMP1. In view of this observation, the most plausible explanation for the BMP1-enhancer activity of the KKN1 fragment is that it shifts a conformational equilibrium of latent myostatin from the closed circular structure of the homodimer to a more open form, making the cleavage sites more accessible to BMP1. On the other hand, the observation that the enhancer activity of KKN1 is superstimulated in the presence of heparin is explained by the fact KKN1, latent myostatin, and BMP1 have affinity for heparin and these interactions with heparin increase the local concentrations of the reactants thereby facilitating the action of BMP1.

ENZYMES

Furin: EC 3.4.21.75; BMP1, bone morphogentic protein 1 or procollagen C-endopeptidase: EC 3.4.24.19.

K153R polymorphism in myostatin gene increases the rate of promyostatin activation by furin

Recent studies demonstrated an association between the K153R polymorphism in the myostatin gene with extreme longevity, lower muscle strength and obesity but the molecular basis of these associations has not been clarified. Here, we show that the K153R mutation significantly increases the rate of proteolysis of promyostatin by furin, but has no effect on the activity of the latent complex or the cleavage of the latent complex by bone morphogenetic protein 1 (BMP-1). The increased rate of activation of K153R mutant promyostatin may explain why this polymorphism is associated with obesity, lower muscle strength and extension of lifespan.

FixPred: a resource for correction of erroneous protein sequences

Protein databases are heavily contaminated with erroneous (mispredicted, abnormal and incomplete) sequences and these erroneous data significantly distort the conclusions drawn from genome-scale protein sequence analyses. In our earlier work we described the MisPred resource that serves to identify erroneous sequences; here we present the FixPred computational pipeline that automatically corrects sequences identified by MisPred as erroneous. The current version of the associated FixPred database contains corrected UniProtKB/Swiss-Prot and NCBI/RefSeq sequences from Homo sapiens, Mus musculus, Rattus norvegicus, Monodelphis domestica, Gallus gallus, Xenopus tropicalis, Danio rerio, Fugu rubripes, Ciona intestinalis, Branchostoma floridae, Drosophila melanogaster and Caenorhabditis elegans; future releases of the FixPred database will include corrected sequences of additional Metazoan species. The FixPred computational pipeline and database (http://www.fixpred.com) are easily accessible through a simple web interface coupled to a powerful query engine and a standard web service. The content is completely or partially downloadable in a variety of formats.

Database URL:http://www.fixpred.com

MisPred: a resource for identification of erroneous protein sequences in public databases

Correct prediction of the structure of protein-coding genes of higher eukaryotes is still a difficult task; therefore, public databases are heavily contaminated with mispredicted sequences. The high rate of misprediction has serious consequences because it significantly affects the conclusions that may be drawn from genome-scale sequence analyses of eukaryotic genomes. Here we present the MisPred database and computational pipeline that provide efficient means for the identification of erroneous sequences in public databases. The MisPred database contains a collection of abnormal, incomplete and mispredicted protein sequences from 19 metazoan species identified as erroneous by MisPred quality control tools in the UniProtKB/Swiss-Prot, UniProtKB/TrEMBL, NCBI/RefSeq and EnsEMBL databases. Major releases of the database are automatically generated and updated regularly. The database (http://www.mispred.com) is easily accessible through a simple web interface coupled to a powerful query engine and a standard web service. The content is completely or partially downloadable in a variety of formats. DATABASE URL: http://www.mispred.com.

Latent myostatin has significant activity and this activity is controlled more efficiently by WFIKKN1 than by WFIKKN2

Myostatin, a negative regulator of skeletal muscle growth, is produced from myostatin precursor by multiple steps of proteolytic processing. After cleavage by a furin-type protease, the propeptide and growth factor domains remain associated, forming a noncovalent complex, the latent myostatin complex. Mature myostatin is liberated from latent myostatin by bone morphogenetic protein 1/tolloid proteases. Here, we show that, in reporter assays, latent myostatin preparations have significant myostatin activity, as the noncovalent complex dissociates at an appreciable rate, and both mature and semilatent myostatin (a complex in which the dimeric growth factor domain interacts with only one molecule of myostatin propeptide) bind to myostatin receptor. The interaction of myostatin receptor with semilatent myostatin is efficiently blocked by WAP, Kazal, immunoglobulin, Kunitz and NTR domain-containing protein 1 or growth and differentiation factor-associated serum protein 2 (WFIKKN1), a large extracellular multidomain protein that binds both mature myostatin and myostatin propeptide [Kondás et al. (2008) J Biol Chem283, 23677–23684]. Interestingly, the paralogous protein WAP, Kazal, immunoglobulin, Kunitz and NTR domain-containing protein 2 or growth and differentiation factor-associated serum protein 1 (WFIKKN2) was less efficient than WFIKKN1 as an antagonist of the interactions of myostatin receptor with semilatent myostatin. Our studies have shown that this difference is attributable to the fact that only WFIKKN1 has affinity for the propeptide domain, and this interaction increases its potency in suppressing the receptor-binding activity of semilatent myostatin. As the interaction of WFIKKN1 with various forms of myostatin permits tighter control of myostatin activity until myostatin is liberated from latent myostatin by bone morphogenetic protein 1/tolloid proteases, WFIKKN1 may have greater potential as an antimyostatic agent than WFIKKN2.

Characterization of a Wnt-binding site of the WIF-domain of Wnt inhibitory factor-1

A Wnt-binding site of the WIF-domain of Wnt inhibitory factor-1 was localized by structure-guided arginine-scanning mutagenesis in combination with surface plasmon resonance assays. Our observation that substitution of some residues of WIF resulted in an increased affinity for Wnt5a, but decreased affinity for Wnt3a, suggests that these residues may define the specificity spectrum of WIF for Wnts. These results hold promise for a more specific targeting of Wnt family members with WIF variants in various forms of cancer.

WFIKKN1 and WFIKKN2 bind growth factors TGFβ1, BMP2 and BMP4 but do not inhibit their signalling activity

WFIKKN1 and WFIKKN2 are large extracellular multidomain proteins consisting of a WAP domain, a follistatin domain, an immunoglobulin domain, two Kunitz-type protease inhibitor domains and an NTR domain. Recent experiments have shown that both proteins have high affinity for growth and differentiation factor (GDF)8 and GDF11. Here we study the interaction of WFIKKN proteins with several additional representatives of the transforming growth factor (TGF)β family using SPR measurements. Analyses of SPR sensorgrams suggested that, in addition to GDF8 and GDF11, both WFIKKN proteins bind TGFβ1, bone morphogenetic protein (BMP)2 and BMP4 with relatively high affinity (K(d) ∼ 10(-6) m). To assess the biological significance of these interactions we studied the effect of WFIKKN proteins on the activity of GDF8, GDF11, TGFβ1, BMP2 and BMP4 using reporter assays. These studies revealed that WFIKKN1 and WFIKKN2 inhibited the biological activity of GDF8 and GDF11 in the nanomolar range, whereas they did not inhibit the activities of TGFβ1, BMP2 and BMP4 even in the micromolar range. Our data indicate that WFIKKN proteins are antagonists of GDF8 and GDF11, but in the case of TGFβ1, BMP2 and BMP4 they function as growth factor binding proteins. It is suggested that the physical association of WFIKKN proteins with these growth factors may localize their action and thus help to establish growth factor gradients in the extracellular space.

Agrin binds BMP2, BMP4 and TGFbeta1

The C-terminal 95 kDa fragment of some isoforms of vertebrate agrins is sufficient to induce clustering of acetylcholine receptors but despite two decades of intense agrin research very little is known about the function of the other isoforms and the function of the larger, N-terminal part of agrins that is common to all isoforms. Since the N-terminal part of agrins contains several follistatin-domains, a domain type that is frequently implicated in binding TGFβs, we have explored the interaction of the N-terminal part of rat agrin (Agrin-Nterm) with members of the TGFβ family using surface plasmon resonance spectroscopy and reporter assays. Here we show that agrin binds BMP2, BMP4 and TGFβ1 with relatively high affinity, the K_D values of the interactions calculated from SPR experiments fall in the 10⁻⁸ M–10⁻⁷ M range. In reporter assays Agrin-Nterm inhibited the activities of BMP2 and BMP4, half maximal inhibition being achieved at ∼5×10⁻⁷ M. Paradoxically, in the case of TGFβ1 Agrin N-term caused a slight increase in activity in reporter assays. Our finding that agrin binds members of the TGFβ family may have important implications for the role of these growth factors in the regulation of synaptogenesis as well as for the role of agrin isoforms that are unable to induce clustering of acetylcholine receptors. We suggest that binding of these TGFβ family members to agrin may have a dual function: agrin may serve as a reservoir for these growth factors and may also inhibit their growth promoting activity. Based on analysis of the evolutionary history of agrin we suggest that agrin's growth factor binding function is more ancient than its involvement in acetylcholine receptor clustering.

[Developments in cancer care with innovative genomics. 2008 report of the National Cancer Consortium]

In the 3rd year of the program 8 new molecular diagnostic services have been introduced to clinic in the management of breast-, lung-, colorectal cancers as well as in GIST and melanoma. Two patents have been filed for innovative modulation of mito/motogenic signaling pathways in cancer cells. In preclinical models of human cancer a functional imaging technique was developed to detect vascular eff ects of erythropoietin. Using a genomic approach, the sequential changes in human melanoma during systemic dissemination were determined revealing several novel potential prognostic factors and some interesting novel targets for therapy.

NMR solution structure of the neurotrypsin Kringle domain

Neurotrypsin is a multidomain protein that serves as a brain-specific serine protease. Here we report the NMR structure of its kringle domain, NT/K. The data analysis was performed with the BACUS (Bayesian analysis of coupled unassigned spins) algorithm. This study presents the first application of BACUS to the structure determination of a 13C unenriched protein for which no prior experimental 3D structure was available. NT/K adopts the kringle fold, consisting of an antiparallel beta-sheet bridged by an overlapping pair of disulfides. The structure reveals the presence of a surface-exposed left-handed polyproline II helix that is closely packed to the core beta-structure. This feature distinguishes NT/K from other members of the kringle fold and points toward a novel functional role for a kringle domain. Functional divergence among kringle domains is discussed on the basis of their surface and electrostatic characteristics.

The second von Willebrand type A domain of cochlin has high affinity for type I, type II and type IV collagens

Cochlin is colocalized with type II collagen in the extracellular matrix of cochlea and has been suggested to interact with this collagen. Here we show that the second von Willebrand type A domain of cochlin has affinity for type II collagen, as well as type I and type IV collagens whereas the LCCL-domain of cochlin has no affinity for these proteins. The implications of these findings for the mechanism whereby cochlin mutations cause the dominant negative DFNA9-type hearing loss are discussed.

Identification and correction of abnormal, incomplete and mispredicted proteins in public databases

Background

Despite significant improvements in computational annotation of genomes, sequences of abnormal, incomplete or incorrectly predicted genes and proteins remain abundant in public databases. Since the majority of incomplete, abnormal or mispredicted entries are not annotated as such, these errors seriously affect the reliability of these databases. Here we describe the MisPred approach that may provide an efficient means for the quality control of databases. The current version of the MisPred approach uses five distinct routines for identifying abnormal, incomplete or mispredicted entries based on the principle that a sequence is likely to be incorrect if some of its features conflict with our current knowledge about protein-coding genes and proteins: (i) conflict between the predicted subcellular localization of proteins and the absence of the corresponding sequence signals; (ii) presence of extracellular and cytoplasmic domains and the absence of transmembrane segments; (iii) co-occurrence of extracellular and nuclear domains; (iv) violation of domain integrity; (v) chimeras encoded by two or more genes located on different chromosomes.

Results

Analyses of predicted EnsEMBL protein sequences of nine deuterostome (Homo sapiens, Mus musculus, Rattus norvegicus, Monodelphis domestica, Gallus gallus, Xenopus tropicalis, Fugu rubripes, Danio rerio and Ciona intestinalis) and two protostome species (Caenorhabditis elegans and Drosophila melanogaster) have revealed that the absence of expected signal peptides and violation of domain integrity account for the majority of mispredictions. Analyses of sequences predicted by NCBI's GNOMON annotation pipeline show that the rates of mispredictions are comparable to those of EnsEMBL. Interestingly, even the manually curated UniProtKB/Swiss-Prot dataset is contaminated with mispredicted or abnormal proteins, although to a much lesser extent than UniProtKB/TrEMBL or the EnsEMBL or GNOMON-predicted entries.

Conclusion

MisPred works efficiently in identifying errors in predictions generated by the most reliable gene prediction tools such as the EnsEMBL and NCBI's GNOMON pipelines and also guides the correction of errors. We suggest that application of the MisPred approach will significantly improve the quality of gene predictions and the associated databases.

Both WFIKKN1 and WFIKKN2 have high affinity for growth and differentiation factors 8 and 11

WFIKKN1 and WFIKKN2 are large extracellular multidomain proteins consisting of a WAP, a follistatin, an immunoglobulin, two Kunitz-type protease inhibitor domains, and an NTR domain. Recent experiments have shown that WFIKKN2 protein binds mature GDF8/myostatin and myostatin propeptide and inhibits the biological activity of myostatin (Hill, J. J., Qiu, Y., Hewick, R. M., and Wolfman, N. M. (2003) Mol. Endocrinol. 17, 1144-1154). Here we show that the paralogue of this protein, WFIKKN1, also binds to both myostatin and myostatin propeptide and that both WFIKKN1 and WFIKKN2 bind GDF11, the growth and differentiation factor most closely related to myostatin, with high affinity. Structure-function studies on WFIKKN1 have revealed that the follistatin domain is primarily responsible for the binding of mature growth factor, whereas the NTR domain contributes most significantly to the interaction with myostatin propeptide. Analysis of the evolutionary histories of WFIKKN1/WFIKKN2 and GDF8/GDF11 proteins indicates that the functional association of an ancestral WFIKKN protein with an ancestor of GDF8/11 may date back to cephalochordates/urochordates. Although duplication of the corresponding genes gave rise to WFIKKN1/WFIKKN2 and GDF8/GDF11 in early vertebrates, the data presented here suggest that there is significant functional overlap of the paralogous proteins.

The implications of alternative splicing in the ENCODE protein complement

Alternative premessenger RNA splicing enables genes to generate more than one gene product. Splicing events that occur within protein coding regions have the potential to alter the biological function of the expressed protein and even to create new protein functions. Alternative splicing has been suggested as one explanation for the discrepancy between the number of human genes and functional complexity. Here, we carry out a detailed study of the alternatively spliced gene products annotated in the ENCODE pilot project. We find that alternative splicing in human genes is more frequent than has commonly been suggested, and we demonstrate that many of the potential alternative gene products will have markedly different structure and function from their constitutively spliced counterparts. For the vast majority of these alternative isoforms, little evidence exists to suggest they have a role as functional proteins, and it seems unlikely that the spectrum of conventional enzymatic or structural functions can be substantially extended through alternative splicing.

[Developments in cancer management by innovative genomics. 2006 report of the National Cancer Consortium]

Research on developing molecular diagnostics for hereditary cancers resulted in establishing diagnostic services for familiar polyposis and non-polyposis patients (mutation determination of APC, MYH, STK11, SMAD4, MLH1, MSH2). In familiar testicular cancers the role of gr/gr gene on Y chromosome was identified. Molecular diagnostic tool was established to monitor the progression of follicular lymphoma using Bcl-2/IgH fusion sequences. Molecular diagnostic tools were developed to monitor circulating endothelial precursor cells (CEP) as well and the technique was tested in lung cancer patients. In malignant melanoma we have tested several potential novel markers among which ryanodine receptor seems to be a promising one, while the functional P2X7 receptor may serve as a therapeutic target. We have determined the tyrosine kinase "kinome" profile of HER-2-amplified breast cancers. Furthermore, the "kinome" profile was found to be characteristic for head and neck cancers of various anatomical location. Based on previous studies on the anti-migratory and antimetastatic potential of low-molecular-weight heparins, we have identified short heparin-derived oligosaccharides with maintained antimetastatic- but non-anticoagulant potentials. Pharmacogenomic studies on the role of polymorphism of the serine-hydroxymethyl-transferase (SHMT) gene in the efficacy of 5-FU and FOLFIRI protocols of colorectal cancer patients revealed a significant effect resulting in altered overall survival as well.

NMR structure of the WIF domain of the human Wnt-inhibitory factor-1

The human Wnt-binding protein Wnt-inhibitory factor-1 (WIF-1) comprises an N-terminal WIF module followed by five EGF-like repeats. Here we report the three-dimensional structure of the WIF domain of WIF-1 determined by NMR spectroscopy. The fold consists of an eight-stranded beta-sandwich reminiscent of the immunoglobulin fold. Residual detergent (Brij-35) used in the refolding protocol was found to bind tightly to the WIF domain. The binding site was identified by intermolecular nuclear Overhauser effects observed between the WIF domain and the alkyl chain of the detergent. The results point to a possible role of WIF domains as a recognition motif of Wnt and Drosophila Hedgehog proteins that are activated by palmitoylation.

Modules, multidomain proteins and organismic complexity

Originally the term 'protein module' was coined to distinguish mobile domains that frequently occur as building blocks of diverse multidomain proteins from 'static' domains that usually exist only as stand-alone units of single-domain proteins. Despite the widespread use of the term 'mobile domain', the distinction between static and mobile domains is rather vague as it is not easy to quantify the mobility of domains. In the present work we show that the most appropriate measure of the mobility of domains is the number of types of local environments in which a given domain is present. Ranking of domains with respect to this parameter in different evolutionary lineages highlighted marked differences in the propensity of domains to form multidomain proteins. Our analyses have also shown that there is a correlation between domain size and domain mobility: smaller domains are more likely to be used in the construction of multidomain proteins, whereas larger domains are more likely to be static, stand-alone domains. It is also shown that shuffling of a limited set of modules was facilitated by intronic recombination in the metazoan lineage and this has contributed significantly to the emergence of novel complex multidomain proteins, novel functions and increased organismic complexity of metazoa.

Hydroxamate-based peptide inhibitors of matrix metalloprotease 2

There is major interest in designing inhibitors for matrix metalloproteinase 2 (MMP-2, gelatinase A) since this enzyme is known to be involved in pathological processes such as tumor invasion or rheumatoid arthritis. The majority of MMP-2 inhibitor candidate drugs block the active site of MMP-2 by binding to its catalytic Zn2+ ion through a chelating (hydroxamate, sulphonate etc.) group. Despite the general interest in designing MMP-2 inhibitors, the results with many of the drug candidates were disappointing, their failure was usually explained by cross-reactions with other MMPs. One way to enhance MMP-2 selectivity is to design inhibitors that interact with both the active site and exosites such as the fibronectin type II (FN2) domains of the enzyme. In the present work, we have examined the inhibitory potential and MMP-2 selectivity of hydroxamates of three groups of peptides known to bind to the collagen-binding FN2 domains of MMP-2. The first type of peptides consisted of collagen-like (Pro-Pro-Gly)(n) repeats, peptides of the second group were identified from a random 15-mer phage display library based on their binding to immobilized FN2 domains of MMP-2. A hydroxamate of peptide p33-42, known to bind to the third FN2 domain of MMP-2 has also been tested. Our studies have shown that these compounds inhibited MMP-2 with IC50 values of 10-100 microM. The fact that their inhibitory potential was nearly identical for MMP-2del, a recombinant version of MMP-2 that lacks the FN2 domains, suggests that inhibition is not mediated by their binding to FN2 domains. It seems likely that the failure to exploit interaction with the FN2 domains is due to the fact that the FN2 domains and the catalytic domain of MMP-2 tumble independently, therefore only a tiny fraction of the conformational isomers can bind peptide hydroxamates via both the active site and the FN2 domain(s).

DoOP: Databases of Orthologous Promoters, collections of clusters of orthologous upstream sequences from chordates and plants

DoOP (http://doop.abc.hu/) is a database of eukaryotic promoter sequences (upstream regions) aiming to facilitate the recognition of regulatory sites conserved between species. The annotated first exons of human and Arabidopsis thaliana genes were used as queries in BLAST searches to collect the most closely related orthologous first exon sequences from Chordata and Viridiplantae species. Up to 3000 bp DNA segments upstream from these first exons constitute the clusters in the chordate and plant sections of the Database of Orthologous Promoters. Release 1.0 of DoOP contains 21 061 chordate clusters from 284 different species and 7548 plant clusters from 269 different species. The database can be used to find and retrieve promoter sequences of a given gene from various species and it is also suitable to see the most trivial conserved sequence blocks in the orthologous upstream regions. Users can search DoOP with either sequence or text (annotation) to find promoter clusters of various genes. In addition to the sequence data, the positions of the conserved sequence blocks derived from multiple alignments, the positions of repetitive elements and the positions of transcription start sites known from the Eukaryotic Promoter Database (EPD) can be viewed graphically.

Insertion of spliceosomal introns in proto-splice sites: the case of secretory signal peptides

Analysis of the exon-intron structures of 2208 human genes has revealed that there is a statistically highly significant excess of phase 1 introns in the vicinity of the signal peptide cleavage sites. It is suggested that amino acid sequences surrounding signal peptide cleavage sites are significantly enriched in phase 1 proto-splice sites and this has favored insertion of spliceosomal introns in these sites.

Modular Autonomy, Ligand Specificity, and Functional Cooperativity of the Three In-tandem Fibronectin Type II Repeats from Human Matrix Metalloproteinase 2

Matrix metalloproteinase 2 (MMP-2) contains three fibronectin type II (col) modules that contribute to its collagen specificity. We observed that the CD spectra of the separate col modules account for the CD and temperature profiles of the in-tandem col-123 construct. Thus, to the extent of not significantly perturbing the secondary structure and thermal stability characteristics of the neighboring units, the domains within col-123 do not interact. Via NMR, we investigated ligand binding properties of the three repeats within col-123: col-123/1 (the col-1 domain within col-123), col-123/2, and col-123/3. Interactions of col-123 with the collagen mimic peptide (Pro-Pro-Gly)6 (PPG6) and propeptide segment PIIKFPGDVA (p33-42) were studied. While col-123/1 and col-123/2 bound PPG6, they interacted more weakly with p33-42. In contrast, col-123/3 exhibited a higher affinity for p33-42 than for PPG6. Thus, despite their structural homology, the col repeats of MMP-2 differ in substrate specificity. Furthermore the binding affinities toward the three in-tandem col repeats were close to those determined for the individual isolated domains or for col-12/1, indicating that vis-a-vis these ligands each module interacts essentially as an autonomous unit. Interestingly the domains within col-123 exhibited enhanced affinities for Hel3, a construct that contains ((Gly-Pro-Pro)12)3 in triple helical configuration. Nevertheless the affinities were significantly higher for col-123/1 and col-123/2 relative to col-123/3 in line with their behaviors toward PPG6. This hints at a cooperative participation toward Hel3, which is a closer mimic of collagen, a hypothesis that is supported by the detected lower affinities of col-12/1, col-12/2, col-2, col-23/2, col-3, and col-23/3 for Hel3.

Low density lipoprotein receptor-related protein mediates endocytic clearance of pro-MMP-2.TIMP-2 complex through a thrombospondin-independent mechanism

The low density lipoprotein receptor-related protein (LRP) mediates the endocytic clearance of various proteinases and proteinase.inhibitor complexes, including thrombospondin (TSP)-dependent endocytosis of matrix metalloproteinase (MMP)-2 (or gelatinase A), a key effector of extracellular matrix remodeling and cancer progression. However, the zymogen of MMP-2 (pro-MMP-2) mostly occurs in tissues as a complex with the tissue inhibitor of MMPs (TIMP-2). Here we show that clearance of the pro-MMP-2.TIMP-2 complex is also mediated by LRP, because addition of receptor-associated protein (RAP), a natural LRP ligand antagonist, inhibited endocytosis and lysosomal degradation of (125)I-pro-MMP-2.TIMP-2. Both TIMP-2 and the pro-MMP-2 collagen-binding domain independently competed for endocytosis of (125)I-pro-MMP-2.TIMP-2 complex. Surface plasmon resonance studies indicated that pro-MMP-2, TIMP-2, and pro-MMP-2.TIMP-2 directly interact with LRP in the absence of TSP. LRP-mediated endocytic clearance of (125)I-pro-MMP-2 was inhibited by anti-TSP antibodies and accelerated upon complexing with TSP-1, but these treatments had no effect on (125)I-pro-MMP-2.TIMP-2 uptake. This implies that mechanisms of clearance by LRP of pro-MMP-2 and pro-MMP-2.TIMP-2 complex are different. Interestingly, RAP did not inhibit binding of (125)I-pro-MMP-2.TIMP-2 to the cell surface. We conclude that clearance of pro-MMP-2.TIMP-2 complex is a TSP-independent two-step process, involving (i) initial binding to the cell membrane in a RAP-insensitive manner and (ii) subsequent LRP-dependent (RAP-sensitive) internalization and degradation.

Evidence that human genes of modular proteins have retained significantly more ancestral introns than their fly or worm orthologues

Comparison of the exon-intron structures of human, fly and worm orthologues of mosaic genes assembled from class 1-1 modules by exon-shuffling has revealed that human genes retained significantly more of the original inter-module introns than their protostome orthologues. It is suggested that the much higher rate of intron loss in the worm- and insect lineages than in the chordate lineage reflects their greater tendency for genome compaction.

Mapping Modifiers Affecting Muscularity of the Myostatin Mutant (MstnCmpt-dl1Abc) Compact Mouse

The hypermuscular Compact phenotype was first noted in a line of mice selected for high body weight and protein content. A new line, based on mice showing the Compact phenotype, was formed and selected for maximum expression of the Compact phenotype. Previously we mapped and identified a 12-bp deletion in the myostatin gene, denoted MstnCmpt-dl1Abc, which can be considered as a major gene responsible for the hypermuscular phenotype. Genetic analysis revealed that full expression of the hypermuscular phenotype requires the action of modifier loci in addition to MstnCmpt-dl1Abc. To map these modifier loci, an interspecific F2 population was generated between Comp9, an inbred line homozygous for MstnCmpt-dl1Abc, and CAST/Ei, an inbred line generated from Mus musculus castaneus. Selective DNA pooling and genotyping, separately by gender, was carried out within a subpopulation of the F2 consisting of individuals homozygous for MstnCmpt-dl1Abc. Significant association with hypermuscularity at a false discovery rate (FDR) of 0.05 was found for markers on chromosomes 3, 5, 7, 11, 16, and X. In all cases, the marker allele derived from the Comp9 parent showed a higher frequency in the hypermuscular group and the CAST/Ei allele in the normal group. The modifier loci apparently exerted their effects on muscularity only in the presence of MstnCmpt-dl1Abc.

Modular assembly of genes and the evolution of new functions

Modular assembly of novel genes from existing genes has long been thought to be an important source of evolutionary novelty. Thanks to major advances in genomic studies it has now become clear that this mechanism contributed significantly to the evolution of novel biological functions in different evolutionary lineages. Analyses of completely sequenced bacterial, archaeal and eukaryotic genomes has revealed that modular assembly of novel constituents of various eukaryotic intracellular signalling pathways played a major role in the evolution of eukaryotes. Comparison of the genomes of single-celled eukaryotes, multicellular plants and animals has also shown that the evolution of multicellularity was accompanied by the assembly of numerous novel extracellular matrix proteins and extracellular signalling proteins that are absolutely essential for multicellularity. There is now strong evidence that exon-shuffling played a general role in the assembly of the modular proteins involved in extracellular communications of metazoa. Although some of these proteins seem to be shared by all major groups of metazoa, others are restricted to certain evolutionary lineages. The genomic features of the chordates appear to have favoured intronic recombination as evidenced by the fact that exon-shuffling continued to be a major source of evolutionary novelty during vertebrate evolution.

Expression, purification and characterization of the second Kunitz-type protease inhibitor domain of the human WFIKKN protein

Recently we have described a novel secreted protein (the WFIKKN protein) that consists of multiple types of protease inhibitory modules, including two tandem Kunitz-type protease inhibitor-domains. On the basis of its homologies we have suggested that the WFIKKN protein is a multivalent protease inhibitor that may control the action of different proteases. In the present work we have expressed the second Kunitz-type protease inhibitor domain of the human protein WFIKKN in Escherichia coli, purified it by affinity chromatography on trypsin-Sepharose and its structure was characterized by CD spectroscopy. The recombinant protein was found to inhibit trypsin (Ki = 9.6 nm), but chymotrypsin, elastase, plasmin, pancreatic kallikrein, lung tryptase, plasma kallikrein, thrombin, urokinase or tissue plasminogen activator were not inhibited by the recombinant protein even at 1 microm concentration. In view of the marked trypsin-specificity of the inhibitor it is suggested that its physiological target may be trypsin.

Peptide ligands for the fibronectin type II modules of matrix metalloproteinase 2 (MMP-2)

The interaction of matrix metalloproteinase 2 (MMP-2) with gelatin is mediated by three repeats homologous to fibronectin type II (FN2) modules, which are inserted in the catalytic domain in proximity of the active site. We screened a random 15-mer phage display library to identify peptides that interact with the FN2 modules of MMP-2. Interestingly, the selected peptides are not gelatin-like and do not share a common, obvious sequence motif. However, they contain a high proportion of aromatic residues. The interactions of two peptides, WHWRH0RIPLQLAAGR and THSHQWRHHQFPAPT, with constructs comprising the in-tandem first and second and second and third FN2 modules of MMP-2 (Col-12 and Col-23, respectively) were characterized by NMR. Both peptides interact with Col-12 and Col-23 with apparent association constants in the mm(-1) range. Peptide binding results in perturbation of signals from residues located in the gelatin-binding pocket and flexible parts of the molecule. Although the former finding suggests that the gelatin-binding site is involved in the contact, the interpretation of the latter is less straightforward and may well reflect both the direct and indirect effects of the interaction.

Expression and characterization of the olfactomedin domain of human myocilin

The olfactomedin-domain has been first identified in olfactomedin, an extracellular matrix protein of the olfactory neuroepithelium. Members of this extracellular domain-family have since been shown to be present in several metazoan proteins, such as latrophilins, myocilins, and noelins, but their biological function is unknown. The olfactomedin-domain of myocilin is of considerable interest, since mutations affecting this domain are associated with primary open angle glaucoma. In order to define structural features of this domain-type we have expressed the olfactomedin-domain of human myocilin in Pichia pastoris. The olfactomedin-domain contains a single disulphide-bond connecting Cys-245 and Cys-433 residues; secondary structure predictions and circular dichroism studies indicate that it consists primarily of beta-strands. It is noteworthy that the majority of mutations associated with severe forms of glaucoma affect residues that reside in conserved secondary structural elements of the olfactomedin-domain or are otherwise critical for the integrity of this protein-fold.

Structural characterization of the second TSP1-module of human thrombospondin

The TSP1-module has been first identified as the type 1 repeat of thrombospondin-1. Members of this extracellular module-family have since been shown to be present in several hundred metazoan proteins as well as in proteins of some protists. Despite the widespread occurrence and biological importance of this module-type, relatively little is known about their three-dimensional structure. To define the structural features of this important module-family, we have expressed the second TSP1-domain of human thrombospondin 1 in Escherichia coli. Amino acid sequencing of proteolytic fragments of the recombinant protein have shown that its disulfide bonds connect the six conserved cysteines in a 1-5, 2-6, 3-4 pattern. Circular dichroism studies on the recombinant protein indicate that the disulfide-bonded TSP1-module consists primarily of distorted beta-strands.

The col-1 module of human matrix metalloproteinase-2 (MMP-2): structural/functional relatedness between gelatin-binding fibronectin type II modules and lysine-binding kringle domains

Human matrix metalloproteinase-2 (MMP-2) contains three in-tandem fibronectin type II (FII) repeats that bind gelatin. Here, we report the NMR solution structure of the first FII module of MMP-2 (col-1). The latter is described as a characteristic, globular FII fold containing two beta-sheets, a stretch of 3(1)-helix, a turn of alpha-helix, and an exposed hydrophobic surface lined with aromatic residues. We show that col-1 binds (Pro-Pro-Gly)6, a mimic of gelatin, with a Ka of approx. 0.42 mm(-1), and that its binding site involves a number of aromatic residues as well as Arg34, as previously found for the second and third homologous repeats. Moreover, the affinity of the in-tandem col-1+2 construct (col-12) toward the longer ligand (Pro-Pro-Gly)12 is twice that for (Pro-Pro-Gly)6, as expected from mass action. A detailed structural comparison between FII and kringle domains indicates that four main conformational features are shared: two antiparallel beta-sheets, a central 3(1)-helix, and the quasiperpendicular orientation of the two proximal Cys-Cys bonds. Structure superposition by optimizing overlap of cystine bridge areas results in close juxtaposition of their main beta-sheets and 31-helices, and reveals that the gelatin binding site of FII modules falls at similar locations and exhibits almost identical topological features to those of the lysine binding site of kringle domains. Thus, despite the minor (<15%) consensus sequence relating FII modules to kringles, there is a strong folding and binding site structural homology between the two domains, enforced by key common conformational determinants.

Distinct expression pattern of two related human proteins containing multiple types of protease-inhibitory modules

We have recently identified a gene (the WFIKKN gene) on human chromosome 16 (16p13.3) that encodes a secreted protein containing WAP-type, Follistatin/ Kazal type, Kunitz-type and NTR-type protease-inhibitory modules and an Immunoglobulin domain [Trexler et al., Proc. Natl. Acad. Sci. USA 98 (2001), 3705 - 3709]. In the present work we show that a gene on chromosome 17 encodes a WFIKKN-related protein (WFIKKNRP) that has the same domain organization as the WFIKKN protein. The exon-intron structure of the two genes is also similar as both genes have a single phase 0 intron that splits their WAP domains in equivalent positions. In view of the presence of several protease inhibitory modules in these proteins it seems likely that they serve to control the action of multiple types of proteases. The tissue expression pattern of the two proteins, however, is markedly different suggesting that they have distinct biological roles. Whereas the WFIKKN gene is expressed primarily in adult pancreas, liver and thymus but not in brain and ovary, significant expression of the WFIKKNRP gene is observed in ovary, testis and brain, but not in liver. Pronounced differences could also be seen in the case of fetal tissues: expression of the WFIKKN gene was highest in the lung, skeletal muscle and liver, whereas the WFIKKNRP gene was expressed primarily in brain, skeletal muscle, thymus and kidney.

Origin of fibronectin type II (FN2) modules: structural analyses of distantly-related members of the kringle family idey the kringle domain of neurotrypsin as a potential link between FN2 domains and kringles

Analysis of complete genome sequences has made it clear that fibronectin type II (FN2) modules are present only in the vertebrate lineage, raising intriguing questions about the origin of this module type. Kringle domains display many similarities to FN2 domains; therefore it was suggested previously that they are highly divergent descendants of the same ancestral protein-fold. Since kringles are present in arthropodes, nematodes, and invertebrate chordates as well as in vertebrates, it is suggested that the FN2 domain arose in the vertebrate lineage through major structural modification of the more ancestral kringle fold. To explore this structural transition, in the present work we compare key structural features of two highly divergent kringle domains (the kringle of Caenorhabditis elegans Ror receptor tyrosine kinase and the kringle of rat neurotrypsin) with those of plasminogen kringles and FN2 domains. Our NMR conformation fingerprinting analysis indicates that characteristic (1)H-NMR markers of kringle or FN2 native folding, such as the dispersion of Trp aromatic connectivities and shifts of the Leu(46)/Thr(16) methyl signals, both decrease in the order kringles > neurotrypsin kringle > FN2 domains. These results suggest that the neurotrypsin kringle may represent an intermediate form between typical kringles and FN2 domains.

NMR structure of the LCCL domain and implications for DFNA9 deafness disorder

The LCCL domain is a recently discovered, conserved protein module named after its presence in Limulus factor C, cochlear protein Coch-5b2 and late gestation lung protein Lgl1. The LCCL domain plays a key role in the autosomal dominant human deafness disorder DFNA9. Here we report the nuclear magnetic resonance (NMR) structure of the LCCL domain from human Coch-5b2, where dominant mutations leading to DFNA9 deafness disorder have been identified. The fold is novel. Four of the five known DFNA9 mutations are shown to involve at least partially solvent-exposed residues. Except for the Trp91Arg mutant, expression of these four LCCL mutants resulted in misfolded proteins. These results suggest that Trp91 participates in the interaction with a binding partner. The unexpected sensitivity of the fold with respect to mutations of solvent-accessible residues might be attributed to interference with the folding pathway of this disulfide-containing domain.

Gelatin-binding region of human matrix metalloproteinase-2: solution structure, dynamics, and function of the COL-23 two-domain construct

Human matrix metalloproteinase-2 (MMP-2) contains an array of three fibronectin type II (FII) modules postulated to interact with gelatin (denatured collagen). Here, we verify that the NMR solution structure of the third FII repeat (COL-3) is similar to that of the second FII repeat (COL-2); characterize its ligand-binding properties; and derive dynamics properties and relative orientation in solution for the two domains of the COL-23 fragment, a construct comprising COL-2 and COL-3 in tandem, with each domain possessing a putative collagen-binding site. Interaction of the synthetic gelatin-like octadecapeptide (Pro-Pro-Gly)(6) (PPG6) with COL-3 is weaker than with COL-2. We found that a synthetic peptide comprising segment 33-42 (peptide 33-42) from the MMP-2 prodomain interacts with COL-3 and, albeit with lower affinity, with COL-2 in a way that mimics PPG6 binding. COL-3 strongly prefers peptide 33-42 over PPG6, which suggests that intramolecular interactions with the prodomain could modulate binding of pro-MMP-2 to its gelatin substrate. In COL-23, the two modules retain their structural individuality and tumble independently. Overall, the NMR data indicate that the relative orientation of the modules in COL-23 is not fixed in solution, that the modules do not interact with one another, and that COL-23 is rather flexible. The binding sites face opposite each other, and their responses to, and normalized affinities for, the longer ligand PPG12 are virtually identical to those of the individual domains for PPG6, thus precluding co- operativity, although they may interact simultaneously with multiple sites of the extracellular matrix.

Involvement of fibronectin type II repeats in the efficient inhibition of gelatinases A and B by long-chain unsaturated fatty acids

The matrix metalloproteinases gelatinase A (MMP-2) and gelatinase B (MMP-9) are implicated in the physiological and pathological breakdown of several extracellular matrix proteins. In the present study, we show that long-chain fatty acids (e.g. oleic acid, elaidic acid, and cis- and trans-parinaric acids) inhibit gelatinase A as well as gelatinase B with K(i) values in the micromolar range but had only weak inhibitory effect on collagenase-1 (MMP-1), as assessed using synthetic or natural substrates. The inhibition of gelatinases depended on fatty acid chain length (with C18 > C16, C14, and C10), and the presence of unsaturations increased their inhibitory capacity on both types of gelatinase. Ex vivo experiments on human skin tissue sections have shown that micromolar concentrations of a long-chain unsaturated fatty acid (elaidic acid) protect collagen and elastin fibers against degradation by gelatinases A and B, respectively. In order to understand why gelatinases are more susceptible than collagenase-1 to inhibition by long-chain fatty acids, the possible role of the fibronectin-like domain (a domain unique to gelatinases) in binding inhibitory fatty acids was investigated. Affinity and kinetic studies with a recombinant fibronectin-like domain of gelatinase A and with a recombinant mutant of gelatinase A from which this domain had been deleted pointed to an interaction of long-chain fatty acids with the fibronectin-like domain of the protease. Surface plasmon resonance studies on the interaction of long-chain fatty acids with the three individual type II modules of the fibronectin-like domain of gelatinase A revealed that the first type II module is primarily responsible for binding these compounds.

Localization of disulfide bonds in the frizzled module of Ror1 receptor tyrosine kinase

The frizzled (FRZ) module is a novel module type that was first identified in G-protein-coupled receptors of the frizzled and smoothened families and has since been shown to be present in several secreted frizzled-related proteins, in some modular proteases, in collagen XVIII, and in various receptor tyrosine kinases of the Ror family. The FRZ modules constitute the extracellular ligand-binding region of frizzled receptors and are known to mediate signals of WNT family members through these receptors. With an eye toward defining the structure of this important module family, we have expressed the FRZ domain of rat Ror1 receptor tyrosine kinase in Pichia pastoris. By proteolytic digestion and amino acid sequencing the disulfide bonds were found to connect the 10 conserved cysteines in a 1-5, 2-4, 3-8, 6-10, and 7-9 pattern. Circular dichroism and differential scanning calorimetry studies on the recombinant protein indicate that the disulfide-bonded FRZ module corresponds to a single, compact, and remarkably stable folding domain possessing both alpha-helices and beta-strands.

A human protein containing multiple types of protease-inhibitory modules

By using sensitive homology-search and gene-finding programs, we have found that a genomic region from the tip of the short arm of human chromosome 16 (16p13.3) encodes a putative secreted protein consisting of a domain related to the whey acidic protein (WAP) domain, a domain homologous with follistatin modules of the Kazal-domain family (FS module), an immunoglobulin-related domain (Ig domain), two tandem domains related to Kunitz-type protease inhibitor modules (KU domains), and a domain belonging to the recently defined NTR-module family (NTR domain). The gene encoding these WAP, FS, Ig, KU, and NTR modules (hereafter referred to as the WFIKKN gene) is intron-depleted--its single 1,157-bp intron splits the WAP module. The validity of our gene prediction was confirmed by sequencing a WFIKKN cDNA cloned from a lung cDNA library. Studies on the tissue-expression pattern of the WFIKKN gene have shown that the gene is expressed primarily in pancreas, kidney, liver, placenta, and lung. As to the function of the WFIKKN protein, it is noteworthy that it contains FS, WAP, and KU modules, i.e., three different module types homologous with domains frequently involved in inhibition of serine proteases. The protein also contains an NTR module, a domain type implicated in inhibition of zinc metalloproteinases of the metzincin family. On the basis of its intriguing homologies, we suggest that the WFIKKN protein is a multivalent protease inhibitor that may control the action of multiple types of serine proteases as well as metalloproteinase(s).

Were protein internal repeats formed by "bricolage"?

Is evolution an engineer, or is it a tinkerer--a "bricoleur"--building up complex molecules in organisms by increasing and adapting the materials at hand? An analysis of completely sequenced genomes suggests the latter, showing that increasing repetition of modules within the proteins encoded by these genomes is correlated with increasing complexity of the organism.

Global patterns of human DNA sequence variation in a 10-kb region on chromosome 1

Human DNA variation is currently a subject of intense research because of its importance for studying human origins, evolution, and demographic history and for association studies of complex diseases. A approximately 10-kb region on chromosome 1, which contains only four small exons (each <155 bp), was sequenced for 61 humans (20 Africans, 20 Asians, and 21 Europeans) and for 1 chimpanzee, 1 gorilla, and 1 orangutan. We found 52 polymorphic sites among the 122 human sequences and 382 variant sites among the human, chimpanzee, gorilla, and orangutan sequences. For the introns sequenced (8,991 bp), the nucleotide diversity (pi) was 0.058% among all sequences, 0.076% among the African sequences, 0.047% among the Asian sequences, and 0.045% among the European sequences. A compilation of data revealed that autosomal regions have, on average, the highest pi value (0.091%), X-linked regions have a somewhat lower pi value (0.079%), and Y-linked regions have a very low pi value (0.008%). The lower polymorphism in the present region may be due to a lower mutation rate and/or selection in the gene containing these introns or in genes linked to this region. The present region and two other 10-kb noncoding regions all show a strong excess of low-frequency variants, indicating a relatively recent population expansion. This region has a low mutation rate, which was estimated to be 0.74 x 10 per nucleotide per year. An average estimate of approximately 12,600 for the long-term effective population size was obtained using various methods; the estimate was not far from the commonly used value of 10,000. Fu and Li's tests rejected the assumption of an equilibrium neutral Wright-Fisher population, largely owing to the high proportion of low-frequency variants. The age of the most recent common ancestor of the sequences in our sample was estimated to be more than 1 Myr. Allowing for some unrealistic assumptions in the model, this estimate would still suggest an age of more than 500,000 years, providing further evidence for a genetic history of humans much more ancient than the emergence of modern humans. The fact that many unique variants exist in Europe and Asia also suggests a fairly long genetic history outside of Africa and argues against a complete replacement of all indigenous populations in Europe and Asia by a small Africa stock. Moreover, the ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. We suggest that both the "Out of Africa" and the multiregional models are too simple to explain the evolution of modern humans.

Worldwide DNA sequence variation in a 10-kilobase noncoding region on human chromosome 22

Human DNA sequence variation data are useful for studying the origin, evolution, and demographic history of modern humans and the mechanisms of maintenance of genetic variability in human populations, and for detecting linkage association of disease. Here, we report worldwide variation data from a approximately 10-kilobase noncoding autosomal region. We identified 75 variant sites in 64 humans (128 sequences) and 463 variant sites among the human, chimpanzee, and orangutan sequences. Statistical tests suggested that the region is selectively neutral. The average nucleotide diversity (pi) across the region was 0.088% among all of the human sequences obtained, 0.085% among African sequences, and 0.082% among non-African sequences, supporting the view of a low nucleotide diversity ( approximately 0.1%) in humans. The comparable pi value in non-Africans to that in Africans indicates no severe bottleneck during the evolution of modern non-Africans; however, the possibility of a mild bottleneck cannot be excluded because non-Africans showed considerably fewer variants than Africans. The present and two previous large data sets all show a strong excess of low frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The mutation rate was estimated to be 1.15 x 10(-9) per nucleotide per year. Estimates of the long-term effective population size N(e) by various statistical methods were similar to those in other studies. The age of the most recent common ancestor was estimated to be approximately 1.29 million years ago among all of the sequences obtained and approximately 634,000 years ago among the non-African sequences, providing the first evidence from a noncoding autosomal region for ancient human histories, even among non-Africans.