Changes in DNA methylation profiles of myalgic encephalomyelitis/chronic fatigue syndrome patients reflect systemic dysfunctions.

BACKGROUND: Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a lifelong debilitating disease with a complex pathology not yet clearly defined. Susceptibility to ME/CFS involves genetic predisposition and exposure to environmental factors, suggesting an epigenetic association. Epigenetic studies with other ME/CFS cohorts have used array-based technology to identify differentially methylated individual sites. Changes in RNA quantities and protein abundance have been documented in our previous investigations with the same ME/CFS cohort used for this study. RESULTS: DNA from a well-characterised New Zealand cohort of 10 ME/CFS patients and 10 age-/sex-matched healthy controls was isolated from peripheral blood mononuclear (PBMC) cells, and used to generate reduced genome-scale DNA methylation maps using reduced representation bisulphite sequencing (RRBS). The sequencing data were analysed utilising the DMAP analysis pipeline to identify differentially methylated fragments, and the MethylKit pipeline was used to quantify methylation differences at individual CpG sites. DMAP identified 76 differentially methylated fragments and Methylkit identified 394 differentially methylated cytosines that included both hyper- and hypo-methylation. Four clusters were identified where differentially methylated DNA fragments overlapped with or were within close proximity to multiple differentially methylated individual cytosines. These clusters identified regulatory regions for 17 protein encoding genes related to metabolic and immune activity. Analysis of differentially methylated gene bodies (exons/introns) identified 122 unique genes. Comparison with other studies on PBMCs from ME/CFS patients and controls with array technology showed 59% of the genes identified in this study were also found in one or more of these studies. Functional pathway enrichment analysis identified 30 associated pathways. These included immune, metabolic and neurological-related functions differentially regulated in ME/CFS patients compared to the matched healthy controls. CONCLUSIONS: Major differences were identified in the DNA methylation patterns of ME/CFS patients that clearly distinguished them from the healthy controls. Over half found in gene bodies with RRBS in this study had been identified in other ME/CFS studies using the same cells but with array technology. Within the enriched functional immune, metabolic and neurological pathways, a number of enriched neurotransmitter and neuropeptide reactome pathways highlighted a disturbed neurological pathophysiology within the patient group.

ToothPrint, a proteomic database for dental tissues.

Increasing demand exists to disseminate and integrate proteomic data as proteome analysis assumes a commanding role in the postgenome era. Databases on the World Wide Web are an effective means to share information obtained from two-dimensional gels and allied proteomic approaches. Here we report the establishment of ToothPrint, a proteomic database for dental tissues accessed at http://toothprint.otago.ac.nz. Using developing rat enamel as a prototype, ToothPrint provides a variety of functionally relevant data (ligand binding, subcellular localisation, developmental regulation) in addition to protein identification maps. Features designed to enhance usability of the website and simplify its computing requirements are also outlined. Customized for mineralizing tissues, ToothPrint should prove to be an effective bioinformatic resource for investigations of dental biology.

Self-programmable, self-assembling two-dimensional genetic matter.

Putative two-dimensional coding systems can be constructed from aqueous solutions of purine and pyrimidine nucleic acid bases evaporated at moderate temperatures on the surfaces of inorganic solids. The resultant structures are monolayers which are formed spontaneously by molecular self-assembly and they have been observed with molecular resolution by scanning tunnelling microscopy (STM). When formed from solutions of a single base, the monolayers of adenine and uracil have crystalline characteristics and the STM images can be interpreted in terms of the geometrical placement of planar arranged molecules that interact laterally by intermolecular hydrogen bonding. When formed from solutions containing a mixture of adenine and uracil, the monolayers have aperiodic structures. Small crystalline domains within these monolayers can be interpreted in terms of the single phase configurations of the molecules and the remaining aperiodic structures can presumably be interpreted, geometrically, in terms of the 21 theoretically possible adenine-adenine, uracil-uracil and adenine-uracil hydrogen bonding interactions. We propose that combinatorial arrangements of planar arranged purine and pyrimidine bases could provide the necessary complexity to act as a primitive genetic mechanism and may have relevance to the origin of life.

Transterm: a database of messenger RNA components and signals.

Transterm facilitates studies of messenger RNAs and translational control signals. Each messenger RNA (mRNA) from GenBank is extracted and broken into its functional components, its coding sequence, initiation context, termination context, flanking sequence representing its 5' UTR (untranslated region), 3' UTR and translational signals. In addition, numerical parameters characterising each coding region in Transterm, including codon and GC bias, are available. For each species in Transterm, the initiation and termination regions are aligned by their start or stop codons and presented as base frequency matrices and tables of the information content of the bases in the alignments. Users can obtain summaries of characteristics of the mRNAs for species of their choice and search for translational signals both in the Transterm database and in their own sequence. The current release contains data from over 10 000 species, including the complete genomes of 20 prokaryotes and three eukaryotes. Both flat-file and relational database forms of Transterm are accessible via the WWW at http://biochem.otago.ac.nz/Transterm/

