Mean Squared Residue Based Biclustering Algorithms for the Analysis of Gene Expression Data

Computational Biology is the research are that contributes to the analysis of biological data through the development of algorithms which will address significant research problems.The data from molecular biology includes DNA,RNA ,Protein and Gene expression data.Gene Expression Data provides the expression level of genes under different conditions.Gene expression is the process of transcribing the DNA sequence of a gene into mRNA sequences which in turn are later translated into proteins.The number of copies of mRNA produced is called the expression level of a gene.Gene expression data is organized in the form of a matrix. Rows in the matrix represent genes and columns in the matrix represent experimental conditions.Experimental conditions can be different tissue types or time points.Entries in the gene expression matrix are real values.Through the analysis of gene expression data it is possible to determine the behavioral patterns of genes such as similarity of their behavior,nature of their interaction,their respective contribution to the same pathways and so on. Similar expression patterns are exhibited by the genes participating in the same biological process.These patterns have immense relevance and application in bioinformatics and clinical research.Theses patterns are used in the medical domain for aid in more accurate diagnosis,prognosis,treatment planning.drug discovery and protein network analysis.To identify various patterns from gene expression data,data mining techniques are essential.Clustering is an important data mining technique for the analysis of gene expression data.To overcome the problems associated with clustering,biclustering is introduced.Biclustering refers to simultaneous clustering of both rows and columns of a data matrix. Clustering is a global whereas biclustering is a local model.Discovering local expression patterns is essential for identfying many genetic pathways that are not apparent otherwise.It is therefore necessary to move beyond the clustering paradigm towards developing approaches which are capable of discovering local patterns in gene expression data.A biclusters is a submatrix of the gene expression data matrix.The rows and columns in the submatrix need not be contiguous as in the gene expression data matrix.Biclusters are not disjoint.Computation of biclusters is costly because one will have to consider all the combinations of columans and rows in order to find out all the biclusters.The search space for the biclustering problem is 2 m+n where m and n are the number of genes and conditions respectively.Usually m+n is more than 3000.The biclustering problem is NP-hard.Biclustering is a powerful analytical tool for the biologist.The research reported in this thesis addresses the problem of biclustering.Ten algorithms are developed for the identification of coherent biclusters from gene expression data.All these algorithms are making use of a measure called mean squared residue to search for biclusters.The objective here is to identify the biclusters of maximum size with the mean squared residue lower than a given threshold. All these algorithms begin the search from tightly coregulated submatrices called the seeds.These seeds are generated by K-Means clustering algorithm.The algorithms developed can be classified as constraint based,greedy and metaheuristic.Constarint based algorithms uses one or more of the various constaints namely the MSR threshold and the MSR difference threshold.The greedy approach makes a locally optimal choice at each stage with the objective of finding the global optimum.In metaheuristic approaches particle Swarm Optimization(PSO) and variants of Greedy Randomized Adaptive Search Procedure(GRASP) are used for the identification of biclusters.These algorithms are implemented on the Yeast and Lymphoma datasets.Biologically relevant and statistically significant biclusters are identified by all these algorithms which are validated by Gene Ontology database.All these algorithms are compared with some other biclustering algorithms.Algorithms developed in this work overcome some of the problems associated with the already existing algorithms.With the help of some of the algorithms which are developed in this work biclusters with very high row variance,which is higher than the row variance of any other algorithm using mean squared residue, are identified from both Yeast and Lymphoma data sets.Such biclusters which make significant change in the expression level are highly relevant biologically.

Personnel Management Practices in the Kerala based Scheduled commercial Banks

