Nucleotide sequence of the 3terminal region of belladonna mottle virus (renamed Physalis mottle virus) RNA and analysis of the evolutionary relationships among tymoviral coat proteins

The 3' terminal 1255 nt sequence of Physalis mottle virus (PhMV) genomic RNA has been determined from a set of overlapping cDNA clones. The open reading frame (ORF) at the 3' terminus corresponds to the amino acid sequence of the coat protein (CP) determined earlier except for the absence of the dipeptide, Lys-Leu, at position 110-111. In addiition, the sequence upstream of the CP gene contains the message coding for 178 amino acid residues of the C-terminus of the putative replicase protein (RP). The sequence downstream of the CP gene contains an untranslated region whose terminal 80 nucleotides can be folded into a characteristic tRNA-like structure. A phylogenetic tree constructed after aligning separately the sequence of the CP, the replicase protein (RP) and the tRNA-like structure determined in this study with the corresponding sequences of other tymoviruses shows that PhMV wrongly named belladonna mottle virus [BDMV(I)] is a separate tymovirus and not another strain of BDMV(E) as originally envisaged. The phylogenetic tree in all the three cases is identical showing that any subset of genomic sequence of sufficient length can be used for establishing evolutionary relationships among tymoviruses.

The primary structure of sheep liver cytosolic serine hydroxymethyltransferase and an analysis of the evolutionary relationships among serine hydroxymethyltransferases

The complete amino-acid sequence of sheep liver cytosolic serine hydroxymethyltransferase was determined from an analysis of tryptic, chymotryptic, CNBr and hydroxylamine peptides. Each subunit of sheep liver serine hydroxymethyltransferase consisted of 483 amino-acid residues. A comparison of this sequence with 8 other serine hydroxymethyltransferases revealed that a possible gene duplication event could have occurred after the divergence of animals and fungi. This analysis also showed independent duplication of SHMT genes in Neurospora crassa. At the secondary structural level, all the serine hydroxymethyltransferases belong to the alpha/beta category of proteins. The predicted secondary structure of sheep liver serine hydroxymethyltransferase was similar to that of the observed structure of tryptophan synthase, another pyridoxal 5'-phosphate containing enzyme, suggesting that sheep liver serine hydroxymethyltransferase might have a similar pyridoxal 5'-phosphate binding domain. In addition, a conserved glycine rich region, G L Q G G P, was identified in all the serine hydroxymethyltransferases and could be important in pyridoxal 5'-phosphate binding. A comparison of the cytosolic serine hydroxymethyltransferases from rabbit and sheep liver with other proteins sequenced from both these sources showed that serine hydroxymethyltransferase was a highly conserved protein. It was slightly less conserved than cytochrome c but better conserved than myoglobin, both of which are well known evolutionary markers. C67 and C203 were specifically protected by pyridoxal 5'-phosphate against modification with [C-14]iodoacetic acid, while C247 and C261 were buried in the native serine hydroxymethyltransferase. However, the cysteines are not conserved among the various serine hydroxymethyltransferases. The exact role of the cysteines in the reaction catalyzed by serine hydroxymethyltransferase remains to be elucidated.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Nucleotide sequence of the 3' terminal region of belladonna mottle virus-Iowa (renamed Physalis mottle virus) RNA and an analysis of the relationships of tymoviral coat proteins

The 3prime terminal 1255nt sequence of Physalis mottle virus (PhMV) genomic RNA has been determined from a set of overlapping cDNA clones. The open reading frame (ORF) at the 3prime terminus corresponds to the amino acid sequence of the coat protein (CP) determined earlier except for the absence of the dipeptide, Lys-Leu, at position 110-111. In addition, the sequence upstream of the CP gene contains the message coding for 178 amino acid residues of the C-terminus of the putative replicase protein (RP). The sequence downstream of the CP gene contains an untranslated region whose terminal 80 nucleotides can be folded into a characteristic tRNA-like structure. A phylogenetic tree constructed after aligning separately the sequence of the CP, the replicase protein (RP) and the tRNA-like structure determined in this study with the corresponding sequences of other tymoviruses shows that PhMV wrongly named belladonna mottle virus [BDMV(I)] is a separate tymovirus and not another strain of BDMV(E) as originally envisaged. The phylogenetic tree in all the three cases is identical showing that any subset of genomic sequence of sufficient length can be used for establishing evolutionary relationships among tymoviruses.

Recognizing drug targets using evolutionary information: implications for repurposing FDA-approved drugs against Mycobacterium tuberculosis H37Rv