TransTerm, the translational signal database, extended to include full coding sequences and untranslated regions.

TransTerm is a database of mRNA sequences and parameters useful for detecting translational control signals in general. TransTerm-98 has been expanded beyond previous years to include full coding sequences and UTRs, while retaining the original small contexts about the coding sequence start- and stop-codons. The database contains more than 130 000 non-redundant coding sequences with associated untranslated regions (UTRs) from over 450 species. This includes the complete genomes of 12 prokaryotic and one eukaryotic organism. Several coding sequence parameters are available: coding sequence length, Nc, GC3 and, when it is computable, Codon Adaptation Index (CAI). Codon usage tables and summaries of start- and stop-codon contexts are also included. TransTerm-98 has both a relational database form with a WWW interface and a flatfile format, also available by Internet browser. TransTerm is available at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml

Is gene deletion in eukaryotes sequence-dependent? A study of nine deletion junctions and nineteen other deletion breakpoints in intron 7 of the human dystrophin gene.

Although large deletions comprise 65% of the mutations that underlie most cases of Duchenne and Becker muscular dystrophies, the DNA sequence characteristics of the deletions and the molecular processes leading to their formation are largely unknown. Intron 7 of the human dystrophin gene is unusually large (110 kb) and a substantial number of deletions have been identified with endpoints within this intron. The distribution of 28 deletion endpoints was mapped to local sequence elements by PCR. The break points were distributed among unique sequence, LINE-1, Alu, MIR, MER and microsatellite sequences with frequencies expected from the frequency of those sequences in the intron. Thus, deletions in this intron are not associated primarily with any one of those sequences in the intron. Nine deletion junctions were amplified and sequenced. Eight were deletions between DNA sequences with minimal homology (0-4 bp) and are therefore unlikely to be products of homologous recombination. In the ninth case, a complex rearrangement was found to be consistent with unequal recombinational exchange between two Alu sequences coupled with a duplication. We have hypothesized that a paucity of matrix attachment regions in this very large intron expanded by the insertion of many mobile elements might provoke a chromatin structure that stimulates deletions (McNaughton et al., 1997, Genomics 40, 294-304). The data presented here are consistent with that idea and demonstrate that the deletion sequences are not usually produced by homologous DNA misalignments.

The translational signal database, TransTerm, is now a relational database.

TransTerm-97 contains more than 97 500 non-redundant coding-sequence initiation and termination contexts compiled from GenBank, release 101 (15-June-1997). In addition, several coding sequence parameters are available: coding sequence length, Nc, GC3, and, when it is computable, codon adaptation index (CAI). Codon usage tables and summaries of start and stop codon contexts are also included. The information covers more than 325 species and organelles, including seven complete bacterial genomes and one complete eukaryotic genome. To promote research in translational control of protein synthesis, TransTerm has been converted into a relational database to ease the process of making queries. The relational database manager, Postgresql, gives access to the database using SQL (Structured Query Language). A World Wide Web interface using forms is being completed to allow the casual user access to the database. Extensions are planned to include the full 5'-UTR, full coding sequence and 3'-UTR. TransTerm-97 is available on the World Wide Web at:http://biochem. otago.ac.nz:800/Transterm/homepage.html

The evolution of an intron: analysis of a long, deletion-prone intron in the human dystrophin gene.

The sequence of a 112-kb region of the human dystrophin (DMD/BMD) gene encompassing the deletion prone intron 7 (110 kb) and the much shorter intron 8 (1.1 kb) has been determined. Recognizable insertion sequences account for approximately 40% of intron 7. LINE-1 and THE-1/LTR sequences occur in intron 7 with significantly higher frequency than would be expected statistically while Alu sequences are underrepresented. Intron 7 also contains numerous mammalian-wide interspersed repeats, a diverse range of medium reiteration repeats of unknown origin, and a sequence derived from a mariner transposon. By contrast, the shorter intron 8 contains no detectable insertion sequences. Dating of the LI and Alu sequences suggests that intron 7 has approximately doubled in size within the past 130 million years, and comparison with the corresponding intron from the pufferfish (Fugu rubripes) suggests that the intron has expanded some 44-fold over a period of 400 million years. The possible contribution of the insertion elements to the instability of intron 7 is discussed.

The translational signal database, TransTerm: more organisms, complete genomes.

