Search PubMed⌕ Search

Biomedical subjects

R Mott

Publications and source records attributed to R Mott.

At least 19 recordsLinked to original sources

Using progenitor strain information to identify quantitative trait nucleotides in outbred mice.

We have developed a fast and economical strategy for dissecting the genetic architecture of quantitative trait loci at a molecular level. The method uses two pieces of information: mapping data from crosses that involve more than two inbred strains and sequence variants in the progenitor strains within the interval containing a quantitative trait locus (QTL). By testing whether the strain distribution pattern in the progenitor strains is consistent with the observed genetic effect of the QTL we can assign a probability that any sequence variant is a quantitative trait nucleotide (QTN). It is not necessary to genotype the animals except at a skeleton of markers; the genotypes at all other polymorphisms are estimated by a multipoint analysis. We apply the method to a 4.8-Mb region on mouse chromosome 1 that contains a QTL influencing anxiety segregating in a heterogeneous stock and show that, under the assumption that a single QTN is present and lies in a region conserved between the human and mouse genomes, it is possible to reduce the number of variants likely to be the quantitative trait nucleotide from many thousands to <20.

Animals↗

Population-specific patterns of linkage disequilibrium in the human 5q31 region.

Linkage disequilibrium across the human genome is generally lower in West Africans than Europeans. However in the 5q31 region, which is rich in immune genes, we find significantly more examples of apparent nonrecombination between distant marker pairs in West Africans. Much of this effect is due to SNPs that are absent in Europeans, possibly reflecting recent positive selection in the West African population.

Black People↗

Unexpected complexity in the haplotypes of commonly used inbred strains of laboratory mice.

Investigation of sequence variation in common inbred mouse strains has revealed a segmented pattern in which regions of high and low variant density are intermixed. Furthermore, it has been suggested that allelic strain distribution patterns also occur in well defined blocks and consequently could be used to map quantitative trait loci (QTL) in comparisons between inbred strains. We report a detailed analysis of polymorphism distribution in multiple inbred mouse strains over a 4.8-megabase region containing a QTL influencing anxiety. Our analysis indicates that it is only partly true that the genomes of inbred strains exist as a patchwork of segments of sequence identity and difference. We show that the definition of haplotype blocks is not robust and that methods for QTL mapping may fail if they assume a simple block-like structure.

Alleles↗

Complex haplotypic structure of the central MHC region flanking TNF in a West African population.

TNF polymorphisms have been associated with susceptibility to malaria and other infectious and inflammatory conditions. We investigated a sample of 150 West African chromosomes to determine linkage disequilibrium (LD) between 25 SNP markers located in an 80 kb segment of the MHC Class III region encompassing TNF and eight neighbouring genes. We observed 45 haplotypes, and 22 of them comprise 80% of the sample. The pattern of LD is remarkably patchy, such that many markers show no LD with adjacent markers but high LD with markers that are much further away. We introduce a method of examining the implications of LD data for disease association studies based on sample size considerations: this shows that certain TNF polymorphisms would be likely to yield positive associations if the true disease allele resided in LTA or BAT1. We conclude that detailed marker maps are needed to resolve the causal origin of disease associations observed at the TNF locus.

Adult↗

Finding the molecular basis of quantitative traits: successes and pitfalls.

Understanding the molecular basis of quantitative genetic variation is a principal goal for biomedicine. Although the complex genetic architecture of quantitative traits has so far largely frustrated attempts to identify genes in humans by standard linkage methodologies, quantitative trait loci (QTL) have been mapped in plants, insects and rodents. However, identifying the molecular bases of QTL remains a challenge. Here, we discuss why this is and how new experimental strategies and analytical techniques, combined with the fruits of the genome projects, are beginning to identify candidate genes for QTL studies in several model organisms.

Animals↗

Automatic analysis of agarose gel images.