The thesis entitled Personnel Management Practices in the Kerala-Based Scheduled Commercial Banks. Personnel management function is of cardinal importance, requiring a sophisticated and scientific approach. In a labour-intensive, service industry like banking. Productivity and ultimate profitability of the entire organization depend considerably on the effectiveness with which personnel management function is executed; and the prudence with which personnel problems are handle. The main objectives of the study are to understand the current status of personnel management functions in the banks and to evaluate the practices in the light of the principles and theories of personnel management so as to identify the strengths and weaknesses. The universe of this study is the eight Scheduled Commercial Banks based in Kerala. The major limitation of the study is that as State Bank of Travancore, the lone public sector bank based in Kerala did not grant permission for collection of data, this study had to be confined to private sector banks only. Almost the entire data used for this study are primary and were collected from the files and other records or the concerned banks. This report has chapters dealing with the functional areas of personnel management such as determination of human resource requirements, recruitment and selection, training and development, performance appraisal, promotions and compensation. Findings reveal that the practice of personnel management in the Kerala-based private sector scheduled commercial banks has not gained a degree of sophistication compatible with its role in modern business management.

A Compiler Integrated Assistance for Optimum Data Allocation in Banked Memory Embedded Processors

Bank switching in embedded processors having partitioned memory architecture results in code size as well as run time overhead. An algorithm and its application to assist the compiler in eliminating the redundant bank switching codes introduced and deciding the optimum data allocation to banked memory is presented in this work. A relation matrix formed for the memory bank state transition corresponding to each bank selection instruction is used for the detection of redundant codes. Data allocation to memory is done by considering all possible permutation of memory banks and combination of data. The compiler output corresponding to each data mapping scheme is subjected to a static machine code analysis which identifies the one with minimum number of bank switching codes. Even though the method is compiler independent, the algorithm utilizes certain architectural features of the target processor. A prototype based on PIC 16F87X microcontrollers is described. This method scales well into larger number of memory blocks and other architectures so that high performance compilers can integrate this technique for efficient code generation. The technique is illustrated with an example

DNA methylation in Dictyostelium discoideum

DNA methyltransferases of type Dnmt2 are a highly conserved protein family with enigmatic function. The aim of this work was to characterize DnmA, the Dnmt2 methyltransferase in Dictyostelium discoideum, and further to investigate its implication in DNA methylation and transcriptional gene silencing. The genome of the social amoeba Dictyostelium encodes DnmA as the sole DNA methyltransferase. The enzyme bears all ten characteristic DNA methyltransferase motifs in its catalytic domain. The DnmA mRNA was found by RT-PCR to be expressed during vegetative growth and down regulated during development. Investigations using fluorescence microscopy showed that both DnmA-myc and DnmA-GFP fusions predominantly localised to the nucleus. The function of DnmA remained initially unclear, but later experiment revealed that the enzyme is an active DNA methyltransferase responsible for all DNA (cytosine) methylation in Dictyostelium. Neither in gel retardation assays, nor by the yeast two hybrid system, clues on the functionality of DnmA could be obtained. However, immunological detection of the methylation mark with an α - 5mC antibody gave initial evidence that the DNA of Dictyostelium was methylated. Furthermore, addition of 5-aza-cytidine as demethylating agent to the Dictyostelium medium and subsequent in vitro incubation of the DNA isolated from these cells with recombinant DnmA showed that the enzyme binds slightly better to this target DNA. In order to investigate further the function of the protein, a gene knock-out for dnmA was generated. The gene was successfully disrupted by homologous recombination, the knock-out strain, however, did not show any obvious phenotype under normal laboratory conditions. To identify specific target sequences for DNA methylation, a microarray analysis was carried out. Setting a threshold of at least 1.5 fold for differences in the strength of gene expression, several such genes in the knock-out strain were chosen for further investigation. Among the up-regulated genes were the ESTs representing the gag and the RT genes respectively of the retrotransposon skipper. In addition Northern blot analysis confirmed the up-regulation of skipper in the DnmA knock-out strain. Bisufite treatment and sequencing of specific DNA stretches from skipper revealed that DnmA is responsible for methylation of mostly asymmetric cytosines. Together with skipper, DIRS-1 retrotransposon was found later also to be methylated but was not present on the microarray. Furthermore, skipper transcription was also up-regulated in strains that had genes disrupted encoding components of the RNA interference pathway. In contrast, DIRS 1 expression was not affected by a loss of DnmA but was strongly increased in the strain that had the RNA directed RNA polymerase gene rrpC disrupted. Strains generated by propagating the usual wild type Ax2 and the DnmA knock-out cells over 16 rounds in development were analyzed for transposon activity. Northern blot analysis revealed activation for skipper expression, but not for DIRS-1. A large number of siRNAs were found to be correspondent to the DIRS-1 sequence, suggesting concerted regulation of DIRS-1 expression by RNAi and DNA methylation. In contrast, no siRNAs corresponding to the standard skipper element were found. The data show that DNA methylation plays a crucial role in epigenetic gene regulation in Dictyostelium and that different, partially overlapping mechanisms control transposon silencing for skipper and DIRS-1. To elucidate the mechanism of targeting the protein to particular genes in the Dictyostelium genome, some more genes which were up-regulated in the DnmA knock-out strain were analyzed by bisulfite sequencing. The chosen genes are involved in the multidrug response in other species, but their function in Dictyostelium is uncertain. Bisulfite data showed that two of these genes were methylated at asymmetrical C-residues in the wild type, but not in DnmA knock-out cells. This suggested that DNA methylation in Dictyostelium is involved not only in transposon regulation but also in transcriptional silencing of specific genes.

