PubMed Health⌕ Search

Biomedical subjects

Terence P Speed

Publications and source records attributed to Terence P Speed.

At least 19 recordsLinked to original sources

Exploration, normalization, and genotype calls of high-density oligonucleotide SNP array data.

In most microarray technologies, a number of critical steps are required to convert raw intensity measurements into the data relied upon by data analysts, biologists, and clinicians. These data manipulations, referred to as preprocessing, can influence the quality of the ultimate measurements. In the last few years, the high-throughput measurement of gene expression is the most popular application of microarray technology. For this application, various groups have demonstrated that the use of modern statistical methodology can substantially improve accuracy and precision of the gene expression measurements, relative to ad hoc procedures introduced by designers and manufacturers of the technology. Currently, other applications of microarrays are becoming more and more popular. In this paper, we describe a preprocessing methodology for a technology designed for the identification of DNA sequence variants in specific genes or regions of the human genome that are associated with phenotypes of interest such as disease. In particular, we describe a methodology useful for preprocessing Affymetrix single-nucleotide polymorphism chips and obtaining genotype calls with the preprocessed data. We demonstrate how our procedure improves existing approaches using data from 3 relatively large studies including the one in which large numbers of independent calls are available. The proposed methods are implemented in the package oligo available from Bioconductor.

Algorithms↗

Evolution and comparative analysis of the MHC Class III inflammatory region.

BACKGROUND: The Major Histocompatibility Complex (MHC) is essential for immune function. Historically, it has been subdivided into three regions (Class I, II, and III), but a cluster of functionally related genes within the Class III region has also been referred to as the Class IV region or "inflammatory region". This group of genes is involved in the inflammatory response, and includes members of the tumour necrosis family. Here we report the sequencing, annotation and comparative analysis of a tammar wallaby BAC containing the inflammatory region. We also discuss the extent of sequence conservation across the entire region and identify elements conserved in evolution. RESULTS: Fourteen Class III genes from the tammar wallaby inflammatory region were characterised and compared to their orthologues in other vertebrates. The organisation and sequence of genes in the inflammatory region of both the wallaby and South American opossum are highly conserved compared to known genes from eutherian ("placental") mammals. Some minor differences separate the two marsupial species. Eight genes within the inflammatory region have remained tightly clustered for at least 360 million years, predating the divergence of the amphibian lineage. Analysis of sequence conservation identified 354 elements that are conserved. These range in size from 7 to 431 bases and cover 15.6% of the inflammatory region, representing approximately a 4-fold increase compared to the average for vertebrate genomes. About 5.5% of this conserved sequence is marsupial-specific, including three cases of marsupial-specific repeats. Highly Conserved Elements were also characterised. CONCLUSION: Using comparative analysis, we show that a cluster of MHC genes involved in inflammation, including TNF, LTA (or its putative teleost homolog TNF-N), APOM, and BAT3 have remained together for over 450 million years, predating the divergence of mammals from fish. The observed enrichment in conserved sequences within the inflammatory region suggests conservation at the transcriptional regulatory level, in addition to the functional level.

Animals↗

Functional genomic analysis of oligodendrocyte differentiation.

To better understand the molecular mechanisms governing oligodendrocyte (OL) differentiation, we have used gene profiling to quantitatively analyze gene expression in synchronously differentiating OLs generated from pure oligodendrocyte precursor cells in vitro. By comparing gene expression in these OLs to OLs generated in vivo, we discovered that the program of OL differentiation can progress normally in the absence of heterologous cell-cell interactions. In addition, we found that OL differentiation was unexpectedly prolonged and occurred in at least two sequential stages, each characterized by changes in distinct complements of transcription factors and myelin proteins. By disrupting the normal dynamic expression patterns of transcription factors regulated during OL differentiation, we demonstrated that these sequential stages of gene expression can be independently controlled. We also uncovered several genes previously uncharacterized in OLs that encode transmembrane, secreted, and cytoskeletal proteins that are as highly upregulated as myelin genes during OL differentiation. Last, by comparing genomic loci associated with inherited increased risk of multiple sclerosis (MS) to genes regulated during OL differentiation, we identified several new positional candidate genes that may contribute to MS susceptibility. These findings reveal a previously unexpected complexity to OL differentiation and suggest that an intrinsic program governs successive phases of OL differentiation as these cells extend and align their processes, ensheathe, and ultimately myelinate axons.

