Phylogenetic Predictions on Grids

A phylogenetic or evolutionary tree is constructed from a set of species or DNA sequences and depicts the relatedness between the sequences. Predictions of future sequences in a phylogenetic tree are important for a variety of applications including drug discovery, pharmaceutical research and disease control. In this work, we predict future DNA sequences in a phylogenetic tree using cellular automata. Cellular automata are used for modeling neighbor-dependent mutations from an ancestor to a progeny in a branch of the phylogenetic tree. Since the number of possible ways of transformations from an ancestor to a progeny is huge, we use computational grids and middleware techniques to explore the large number of cellular automata rules used for the mutations. We use the popular and recurring neighbor-based transitions or mutations to predict the progeny sequences in the phylogenetic tree. We performed predictions for three types of sequences, namely, triose phosphate isomerase, pyruvate kinase, and polyketide synthase sequences, by obtaining cellular automata rules on a grid consisting of 29 machines in 4 clusters located in 4 countries, and compared the predictions of the sequences using our method with predictions by random methods. We found that in all cases, our method gave about 40% better predictions than the random methods.

Testing concordance in species boundaries using acoustic, morphological, and molecular data in the field cricket genus Itaropsis (Orthoptera: Grylloidea, Gryllidae: Gryllinae)

In most taxa, species boundaries are inferred based on differences in morphology or DNA sequences revealed by taxonomic or phylogenetic analyses. In crickets, acoustic mating signals or calling songs have species-specific structures and provide a third data set to infer species boundaries. We examined the concordance in species boundaries obtained using acoustic, morphological, and molecular data sets in the field cricket genus Itaropsis. This genus is currently described by only one valid species, Itaropsis tenella, with a broad distribution in western peninsular India and Sri Lanka. Calling songs of males sampled from four sites in peninsular India exhibited significant differences in a number of call features, suggesting the existence of multiple species. Cluster analysis of the acoustic data, molecular phylogenetic analyses, and phylogenetic analyses combining all data sets suggested the existence of three clades. Whatever the differences in calling signals, no full congruence was obtained between all the data sets, even though the resultant lineages were largely concordant with the acoustic clusters. The genus Itaropsis could thus be represented by three morphologically cryptic incipient species in peninsular India; their distributions are congruent with usual patterns of endemism in the Western Ghats, India. Song evolution is analysed through the divergence in syllable period, syllable and call duration, and dominant frequency.

Mechanistic insights into type III restriction enzymes

Type III restriction-modification (R-M) enzymes need to interact with two separate unmethylated DNA sequences in indirectly repeated, head-to-head orientations for efficient cleavage to occur at a defined location next to only one of the two sites. However, cleavage of sites that are not in head-to-head orientation have been observed to occur under certain reaction conditions in vitro. ATP hydrolysis is required for the long-distance communication between the sites prior to cleavage. Type III R-M enzymes comprise two subunits, Res and Mod that form a homodimeric Mod(2) and a heterotetrameric Res(2)Mod(2) complex. The Mod subunit in M-2 or R2M2 complex recognizes and methylates DNA while the Res subunit in R2M2 complex is responsible for ATP hydrolysis, DNA translocation and cleavage. A vast majority of biochemical studies on Type III R-M enzymes have been undertaken using two closely related enzymes, EcoP1I and EcoP15I. Divergent opinions about how the long-distance interaction between the recognition sites exist and at least three mechanistic models based on 1D- diffusion and/or 3D-DNA looping have been proposed.

A Comparative Kinetic and Thermodynamic Perspective of the sigma-Competition Model in Escherichia coli

