Construction of phylogenetic trees for proteins and nucleic acids: empirical evaluation of alternative matrix methods

Visualization of protein sequence space with force-directed graphs, and their application to the choice of target-template pairs for homology modelling

The protein sequence-structure gap results from the contrast between rapid, low-cost deep sequencing, and slow, expensive experimental structure determination techniques. Comparative homology modelling may have the potential to close this gap by predicting protein structure in target sequences using existing experimentally solved structures as templates. This paper presents the first use of force-directed graphs for the visualization of sequence space in two dimensions, and applies them to the choice of suitable RNA-dependent RNA polymerase (RdRP) target-template pairs within human-infective RNA virus genera. Measures of centrality in protein sequence space for each genus were also derived and used to identify centroid nearest-neighbour sequences (CNNs) potentially useful for production of homology models most representative of their genera. Homology modelling was then carried out for target-template pairs in different species, different genera and different families, and model quality assessed using several metrics. Reconstructed ancestral RdRP sequences for individual genera were also used as templates for the production of ancestral RdRP homology models. High quality ancestral RdRP models were consistently produced, as were good quality models for target-template pairs in the same genus. Homology modelling between genera in the same family produced mixed results and inter-family modelling was unreliable. We present a protocol for the production of optimal RdRP homology models for use in further experiments, e.g. docking to discover novel anti-viral compounds. (219 words)

•

The first use of force-directed graphs for the visualization of multidimensional protein sequence space in two dimensions.

•

Measures of centrality in protein sequence space to identify sequences for production of homology models.

•

Homology modelling for RNA-dependent RNA polymerase (RdRP) target-template pairs in different species, genera and families.

•

A protocol for the production of optimal RdRP homology models for use in further experiments.

Fosfomycin Protects Mice From Staphylococcus aureus Pneumonia Caused by α-Hemolysin in Extracellular Vesicles by Inhibiting MAPK-Regulated NLRP3 Inflammasomes

α-Hemolysin (Hla) is a significant virulence factor in Staphylococcus aureus (S. aureus)-caused infectious diseases such as pneumonia. Thus, to prevent the production of Hla when treating S. aureus infection, it is necessary to choose an antibiotic with good antibacterial activity and effect. In our study, we observed that Fosfomycin (FOM) at a sub-inhibitory concentration inhibited expression of Hla. Molecular dynamics demonstrated that FOM bound to the binding sites LYS 154 and ASP 108 of Hla, potentially inhibiting Hla. Furthermore, we verified that staphylococcal membrane-derived vesicles (SMVs) contain Hla and that FOM treatment significantly reduced the production of SMVs and Hla. Based on our pharmacological inhibition analysis, ERK and p38 activated NLRP3 inflammasomes. Moreover, FOM inhibited expression of MAPKs and NLRP3 inflammasome-related proteins in S. aureus as well as SMV-infected human macrophages (MΦ) and alveolar epithelial cells. In vivo, SMVs isolated from S. aureus DU1090 (an isogenic Hla deletion mutant) or the strain itself caused weaker inflammation than that of its parent strain 8325-4. FOM also significantly reduced the phosphorylation levels of ERK and P38 and expression of NLRP3 inflammasome-related proteins. In addition, FOM decreased MPO activity, pulmonary vascular permeability and edema formation in the lungs of mice with S. aureus-caused pneumonia. Taken together, these data indicate that FOM exerts protective effects against S. aureus infection in vitro and in vivo by inhibiting Hla in SMVs and blocking ERK/P38-mediated NLRP3 inflammasome activation by Hla.

Characterization and validation of somatic mutation spectrum to reveal heterogeneity in gastric cancer by single cell sequencing