Animals↗

Proximal genomic localization of STAT1 binding and regulated transcriptional activity.

BACKGROUND: Signal transducer and activator of transcription (STAT) proteins are key regulators of gene expression in response to the interferon (IFN) family of anti-viral and anti-microbial cytokines. We have examined the genomic relationship between STAT1 binding and regulated transcription using multiple tiling microarray and chromatin immunoprecipitation microarray (ChIP-chip) experiments from public repositories. RESULTS: In response to IFN-gamma, STAT1 bound proximally to regions of the genome that exhibit regulated transcriptional activity. This finding was consistent between different tiling microarray platforms, and between different measures of transcriptional activity, including differential binding of RNA polymerase II, and differential mRNA transcription. Re-analysis of tiling microarray data from a recent study of IFN-gamma-induced STAT1 ChIP-chip and mRNA expression revealed that STAT1 binding is tightly associated with localized mRNA transcription in response to IFN-gamma. Close relationships were also apparent between STAT1 binding, STAT2 binding, and mRNA transcription in response to IFN-alpha. Furthermore, we found that sites of STAT1 binding within the Encyclopedia of DNA Elements (ENCODE) region are precisely correlated with sites of either enhanced or diminished binding by the RNA polymerase II complex. CONCLUSION: Together, our results indicate that STAT1 binds proximally to regions of the genome that exhibit regulated transcriptional activity. This finding establishes a generalized basis for the positioning of STAT1 binding sites within the genome, and supports a role for STAT1 in the direct recruitment of the RNA polymerase II complex to the promoters of IFN-gamma-responsive genes.

Binding Sites↗

SNPs on chips: the hidden genetic code in expression arrays.

Gene expression microarray analysis in postmortem brains is one of the fastest growing fields of psychiatric research. Here we show that common polymorphisms (SNPs) present on probe sets can masquerade as significant "gene expression" differences. After first observing this artifact in the Catechol-O-methyl transferase (COMT) gene, we replicate the finding in two additional genes predicted to show this artifact. Many Affymetrix chips contain thousands of SNPs that are both common and in the central probe region affecting hybridization, and thus have the potential to confound expression analysis.

Analysis of Variance↗

Identification and stoichiometry of glycosylphosphatidylinositol-anchored membrane proteins of the human malaria parasite Plasmodium falciparum.

Most proteins that coat the surface of the extracellular forms of the human malaria parasite Plasmodium falciparum are attached to the plasma membrane via glycosylphosphatidylinositol (GPI) anchors. These proteins are exposed to neutralizing antibodies, and several are advanced vaccine candidates. To identify the GPI-anchored proteome of P. falciparum we used a combination of proteomic and computational approaches. Focusing on the clinically relevant blood stage of the life cycle, proteomic analysis of proteins labeled with radioactive glucosamine identified GPI anchoring on 11 proteins (merozoite surface protein (MSP)-1, -2, -4, -5, -10, rhoptry-associated membrane antigen, apical sushi protein, Pf92, Pf38, Pf12, and Pf34). These proteins represent approximately 94% of the GPI-anchored schizont/merozoite proteome and constitute by far the largest validated set of GPI-anchored proteins in this organism. Moreover MSP-1 and MSP-2 were present in similar copy number, and we estimated that together these proteins comprise approximately two-thirds of the total membrane-associated surface coat. This is the first time the stoichiometry of MSPs has been examined. We observed that available software performed poorly in predicting GPI anchoring on P. falciparum proteins where such modification had been validated by proteomics. Therefore, we developed a hidden Markov model (GPI-HMM) trained on P. falciparum sequences and used this to rank all proteins encoded in the completed P. falciparum genome according to their likelihood of being GPI-anchored. GPI-HMM predicted GPI modification on all validated proteins, on several known membrane proteins, and on a number of novel, presumably surface, proteins expressed in the blood, insect, and/or pre-erythrocytic stages of the life cycle. Together this work identified 11 and predicted a further 19 GPI-anchored proteins in P. falciparum.