Transcription is the most fundamental step in gene expression in any living organism. Various environmental cues help in the maturation of core RNA polymerase (RNAP; alpha(2)beta beta'omega) with different sigma-factors, leading to the directed recruitment of RNAP to different promoter DNA sequences. Thus it is essential to determine the sigma-factors that affect the preferential partitioning of core RNAP among various a-actors, and the role of sigma-switching in transcriptional gene regulation. Further, the macromolecular assembly of holo RNAP takes place in an extremely crowded environment within a cell, and thus far the kinetics and thermodynamics of this molecular recognition process have not been well addressed. In this study we used a site-directed bioaffinity immobilization method to evaluate the relative binding affinities of three different Escherichia coli sigma-factors to the same core RNAP with variations in temperature and ionic strength while emulating the crowded cellular milieu. Our data indicate that the interaction of core RNAP-sigma is susceptible to changes in external stimuli such as osmolytic and thermal stress, and the degree of susceptibility varies among different sigma-factors. This allows for a reversible sigma-switching from housekeeping factors to alternate sigma-factors when the organism senses a change in its physiological conditions.

Sequence and structural basis for chromosomal fragility during translocations in cancer

Chromosomal aberration is considered to be one of the major characteristic features in many cancers. Chromosomal translocation, one type of genomic abnormality, can lead to deregulation of critical genes involved in regulating important physiological functions such as cell proliferation and DNA repair. Although chromosomal translocations were thought to be random events, recent findings suggest that certain regions in the human genome are more susceptible to breakage than others. The possibility of deviation from the usual B-DNA conformation in such fragile regions has been an active area of investigation. This review summarizes the factors that contribute towards the fragility of these regions in the chromosomes, such as DNA sequences and the role of different forms of DNA structures. Proteins responsible for chromosomal fragility, and their mechanism of action are also discussed. The effect of positioning of chromosomes within the nucleus favoring chromosomal translocations and the role of repair mechanisms are also addressed.

Type III restriction-modification enzymes: a historical perspective

Restriction endonucleases interact with DNA at specific sites leading to cleavage of DNA. Bacterial DNA is protected from restriction endonuclease cleavage by modifying the DNA using a DNA methyltransferase. Based on their molecular structure, sequence recognition, cleavage position and cofactor requirements, restriction-modification (R-M) systems are classified into four groups. Type III R-M enzymes need to interact with two separate unmethylated DNA sequences in inversely repeated head-to-head orientations for efficient cleavage to occur at a defined location (25-27 bp downstream of one of the recognition sites). Like the Type I R-M enzymes, Type III R-M enzymes possess a sequence-specific ATPase activity for DNA cleavage. ATP hydrolysis is required for the long-distance communication between the sites before cleavage. Different models, based on 1D diffusion and/or 3D-DNA looping, exist to explain how the long-distance interaction between the two recognition sites takes place. Type III R-M systems are found in most sequenced bacteria. Genome sequencing of many pathogenic bacteria also shows the presence of a number of phase-variable Type III R-M systems, which play a role in virulence. A growing number of these enzymes are being subjected to biochemical and genetic studies, which, when combined with ongoing structural analyses, promise to provide details for mechanisms of DNA recognition and catalysis.

Mud-packing frog: A novel breeding behaviour and parental care in a stream dwelling new species of Nyctibatrachus (Amphibia, Anura, Nyctibatrachidae)

Reproductive modes are diverse and unique in anurans. Selective pressures of evolution, ecology and environment are attributed to such diverse reproductive modes. Globally forty different reproductive modes in anurans have been described to date. The genus Nyctibatrachus has been recently revised and belongs to an ancient lineage of frog families in the Western Ghats of India. Species of this genus are known to exhibit mountain associated clade endemism and novel breeding behaviours. The purpose of this study is to present unique reproductive behaviour, oviposition and parental care in a new species Nyctibatrachus kumbara sp. nov. which is described in the paper. Nyctibatrachus kumbara sp. nov. is a medium sized stream dwelling frog. It is distinct from the congeners based on a suite of morphological characters and substantially divergent in DNA sequences of the mitochondrial 16S rRNA gene. Males exhibit parental care by mud packing the egg clutch. Such parental care has so far not been described from any other frog species worldwide. Besides this, we emphasize that three co-occurring congeneric species of Nyctibatrachus, namely N. jog, N. kempholeyensis and Nyctibatrachus kumbara sp. nov. from the study site differ in breeding behaviour, which could represent a case of reproductive character displacement. These three species are distinct in their size, call pattern, reproductive behaviour, maximum number of eggs in a clutch, oviposition and parental care, which was evident from the statistical analysis. The study throws light on the reproductive behaviour of Nyctibatrachus kumbara sp. nov. and associated species to understand the evolution and adaptation of reproductive modes of anurans in general, and Nyctibatrachus in particular from the Western Ghats.

Analysis of the Genome and Transcriptome of Cryptococcus neoformans var. grubii Reveals Complex RNA Expression and Microevolution Leading to Virulence Attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Distribution of simple repetitive (TG/CA)n and (CT/AG)n sequences in human and rodent genomes

Sixteen million nucleotide sequence of genome of various organisms have been analysed to detect and study the extent of occurrence of simple repetitive sequences. Two sequence motifs (TG/CA)n and (CT/AG)n capable of adopting unusual DNA structures, left handed Z-conformation and triple-helical conformation respectively, are found to be abundant in rodent and human genomes, but almost completely absent in bacterial genome. (TG/CA)n and (CT/AG)n sequences are present mostly in the intron or 5'/3' flanking regions of the genes. The presence of such repeat motifs in genomic sequence of higher eukaryotes has been correlated with their possible functional significance in nucleosome organization, recombination and gene expression.

Solution conformations of penta and heptapeptides containing repetitive α-aminoisobutyryl-Image -alanyl and α -aminoisobutyryl-Image -valyl sequences

The presence of folded solution conformations in the peptides Boc-Ala-(Aib-Ala)2-OMe, Boc-Val-(Aib-Val) 2-OMe, Boc-Ala-(Aib-Ala)3-OMe and Boc-Val-(Aib-Val)3-OMe has been established by 270MHz 1H NMR. Intramolecularly H-bonded NH groups have been identified using temperature and solvent dependence of NH chemical shifts and paramagnetic radical induced broadening of NH resonances. Both pentapeptides adopt 310 helical conformations possessing 3 intramolecular H-bonds in CDCl3 and (CD3)2SO. The heptapeptides favour helical structures with 5 H-bonds in CDCl3. In (CD3)2SO only 4 H-bonds are readily detected.

Sequences that adopt non-B-DNA conformation in form V DNA as probed by enzymic methylation

pBR322 form V DNA is a highly torsionally strained molecule with a linking number of zero. We have used sequence- specific DNA methylases as probes for B-DNA in this molecule, exploiting the inability of methylases to methylate single-stranded DNA and Z-DNA, both of which are known to occur in form V DNA. Some sequences in form V DNA were shown to be totally in the B-form, others were totally in an altered, unmethylatable conformation, while still other sites appeared to exist partly in altered and partly in normal B-conformation. Some potential Z-forming sequences (alternating pyrimidine/purine) of less than seven base-pairs were not in the Z conformation in form V DNA, whereas others did adopt an altered structure, indicating a modulating influence of flanking sequences. Furthermore, regions of imperfect alternating pyrimidine/purine structure were sometimes capable of adopting an altered structure. In addition, some regions of altered structure had no apparent Z-forming sequences, nor were they in polypurine stretches, which have also been proposed to form left-handed DNA. These non-B-DNA conformations may represent novel left-handed helical structures or sequences that become single stranded under torsional strain. Long regions of either altered (unmethylatable) DNA or B-DNA were not always observed. In fact, one region showed three transitions between B-like DNA and altered structure within 26 base-pairs.

Polypurine-polypyrimidine sequences adopt unwound structure in pBR322 form V DNA as probed by single-hit analysis of HpaII sites.

Structure at the polypurine-polypyrimidine sequences flanking the HpaII sites (CCGG) in pBR322 form V DNA was probed employing single-hit analysis using HpaII restriction endonuclease. Reduced cleavage efficiency of HpaII sites flanked by polypurine-polypyrimidine sequences suggested that under high torsional stress these sequences adopt unwound structures rendering these sites insensitive to restriction enzyme cleavage. In addition to polypurine-polypyrimidine sequences. HpaII sites flanked by alternating purine-pyrimidine sequence, a potential motif of left handed Z-DNA, were also found to be resistant to HpaII cleavage. Results obtained from various studies implicating structure sensitivity of restriction endonucleases and methylases were compiled and a direct correlation was observed between the occurrence of altered sites in a domain and its G/C content in pBR322 form V DNA.

Unwinding of heterologous DNA by RecA protein during the search for homologous sequences

The search for homologous sequences promoted by RecA protein in vitro involves a presynaptic filament and naked duplex DNA, the multiple contacts of which produce nucleoprotein networks or coaggregates. The single-stranded DNA within the presynaptic filaments, however, is extended to an axial spacing 1.5 times that of B-form DNA. To investigate this paradoxical difference between the spacing of bases in the RecA presynaptic filament versus the target duplex DNA, we explored the effect of heterologous contacts on the conformation of DNA, and vice versa. In the presence of wheat germ topoisomerase I, RecA presynaptic filaments induced a rapid, limited reduction in the linking number of heterologous circular duplex DNA. This limited unwinding of heterologous duplex DNA, termed heterologous unwinding, was detected within 30 seconds and reached a steady state within a few minutes. Presynaptic filaments that were formed in the presence of ATP?S and separated from free RecA protein by gel filtration also generated a ladder of topoisomers upon incubation with relaxed duplex DNA and topoisomerase. The inhibition of heterologous contacts by 60 mImage -NaCl or 5 mImage -ADP resulted in a corresponding decrease in heterologous unwinding. In reciprocal fashion, the stability or number of heterologous contacts with presynaptic filaments was inversely related to the linking number of circular duplex DNA. These observations show that heterologous contacts with the presynaptic filament cause a limited unwinding of the duplex DNA, and conversely that the ability of the DNA to unwind stabilizes transient heterologous contacts.

Left-Handed DNA in Synthetic and Topologically Constrained Form V-DNA

We have investigated structural transitions in Poly(dG-dC) and Poly(dG-Me5dC) in order to understand the exact role of cations in stabilizing left-handed helical structures in specific sequences andthe biological role, if any, of these structures. From a novel temperature dependent transition it has been shown that a minor fluctuation in Na+ concentration at ambient temperature can bring about Β to Ζ transition. Forthe first time, wehave observed a novel double transition in poly(dG-Me5dC) as the Na+ concentration is gradually increased. This suggests that a minor fluctuation in Na+ concentration in conjunction with methylation may transform small stretches of CG sequences from one conformational state to another. These stretches could probably serve as sites for regulation. Supercoiled formV DNA reconstituted from pBR322 and pßG plasmids have been studied as model systems, in order to understand the nature and role of left-handed helical conformation in natural sequences. A large portion of DNA in form V, obtained by reannealing the two complementary singlestranded circles is forced to adopt left-handed double helical structure due to topological constraints (Lk = 0). Binding studies with Z-DNA specific antibody and spectroscopic studies confirm the presence of left-handed Z-structure in the pßG and pßR322 form V DNA. Cobalt hexamine chloride, which induces Z-form in Poly(dG-dC) stabilizes the Z-conformation in form V DNA even in the non-alternating purine-pyrimidine sequences. A reverse effect is observed with ethidium bromide. Interestingly, both topoisomerase I and II (from wheat germ) act effectively on form V DNA to give rise to a species having an electrophoretic mobility on agarose gel similar to that of open circular (form II) DNA. Whether this molecule is formed as a result of the left-handed helical segments of form V DNA undergoing a transition to the right-handed B-form during the topoisomerase action remains to be solved.

Effects Of 5' Flanking Sequences And Changes In The 5' Internal Control Region On The Transcription Of Rice Transfer-Rna Gcc Cly Gene

A stretch of 71 nucleotides in a 1.2 kilobase pair Pst I fragment of rice DNA was identified as tRNA~ gene by hybridization and nucleotide sequence analyses. The hybridization of genomic DNA with the tRNA gene showed that there are about 10 glycine tRNA genes per diploid rice genome. The 3' and 5' internal control regions, where RNA polymerase III and transcription factors bind, were found to be present in the coding sequence. The gene was transcribed into a 4S product in an yeast cell-free extract. The substitution of 5' internal control region with analogous sequences from either M13mpl9 or M13mpl8 DNA did not affect the transcription of the gene in vitro. The changes in three highly conserved nucleotides in the consensus 5' internal control region (RGYNNARYGG; R = purine, Y = pyrimidine, N = any nucleotide) did not affect transcription showing that these nucleotides are not essential for promotion of transcription. There were two 16 base pair repeats, 'TGTTTGTTTCAGCTTA' at - 130 and - 375 positions upstream from the start of the gene. Deletion of 5' flanking sequences including the 16 base pair repeat at - 375 showed increased transcription indicating that these sequences negatively modulate the expression of the gene.