Spectral Analysis of Guanine and Cytosine Fluctuations of Mouse Genomic DNA

We study global fluctuations of the guanine and cytosine base content (GC%) in mouse genomic DNA using spectral analyses. Power spectra S(f) of GC% fluctuations in all nineteen autosomal and two sex chromosomes are observed to have the universal functional form S(f) \sim 1/f^alpha (alpha \approx 1) over several orders of magnitude in the frequency range 10^-7< f < 10^-5 cycle/base, corresponding to long-ranging GC% correlations at distances between 100 kb and 10 Mb. S(f) for higher frequencies (f > 10^-5 cycle/base) shows a flattened power-law function with alpha < 1 across all twenty-one chromosomes. The substitution of about 38% interspersed repeats does not affect the functional form of S(f), indicating that these are not predominantly responsible for the long-ranged multi-scale GC% fluctuations in mammalian genomes. Several biological implications of the large-scale GC% fluctuation are discussed, including neutral evolutionary history by DNA duplication, chromosomal bands, spatial distribution of transcription units (genes), replication timing, and recombination hot spots.Comment: 15 pages (figures included), 2 figure

Aubergene - a sensitive genome alignment tool.

Motivation: The accumulation of genome sequences will only accelerate in the coming years. We aim to use this abundance of data to improve the quality of genomic alignments and devise a method which is capable of detecting regions evolving under weak or no evolutionary constraints. Results: We describe a genome alignment program AuberGene, which explores the idea of transitivity of local alignments. Assessment of the program was done based on a 2 Mbp genomic region containing the CFTR gene of 13 species. In this region, we can identify 53% of human sequence sharing common ancestry with mouse, as compared with 44% found using the usual pairwise alignment. Between human and tetraodon 93 orthologous exons are found, as compared with 77 detected by the pairwise human-tetraodon comparison. AuberGene allows the user to (1) identify distant, previously undetected, conserved orthogonal regions such as ORFs or regulatory regions; (2) identify neutrally evolving regions in related species which are often overlooked by other alignment programs; (3) recognize false orthologous genomic regions. The increased sensitivity of the method is not obtained at the cost of reduced specificity. Our results suggest that, over the CFTR region, human shares 10% more sequence with mouse than previously thought (∼50%, instead of 40% found with the pairwise alignment). © 2006 Oxford University Press

Copy number variant detection in inbred strains from short read sequence data

Summary: We have developed an algorithm to detect copy number variants (CNVs) in homozygous organisms, such as inbred laboratory strains of mice, from short read sequence data. Our novel approach exploits the fact that inbred mice are homozygous at virtually every position in the genome to detect CNVs using a hidden Markov model (HMM). This HMM uses both the density of sequence reads mapped to the genome, and the rate of apparent heterozygous single nucleotide polymorphisms, to determine genomic copy number. We tested our algorithm on short read sequence data generated from re-sequencing chromosome 17 of the mouse strains A/J and CAST/EiJ with the Illumina platform. In total, we identified 118 copy number variants (43 for A/J and 75 for CAST/EiJ). We investigated the performance of our algorithm through comparison to CNVs previously identified by array-comparative genomic hybridization (array CGH). We performed quantitative-PCR validation on a subset of the calls that differed from the array CGH data sets

How useful are post consultation letters to patients?

BACKGROUND: As part of the NHS plan it was suggested that all patients receive copies of letters sent to their General Practitioner following outpatient consultations. The former Secretary of State for Health extended this proposal, suggesting that patients have a specific letter to themselves after a hospital consultation. METHODS: The aim of this study was to send cardiorespiratory patients attending Charing Cross Hospital, a copy of the letter sent to their G.P. plus a specific letter to themselves and to assess the usefulness and comprehensibility of each. The letters were analysed for dictation time, Flesch Reading Ease Score, Flesch-Kincaid Grade Level and word count. Eighty-four out of 105 sequential patients (80%) consented and were sent both types of letter after their attendance. Patients returned both letters circling any items they did not understand and stated a preference for the GP letter, patient letter, or both. The patients' GPs were subsequently also asked for their views on each letter. RESULTS: GP letters took significantly longer to dictate than patient letters. The Flesch Reading Ease Score was significantly higher in the patient letters, indicating that the patient letters were easier to read. The GP letters were significantly longer than the patient letters and patients were significantly more likely to circle more items in the GP letters (p < 0.001). The content of letters is sometimes inaccurate. Thirty-six out of 62 patients (58%) would like to receive both letters, 13/62 (21.6%) would prefer the GP letter and 13/62 (20%) wanted only the patient letter. 45 GPs replied (62.5%), 28/45 (62.5%) wanted the GP letter, 14 GPs (31.1%) wanted both letters and 3/45 (6.7%) wanted the patient letter only. General themes concerned insufficient clinical details and the GPs preferred the structure of the letters written to them. CONCLUSION: Patients appreciate copies of the letter being sent to their GP but comprehension is less good than with a shorter letter written especially to the patient. More attention needs to be paid to making letters to GPs simpler to read without losing the structure and detail liked by GPs. A compromise might be to dictate the letter in front of the patient and to provide a speciality-specific glossary to accompany each letter

MITE-Hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences

Miniature inverted-repeat transposable elements (MITEs) are a special type of Class 2 non-autonomous transposable element (TE) that are abundant in the non-coding regions of the genes of many plant and animal species. The accurate identification of MITEs has been a challenge for existing programs because they lack coding sequences and, as such, evolve very rapidly. Because of their importance to gene and genome evolution, we developed MITE-Hunter, a program pipeline that can identify MITEs as well as other small Class 2 non-autonomous TEs from genomic DNA data sets. The output of MITE-Hunter is composed of consensus TE sequences grouped into families that can be used as a library file for homology-based TE detection programs such as RepeatMasker. MITE-Hunter was evaluated by searching the rice genomic database and comparing the output with known rice TEs. It discovered most of the previously reported rice MITEs (97.6%), and found sixteen new elements. MITE-Hunter was also compared with two other MITE discovery programs, FINDMITE and MUST. Unlike MITE-Hunter, neither of these programs can search large genomic data sets including whole genome sequences. More importantly, MITE-Hunter is significantly more accurate than either FINDMITE or MUST as the vast majority of their outputs are false-positives

Non-perturbative dynamics of hot non-Abelian gauge fields: beyond leading log approximation

Many aspects of high-temperature gauge theories, such as the electroweak baryon number violation rate, color conductivity, and the hard gluon damping rate, have previously been understood only at leading logarithmic order (that is, neglecting effects suppressed only by an inverse logarithm of the gauge coupling). We discuss how to systematically go beyond leading logarithmic order in the analysis of physical quantities. Specifically, we extend to next-to-leading-log order (NLLO) the simple leading-log effective theory due to Bodeker that describes non-perturbative color physics in hot non-Abelian plasmas. A suitable scaling analysis is used to show that no new operators enter the effective theory at next-to-leading-log order. However, a NLLO calculation of the color conductivity is required, and we report the resulting value. Our NLLO result for the color conductivity can be trivially combined with previous numerical work by G. Moore to yield a NLLO result for the hot electroweak baryon number violation rate.Comment: 20 pages, 1 figur