Animals↗

Mind the gap: analysis of marker-assisted breeding strategies for inbred mouse strains.

The development of congenic mouse strains is the principal approach for confirming and fine mapping quantitative trait loci, as well as for comparing the phenotypic effect of a transgene or gene-targeted disruption between different inbred mouse strains. The traditional breeding scheme calls for at least nine consecutive backcrosses before establishing a congenic mouse strain. Recent availability of genome sequence and high-throughput genotyping now permit the use of polymorphic DNA markers to reduce this number of backcrosses, and empirical data suggest that marker-assisted breeding may require as few as four backcrosses. We used simulation studies to investigate the efficiency of different marker-assisted breeding schemes by examining the trade-off between the number of backcrosses, the number of mice produced per generation, and the number of genotypes per mouse required to achieve a quality congenic mouse strain. An established model of crossover interference was also incorporated into these simulations. The quality of the strain produced was assessed by the probability of an undetected region of heterozygosity (i.e., "gaps") in the recipient genetic background, while maintaining the desired donor-derived interval. Somewhat surprisingly, we found that there is a relatively high probability for undetected gaps in potential breeders for establishing a congenic mouse strain. Marker-assisted breeding may decrease the number of backcross generations required to generate a congenic strain, but only additional backcrossing will guarantee a reduction in the number and length of undetected gaps harboring contaminating donor alleles.

Animals↗

Expression profiling in primates reveals a rapid evolution of human transcription factors.

Although it has been hypothesized for thirty years that many human adaptations are likely to be due to changes in gene regulation, almost nothing is known about the modes of natural selection acting on regulation in primates. Here we identify a set of genes for which expression is evolving under natural selection. We use a new multi-species complementary DNA array to compare steady-state messenger RNA levels in liver tissues within and between humans, chimpanzees, orangutans and rhesus macaques. Using estimates from a linear mixed model, we identify a set of genes for which expression levels have remained constant across the entire phylogeny (approximately 70 million years), and are therefore likely to be under stabilizing selection. Among the top candidates are five genes with expression levels that have previously been shown to be altered in liver carcinoma. We also find a number of genes with similar expression levels among non-human primates but significantly elevated or reduced expression in the human lineage, features that point to the action of directional selection. Among the gene set with a human-specific increase in expression, there is an excess of transcription factors; the same is not true for genes with increased expression in chimpanzee.

Animals↗

Lineage-specific expansion of proteins exported to erythrocytes in malaria parasites.

BACKGROUND: The apicomplexan parasite Plasmodium falciparum causes the most severe form of malaria in humans. After invasion into erythrocytes, asexual parasite stages drastically alter their host cell and export remodeling and virulence proteins. Previously, we have reported identification and functional analysis of a short motif necessary for export of proteins out of the parasite and into the red blood cell. RESULTS: We have developed software for the prediction of exported proteins in the genus Plasmodium, and identified exported proteins conserved between malaria parasites infecting rodents and the two major causes of human malaria, P. falciparum and P. vivax. This conserved 'exportome' is confined to a few subtelomeric chromosomal regions in P. falciparum and the synteny of these and surrounding regions is conserved in P. vivax. We have identified a novel gene family PHIST (for Plasmodium helical interspersed subtelomeric family) that shares a unique domain with 72 paralogs in P. falciparum and 39 in P. vivax; however, there is only one member in each of the three species studied from the P. berghei lineage. CONCLUSION: These data suggest radiation of genes encoding remodeling and virulence factors from a small number of loci in a common Plasmodium ancestor, and imply a closer phylogenetic relationship between the P. vivax and P. falciparum lineages than previously believed. The presence of a conserved 'exportome' in the genus Plasmodium has important implications for our understanding of both common mechanisms and species-specific differences in host-parasite interactions, and may be crucial in developing novel antimalarial drugs to this infectious disease.