MOTIVATION: Automatic tools to speed up routine biological processes are very much sought after in bio-medical research. Much repetitive work in molecular biology, such as allele calling in genetic analysis, can be made semi-automatic or task specific automatic by using existing techniques from computer science and signal processing. Computerized analysis is reproducible and avoids various forms of human error. Semi-automatic techniques with an interactive check on the results speed up the analysis and reduce the error. RESULTS: We have successfully implemented an image processing software package to automatically analyze agarose gel images of polymorphic DNA markers. We have obtained up to 90% accuracy for the classification of alleles in good quality images and up to 70% accuracy in average quality images. These results are obtained within a few seconds. Even after subsequent interactive checking to increase the accuracy of allele classification to 100%, the overall speed with which the data can be processed is greatly increased, compared to manual allele classification. AVAILABILITY: The IDL source code of the software is available on request from jonathan.flint@well.ox.ac.uk

Animals↗

Novel protein domains and repeats in Drosophila melanogaster: insights into structure, function, and evolution.

Sequence database searching methods such as BLAST, are invaluable for predicting molecular function on the basis of sequence similarities among single regions of proteins. Searches of whole databases however, are not optimized to detect multiple homologous regions within a single polypeptide. Here we have used the prospero algorithm to perform self-comparisons of all predicted Drosophila melanogaster gene products. Predicted repeats, and their homologs from all species, were analyzed further to detect hitherto unappreciated evolutionary relationships. Results included the identification of novel tandem repeats in the human X-linked retinitis pigmentosa type-2 gene product, repeated segments in cystinosin, associated with a defect in cystine transport, and 'nested' homologous domains in dysferlin, whose gene is mutated in limb girdle muscular dystrophy. Novel signaling domain families were found that may regulate the microtubule-based cytoskeleton and ubiquitin-mediated proteolysis, respectively. Two families of glycosyl hydrolases were shown to contain internal repetitions that hint at their evolution via a piecemeal, modular approach. In addition, three examples of fruit fly genes were detected with tandem exons that appear to have arisen via internal duplication. These findings demonstrate how completely sequenced genomes can be exploited to further understand the relationships between molecular structure, function, and evolution.

Amino Acid Sequence↗

A method for fine mapping quantitative trait loci in outbred animal stocks.

High-resolution mapping of quantitative trait loci (QTL) in animals has proved to be difficult because the large effect sizes detected in crosses between inbred strains are often caused by numerous linked QTLs, each of small effect. In a study of fearfulness in mice, we have shown it is possible to fine map small-effect QTLs in a genetically heterogeneous stock (HS). This strategy is a powerful general method of fine mapping QTLs, provided QTLs detected in crosses between inbred strains that formed the HS can be reliably detected in the HS. We show here that single-marker association analysis identifies only two of five QTLs expected to be segregating in the HS and apparently limits the strategy's usefulness for fine mapping. We solve this problem with a multipoint analysis that assigns the probability that an allele descends from each progenitor in the HS. The analysis does not use pedigrees but instead requires information about the HS founder haplotypes. With this method we mapped all three previously undetected loci [chromosome (Chr.) 1 logP 4.9, Chr. 10 logP 6.0, Chr. 15 logP 4.0]. We show that the reason for the failure of single-marker association to detect QTLs is its inability to distinguish opposing phenotypic effects when they occur on the same marker allele. We have developed a robust method of fine mapping QTLs in genetically heterogeneous animals and suggest it is now cost effective to undertake genomewide high-resolution analysis of complex traits in parallel on the same set of mice.

Animals↗

Accurate formula for P-values of gapped local sequence and profile alignments.

A simple general approximation for the distribution of gapped local alignment scores is presented, suitable for assessing significance of comparisons between two protein sequences or a sequence and a profile. The approximation takes account of the scoring scheme (i.e. gap penalty and substitution matrix or profile), sequence composition and length. Use of this formula means it is unnecessary to fit an extreme-value distribution to simulations or to the results of databank searches. The method is based on the theoretical ideas introduced by R. Mott and R. Tribe in 1999. Extensive simulation studies show that score-thresholds produced by the method are accurate to within +/-5 % 95 % of the time. We also investigate factors which effect the accuracy of alignment statistics, and show that any method based on asymptotic theory is limited because asymptotic behaviour is not strictly achieved for many real protein sequences, due to extreme composition effects. Consequently, it may not be practicable to find a general formula that is significantly more accurate until the sub-asymptotic behaviour of alignments is better understood.