Drug repurposing to explore target space has been gaining pace over the past decade with the upsurge in the use of systematic approaches for computational drug discovery. Such a cost and time-saving approach gains immense importance for pathogens of special interest, such as Mycobacterium tuberculosis H37Rv. We report a comprehensive approach to repurpose drugs, based on the exploration of evolutionary relationships inferred from the comparative sequence and structural analyses between targets of FDA-approved drugs and the proteins of M. tuberculosis. This approach has facilitated the identification of several polypharmacological drugs that could potentially target unexploited M. tuberculosis proteins. A total of 130 FDA-approved drugs, originally intended against other diseases, could be repurposed against 78 potential targets in M. tuberculosis. Additionally, we have also made an attempt to augment the chemical space by recognizing compounds structurally similar to FDA-approved drugs. For three of the attractive cases we have investigated the probable binding modes of the drugs in their corresponding M. tuberculosis targets by means of structural modelling. Such prospective targets and small molecules could be prioritized for experimental endeavours, and could significantly influence drug-discovery and drug-development programmes for tuberculosis.

Mitochondrial and nuclear markers suggest Hanuman langur (Primates: Colobinae) polyphyly: Implications for their species status

Recent molecular studies on langurs of the Indian subcontinent suggest that the widely-distributed and morphologically variable Hanuman langurs (Semnopithecus entellus) are polyphyletic with respect to Nilgiri and urple-faced langurs. To further investigate this scenario, we have analyzed additional sequences of mitochondrial cytochrome b as well as nuclear protamine P1 genes from these species. The results confirm Hanuman langur polyphyly in the mitochondrial tree and the nuclear markers suggest that the Hanuman langurs share protamine P1 alleles with Nilgiri and purple-faced langurs. We recommend provisional splitting of the so-called Hanuman langurs into three species such that the taxonomy is consistent with their evolutionary relationships.

Analysis of DNA sequence transformations on grids

Study of the evolution of species or organisms is essential for various biological applications. Evolution is typically studied at the molecular level by analyzing the mutations of DNA sequences of organisms. Techniques have been developed for building phylogenetic or evolutionary trees for a set of sequences. Though phylogenetic trees capture the overall evolutionary relationships among the sequences, they do not reveal fine-level details of the evolution. In this work, we attempt to resolve various fine-level sequence transformation details associated with a phylogenetic tree using cellular automata. In particular, our work tries to determine the cellular automata rules for neighbor-dependent mutations of segments of DNA sequences. We also determine the number of time steps needed for evolution of a progeny from an ancestor and the unknown segments of the intermediate sequences in the phylogenetic tree. Due to the existence of vast number of cellular automata rules, we have developed a grid system that performs parallel guided explorations of the rules on grid resources. We demonstrate our techniques by conducting experiments on a grid comprising machines in three countries and obtaining potentially useful statistics regarding evolutions in three HIV sequences. In particular, our work is able to verify the phenomenon of neighbor-dependent mutations and find that certain combinations of neighbor-dependent mutations, defined by a cellular automata rule, occur with greater than 90% probability. We also find the average number of time steps for mutations for some branches of phylogenetic tree over a large number of possible transformations with standard deviations less than 2.

Identification of Local Conformational Similarity in Structurally Variable Regions of Homologous Proteins Using Protein Blocks