Amino Acid Sequence↗

The comparative roles of suppressor of cytokine signaling-1 and -3 in the inhibition and desensitization of cytokine signaling.

Negative feedback is a mechanism commonly employed in biological processes as a means of maintaining homeostasis. We have investigated the roles of suppressor of cytokine signaling (SOCS) proteins in regulating the kinetics of negative feedback in response to cytokine signaling. In mouse livers and bone marrow-derived macrophages, both interferon-gamma (IFNgamma) and interleukin-6 (IL-6) rapidly induced the tyrosine phosphorylation of signal transducer and activator of transcription-1 (STAT1) and STAT3. STAT3 tyrosine phosphorylation was bi-phasic in response to continuous IL-6 signaling. In macrophages lacking Socs3, however, continuous IL-6 signaling induced uniformly high levels of STAT3 tyrosine phosphorylation, and early IL-6-inducible genes were inappropriately expressed at intermediate time points. SOCS3 therefore imposes bi-phasic kinetics upon IL-6 signaling. Compared with Socs3 mRNA, Socs1 mRNA was induced relatively slowly, and SOCS1 simply attenuated the duration of IFNgamma signaling. Surprisingly, heightened Socs1 mRNA expression but minimal STAT1 tyrosine phosphorylation was observed after prolonged stimulation with IFNgamma, indicating that STAT1 may not play a large role in inducing Socs1 mRNA during steady-state IFNgamma signaling. We also demonstrate that both SOCS1 and SOCS3 can desensitize primary bone marrow-derived macrophages to IFNgamma and IL-6 signaling, respectively. Consistent with the kinetics with which Socs1 and Socs3 mRNAs were induced, SOCS3 desensitized cells to IL-6 rapidly, whereas SOCS1-mediated desensitization to IFNgamma occurred at later time points. The kinetics with which SOCS proteins are induced by cytokine may therefore be a parameter that is "hard-wired" into specific cytokine signaling pathways as a means of tailoring the kinetics with which cells become desensitized.

Animals↗

Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex.

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used the Monodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral "immune supercomplex" that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.

Animals↗

On the utility of data from the International HapMap Project for Australian association studies.

We compare patterns of linkage disequilibrium (LD) for 633 SNPs in two regions between samples collected in two Australian states and HapMap samples collected from Utah residents of Northern and Western (NW) European ancestry (CEU). Patterns of LD in the Australian and HapMap samples are similar, and tag SNPs chosen using HapMap genotypes perform almost as well on Australian samples as tags chosen using Australian genotypes.

Europe↗

Regulation of apicomplexan actin-based motility.

Apicomplexan parasites are an ancient group of protozoan parasites that includes several significant pathogens of humans and animals. To target and invade host cells they use a unique form of actin-based motility, called gliding motility. At the centre of the molecular motor that underlies this unique mode of locomotion are short, highly dynamic actin filaments. Recent molecular work, along with the availability of completed genomes for several Apicomplexa, has highlighted unique features of parasite actin and its regulation - features that might provide new ways to block motility and, consequently, prevent infection and disease.

Actins↗

Mapping of the Plasmodium chabaudi resistance locus char2.

Animals congenic for the char2 host response locus to the murine malarial parasite Plasmodium chabaudi have been bred, and they demonstrated a phenotypic difference from the parental lines. These congenic lines have been crossed back to the parental line to generate recombinants across the congenic intervals. The recombinants were inbred, and the subcongenic intervals were fixed. These lines were then challenged with parasites and assessed as being either resistant or susceptible. From the analysis of many subcongenic lines, it has become obvious that there are at least two loci underlying the char2 locus and that both of these mediate resistance when the haplotype derives from the resistant C57BL/6 strain.

Animals↗

Detecting genome wide haplotype sharing using SNP or microsatellite haplotype data.