Computer Simulation↗

Approximate statistics of gapped alignments.

A heuristic approximation to the score distribution of gapped alignments in the logarithmic domain is presented. The method applies to comparisons between random, unrelated protein sequences, using standard score matrices and arbitrary gap penalties. It is shown that gapped alignment behavior is essentially governed by a single parameter, alpha, depending on the penalty scheme and sequence composition. This treatment also predicts the position of the transition point between logarithmic and linear behavior. The approximation is tested by simulation and shown to be accurate over a range of commonly used substitution matrices and gap-penalties.

Amino Acid Sequence↗

Local sequence alignments with monotonic gap penalties.

MOTIVATION: Sequence alignments obtained using affine gap penalties are not always biologically correct, because the insertion of long gaps is over-penalised. There is a need for an efficient algorithm which can find local alignments using non-linear gap penalties. RESULTS: A dynamic programming algorithm is described which computes optimal local sequence alignments for arbitrary, monotonically increasing gap penalties, i.e. where the cost g(k) of inserting a gap of k symbols is such that g(k) >/= g(k-1). The running time of the algorithm is dependent on the scoring scheme; if the expected score of an alignment between random, unrelated sequences of lengths m, n is proportional to log mn, then with one exception, the algorithm has expected running time O(mn). Elsewhere, the running time is no greater than O(mn(m+n)). Optimisations are described which appear to reduce the worst-case run-time to O(mn) in many cases. We show how using a non-affine gap penalty can dramatically increase the probability of detecting a similarity containing a long gap. AVAILABILITY: The source code is available to academic collaborators under licence.

Algorithms↗

Comparative gene expression profiling by oligonucleotide fingerprinting.

The use of hybridisation of synthetic oligonucleotides to cDNAs under high stringency to characterise gene sequences has been demonstrated by a number of groups. We have used two cDNA libraries of 9 and 12 day mouse embryos (24 133 and 34 783 clones respectively) in a pilot study to characterise expressed genes by hybridisation with 110 hybridisation probes. We have identified 33 369 clusters of cDNA clones, that ranged in representation from 1 to 487 copies (0.7%). 737 were assigned to known rodent genes, and a further 13 845 showed significant homologies. A total of 404 clusters were identified as significantly differentially represented (P < 0.01) between the two cDNA libraries. This study demonstrates the utility of the fingerprinting approach for the generation of comparative gene expression profiles through the analysis of cDNAs derived from different biological materials.

Animals↗

Trace alignment and some of its applications.

MOTIVATION: Extra useful information can be extracted from a DNA chromatogram trace, over that contained in the base-called DNA sequence. Many sequencing applications can benefit from examination of these traces. RESULTS: An algorithm, based on dynamic programming, for aligning a DNA chromatogram to a DNA sequence is described and implemented. Its applications to vector clipping, EST alignment and mutation detection are discussed.

Algorithms↗

Sequence assembly with CAFTOOLS.

Large-scale genomic sequencing requires a software infrastructure to support and integrate applications that are not directly compatible. We describe a suite of software tools built around the Common Assembly Format (CAF), a comprehensive representation of a sequence assembly as a text file. These tools form the backbone of sequencing informatics at the Sanger Centre and the Genome Sequencing Center. The CAF format is intentionally flexible, and our Perl and C libraries, which parse and manipulate it, provide powerful tools for creating new applications as well as wrappers to incorporate other software. The tools are available free by anonymous FTP from ftp://ftp.sanger.ac.uk/pub/badger/.

Algorithms↗

Instability of highly expanded CAG repeats in mice transgenic for the Huntington's disease mutation.