TransTerm is a database of initiation and termination sequence contexts from more than 250 organisms listed in GenBank, including the four complete genomes:Haemophilus influenzae, Methanococcus jannaschii, Mycoplasma genitalium,and Saccharomyces cerevisiae. For the current release, more than 60 000 coding sequences were analysed. The tabulated data include initiation and termination contexts organised by species along with quantitative parameters about individual coding sequences (length, %GC, GC3, Nc and CAI). There are also tables of initiation- and termination-region nucleotide-frequencies, codon usage tables and summaries of stop signal usage. TransTerm is available on the World Wide Web at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml

TransTerm: a database of translational signals.

The TransTerm database of sequence contexts of stop and start codons has been expanded to include approximately 50% more species than last year's release. It now contains 148 organisms and >39 500 coding sequences; it is now available on the World Wide Web. The database includes: (i) initiation and termination sequence contexts organized by species; (ii) summary parameters about the individual sequences (sequence length, GC%, GC3, Nc, CAI) in addition to tables of base frequencies for each species' stop and start codon sequence context; (iii) species codon usage tables; and (iv) summary tables of stop signal frequency.

Phylogenetic relationships among transposon-like elements in human and primate DNA.

A THE-1 sequence in intron 7 of the human dystrophin gene has been found to represent a new subfamily of THE-1 elements. The sequence is closely related to the MstII family of repetitive sequences and is more like single-copy sequences found in the galago genome than any other THE-1 sequence previously reported. This new THE-1 sequence has been compared with two other complete THE-1 sequences and three related long-terminal repeat elements that we have previously found in intron 7 of the dystrophin gene, and with members of the same family from elsewhere in the primate genome. Parsimony and deletion analysis show that the cluster of THE-1 sequences in intron 7 of the dystrophin gene has arisen from at least three individual insertion events, rather than from the insertion and duplication of a single progenitor sequence.

Sequence of a transposon identified as Tn1000 (gamma delta).

We report the complete sequence of a transposon found in a cosmid clone of a human DNA sequence. The transposon is identified as the Escherichia coli transposon Tn1000 (also known as gamma delta) on the basis of the identity of the restriction map of the new sequence with that previously recorded for Tn1000 and homology between parts of the new sequence and that of published fragments of Tn1000 sequence. The transposon, which comprises 5,981 nucleotides including two 35 bp inverted terminal repeat sequences (ITRs), contains three open reading frames. The sequence of the resolvase coding region (tnpR) is identical to that published by others. A second reading frame can be identified as the tnpA gene, coding for the transposase, on the grounds of its strong homology with the corresponding gene from transposon Tn3. The third reading frame has the potential to code for a protein of unknown function containing 698 amino acids.

The translational termination signal database (TransTerm) now also includes initiation contexts.

The TransTerm database of termination codon contexts has been extended to include sense codon usage, and initiation codon contexts. The database was constructed from 23,721 coding sequences from 93 organisms. The database contains: a) the sequence around the termination codon (-10, +10); b) the sequence around the initiation codon (-20, +10); c) the length, 'G+C%' of the third position of codons (GC3), the 'codon adaptation index' (CAI) and the 'effective number of codons' statistic (Nc); d) summary tables for each organism including total codon usage, stop codon and tetranucleotide stop-signal usage, and matrices tallying base frequencies at each position around the initiation and termination codons. The data are arranged to facilitate investigation of the relationships between the three phases of protein synthesis. The database is available electronically from EMBL.

The translational termination signal database.

The Translational Termination Database (TransTerm) consists of the immediate context sequences around the natural termination codons from 45 organisms, and summary tables. The influence of termination codon context on their effectivness as stop signals has been widely documented. The SPECIES--TRI.DAT table shows trinucleotide stop codon usage in each organism and for comparison the occurrence of these sequences in the noncoding region. The SPECIES--TETRA.DAT table contains is a similar table of tetranucleotide stop signal usage. The database is available from EMBL.

A cluster of transposon-like repetitive sequences in intron 7 of the human dystrophin gene.

A 32 kilobase-pair fragment of intron 7 of the human dystrophin gene has been sequenced and analysed for the presence of repetitive elements and open reading frames. Two transposon-like human elements (THE-1 sequences), and three intervening, and related, long terminal repeat elements, together with three Alu sequences and a LINE sequence have been identified. These represent an unexpected clustering of highly-repetitive sequences within this single segment of intron DNA. Amplification of a region of chimpanzee genomic DNA by the polymerase chain reaction has provided evidence that at least one of the THE-1 sequences is present in the same position in the chimpanzee genome and the high homology between the human and chimpanzee sequences indicates that this element was fixed within the ancestral genome before the divergence of the two species. The possible role of repetitive, transposon-like sequences in natural mutagenesis of the dystrophin gene is discussed.

The putative single-stranded DNA-binding protein of the filamentous bacteriophage, Ifl. Amino acid sequence of the protein and structure of the gene.