Genome wide association studies using high throughput technology are already being conducted despite the significant hurdles that need to be overcome (Nat Rev Genet 6:95-108, 2005; Nat Rev Genet 6:109-118, 2005). Methods for detecting haplotype association signals in genome wide haplotype datasets are as yet very limited. Much methodological research has already been devoted to linkage disequilibrium (LD) fine mapping where the focus is the identification of the disease locus rather than the detection of a disease signal. Applications of these approaches to genome wide scanning are limited by the strong model assumptions of the sharing process, which lead to computational complexity. We describe a new algorithm for the initial identification of disease susceptibility loci in genome wide haplotype association studies. Excess sharing of ancestral haplotypes, which indicates the presence of a disease locus, is detected with a simple, easy to interpret, chi2 based statistic. The method allows genome wide scanning for qualitative traits within reasonable computational timeframes and can serve as a first pass analysis prior to the usage of likelihood based methods, providing candidate regions and inferred susceptibility haplotypes. Our method makes no assumptions regarding the population history or the pattern of background LD. Statistical significance is evaluated with permutation tests. The method is illustrated on simulated and real data where it is applied to simple (cystic fibrosis) and complex disease (multiple sclerosis) examples. The statistic has low type I error and greater power to map disease loci over conventional single marker tests for low to moderate levels of LD.

Algorithms↗

Incorporating interference into linkage analysis for experimental crosses.

The phenomenon of interference in genetic recombination is well-known and studied in a wide variety of organisms. Multilocus linkage analysis, which makes use of recombination patterns among all genetic markers simultaneously, is routinely used with data on humans and experimental organisms to build genetic maps. It is also used to try to determine the genes involved in traits of interest, such as common diseases. Most linkage analyses performed today ignore the occurrence of genetical interference. We present an extension to the Lander-Green algorithm for experimental crosses (backcross and intercross) to incorporate crossover interference according to the chi2 model. Simulation results show the impact of using this model on the accuracy of estimated genetic maps.

Algorithms↗

Long-term effect of therapeutic laser photocoagulation on gene expression in the eye.

Microarray-based gene expression analysis demonstrated that laser photocoagulation (LPC) of mouse eyes had a long-term effect on the expression of genes functionally related to tissue repair, cell migration, proliferation, ion, protein and nucleic acid metabolism, cell signaling, and angiogenesis. Six structural genes, including five crystallins (Cryaa, Cryba1, Crybb2, Crygc, Crygs) and keratin 1-12 (Krt1-12), the anti-angiogenic factor thrombospondin 1 (Tsp1), the retina- and brain-specific putative transcription factor tubby-like protein 1 (Tulp1), and transketolase (Tkt), a key enzyme in the pentose-phosphate pathway, were all shown to be up-regulated by real-time PCR and/or Western blotting. Immunohistochemistry localized five of these proteins to the laser lesions and surrounding tissue within the retina and pigmented epithelium. This is the first study demonstrating long-term changes in the expression of these genes associated with LPC. Therefore, it suggests that modulated gene expression might contribute to the long-term inhibitory effect of LPC. In addition, these genes present novel targets for gene-based therapies aimed at treating microangiopathies, especially diabetic retinopathy, a disease currently only treatable with LPC.

Animals↗

A genotype calling algorithm for affymetrix SNP arrays.

MOTIVATION: A classification algorithm, based on a multi-chip, multi-SNP approach is proposed for Affymetrix SNP arrays. Current procedures for calling genotypes on SNP arrays process all the features associated with one chip and one SNP at a time. Using a large training sample where the genotype labels are known, we develop a supervised learning algorithm to obtain more accurate classification results on new data. The method we propose, RLMM, is based on a robustly fitted, linear model and uses the Mahalanobis distance for classification. The chip-to-chip non-biological variance is reduced through normalization. This model-based algorithm captures the similarities across genotype groups and probes, as well as across thousands of SNPs for accurate classification. In this paper, we apply RLMM to Affymetrix 100 K SNP array data, present classification results and compare them with genotype calls obtained from the Affymetrix procedure DM, as well as to the publicly available genotype calls from the HapMap project.

Algorithms↗