Land plants and DNA barcodes: short-term and long-term goals

Land plants have had the reputation of being problematic for DNA barcoding for two general reasons: (i) the standard DNA regions used in algae, animals and fungi have exceedingly low levels of variability and (ii) the typically used land plant plastid phylogenetic markers (e.g. rbcL, trnL-F, etc.) appear to have too little variation. However, no one has assessed how well current phylogenetic resources might work in the context of identification (versus phylogeny reconstruction). In this paper, we make such an assessment, particularly with two of the markers commonly sequenced in land plant phylogenetic studies, plastid rbcL and internal transcribed spacers of the large subunits of nuclear ribosomal DNA (ITS), and find that both of these DNA regions perform well even though the data currently available in GenBank/EBI were not produced to be used as barcodes and BLAST searches are not an ideal tool for this purpose. These results bode well for the use of even more variable regions of plastid DNA (such as, for example, psbA-trnH) as barcodes, once they have been widely sequenced. In the short term, efforts to bring land plant barcoding up to the standards being used now in other organisms should make swift progress. There are two categories of DNA barcode users, scientists in fields other than taxonomy and taxonomists. For the former, the use of mitochondrial and plastid DNA, the two most easily assessed genomes, is at least in the short term a useful tool that permits them to get on with their studies, which depend on knowing roughly which species or species groups they are dealing with, but these same DNA regions have important drawbacks for use in taxonomic studies (i.e. studies designed to elucidate species limits). For these purposes, DNA markers from uniparentally (usually maternally) inherited genomes can only provide half of the story required to improve taxonomic standards being used in DNA barcoding. In the long term, we will need to develop more sophisticated barcoding tools, which would be multiple, low-copy nuclear markers with sufficient genetic variability and PCR-reliability; these would permit the detection of hybrids and permit researchers to identify the 'genetic gaps' that are useful in assessing species limits.

Phylogeographic analysis of the chloroplast DNA variation in wild common bean (Phaseolus vulgaris L.) in the Americas

The wild common bean (Phaseolus vulgaris) is widely but discontinuously distributed from northern Mexico to northern Argentina on both sides of the Isthmus of Panama. Little is known on how the species has reached its current disjunct distribution. In this research, chloroplast DNA polymorphisms in seven non-coding regions were used to study the history of migration of wild P. vulgaris between Mesoamerica and South America. A penalized likelihood analysis was applied to previously published Leguminosae ITS data to estimate divergence times between P. vulgaris and its sister taxa from Mesoamerica, and divergence times of populations within P. vulgaris. Fourteen chloroplast haplotypes were identified by PCR-RFLP and their geographical associations were studied by means of a Nested Clade Analysis and Mantel Tests. The results suggest that the haplotypes are not randomly distributed but occupy discrete parts of the geographic range of the species. The current distribution of haplotypes may be explained by isolation by distance and by at least two migration events between Mesoamerica and South America: one from Mesoamerica to South America and another one from northern South America to Mesoamerica. Age estimates place the divergence of P. vulgaris from its sister taxa from Mesoamerica at or before 1.3 Ma, and divergence of populations from Ecuador-northern Peru at or before 0.6 Ma. As these ages are taken as minimum divergence times, the influence of past events, such as the closure of the Isthmus of Panama and the final uplift of the Andes, on the migration history and population structure of this species cannot be disregarded.

