Biology is encoded in molecular sequences: deciphering this encoding remains a grand scientific challenge. Functional regions of DNA, RNA, and protein sequences often exhibit characteristic but subtle motifs; thus, computational discovery of motifs in sequences is a fundamental and much-studied problem. However, most current algorithms do not allow for insertions or deletions (indels) within motifs, and the few that do have other limitations. We present a method, GLAM2 (Gapped Local Alignment of Motifs), for discovering motifs allowing indels in a fully general manner, and a companion method GLAM2SCAN for searching sequence databases using such motifs. glam2 is a generalization of the gapless Gibbs sampling algorithm. It re-discovers variable-width protein motifs from the PROSITE database significantly more accurately than the alternative methods PRATT and SAM-T2K. Furthermore, it usefully refines protein motifs from the ELM database: in some cases, the refined motifs make orders of magnitude fewer overpredictions than the original ELM regular expressions. GLAM2 performs respectably on the BAliBASE multiple alignment benchmark, and may be superior to leading multiple alignment methods for “motif-like” alignments with N- and C-terminal extensions. Finally, we demonstrate the use of GLAM2 to discover protein kinase substrate motifs and a gapped DNA motif for the LIM-only transcriptional regulatory complex: using GLAM2SCAN, we identify promising targets for the latter. GLAM2 is especially promising for short protein motifs, and it should improve our ability to identify the protein cleavage sites, interaction sites, post-translational modification attachment sites, etc., that underlie much of biology. It may be equally useful for arbitrarily gapped motifs in DNA and RNA, although fewer examples of such motifs are known at present. GLAM2 is public domain software, available for download at http://bioinformatics.org.au/glam2. In recent decades, scientists have extracted genetic sequences—DNA, RNA, and protein sequences—from numerous organisms. These sequences hold the information for the construction and functioning of these organisms, but as yet we are mostly unable to read them. It has long been known that these sequences contain many kinds of “motifs”, i.e. re-occurring patterns, associated with specific biological functions. Thus, much research has been devoted to computer algorithms for automatically discovering subtle, recurring motifs in sequences. However, previous algorithms search for rigid motifs whose instances vary only by substitutions, and not by insertions or deletions. Real motifs are flexible, and do vary by insertions and deletions. This study describes a new computer algorithm for discovering motifs, which allows for arbitrary insertions and deletions. This algorithm can discover real, flexible motifs, and should be able to help us determine the functions of many biological molecules.

We generated draft genome sequences for two cold-adaptedArchaea, Methanogenium frigidumandMethanococcoides burtonii, to identify genotypic characteristics that distinguish them fromArchaeawith a higher optimal growth temperature (OGT). Comparative genomics revealed trends in amino acid and tRNA composition, and structural features of proteins. Proteins from the cold-adaptedArchaeaare characterized by a higher content of noncharged polar amino acids, particularly Gln and Thr and a lower content of hydrophobic amino acids, particularly Leu. Sequence data from nine methanogen genomes (OGT 15°–98°C) were used to generate 1111 modeled protein structures. Analysis of the models from the cold-adaptedArchaeashowed a strong tendency in the solvent-accessible area for more Gln, Thr, and hydrophobic residues and fewer charged residues. A cold shock domain (CSD) protein (CspA homolog) was identified inM. frigidum, two hypothetical proteins with CSD-folds inM. burtonii, and a unique winged helix DNA-binding domain protein inM. burtonii. This suggests that these types of nucleic acid binding proteins have a critical role in cold-adaptedArchaea. Structural analysis of tRNA sequences from theArchaeaindicated that GC content is the major factor influencing tRNA stability in hyperthermophiles, but not in the psychrophiles, mesophiles or moderate thermophiles. Below an OGT of 60°C, the GC content in tRNA was largely unchanged, indicating that any requirement for flexibility of tRNA in psychrophiles is mediated by other means. This is the first time that comparisons have been performed with genome data fromArchaeaspanning the growth temperature extremes from psychrophiles to hyperthermophiles.

Background: The development of colorectal cancer (CRC) is accompanied by extensive epigenetic changes, including frequent regional hypermethylation particularly of gene promoter regions. Specific genes, including SEPT9, VIM1 and TMEFF2 become methylated in a high fraction of cancers and diagnostic assays for detection of cancer-derived methylated DNA sequences in blood and/or fecal samples are being developed. There is considerable potential for the development of new DNA methylation biomarkers or panels to improve the sensitivity and specificity of current cancer detection tests. Methods: Combined epigenomic methods – activation of gene expression in CRC cell lines following DNA demethylating treatment, and two novel methods of genome-wide methylation assessment – were used to identify candidate genes methylated in a high fraction of CRCs. Multiplexed amplicon sequencing of PCR products from bisulfite-treated DNA of matched CRC and non-neoplastic tissue as well as healthy donor peripheral blood was performed using Roche 454 sequencing. Levels of DNA methylation in colorectal tissues and blood were determined by quantitative methylation specific PCR (qMSP). Results: Combined analyses identified 42 candidate genes for evaluation as DNA methylation biomarkers. DNA methylation profiles of 24 of these genes were characterised by multiplexed bisulfite-sequencing in ten matched tumor/normal tissue samples; differential methylation in CRC was confirmed for 23 of these genes. qMSP assays were developed for 32 genes, including 15 of the sequenced genes, and used to quantify methylation in tumor, adenoma and non-neoplastic colorectal tissue and from healthy donor peripheral blood. 24 of the 32 genes were methylated in >50% of neoplastic samples, including 11 genes that were methylated in 80% or more CRCs and a similar fraction of adenomas. Conclusions: This study has characterised a panel of 23 genes that show elevated DNA methylation in >50% of CRC tissue relative to non-neoplastic tissue. Six of these genes (SOX21, SLC6A15, NPY, GRASP, ST8SIA1 and ZSCAN18) show very low methylation in non-neoplastic colorectal tissue and are candidate biomarkers for stool-based assays, while 11 genes (BCAT1, COL4A2, DLX5, FGF5, FOXF1, FOXI2, GRASP, IKZF1, IRF4, SDC2 and SOX21) have very low methylation in peripheral blood DNA and are suitable for further evaluation as blood-based diagnostic markers.

Summary: A global view of the biology of the cold‐adapted archaeon Methanococcoides burtonii was achieved using proteomics. Proteins specific to growth at 4°C versus T_opt (23°C) were identified by mass spectrometry using the draft genome sequence of M. burtonii. mRNA levels were determined for all genes identified by proteomics, and specific enzyme assays confirmed the protein expression results. Key aspects of cold adaptation related to transcription, protein folding and metabolism, including specific roles for RNA polymerase subunit E, a response regulator and peptidyl prolyl cis/trans isomerase. Heat shock protein DnaK was expressed during growth at T_opt, indicating that growth at ‘optimal’ temperatures was stressful for this cold‐adapted organism. Expression of trimethylamine methyltransferase involves contiguous translation of two open reading frames, which is likely to result from incorporation of pyrrolysine at an amber stop codon. Thermal regulation in M. burtonii is achieved through complex gene expression events involving gene clusters and operons, through to protein modifications.