Gastric cancer (GC) is a highly heterogeneous disease with multiple cellular types and poor prognosis. However, the cellular evolution and molecular basis of GC at the individual intra-tumor level has not been well demonstrated. We performed single-cell whole exome sequencing to detect somatic single-nucleotide variants (SNVs) and significantly mutated genes (SMGs) among 34 tumor cells and 9 normal cells from a patient with GC. The Complete Prediction for Protein Conformation (CPPC) approach directly predicting the folding conformation of the protein 3D structure with Protein Folding Shape Code, combined with functional experiments were used to confirm the characterization of mutated SMGs in GC cells. We identified 201 somatic SNVs, including 117 non-synonymous mutations in GC cells. Further analysis identified 24 significant mutated genes (SMGs) in single cells, for which a single amino acid change might affect protein conformation. Among them, two genes (CDC27 and FLG) that were mutated only in single cells but not in the corresponding tumor tissue, were recurrently present in another GC tissue cohort, and may play a potential role to promote carcinogenesis, as confirmed by functional characterization. Our findings showed a mutational landscape of GC at intra-tumor level for the first time and provided opportunities for understanding the heterogeneity and individualized target therapy for this disease.

Exploring the interaction between human focal adhesion kinase and inhibitors: a molecular dynamic simulation and free energy calculations

Focal adhesion kinase is an important target for the treatment of many kinds of cancers. Inhibitors of FAK are proposed to be the anticancer agents for multiple tumors. The interaction characteristic between FAK and its inhibitors is crucial to develop new inhibitors. In the present article, we used Molecular Dynamic (MD) simulation method to explore the characteristic of interaction between FAK and three inhibitors (PHM16, TAE226, and ligand3). The MD simulation results together with MM-GB/SA calculations show that the combinations are enthalpy-driven process. Cys502 and Asp564 are both essential residues due to the hydrogen bond interactions with inhibitors, which was in good agreement with experimental data. Glu500 can form a non-classical hydrogen bond with each inhibitor. Arg426 can form electrostatic interactions with PHM16 and ligand3, while weaker with TAE226. The electronic static potential was employed, and we found that the ortho-position methoxy of TAE226 has a weaker negative charge than the meta-position one in PHM16 or ligand3. Ile428, Val436, Ala452, Val484, Leu501, Glu505, Glu506, Leu553, Gly563 Leu567, Ser568 are all crucial residues in hydrophobic interactions. The key residues in this work will be available for further inhibitor design of FAK and also give assistance to further research of cancer.

Structural analysis of in silico mutant experiments of human inner-kinetochore structure

Large multi-molecular complexes like the kinetochore are lacking of suitable methods to determine their spatial structure. Here, we use and evaluate a novel modeling approach that combines rule-bases reaction network models with spatial molecular geometries. In particular, we introduce a method that allows to study in silico the influence of single interactions (e.g. bonds) on the spatial organization of large multi-molecular complexes and apply this method to an extended model of the human inner-kinetochore. Our computational analysis method encompasses determination of bond frequency, geometrical distances, statistical moments, and inter-dependencies between bonds using mutual information. For the analysis we have extend our previously reported human inner-kinetochore model by adding 13 new protein interactions and three protein geometry details. The model is validated by comparing the results of in silico with reported in vitro single protein deletion experiments. Our studies revealed that most simulations mimic the in vitro behavior of the kinetochore complex as expected. To identify the most important bonds in this model, we have created 39 mutants in silico by selectively disabling single protein interactions. In a total of 11,800 simulation runs we have compared the resulting structures to the wild-type. In particular, this allowed us to identify the interaction Cenp-W-H3 and Cenp-S-Cenp-X as having the strongest influence on the inner-kinetochore's structure. We conclude that our approach can become a useful tool for the in silico dynamical study of large, multi-molecular complexes.

Discovery of biclonal origin and a novel oncogene SLC12A5 in colon cancer by single-cell sequencing

Single-cell sequencing is a powerful tool for delineating clonal relationship and identifying key driver genes for personalized cancer management. Here we performed single-cell sequencing analysis of a case of colon cancer. Population genetics analyses identified two independent clones in tumor cell population. The major tumor clone harbored APC and TP53 mutations as early oncogenic events, whereas the minor clone contained preponderant CDC27 and PABPC1 mutations. The absence of APC and TP53 mutations in the minor clone supports that these two clones were derived from two cellular origins. Examination of somatic mutation allele frequency spectra of additional 21 whole-tissue exome-sequenced cases revealed the heterogeneity of clonal origins in colon cancer. Next, we identified a mutated gene SLC12A5 that showed a high frequency of mutation at the single-cell level but exhibited low prevalence at the population level. Functional characterization of mutant SLC12A5 revealed its potential oncogenic effect in colon cancer. Our study provides the first exome-wide evidence at single-cell level supporting that colon cancer could be of a biclonal origin, and suggests that low-prevalence mutations in a cohort may also play important protumorigenic roles at the individual level.