The protein product corresponding to the gene located in the region of the coliphage Ifl genome shown to contain the code for the single-stranded DNA (ssDNA)-binding proteins of all filamentous phages so far studied has been isolated from infected bacterial cells and its amino acid sequence determined. The mature protein contains 95 amino acids (calculated molecular mass 10553 Da). Its sequence corresponds to that predicted from the DNA sequence but lacks the initiating methionine residue. Although there is little direct sequence homology between the phage Ifl protein and the ssDNA-binding proteins of the other filamentous phages that have been studied, computer-based comparisons of various physical and structural parameters showed that the phage Ifl protein contains a domain that is closely related to domains in the coliphage T4 gene 32 protein and the Pseudomonas phage Pfl ssDNA-binding protein and suggest that the Ifl protein does have a ssDNA-binding function although we were unable to show this directly.

Sequence analysis suggests that tetra-nucleotides signal the termination of protein synthesis in eukaryotes.

An increasing number of cases where tri-nucleotide stop codons do not signal the termination of protein synthesis are being reported. In order to identify what constitutes an efficient stop signal, we analysed the region around natural stop codons in genes from a wide variety of eukaryotic species and gene families. Certain stop codons and nucleotides following stop codons are over-represented, and this pattern is accentuated in highly expressed genes. For example, the preferred signal for Saccharomyces cerevisiae and Drosophila melanogaster highly expressed genes is UAAG, and generally the signals UAA(A/G) and UGA(A/G) are preferred in eukaryotes. The GC% of the organism or DNA region can affect whether there is A or G in the second or fourth positions. We suggest therefore, that the stop codon and the nucleotide following it comprise a tetra-nucleotide stop signal. A model is proposed in which the polypeptide chain release factor, a protein, recognises this sequence, but will tolerate some substitution, particularly A to G in the second or third positions.

The signal for the termination of protein synthesis in procaryotes.

The sequences around the stop codons of 862 Escherichia coli genes have been analysed to identify any additional features which contribute to the signal for the termination of protein synthesis. Highly significant deviations from the expected nucleotide distribution were observed, both before and after the stop codon. Immediately prior to UAA stop codons in E. coli there is a preference for codons of the form NAR (any base, adenine, purine), and in particular those that code for glutamine or the basic amino acids. In contrast, codons for threonine or branched nonpolar amino acids were under-represented. Uridine was over-represented in the nucleotide position immediately following all three stop codons, whereas adenine and cytosine were under-represented. This pattern is accentuated in highly expressed genes, but is not as marked in either lowly expressed genes or those that terminate in UAG, the codon specifically recognised by polypeptide chain release factor-1. These observations suggest that for the efficient termination of protein synthesis in E. coli, the 'stop signal' may be a tetranucleotide, rather than simply a tri-nucleotide codon, and that polypeptide chain release factor-2 recognises this extended signal. The sequence following stop codons was analysed in genes from several other procaryotes and bacteriophages. Salmonella typhimurium, Bacillus subtilis, bacteriophages and the methanogenic archaebacteria showed a similar bias to E. coli.

A homologue of retroviral pseudoproteases in the parapoxvirus, orf virus.

The nucleotide sequence of a near-terminal region of orf virus DNA was determined. Examination of the sequence revealed an open reading frame encoding a peptide with significant amino acid homology to the pseudoprotease domains recently identified in a number of retroviruses including mouse mammary tumor virus, simian Mason-Pfizer virus, maedi-visna virus, and equine infectious anaemia virus. The orf virus pseudoprotease shares up to 28% amino acid homology with retroviral pseudoproteases and appears to be a discrete transcriptional unit rather than a subunit of a larger polypeptide as is the case in retroviruses. The sharing of amino acid composition across such wide taxonomic boundaries suggests that this polypeptide has a functional significance in both retroviruses and poxviruses.

Laser excitation of fluorescent-labeled polypeptides in polyacrylamide gels.

A laser beam at 488 nm, converted into a fan of light by a surface-coated mirror oscillated in response to a triangular wave, was inserted into the base of a polyacrylamide gel. The laser light was trapped by internal reflection and gave uniform illumination throughout the entire gel slab. Photography with color film detected 50 fmol of fluorescein covalently coupled to ovalbumin, gave 80-fold greater sensitivity than transillumination in detection of fluorescein-labeled polypeptides, and was about 25-fold more sensitive than protein staining with silver. Laser illumination visualized end-labeled beta-galactosidase, afforded quality control of such preparations, and demonstrated that the end-labeled derivative contained about 25-fold less fluorescein than uniformly labeled beta-galactosidase. The latter result was confirmed by dot-blot analysis using a polyclonal antibody specific for fluorescein. The application of end-labeling to the location of features of protein primary structure is discussed.