Effect of corn silage from an herbicide-tolerant genetically modified variety on milk production and absence of transgenic DNA in milk

Data from 60 multiparous Holstein cows were used in a 12-wk continuous design feeding trial. Cows were allocated to 1 of 4 experimental treatments (T1 to T4). In T1 and T2, the total mixed ration (TMR) contained either corn silage from the genetically modified (GM) variety Chardon Liberty Link, which is tolerant to the herbicide glufosinate ammonium, or its near isogenic nonGM counterpart, whereas the TMR used in T3 and T4 contained corn silage from the commercially available nonGM varieties Fabius and Antares, respectively. The objectives of the study were to determine if the inserted gene produced a marked effect on chemical composition, nutritive value, feed intake, and milk production, and to determine if transgenic DNA and the protein expressed by the inserted gene could be detected in bovine milk. The nutritive value, fermentation characteristics, mineral content, and amino acid composition of all 4 silages were similar. There were no significant treatment effects on milk yield, milk composition, and yield of milk constituents, and the dry matter (DM) intake of the GM variety was not significantly different from the 2 commercial varieties. However, although the DM intake noted for the nonGM near-isogenic variety was similar to the commercial varieties, it was significantly lower when compared with the GM variety. Polymerase chain reaction analyses of milk samples collected at wk 1, 6, and 12 of the study showed that none of the 90 milk samples tested positive, above a detection limit of 2.5 ng of total genomic DNA/mL of milk, for either tDNA (event T25) or the single-copy endogenous Zea mays gene, alcohol dehydrogenase. Using ELISA assays, the protein expressed by the T25 gene was not detected in milk.

DNA barcoding of a large genus, Aspalathus L. (Fabaceae)

The DNA barcode potential of three regions (the nuclear ribosomal ITS and the plastid psbA-trnH and trnT-trnL intergenic spacers) was investigated for the plant genus Aspalathus L. (Fabaceac: Crotalarieae). Aspalathus is a large genus (278 species) that revealed low levels of DNA variation in phylogenetic studies. In a 51-species dataset for the psbA-trnH and ITS regions, 45%, and 16% of sequences respectively were identical to the sequence of at least one other species, with two species undiscriminated even when the two regions were combined. In contrast, trnT-trnL, discriminated between all species in this dataset. In a larger ITS and trnT-trnL dataset. including a further 82 species. 7 species in five pairwise comparisons remained Undiscriminated when the two regions were combined. Four of the five pairs of species not discriminated by sequence data were readily distinguished using a combination of qualitative and quantitative morphological data. The difficulty of barcoding in this group is increased by the presence of intraspecific variation in all three regions studied. In the case of psbA-trnH, three intraspecific samples had a sequence identical to at least one other species. Overall, psbA-trnH. currently a candidate for plant barcoding, was the least discriminatory region in our study.

Novel high-resolution characterization of ancient DNA reveals C > U-type base modification events as the sole cause of post mortem miscoding lesions

