Accurate and scalable identification of functional sites by evolutionary tracing

A common difficulty in post genomics biology is that large-scale techniques of data collection often strip away information on the biological context of these data. The result is a massive number of disconnected observations on sequence, structure, and function from which underlying patterns and biological meaning are obscured. One solution is to build computational filters that pick out sufficiently few facts, relevant to a query, that their relationship is immediately apparent and experimentally testable. Typically, these filters rely on mathematics and statistics, and on first principles from physics and chemistry. We show here that evolution itself can be used to filter sequence and structure data in order to identify evolutionarily important amino acids. A general property of these residues is that they form clusters in native protein structures and point to regions where mutations have the greatest biological impact. The result is an accurate method of functional site annotation that is scalable for structural proteomics.     

Genomics of Actinobacteria: tracing the evolutionary history of an ancient phylum.

Actinobacteria constitute one of the largest phyla among bacteria and represent gram-positive bacteria with a high G+C content in their DNA. This bacterial group includes microorganisms exhibiting a wide spectrum of morphologies, from coccoid to fragmenting hyphal forms, as well as possessing highly variable physiological and metabolic properties. Furthermore, Actinobacteria members have adopted different lifestyles, and can be pathogens (e.g., Corynebacterium, Mycobacterium, Nocardia, Tropheryma, and Propionibacterium), soil inhabitants (Streptomyces), plant commensals (Leifsonia), or gastrointestinal commensals (Bifidobacterium). The divergence of Actinobacteria from other bacteria is ancient, making it impossible to identify the phylogenetically closest bacterial group to Actinobacteria. Genome sequence analysis has revolutionized every aspect of bacterial biology by enhancing the understanding of the genetics, physiology, and evolutionary development of bacteria. Various actinobacterial genomes have been sequenced, revealing a wide genomic heterogeneity probably as a reflection of their biodiversity. This review provides an account of the recent explosion of actinobacterial genomics data and an attempt to place this in a biological and evolutionary context.     

The ribonuclease A superfamily of mammals and birds: identifying new members and tracing evolutionary histories

The RNase A superfamily has been important in biochemical, structural, and evolutionary studies and is believed to be the sole vertebrate-specific enzyme family. To understand the origin and diversification of the superfamily, we here determine its entire repertoire in the sequenced genomes of human, mouse, rat, and chicken. We report a previously unnoticed gene cluster in mouse chromosome 10 and a number of new genes, including mammalian RNases 11-13, which are close relatives of the recently identified RNases 9 and 10. Gene expression data imply male-reproductive functions for RNases 9-13, although their sequences suggest the lack of ribonucleolytic activities. In contrast to the presence of 13-20 functional genes in mammals, chicken has only 3 RNase genes, which are evolutionarily close to mammalian RNase 5, like other nonmammalian RNases. This and other evidence suggests that the RNase A superfamily originated from an RNase 5-like gene and expanded in mammals. Together with the fact that multiple lineages of the superfamily, including RNases 2, 3, 5, and 7, have antipathogenic activities, we suggest that the superfamily started off as a host-defense mechanism in vertebrates. Consistent with this hypothesis, all members of the superfamily exhibit high rates of amino acid substitution as is commonly observed in immunity genes.     

Ecological,historical and evolutionary determinants of modularity in weighted seed‐dispersal networks

Modularity is a recurrent and important property of bipartite ecological networks. Although well‐resolved ecological networks describe interaction frequencies between species pairs, modularity of bipartite networks has been analysed only on the basis of binary presence–absence data. We employ a new algorithm to detect modularity in weighted bipartite networks in a global analysis of avian seed‐dispersal networks. We define roles of species, such as connector values, for weighted and binary networks and associate them with avian species traits and phylogeny. The weighted, but not binary, analysis identified a positive relationship between climatic seasonality and modularity, whereas past climate stability and phylogenetic signal were only weakly related to modularity. Connector values were associated with foraging behaviour and were phylogenetically conserved. The weighted modularity analysis demonstrates the dominating impact of ecological factors on the structure of seed‐dispersal networks, but also underscores the relevance of evolutionary history in shaping species roles in ecological communities.     

Identification of monomorphic and divergent haplotypes in the 2006-2007 norovirus GII/4 epidemic population by genomewide tracing of evolutionary history