Phylogeographic population structure of Red-winged Blackbirds assessed by mitochondrial DNA

A continent-wide survey of restriction-site variation in mitochondrial DNA (mtDNA) of the Red-winged Blackbird (Agelaius phoeniceus) was conducted to assess the magnitude of phylogeographic population structure in an avian species. A total of 34 mtDNA genotypes was observed among the 127 specimens assayed by 18 restriction endonucleases. Nonetheless, population differentiation was minor, as indicated by (i) small genetic distances in terms of base substitutions per nucleotide site between mtDNA genotypes (maximum P approximately 0.008) and by (ii) the widespread geographic distributions of particular mtDNA clones and phylogenetic arrays of clones. Extensive morphological differentiation among redwing populations apparently has occurred in the context of relatively little phylogenetic separation. A comparison between mtDNA data sets for Red-winged Blackbirds and deermice (Peromyscus maniculatus) also sampled from across North America shows that intraspecific population structures of these two species differ dramatically. The lower phylogeographic differentiation in redwings is probably due to historically higher levels of gene flow.

GIGA: a simple, efficient algorithm for gene tree inference in the genomic age

Background

Phylogenetic relationships between genes are not only of theoretical interest: they enable us to learn about human genes through the experimental work on their relatives in numerous model organisms from bacteria to fruit flies and mice. Yet the most commonly used computational algorithms for reconstructing gene trees can be inaccurate for numerous reasons, both algorithmic and biological. Additional information beyond gene sequence data has been shown to improve the accuracy of reconstructions, though at great computational cost.

Results

We describe a simple, fast algorithm for inferring gene phylogenies, which makes use of information that was not available prior to the genomic age: namely, a reliable species tree spanning much of the tree of life, and knowledge of the complete complement of genes in a species' genome. The algorithm, called GIGA, constructs trees agglomeratively from a distance matrix representation of sequences, using simple rules to incorporate this genomic age information. GIGA makes use of a novel conceptualization of gene trees as being composed of orthologous subtrees (containing only speciation events), which are joined by other evolutionary events such as gene duplication or horizontal gene transfer. An important innovation in GIGA is that, at every step in the agglomeration process, the tree is interpreted/reinterpreted in terms of the evolutionary events that created it. Remarkably, GIGA performs well even when using a very simple distance metric (pairwise sequence differences) and no distance averaging over clades during the tree construction process.

Conclusions

GIGA is efficient, allowing phylogenetic reconstruction of very large gene families and determination of orthologs on a large scale. It is exceptionally robust to adding more gene sequences, opening up the possibility of creating stable identifiers for referring to not only extant genes, but also their common ancestors. We compared trees produced by GIGA to those in the TreeFam database, and they were very similar in general, with most differences likely due to poor alignment quality. However, some remaining differences are algorithmic, and can be explained by the fact that GIGA tends to put a larger emphasis on minimizing gene duplication and deletion events.

DETECTING EVOLUTIONARY TRANSFER OF GENES USING PhIGs(1)