Ancient DNA (aDNA) research has long depended on the power of PCR to amplify trace amounts of surviving genetic material from preserved specimens. While PCR permits specific loci to be targeted and amplified, in many ways it can be intrinsically unsuited to damaged and degraded aDNA templates. PCR amplification of aDNA can produce highly-skewed distributions with significant contributions from miscoding lesion damage and non-authentic sequence artefacts. As traditional PCR-based approaches have been unable to fully resolve the molecular nature of aDNA damage over many years, we have developed a novel single primer extension (SPEX)-based approach to generate more accurate sequence information. SPEX targets selected template strands at defined loci and can generate a quantifiable redundancy of coverage; providing new insights into the molecular nature of aDNA damage and fragmentation. SPEX sequence data reveals inherent limitations in both traditional and metagenomic PCR-based approaches to aDNA, which can make current damage analyses and correct genotyping of ancient specimens problematic. In contrast to previous aDNA studies, SPEX provides strong quantitative evidence that C U-type base modifications are the sole cause of authentic endogenous damage-derived miscoding lesions. This new approach could allow ancient specimens to be genotyped with unprecedented accuracy.

Effect of corn silage from an herbicide-tolerant genetically modified variety on milk production and absence of transgenic DNA in milk

Data from 60 multiparous Holstein cows were used in a 12-wk continuous design feeding trial. Cows were allocated to 1 of 4 experimental treatments (T1 to T4). In T1 and T2, the total mixed ration (TMR) contained either corn silage from the genetically modified (GM) variety Chardon Liberty Link, which is tolerant to the herbicide glufosinate ammonium, or its near isogenic nonGM counterpart, whereas the TMR used in T3 and T4 contained corn silage from the commercially available nonGM varieties Fabius and Antares, respectively. The objectives of the study were to determine if the inserted gene produced a marked effect on chemical composition, nutritive value, feed intake, and milk production, and to determine if transgenic DNA and the protein expressed by the inserted gene could be detected in bovine milk. The nutritive value, fermentation characteristics, mineral content, and amino acid composition of all 4 silages were similar. There were no significant treatment effects on milk yield, milk composition, and yield of milk constituents, and the dry matter (DM) intake of the GM variety was not significantly different from the 2 commercial varieties. However, although the DM intake noted for the nonGM near-isogenic variety was similar to the commercial varieties, it was significantly lower when compared with the GM variety. Polymerase chain reaction analyses of milk samples collected at wk 1, 6, and 12 of the study showed that none of the 90 milk samples tested positive, above a detection limit of 2.5 ng of total genomic DNA/mL of milk, for either tDNA (event T25) or the single-copy endogenous Zea mays gene, alcohol dehydrogenase. Using ELISA assays, the protein expressed by the T25 gene was not detected in milk.

Phylogeny of Pelargonium (Geraniaceae) based on DNA sequences from three genomes

Phylogenetic hypotheses for the largely South African genus Pelargonium L'Hér. (Geraniaceae) were derived based on DNA sequence data from nuclear, chloroplast and mitochondrial encoded regions. The datasets were unequally represented and comprised cpDNA trnL-F sequences for 152 taxa, nrDNA ITS sequences for 55 taxa, and mtDNA nad1 b/c exons for 51 taxa. Phylogenetic hypotheses derived from the separate three datasets were overall congruent. A single hypothesis synthesising the information in the three datasets was constructed following a total evidence approach and implementing dataset specific stepmatrices in order to correct for substitution biases. Pelargonium was found to consist of five main clades, some with contrasting evolutionary patterns with respect to biogeographic distributions, dispersal capacity, pollination biology and karyological diversification. The five main clades are structured in two (subgeneric) clades that correlate with chromosome size. One of these clades includes a "winter rainfall clade" containing more than 70% of all currently described Pelargonium species, and all restricted to the South African Cape winter rainfall region. Apart from (woody) shrubs and small herbaceous rosette subshrubs, this clade comprises a large "xerophytic" clade including geophytes, stem and leaf succulents, harbouring in total almost half of the genus. This clade is considered to be the result of in situ proliferation, possibly in response to late-Miocene and Pliocene aridification events. Nested within it is a radiation comprising c. 80 species from the geophytic Pelargonium section Hoarea, all characterised by the possession of (a series of) tunicate tubers.

Crops and cultures in the Pacific: new data and new techniques for the investigation of old questions