Our norovirus (NoV) surveillance group reported a >4-fold increase in NoV infection in Japan during the winter of 2006-2007 compared to the previous winter. Because the increase was not linked to changes in the surveillance system, we suspected the emergence of new NoV GII/4 epidemic variants. To obtain information on viral changes, we conducted full-length genomic analysis. Stool specimens from 55 acute gastroenteritis patients of various ages were collected at 11 sites in Japan between May 2006 and January 2007. Direct sequencing of long PCR products revealed 37 GII/4 genome sequences. Phylogenetic study of viral genome and partial sequences showed that the two new GII/4 variants in Europe, termed 2006a and 2006b, initially coexisted as minorities in early 2006 in Japan and that 2006b alone had dominated over the resident GII/4 variants during 2006. A combination of phylogenetic and entropy analyses revealed for the first time the unique amino acid substitutions in all eight proteins of the new epidemic strains. These data and computer-assisted structural study of the NoV capsid protein are compatible with a model of antigenic drift with tuning of the structure and functions of multiple proteins for the global outgrowth of new GII/4 variants. The availability of comprehensive information on genome sequences and unique protein changes of the recent global epidemic variants will allow studies of diagnostic assays, molecular epidemiology, molecular biology, and adaptive changes of NoV in nature.     

Haseman-Elston weighted by marker informativity

In the Haseman-Elston approach the squared phenotypic difference is regressed on the proportion of alleles shared identical by descent (IBD) to map a quantitative trait to a genetic marker. In applications the IBD distribution is estimated and usually cannot be determined uniquely owing to incomplete marker information. At Genetic Analysis Workshop (GAW) 13, Jacobs et al. [BMC Genet 2003, 4(Suppl 1):S82] proposed to improve the power of the Haseman-Elston algorithm by weighting for information available from marker genotypes. The authors did not show, however, the validity of the employed asymptotic distribution. In this paper, we use the simulated data provided for GAW 14 and show that weighting Haseman-Elston by marker information results in increased type I error rates. Specifically, we demonstrate that the number of significant findings throughout the chromosome is significantly increased with weighting schemes. Furthermore, we show that the classical Haseman-Elston method keeps its nominal significance level when applied to the same data. We therefore recommend to use Haseman-Elston with marker informativity weights only in conjunction with empirical p-values. Whether this approach in fact yields an increase in power needs to be investigated further.     

Oligomers of ERBB3 have two distinct interfaces that differ in their sensitivity to disruption by heregulin

ErbB receptors associate in a ligand-dependent or -independent manner, and overexpression of epidermal growth factor receptor (ErbB1) or ErbB2 results in ligand-independent activation. Ligand-independent activation is poorly understood, and dimerization alone is not sufficient for activation. ErbB receptors also form higher order oligomers, but the mechanism of oligomer formation and their contribution to signaling are not known. The kinase-deficient ErbB3 as well as its extracellular domains are particularly prone to ligand-independent oligomerization, and oligomers are destabilized by binding of the ligand heregulin. In contrast, ligand binding facilitates heterodimerization with ErbB2 and is expected to stabilize an extended conformation of the ErbB3 extracellular domain (ECD) in which the dimerization interface is exposed. In the absence of ligand, ErbB3 can adopt a closed conformation that is held together by an intramolecular tether. We used a constitutively extended form of the ErbB3-ECD to analyze the conformation of the ECD in oligomers and the mechanism of oligomer disruption by heregulin. The extended conformation of the ECD forms oligomers more readily, suggesting the crystallographically defined dimer interface is one of the interfaces involved in oligomerization. Heregulin destabilizes oligomeric complexes but not dimers, which are neither stabilized nor disrupted by ligand binding, indicating a distinct second interface in oligomers of ErbB3. Cross-linking and activation studies on membrane-embedded ErbB3/ErbB2 chimeras confirm this dual effect of heregulin. Most of the ErbB3-ECD on the cell surface is apparently kept in an open conformation through oligomerization, and the resulting oligomers adopt a conformation representing a state of reduced activity.     

Identification of a heregulin binding site in HER3 extracellular domain