Organisms have acquired plastids by convoluted paths that have provided multiple opportunities for gene transfer into a host nucleus from intracellular organisms, including the cyanobacterial ancestor of plastids, the proteobacterial ancestor of mitochondria, and both green and red algae whose engulfment has led to secondary acquisition of plastids. These gene movements are most accurately demonstrated by building phylogenetic trees that identify the evolutionary origin of each gene, and one effective tool for this is "PhIGs" (Phylogenetically Inferred Groups; http://PhIGs.org), a set of databases and computer tools with a Web interface for whole-genome evolutionary analysis. PhIGs takes as input gene sets of completely sequenced genomes, builds clusters of genes using a novel, graph-based approach, and reconstructs the evolutionary relationships among all gene families. The user can view and download the sequence alignments, compare intron-exon structures, and follow links to functional genomic databases. Currently, PhIGs contains 652,756 genes from 45 genomes grouped into 61,059 gene families. Graphical displays show the relative positions of these genes among genomes. PhIGs has been used to detect the evolutionary transfer of hundreds of genes from cyanobacteria and red algae into oömycete nuclear genomes, revealing that even though they have no plastids, their ancestors did, having secondarily acquired them from an intracellular red alga. A great number of genomes are soon to become available that are relevant to our broader understanding of the movement of genes among intracellular compartments after engulfing other organisms, and PhIGs will be an effective tool to interpret these gene movements.

Full-Length Sequences of Two Hepatitis E Virus Isolates Representing an Eastern China-Indigenous Subgroup of Genotype 4

Although the majority of hepatitis E virus (HEV) infections in Japan are 'domestic' due to the presence of indigenous strains, there are still 'imported' cases as well. Among 83 patients with non-A, non-B and non-C acute liver diseases admitted to Saitama Medical University Hospital, 7 (8.4%) were positive for serum HEV-RNA, of whom 2 had a recent history of traveling to China, one to Xian and another to Shanghai. We determined the full-genome sequences of HEV from these 2 patients (isolate names are JKO-ChiSai98c and JYI-ChiSai01c, genotype 4 in both) for phylogenetic analyses. Initially, when compared only to the 13 full-genome sequences of genotype 4 so far reported, our 2 isolates were thought to be novel strains because they showed a significant genetic difference from the sequences known to date. However, when we included a set of short sequences (150 nt) recently reported from China in the comparison, we found that our 2 isolates represent a subgroup of genotype 4, which seems to be restricted to eastern China. In conclusion, the 2 HEV isolates reported here could serve as full-genome prototypes for an eastern China-indigenous subgroup of the genotype 4 HEV.

Multilocus sequence typing of bacteria

Multilocus sequence typing (MLST) was proposed in 1998 as a portable, universal, and definitive method for characterizing bacteria, using the human pathogen Neisseria meningitidis as an example. In addition to providing a standardized approach to data collection, by examining the nucleotide sequences of multiple loci encoding housekeeping genes, or fragments of them, MLST data are made freely available over the Internet to ensure that a uniform nomenclature is readily available to all those interested in categorizing bacteria. At the time of writing, over thirty MLST schemes have been published and made available on the Internet, mostly for pathogenic bacteria, although there are schemes for pathogenic fungi and some nonpathogenic bacteria. MLST data have been employed in epidemiological investigations of various scales and in studies of the population biology, pathogenicity, and evolution of bacteria. The increasing speed and reduced cost of nucleotide sequence determination, together with improved web-based databases and analysis tools, present the prospect of increasingly wide application of MLST.

Evidence from nuclear DNA sequences sheds light on the phylogenetic relationships of Pinnipedia: single origin with affinity to Musteloidea

Considerable long-standing controversy and confusion surround the phylogenetic affinities of pinnipeds, the largely marine group of "fin-footed" members of the placental mammalian order Carnivora. Until most recently, the two major competing hypotheses were that the pinnipeds have a single (monophyletic) origin from a bear-like ancestor, or that they have a dual (diphyletic) origin, with sea lions (Otariidae) derived from a bear-like ancestor, and seals (Phocidae) derived from an otter-, mustelid-, or musteloid-like ancestor. We examined phylogenetic relationships among 29 species of arctoid carnivorans using a concatenated sequence of 3228 bp from three nuclear loci (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). The species represented Pinnipedia (Otariidae: Callorhinus, Eumetopias; Phocidae: Phoca), bears (Ursidae: Ursus, Melursus), and Musteloidea (Mustelidae: Mustela, Enhydra, Melogale, Martes, Gulo, Meles; Procyonidae: Procyon; Ailuridae: Ailurus; Mephitidae: Mephitis). Maximum parsimony, maximum likelihood, and Bayesian inference phylogenetic analyses of separate and combined datasets produced trees with largely congruent topologies. The analyses of the combined dataset resulted in well-resolved and well-supported phylogeny reconstructions. Evidence from nuclear DNA evolution presented here contradicts the two major hypotheses of pinniped relationships and strongly suggests a single origin of the pinnipeds from an arctoid ancestor shared with Musteloidea to the exclusion of Ursidae.

Identification and characterization of subfamily-specific signatures in a large protein superfamily by a hidden Markov model approach

BACKGROUND

Most profile and motif databases strive to classify protein sequences into a broad spectrum of protein families. The next step of such database studies should include the development of classification systems capable of distinguishing between subfamilies within a structurally and functionally diverse superfamily. This would be helpful in elucidating sequence-structure-function relationships of proteins.

RESULTS

Here, we present a method to diagnose sequences into subfamilies by employing hidden Markov models (HMMs) to find windows of residues that are distinct among subfamilies (called signatures). The method starts with a multiple sequence alignment (MSA) of the subfamily. Then, we build a HMM database representing all sliding windows of the MSA of a fixed size. Finally, we construct a HMM histogram of the matches of each sliding window in the entire superfamily. To illustrate the efficacy of the method, we have applied the analysis to find subfamily signatures in two well-studied superfamilies: the cadherin and the EF-hand protein superfamilies. As a corollary, the HMM histograms of the analyzed subfamilies revealed information about their Ca2+ binding sites and loops.

CONCLUSIONS

The method is used to create HMM databases to diagnose subfamilies of protein superfamilies that complement broad profile and motif databases such as BLOCKS, PROSITE, Pfam, SMART, PRINTS and InterPro.

Sequence of mdr3 cDNA encoding a human P-glycoprotein

We have determined the sequence of the human mdr3 gene using cDNA derived from liver RNA. The mdr3 gene codes for a member of a family of membrane proteins, the P-glycoproteins, overproduced in many multi-drug-resistant (MDR) cell lines. Like its relatives, the protein encoded by mdr3 has a deduced Mr of 140,000, which is presumably increased by glycosylation after synthesis. The sequence consists of two similar halves, each with a series of six hydrophobic segments that may form a membrane channel. The halves also possess nucleotide-binding consensus sequences, which presumably act as ATPases and drive drug transport. The presumed ATPase domains are all but identical to those of the human mdr1 gene product [Chen et al., Cell 47 (1986) 381-389]. We attribute this high level of sequence conservation to the repeated gene conversion that is evident from segments in which mdr1 and mdr3 differ only in a few silent mutations. Divergence between P-glycoprotein family members is greatest at the N terminus and in the 60 amino acid linker connecting the two halves. In the putative trans-membrane domains approx. 80% of the amino acids are conserved between the products of mdr1 and mdr3. Although the function of mdr3 is not yet known, its high homology with mdr1 suggests that it also encodes an efflux pump with broad specificity.

Determining evolutionary distances from highly diverged nucleic acid sequences: operator metrics

Operator metrics are explicitly designed to measure evolutionary distances from nucleic acid sequences when substitution rates differ greatly among the organisms being compared, or when substitutions have been extensive. Unlike lengths calculated by the distance matrix and parsimony methods, in which substitutions in one branch of a tree can alter the measured length of another branch, lengths determined by operator metrics are not affected by substitutions outside the branch. In the method, lengths (operator metrics) corresponding to each of the branches of an unrooted tree are calculated. The metric length of a branch reconstructs the number of (transversion) differences between sequences at a tip and a node (or between nodes) of a tree. The theory is general and is fundamentally independent of differences in substitution rates among the organisms being compared. Mathematically, the independence has been obtained because the metrics are eigenvectors of fundamental equations which describe the evolution of all unrooted trees. Even under conditions when both the distance matrix method or a simple parismony length method are shown to indicate lengths that are an order of magnitude too large or too small, the operator metrics are accurate. Examples, using data calculated with evolutionary rates and branchings designed to confuse the measurement of branch lengths and to camouflage the topology of the true tree, demonstrate the validity of operator metrics. The method is robust. Operator metric distances are easy to calculate, can be extended to any number of taxa, and provide a statistical estimate of their variances. The utility of the method is demonstrated by using it to analyze the origins and evolution of chloroplasts, mitochondria, and eubacteria.

A structural comparison of type c lysozymes based on their hydropathic profiles

Hydropathic profiles can be considered as an approach to the three-dimensional structure of a protein and so their use for comparison of homologous proteins is proposed, as they provide information on relative structural conservativeness. A simple approach was developed for comparison of hydropathic profiles and applied to 19 lysozymes c of known primary structure. Trees were constructed in order to discover which method yielded the best estimation of the phenotypic differences between the proteins considered, by means of the goodness-of-fit criterion. Iterative methods, such as the Fitch-and-Margoliash and the unweighted-pair-group methods, gave a better fit than did a non-iterative method. When the hydropathic approach is used for comparison of lysozymes c, the enzyme obtained from chachalaca egg-white is placed closer to those from pheasant-like birds than to those of ducks; this result agrees with the morphological resemblance of the chachalaca to pheasant-like birds. Pigeon egg-white and equine milk lysozymes differ greatly in sequence from other lysozymes c and their hydropathic analysis shows important differences with respect to the other homologous enzymes.

Evolutionary relationships of a "primitive" shark (Heterodontus) assessed by micro-complement fixation of serum transferrin

The evolutionary relationships of six sharks were investigated by comparing their transferrins using the micro-complement fixation method. The immunological distances observed were used to build a tree that confirms that the squaloid and galeoid species examined belong to two separate groups and that Heterodontus, a genus of hitherto uncertain position, belongs with the galeoids. The divergence time estimated from the transferrin comparisons is roughly 240 +/- 65 million years between Heterodontus and galeoids.

Structure of the pigeon lysozyme and its relationship with other type c lysozymes

1. The secondary structure of the pigeon egg-white lysozyme shows important differences when compared to other type c lysozymes. These differences are mainly located at the region comprising residues 77-84. This segment contains one alpha-helix in the lysozymes c studied by means of an X-ray analysis, while the residues at such positions in pigeon lysozyme would form two beta-bends. 2. Analysis of the tertiary structure of the pigeon lysozyme by means of hydropathy profiles reveals that the above segment seems to be more hydrophilic in the pigeon enzyme than in other type c lysozymes. 3. Though a certain similarity to the calcium-binding loop of alpha-lactalbumins is detected in pigeon lysozyme, the circular dichroism spectra of the protein at neutral pH do not change in the presence of Ca2+ ions. 4. The presented structural analysis is discussed in terms of function-structure and antigenicity relationships between the type c lysozymes.

Evolution of the genus Leishmania as revealed by comparisons of nuclear DNA restriction fragment patterns

Restriction endonuclease DNA fragment patterns have been used to examine the relationships among 28 isolates of Leishmania as well as Crithidia, Endotrypanum, and Trypanosoma cruzi. Fragments of nuclear DNA were generated with six restriction enzymes, and blots were hybridized with probes from three loci. Among the major lineages the fragment patterns are essentially completely different, while within the major lineages various degrees of divergence are found. Molecular evolutionary trees were constructed using the method of Nei and Li to estimate the percent nucleotide sequence divergence among strains from the fraction of fragments shared. Defined groups, such as species or subspecies within the major lineages, are also grouped by nuclear DNA comparisons. Within the donovani complex, we find Leishmania donovani chagasi and Leishmania donovani infantum to be as similar as strains within Leishmania donovani donovani, consistent with the proposal by other workers that New World visceral leishmaniasis originated quite recently.

Phylogenetic trees constructed from hydrophobicity values of protein sequences

Information about conformational properties of a protein is contained in the hydrophobicity values of the amino acids in its primary sequence. We have investigated the possibility of extracting meaningful evolutionary information from the comparison of the hydrophobicity values of the corresponding amino acids in the sequences of homologous proteins. Distance matrices for six families of homologous proteins were made on the basis of the differences in hydrophobicity values of the amino acids. The phylogenetic trees constructed from such matrices were at least as good (as judged from their faithful reflection of evolutionary relationships), as trees constructed from the usual minimum mutation distance matrix.