Fifty years ago Carl Sauer suggested, controversially and on the basis of theory rather than evidence, that Southeast Asia was the source area for agriculture throughout the Old World, including the Pacific. Since then, the archaeobotanical record (macroscopic and microscopic) from the Pacific islands has increased, leading to suggestions, also still controversial, that Melanesia was a center of origin of agriculture independent of South-east Asia, based on tree fruits and nuts and vegetatively propagated starchy staples. Such crops generally lack morphological markers of domestication, so exploitation, cultivation and domestication cannot easily be distinguished in the archaeological record. Molecular studies involving techniques such as chromosome painting, DNA fingerprinting and DNA sequencing, can potentially complement the archaeological record by suggesting where species which were spread through the Pacific by man originated and by what routes they attained their present distributions. A combination of archaeobotanical and molecular studies should therefore eventually enable the rival claims of Melanesia versus South-east Asia as independent centers of invention of agriculture to be assessed.

Phylogenetic relationships within the moss family Bryaceae based on chloroplast DNA evidence

The Bryaceae are a large cosmopolitan family of mosses containing genera of considerable taxonomic difficulty. Phylogenetic relationships within the family were inferred using data from chloroplast DNA sequences (rps4 and trnL-trnF region). Parsimony and maximum likelihood optimality criteria, and Bayesian phylogenetic inference procedures were employed to reconstruct relationships. The genera Bryum and Brachymenium are not monophyletic groups. A clade comprising Plagiobryum, Acidodontium, Mielichhoferia macrocarpa, Bryum sects. Bryum, Apalodictyon, Limbata, Leucodontium, Caespiticia, Capillaria (in part: sect. Capillaria), and Brachymenium sect. Dicranobryum, is well supported in all analyses and represents a major lineage within the family. Section Dicranobryum of Brachymenium is more closely related to section Bryum than to the other sections of Brachymenium, as are Mielichhoferia macrocarpa and M. himalayana. Species of Acidodontium form a clade with Anomobryum julaceum. The grouping of species with a rosulate gametophytic growth form suggests the presence of a 'rosulate' clade similar in circumscription to the genus Rosulabryum. Mielichhoferia macrocarpa and M. himalayana are transferred to Bryum as B. porsildii and B. caucasicum, respectively.

Structural studies on acridine derivatives binding to telomeric DNA

Acridine derivatives can inhibit a variety of nuclear enzymes by binding or intercalating to DNA. This class of compounds is of great interest in the development of novel anticancer agents. Despite the availability of crystallographic data for some of the compounds complexed with DNA, uncertainties remain about the mechanisms of action, binding preferences and biological targets. To investigate the intercalation of several acridine derivatives, a variety of techniques are being employed. Single-crystal X-ray diffraction is being used to determine the high resolution three-dimensional structure of short sequences of quadruplex telomeric DNA with bound drug. This will be compared to the effect of drug binding to long segments of double-stranded DNA using fibre diffraction, with neutron diffraction studies planned to analyse the hydrogen bonding patterns of the DNA-drug complexes. Small-angle neutron scattering (SANS) will also be applied to study drug binding to both short and long sequences of quadruplex and double-stranded DNA in solution. Initial SANS measurements of the telomeric repeat d(TGGGGT) imply that this hexamer is present as a quadruplex. (c) 2006 Elsevier B.V. All rights reserved.

Implications for DNA identification arising from an analysis of Australian forensic databases

Previous analyses of Australian samples have suggested that populations of the same broad racial group (Caucasian, Asian, Aboriginal) tend to be genetically similar across states. This suggests that a single national Australian database for each such group may be feasible, which would greatly facilitate casework. We have investigated samples drawn from each of these groups in different Australian states, and have quantified the genetic homogeneity across states within each racial group in terms of the "coancestry coefficient" F(ST). In accord with earlier results, we find that F(ST) values, as estimated from these data, are very small for Caucasians and Asians, usually <0.5%. We find that "declared" Aborigines (which includes many with partly Aboriginal genetic heritage) are also genetically similar across states, although they display some differentiation from a "pure" Aboriginal population (almost entirely of Aboriginal genetic heritage).