HER3 (also known as c-Erb-b3) is a type I receptor tyrosine kinase similar in sequence to the epidermal growth factor (EGF) receptor. The extracellular segment of this transmembrane receptor contains four domains. Domains I and II are similar in sequence to domains III and IV, respectively, and domains II and IV are cysteine-rich. We show that the EGF-like domain of heregulin (hrg) binds to domains I and II of HER3, in contrast to the EGF receptor, for which prior studies have shown that a construct consisting of domains III and portions of domain IV binds EGF. Next, we identified a putative hrg binding site by limited proteolysis of the recombinant extracellular domains of HER3 (HER3-ECD(I-IV)) in both the presence and absence of hrg. In the absence of hrg, HER3-ECD(I-IV) is cleaved after position Tyr(50), near the beginning of domain I. Binding of hrg to HER3-ECD(I-IV) fully protects position Tyr(50) from proteolysis. To confirm that domain I contains a hrg binding site, we expressed domains I and II (HER3-ECD(I-II)) and find that it binds hrg with 68 nm affinity. These data suggest that domains I and II of HER3-ECD(I-IV) act as a functional unit in folding and binding of hrg. Thus, our biochemical findings reinforce the structural hypothesis of others that HER3-ECD(I-IV) is similar to the insulin-like growth factor-1 receptor (IGF-1R), as follows: 1) The protected cleavage site in HER3-ECD(I-IV) corresponds to a binding footprint in domain I of IGF-1R; 2) HER3-ECD(I-II) binds hrg with a 68 nm dissociation constant, supporting the hypothesis that domain I is involved in ligand binding; and 3) the large accessible surface area (1749 A) of domain L1 of IGF-1R that is buried by domain S1, as well as the presence of conserved contacts in this interface of type 1 RTKs, suggests that domains L1 and S1 of IGF-1R function as a unit as observed for HER3-ECD(I-II). Our results are consistent with the proposal that HER3 has a structure similar to IGF-1R and binds ligand at a site in corresponding domains.     

A comparison of two methods for constructing evolutionary distances from a weighted contribution of transition and transversion differences

Since the initial work of Jukes and Cantor (1969), a number of procedures have been developed to estimate the expected number of nucleotide substitutions corresponding to a given observed level of nucleotide differentiation assuming particular evolutionary models. Unlike the proportion of different sites, the expected number of substitutions that would have occurred grows linearly with time and therefore has had great appeal as an evolutionary distance. Recently, however, a number of authors have tried to develop improved statistical approaches for generating and evaluating evolutionary distances (Schoniger and von Haeseler 1993; Goldstein and Polock 1994; Tajima and Takezaki 1994). These studies clearly show that the estimated number of nucleotide substitutions is generally not the best estimator for use in reconstruction of phylogenetic relationships. The reason for this is that there is often a large error associated with the estimation of this number. Therefore, even though its expectation is correct (i.e., on average the expected number of substitutions is proportional to time- -but see Tajima 1993), it is not expected to be as useful as estimators designed to have a lower variance.      

Differential regulation of components of the focal adhesion complex by heregulin: role of phosphatase SHP-2.

Heregulin (HRG) has been implicated in the progression of breast cancer cells to a malignant phenotype, a process that involves changes in cell motility and adhesion. Here we demonstrate that HRG differentially regulates the site-specific phosphorylation of the focal adhesion components focal adhesion kinase (FAK) and paxilin in a dose-dependent manner. HRG at suboptimal doses (0.01 and 0.1 nM) increased adhesion of cells to the substratum, induced phosphorylation of FAK at Tyr-577, -925, and induced formation of well-defined focal points in breast cancer cell line MCF-7. HRG at a dose of 1 nM, increased migratory potential of breast cancer cells, selectively dephosphorylated FAK at Tyr-577, -925, and paxillin at Tyr-31. Tyrosine phosphorylation of FAK at Tyr-397 remained unaffected by HRG stimulation. FAK associated with HER2 only in response to 0.01 nM HRG. In contrast, 1 nM HRG induced activation and increased association of tyrosine phosphatase SHP-2 with HER2 but decreased association of HER2 with FAK. Expression of dominant-negative SHP-2 blocked HRG-mediated dephosphorylation of FAK and paxillin, leading to persistent accumulation of mature focal points. Our results suggest that HRG differentially regulates signaling from focal adhesion complexes through selective phosphorylation and dephosphorylation and that tyrosine phosphatase SHP-2 has a role in the HRG signaling.     

Analysis of bilateral inverse symmetry in whole bacterial chromosomes

The positions of the 64 DNA tri-nucleotides (triplets) along the Borrelia burgdorferi chromosome were determined and cumulative position plots (CPP) were obtained. Analysis of CPP for complementary triplets revealed close correlations in complementary triplet frequencies (CTF) between opposing leading and lagging strands. Such bilateral inverse symmetry (BIS) applied also to complementary mono- and di-nucleotides and to some >3 n-tuples. At the level of individual bases BIS explains Chargaff's second parity rule for whole bacterial chromosomes. Using shuffled control sequences we show that single-base BIS was not the source of higher-order BIS. Analysis of CTF in 45 other chromosomes suggests that BIS is a general property of eubacteria. BIS at the various levels may be due to the very similar numbers of codons used in chromosomal halves. Evolutionarily, BIS could have resulted from asymmetric substitution of bases combined with genetic rearrangements. However, the provocative theoretical alternative of whole-genome inverse duplication is here considered.     