Structure comparison tools can be used to align related protein structures to identify structurally conserved and variable regions and to infer functional and evolutionary relationships. While the conserved regions often superimpose well, the variable regions appear non superimposable. Differences in homologous protein structures are thought to be due to evolutionary plasticity to accommodate diverged sequences during evolution. One of the kinds of differences between 3-D structures of homologous proteins is rigid body displacement. A glaring example is not well superimposed equivalent regions of homologous proteins corresponding to a-helical conformation with different spatial orientations. In a rigid body superimposition, these regions would appear variable although they may contain local similarity. Also, due to high spatial deviation in the variable region, one-to-one correspondence at the residue level cannot be determined accurately. Another kind of difference is conformational variability and the most common example is topologically equivalent loops of two homologues but with different conformations. In the current study, we present a refined view of the ``structurally variable'' regions which may contain local similarity obscured in global alignment of homologous protein structures. As structural alphabet is able to describe local structures of proteins precisely through Protein Blocks approach, conformational similarity has been identified in a substantial number of `variable' regions in a large data set of protein structural alignments; optimal residue-residue equivalences could be achieved on the basis of Protein Blocks which led to improved local alignments. Also, through an example, we have demonstrated how the additional information on local backbone structures through protein blocks can aid in comparative modeling of a loop region. In addition, understanding on sequence-structure relationships can be enhanced through our approach. This has been illustrated through examples where the equivalent regions in homologous protein structures share sequence similarity to varied extent but do not preserve local structure.

AlignHUSH: Alignment of HMMs using structure and hydrophobicity information

Background: Sensitive remote homology detection and accurate alignments especially in the midnight zone of sequence similarity are needed for better function annotation and structural modeling of proteins. An algorithm, AlignHUSH for HMM-HMM alignment has been developed which is capable of recognizing distantly related domain families The method uses structural information, in the form of predicted secondary structure probabilities, and hydrophobicity of amino acids to align HMMs of two sets of aligned sequences. The effect of using adjoining column(s) information has also been investigated and is found to increase the sensitivity of HMM-HMM alignments and remote homology detection. Results: We have assessed the performance of AlignHUSH using known evolutionary relationships available in SCOP. AlignHUSH performs better than the best HMM-HMM alignment methods and is observed to be even more sensitive at higher error rates. Accuracy of the alignments obtained using AlignHUSH has been assessed using the structure-based alignments available in BaliBASE. The alignment length and the alignment quality are found to be appropriate for homology modeling and function annotation. The alignment accuracy is found to be comparable to existing methods for profile-profile alignments. Conclusions: A new method to align HMMs has been developed and is shown to have better sensitivity at error rates of 10% and above when compared to other available programs. The proposed method could effectively aid obtaining clues to functions of proteins of yet unknown function. A web-server incorporating the AlignHUSH method is available at http://crick.mbu.iisc.ernet.in/similar to alignhush/

On the systematic position of some Asian enigmatic genera of Asclepiadoideae (Apocynaceae)

The phylogenetic structure of Asclepiadoideae (Apocynaceae) has been elucidated at the tribal and subtribal levels in the last two decades. However, to date, the systematic positions of seven Asian genera, Cosmostigma, Graphistemma, Holostemma, Pentasachme, Raphistemma, Seshagiria and Treutlera, have not been investigated. In this study, we examine the evolutionary relationships among these seven small enigmatic Asian genera and clarify their positions in Asclepiadoideae, using a combination of plastid sequences of rbcL, rps16, trnL and trnL- F regions. Cosmostigma and Treutlera are resolved as members of the non-Hoya clade of Marsdenieae with strong support (maximum parsimony bootstrap support value BSMP = 96, maximum likelihood bootstrap support value BSML = 98, Bayesian-inferred posterior probability PP = 1.0). Pentasachme is resolved as sister of Stapeliinae to Ceropegieae with moderate support (BSMP = 64, BSML = 66, PP = 0.94). Graphistemma, Holostemma, Raphistemma and Seshagiria are all nested in the Asclepiadeae-Cynanchinae clade (BSMP = 97, BSML = 100, PP = 1.0). The study confirms the generally accepted tribal and subtribal structure of the subfamily. One exception is Eustegia minuta, which is placed here as sister to all Asclepiadeae (BSMP = 58, BSML = 76, PP = 0.99) and not as sister to the Marsdenieae + Ceropegieae clade. The weak support and conflicting position indicate the need for a placement of Eustegia as an independent tribe. In Asclepiadeae, a sister group position of Cynanchinae to the Asclepiadinae + Tylophorinae clade is favoured (BSMP = 84, BSML = 88, PP = 1.0), whereas Schizostephanus is retrieved as unresolved. Oxystelma appears as an early-branching member of Asclepiadinae with weak support (BSMP = 52, BSML = 74, PP = 0.69). Calciphila and Solenostemma are also associated with Asclepiadinae with weak support (BSMP = 37, BSML = 45, PP = 0.79), but all alternative positions are essentially without support. The position of Indian Asclepiadoideae in the family phylogeny is discussed. (c) 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 174, 601-619.

Differences in host species relationships and biogeographic influences produce contrasting patterns of prevalence, community composition and genetic structure in two genera of avian malaria parasites in southern Melanesia

1. Host-parasite interactions have the potential to influence broadscale ecological and evolutionary processes, levels of endemism, divergence patterns and distributions in host populations. Understanding the mechanisms involved requires identification of the factors that shape parasite distribution and prevalence. 2. A lack of comparative information on community-level host-parasite associations limits our understanding of the role of parasites in host population divergence processes. Avian malaria (haemosporidian) parasites in bird communities offer a tractable model system to examine the potential for pathogens to influence evolutionary processes in natural host populations. 3. Using cytochrome b variation, we characterized phylogenetic diversity and prevalence of two genera of avian haemosporidian parasites, Plasmodium and Haemoproteus, and analysed biogeographic patterns of lineages across islands and avian hosts, in southern Melanesian bird communities to identify factors that explain patterns of infection. 4. Plasmodium spp. displayed isolation-by-distance effects, a significant amount of genetic variation distributed among islands but insignificant amounts among host species and families, and strong local island effects with respect to prevalence. Haemoproteus spp. did not display isolation-by-distance patterns, showed marked structuring of genetic variation among avian host species and families, and significant host species prevalence patterns. 5. These differences suggest that Plasmodium spp. infection patterns were shaped by geography and the abiotic environment, whereas Haemoproteus spp. infection patterns were shaped predominantly by host associations. Heterogeneity in the complement and prevalence of parasite lineages infecting local bird communities likely exposes host species to a mosaic of spatially divergent disease selection pressures across their naturally fragmented distributions in southern Melanesia. Host associations for Haemoproteus spp. indicate a capacity for the formation of locally co-adapted host-parasite relationships, a feature that may limit intraspecific gene flow or range expansions of closely related host species.

Stress-state permeability relationships for clays with coarse fraction

In standard laboratory consolidation tests, only the fraction of soil passing through a particular size of the sieve, called the matrix material, is used. This size is usually restricted to 1/10 of the height of the consolidation ring. Particles larger than this size that are removed before the test may consist of gravel, fragments of rock, or other coarse materials. Hence, it is not possible to estimate the compressibility and permeability of the total material based on the compressibility and permeability behavior obtained from laboratory consolidation tests on the matrix material. In the present investigation an attempt has been made to estimate the compressibility and permeability behavior of the total material based on the compressibility and permeability behavior of the matrix material. The results indicate that the presence of coarse particles will reduce the compressibility of the soil in proportion to the coarse fraction present in the whole soil and will not affect the permeability of the soil for the range investigated. If the coarse fraction exceeds the Limiting percentage, the void ratio-vertical effective stress path will also start to deviate from the predicted path. An expression has been developed to estimate approximately the deviating pressure, and it is found to depend on the soil type as well as the percent clay fraction.

Evolutionary games with two timescales

We consider a two timescale model of learning by economic agents wherein active or 'ontogenetic' learning by individuals takes place on a fast scale and passive or 'phylogenetic' learning by society as a whole on a slow scale, each affecting the evolution of the other. The former is modelled by the Monte Carlo dynamics of physics, while the latter is modelled by the replicator dynamics of evolutionary biology. Various qualitative aspects of the dynamics are studied in some simple cases, both analytically and numerically, and its role as a useful modelling device is emphasized.

Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots

1 Species-accumulation curves for woody plants were calculated in three tropical forests, based on fully mapped 50-ha plots in wet, old-growth forest in Peninsular Malaysia, in moist, old-growth forest in central Panama, and in dry, previously logged forest in southern India. A total of 610 000 stems were identified to species and mapped to < Im accuracy. Mean species number and stem number were calculated in quadrats as small as 5 m x 5 m to as large as 1000 m x 500 m, for a variety of stem sizes above 10 mm in diameter. Species-area curves were generated by plotting species number as a function of quadrat size; species-individual curves were generated from the same data, but using stem number as the independent variable rather than area. 2 Species-area curves had different forms for stems of different diameters, but species-individual curves were nearly independent of diameter class. With < 10(4) stems, species-individual curves were concave downward on log-log plots, with curves from different forests diverging, but beyond about 104 stems, the log-log curves became nearly linear, with all three sites having a similar slope. This indicates an asymptotic difference in richness between forests: the Malaysian site had 2.7 times as many species as Panama, which in turn was 3.3 times as rich as India. 3 Other details of the species-accumulation relationship were remarkably similar between the three sites. Rectangular quadrats had 5-27% more species than square quadrats of the same area, with longer and narrower quadrats increasingly diverse. Random samples of stems drawn from the entire 50 ha had 10-30% more species than square quadrats with the same number of stems. At both Pasoh and BCI, but not Mudumalai. species richness was slightly higher among intermediate-sized stems (50-100mm in diameter) than in either smaller or larger sizes, These patterns reflect aggregated distributions of individual species, plus weak density-dependent forces that tend to smooth the species abundance distribution and 'loosen' aggregations as stems grow. 4 The results provide support for the view that within each tree community, many species have their abundance and distribution guided more by random drift than deterministic interactions. The drift model predicts that the species-accumulation curve will have a declining slope on a log-log plot, reaching a slope of O.1 in about 50 ha. No other model of community structure can make such a precise prediction. 5 The results demonstrate that diversity studies based on different stem diameters can be compared by sampling identical numbers of stems. Moreover, they indicate that stem counts < 1000 in tropical forests will underestimate the percentage difference in species richness between two diverse sites. Fortunately, standard diversity indices (Fisher's sc, Shannon-Wiener) captured diversity differences in small stem samples more effectively than raw species richness, but both were sample size dependent. Two nonparametric richness estimators (Chao. jackknife) performed poorly, greatly underestimating true species richness.

A fuzzy multistage evolutionary (FUME) clustering technique

n this paper, a multistage evolutionary scheme is proposed for clustering in a large data base, like speech data. This is achieved by clustering a small subset of the entire sample set in each stage and treating the cluster centroids so obtained as samples, together with another subset of samples not considered previously, as input data to the next stage. This is continued till the whole sample set is exhausted. The clustering is accomplished by constructing a fuzzy similarity matrix and using the fuzzy techniques proposed here. The technique is illustrated by an efficient scheme for voiced-unvoiced-silence classification of speech.