Six inherited neurodegenerative diseases are caused by a CAG/polyglutamine expansion, including spinal and bulbar muscular atrophy (SBMA), Huntington's disease (HD), spinocerebellar ataxia type 1 (SCA1), dentatorubral pallidoluysian atrophy (DRPLA) Machado-Joseph disease (MJD or SCA3) and SCA2. Normal and expanded HD allele sizes of 6-39 and 35-121 repeats have been reported, and the allele distributions for the other diseases are comparable. Intergenerational instability has been described in all cases, and repeats tend to be more unstable on paternal transmission. This may present as larger increases on paternal inheritance as in HD, or as a tendency to increase on male and decrease on female transmission as in SCA1 (ref. 15). Somatic repeat instability is also apparent and appears most pronounced in the CNS. The major exception is the cerebellum, which in HD, DRPLA, SCA1 and MJD has a smaller repeat relative to the other brain regions tested. Of non-CNS tissues, instability was observed in blood, liver, kidney and colon. A mouse model of CAG repeat instability would be helpful in unravelling its molecular basis although an absence of CAG repeat instability in transgenic mice has so far been reported. These studies include (CAG) in the androgen receptor cDNA, (CAG) in the HD cDNA, (CAG) in the SCA1 cDNA, (CAG) in the SCA3 cDNA and as an isolated (CAG) tract.

Animals↗

FPC: a system for building contigs from restriction fingerprinted clones.

MOTIVATION: To meet the demands of large-scale sequencing, thousands of clones must be fingerprinted and assembled into contigs. To determine the order of clones, a typical experiment is to digest the clones with one or more restriction enzymes and measure the resulting fragments. The probability of two clones overlapping is based on the similarity of their fragments. A contig contains two or more overlapping clones and a minimal tiling path of clones is selected to be sequenced. Interactive software with algorithmic support is necessary to assemble the clones into contigs quickly. RESULTS: FPC (fingerprinted contigs) is an interactive program for building contigs from restriction fingerprinted clones. FPC uses an algorithm to cluster clones into contigs based on their probability of coincidence score. For each contig, it builds a consensus band (CB) map which is similar to a restriction map; but it does not try to resolve all the errors. The CB map is used to assign coordinates to the clones based on their alignment to the map and to provide a detailed visualization of the clone overlap. FPC has editing facilities for the user to refine the coordinates and to remove poorly fingerprinted clones. Functions are available for updating an FPC database with new clones. Contigs can easily be merged, split or deleted. Markers can be added to clones and are displayed with the appropriate contig. Sequence-ready clones can be selected and their sequencing status displayed. As such, FPC is an integrated program for the assembly of sequence-ready clones for large-scale sequencing projects.

Algorithms↗

Cloning and expression of cystolic phospholipase A2 (cPLA2) and a naturally occurring variant. Phosphorylation of Ser505 of recombinant cPLA2 by p42 mitogen-activated protein kinase results in an increase in specific activity.

Full-length cytosolic phospholipase A2 (cPLA2) was cloned from U937 cells and polymorphonuclear leukocytes (PMNLs) while a naturally occurring variant of cPLA2, which lacks residues Val473-Ala749 but has a C-terminal extension of ILMNLSEYMLWMSKVKRFM (DcPLA2) was cloned from PMNLs and mononuclear leukocytes. We were unable to clone DcPLA2 from U937 cells. When cPLA2 and DcPLA2 were expressed in insect cells, both proteins were detected in cell lysates by SDS/PAGE as single bands of apparent molecular masses 100 kDa and 57 kDa, respectively. Full-length cPLA2 was phosphorylated stoichiometrically by p42 mitogen-activated protein (MAP) kinase in vitro at a similar rate to other physiological substrates of this protein kinase and the major site of phosphorylation was identified by amino acid sequencing as Ser505. [32P]Ser(P)505 in cPLA2 was only dephosphorylated at a slow rate by mammalian tissue homogenates. Protein phosphatases 2A, 2B and 2C all contributed significantly to the overall dephosphorylation of cPLA2. The phosphorylation of cPLA2 by p42 MAP kinase correlated with an approximately 1.5-fold increase in specific enzyme activity which was reversed by dephosphorylation.

Amino Acid Sequence↗