Right-left symmetry preservation by flowers of malvaceae family

The flowers of malvaceae family preserves the symmetry between right and left in a peculiar manner. Plots belonging to this family bear two kinds of flowers, right-handed flowers with anticlockwise twisted petals and left-handed flowers with clockwise twisted petals. The branches of the plant prefers production of one type of flowers in excess of the other. There are two distinct types of branches, dextral branches and sinistral branches. Dextral (sinistral) branches produce more right-handed (left-handed) flowers than left-handed (right-handed) flowers. The average percentage of right-handed flowers in a dextral branch is same as that of left-handed flowers in a sinistral branch.     

Scaffold-mediated symmetry breaking by Cdc42p

Cell polarization generally occurs along a single well-defined axis that is frequently determined by environmental cues such as chemoattractant gradients or cell-cell contacts, but polarization can also occur spontaneously in the apparent absence of such cues, through a process called symmetry breaking. In Saccharomyces cerevisiae, cells are born with positional landmarks that mark the poles of the cell and guide subsequent polarization and bud emergence to those sites, but cells lacking such landmarks polarize towards a random cortical site and proliferate normally. The landmarks employ a Ras-family GTPase, Rsr1p, to communicate with the conserved Rho-family GTPase Cdc42p, which is itself polarized and essential for cytoskeletal polarization. We found that yeast Cdc42p was effectively polarized to a single random cortical site even in the combined absence of landmarks, microtubules and microfilaments. Among a panel of Cdc42p effectors and interacting proteins, we found that the scaffold protein Bem1p was uniquely required for this symmetry-breaking behaviour. Moreover, polarization was dependent on GTP hydrolysis by Cdc42p, suggesting that assembly of a polarization site involves cycling of Cdc42p between GTP- and GDP-bound forms, rather than functioning as a simple on/off switch.     

环节动物和软体动物之间进化关系分析

根据国内外相关文献,分析了环节动物和软体动物之间的进化关系的不同观点及其立论依据,提出就现有证据而言,软体动物起源于环节动物更具说服力,并对软体动物和环节动物进化关系的进一步研究提出设想.     

Haplotype fine mapping by evolutionary trees

To refine the location of a disease gene within the bounds provided by linkage analysis, many scientists use the pattern of linkage disequilibrium between the disease allele and alleles at nearby markers. We describe a method that seeks to refine location by analysis of "disease" and "normal" haplotypes, thereby using multivariate information about linkage disequilibrium. Under the assumption that the disease mutation occurs in a specific gap between adjacent markers, the method first combines parsimony and likelihood to build an evolutionary tree of disease haplotypes, with each node (haplotype) separated, by a single mutational or recombinational step, from its parent. If required, latent nodes (unobserved haplotypes) are incorporated to complete the tree. Once the tree is built, its likelihood is computed from probabilities of mutation and recombination. When each gap between adjacent markers is evaluated in this fashion and these results are combined with prior information, they yield a posterior probability distribution to guide the search for the disease mutation. We show, by evolutionary simulations, that an implementation of these methods, called "FineMap," yields substantial refinement and excellent coverage for the true location of the disease mutation. Moreover, by analysis of hereditary hemochromatosis haplotypes, we show that FineMap can be robust to genetic heterogeneity.     

Generation of evolutionary novelty by functional shift

That biological features may change their function during evolution has long been recognized. Particularly, the acquisition of new functions by molecules involved in developmental pathways is suspected to cause important morphologic novelties. However, the current terminology describing functional changes during evolution (co-option or recruitment) fails to recognize important biologic distinctions between diverse evolutionary routes involving functional shifts. The main goal of our work is to stress the importance of an apparently trivial distinction: Whether or not the element that adopts a new function (anything from a morphologic structure to a protein domain) is a single or a duplicated element. We propose that natural selection must act in a radically different way, depending on the historic succession of co-option and duplication events; that is, co-option may provide the selective pressure for a subsequent gene duplication or could be a stabilizing factor that helps maintain redundancy after gene duplication. We review the evidence available on functional changes, focusing whenever possible on developmental molecules, and we propose a conceptual framework for the study of functional shifts during evolution with a level of resolution appropriate to the power of our current methodologies. BioEssays 21:432–439, 1999. © 1999 John Wiley & Sons, Inc.