Search

Using DNA methylation profiles (n = 15,456) from 348 mammalian species, we constructed phyloepigenetic trees that bear marked similarities to traditional phylogenetic ones. Using unsupervised clustering across all samples, we identified 55 distinct cytosine modules, of which 30 are related to traits such as maximum life span, adult weight, age, sex, and human mortality risk. Maximum life span is associated with methylation levels in HOXL subclass homeobox genes and developmental processes and is potentially regulated by pluripotency transcription factors. The methylation state of some modules responds to perturbations such as caloric restriction, ablation of growth hormone receptors, consumption of high-fat diets, and expression of Yamanaka factors. This study reveals an intertwined evolution of the genome and epigenome that mediates the biological characteristics and traits of different mammalian species.

DNA methylation-based biomarkers of aging have been developed for many mammals but not yet for the vervet monkey (Chlorocebus sabaeus), which is a valuable non-human primate model for biomedical studies. We generated novel DNA methylation data from vervet cerebral cortex, blood, and liver using highly conserved mammalian CpGs represented on a custom array (HorvathMammalMethylChip40). We present six DNA methylation-based estimators of age: vervet multi-tissue epigenetic clock and tissue-specific clocks for brain cortex, blood, and liver. In addition, we developed two dual species clocks (human-vervet clocks) for measuring chronological age and relative age, respectively. Relative age was defined as ratio of chronological age to maximum lifespan to address the species differences in maximum lifespan. The high accuracy of the human-vervet clocks demonstrates that epigenetic aging processes are evolutionary conserved in primates. When applying these vervet clocks to tissue samples from another primate species, rhesus macaque, we observed high age correlations but strong offsets. We characterized CpGs that correlate significantly with age in the vervet. CpG probes that gain methylation with age across tissues were located near the targets of Polycomb proteins SUZ12 and EED and genes possessing the trimethylated H3K27 mark in their promoters. The epigenetic clocks are expected to be useful for anti-aging studies in vervets.

By analyzing multitissue gene expression and genome-wide genetic variation data in samples from a vervet monkey pedigree, we generated a transcriptome resource and produced the first catalog of expression quantitative trait loci (eQTLs) in a nonhuman primate model. This catalog contains more genome-wide significant eQTLs per sample than comparable human resources and identifies sex- and age-related expression patterns. Findings include a master regulatory locus that likely has a role in immune function and a locus regulating hippocampal long noncoding RNAs (lncRNAs), whose expression correlates with hippocampal volume. This resource will facilitate genetic investigation of quantitative traits, including brain and behavioral phenotypes relevant to neuropsychiatric disorders.

Schizophrenia is a common complex disorder with polygenic inheritance. Here we show that by using an approach that compares the individual loads of rare variants in 1,042 schizophrenia cases and 961 controls, schizophrenia cases carry an increased burden of deleterious mutations. At a genome-wide level, our results implicate non-synonymous, splice site as well as stop-altering single-nucleotide variations occurring at minor allele frequency of ≥ 0.01% in the population. In an independent replication sample of 5,585 schizophrenia cases and 8,103 controls of European ancestry we confirm an enrichment in cases of the alleles identified in our study. In addition, the genes implicated by the increased burden of rare coding variants highlight the involvement of neurodevelopment in the aetiology of schizophrenia.

Abstract Background High density genetic maps of plants have, nearly without exception, made use of marker datasets containing missing or questionable genotype calls derived from a variety of genic and non-genic or anonymous markers, and been presented as a single linear order of genetic loci for each linkage group. The consequences of missing or erroneous data include falsely separated markers, expansion of cM distances and incorrect marker order. These imperfections are amplified in consensus maps and problematic when fine resolution is critical including comparative genome analyses and map-based cloning. Here we provide a new paradigm, a high-density consensus genetic map of barley based only on complete and error-free datasets and genic markers, represented accurately by graphs and approximately by a best-fit linear order, and supported by a readily available SNP genotyping resource. Results Approximately 22,000 SNPs were identified from barley ESTs and sequenced amplicons; 4,596 of them were tested for performance in three pilot phase Illumina GoldenGate assays. Data from three barley doubled haploid mapping populations supported the production of an initial consensus map. Over 200 germplasm selections, principally European and US breeding material, were used to estimate minor allele frequency (MAF) for each SNP. We selected 3,072 of these tested SNPs based on technical performance, map location, MAF and biological interest to fill two 1536-SNP "production" assays (BOPA1 and BOPA2), which were made available to the barley genetics community. Data were added using BOPA1 from a fourth mapping population to yield a consensus map containing 2,943 SNP loci in 975 marker bins covering a genetic distance of 1099 cM. Conclusion The unprecedented density of genic markers and marker bins enabled a high resolution comparison of the genomes of barley and rice. Low recombination in pericentric regions is evident from bins containing many more than the average number of markers, meaning that a large number of genes are recombinationally locked into the genetic centromeric regions of several barley chromosomes. Examination of US breeding germplasm illustrated the usefulness of BOPA1 and BOPA2 in that they provide excellent marker density and sensitivity for detection of minor alleles in this genetically narrow material.

Genome-wide gene expression studies may provide substantial insight into gene activities and biological pathways differing between tissues and individuals. We investigated such gene expression variation by analyzing expression profiles in brain tissues derived from eight different brain regions and from blood in 12 monkeys from a biomedically important non-human primate model, the vervet (Chlorocebus aethiops sabaeus). We characterized brain regional differences in gene expression, focusing on transcripts for which inter-individual variation of expression in brain correlates well with variation in blood from the same individuals. Using stringent criteria, we identified 29 transcripts whose expression is measurable, stable, replicable, variable between individuals, relevant to brain function and heritable. Polymorphisms identified in probe regions could, in a minority of transcripts, confound the interpretation of the observed inter-individual variation. The high heritability of levels of these transcripts in a large vervet pedigree validated our approach of focusing on transcripts that showed higher inter-individual compared with intra-individual variation. These selected transcripts are candidate expression Quantitative Trait Loci, differentially regulating transcript levels in the brain among individuals. Given the high degree of conservation of tissue expression profiles between vervets and humans, our findings may facilitate the understanding of regional and individual transcriptional variation and its genetic mechanisms in humans. The approach employed here-utilizing higher quality tissue and more precise dissection of brain regions than is usually possible in humans-may therefore provide a powerful means to investigate variation in gene expression relevant to complex brain related traits, including human neuropsychiatric diseases.

Abstract Background Amyotrophic Lateral Sclerosis (ALS) is a lethal disorder characterized by progressive degeneration of motor neurons in the brain and spinal cord. Diagnosis is mainly based on clinical symptoms, and there is currently no therapy to stop the disease or slow its progression. Since access to spinal cord tissue is not possible at disease onset, we investigated changes in gene expression profiles in whole blood of ALS patients. Results Our transcriptional study showed dramatic changes in blood of ALS patients; 2,300 probes (9.4%) showed significant differential expression in a discovery dataset consisting of 30 ALS patients and 30 healthy controls. Weighted gene co-expression network analysis (WGCNA) was used to find disease-related networks (modules) and disease related hub genes. Two large co-expression modules were found to be associated with ALS. Our findings were replicated in a second (30 patients and 30 controls) and third dataset (63 patients and 63 controls), thereby demonstrating a highly significant and consistent association of two large co-expression modules with ALS disease status. Ingenuity Pathway Analysis of the ALS related module genes implicates enrichment of functional categories related to genetic disorders, neurodegeneration of the nervous system and inflammatory disease. The ALS related modules contain a number of candidate genes possibly involved in pathogenesis of ALS. Conclusion This first large-scale blood gene expression study in ALS observed distinct patterns between cases and controls which may provide opportunities for biomarker development as well as new insights into the molecular mechanisms of the disease.

Abstract Background High density genetic maps of plants have, nearly without exception, made use of marker datasets containing missing or questionable genotype calls derived from a variety of genic and non-genic or anonymous markers, and been presented as a single linear order of genetic loci for each linkage group. The consequences of missing or erroneous data include falsely separated markers, expansion of cM distances and incorrect marker order. These imperfections are amplified in consensus maps and problematic when fine resolution is critical including comparative genome analyses and map-based cloning. Here we provide a new paradigm, a high-density consensus genetic map of barley based only on complete and error-free datasets and genic markers, represented accurately by graphs and approximately by a best-fit linear order, and supported by a readily available SNP genotyping resource. Results Approximately 22,000 SNPs were identified from barley ESTs and sequenced amplicons; 4,596 of them were tested for performance in three pilot phase Illumina GoldenGate assays. Data from three barley doubled haploid mapping populations supported the production of an initial consensus map. Over 200 germplasm selections, principally European and US breeding material, were used to estimate minor allele frequency (MAF) for each SNP. We selected 3,072 of these tested SNPs based on technical performance, map location, MAF and biological interest to fill two 1536-SNP "production" assays (BOPA1 and BOPA2), which were made available to the barley genetics community. Data were added using BOPA1 from a fourth mapping population to yield a consensus map containing 2,943 SNP loci in 975 marker bins covering a genetic distance of 1099 cM. Conclusion The unprecedented density of genic markers and marker bins enabled a high resolution comparison of the genomes of barley and rice. Low recombination in pericentric regions is evident from bins containing many more than the average number of markers, meaning that a large number of genes are recombinationally locked into the genetic centromeric regions of several barley chromosomes. Examination of US breeding germplasm illustrated the usefulness of BOPA1 and BOPA2 in that they provide excellent marker density and sensitivity for detection of minor alleles in this genetically narrow material.

Abstract Background Amyotrophic Lateral Sclerosis (ALS) is a lethal disorder characterized by progressive degeneration of motor neurons in the brain and spinal cord. Diagnosis is mainly based on clinical symptoms, and there is currently no therapy to stop the disease or slow its progression. Since access to spinal cord tissue is not possible at disease onset, we investigated changes in gene expression profiles in whole blood of ALS patients. Results Our transcriptional study showed dramatic changes in blood of ALS patients; 2,300 probes (9.4%) showed significant differential expression in a discovery dataset consisting of 30 ALS patients and 30 healthy controls. Weighted gene co-expression network analysis (WGCNA) was used to find disease-related networks (modules) and disease related hub genes. Two large co-expression modules were found to be associated with ALS. Our findings were replicated in a second (30 patients and 30 controls) and third dataset (63 patients and 63 controls), thereby demonstrating a highly significant and consistent association of two large co-expression modules with ALS disease status. Ingenuity Pathway Analysis of the ALS related module genes implicates enrichment of functional categories related to genetic disorders, neurodegeneration of the nervous system and inflammatory disease. The ALS related modules contain a number of candidate genes possibly involved in pathogenesis of ALS. Conclusion This first large-scale blood gene expression study in ALS observed distinct patterns between cases and controls which may provide opportunities for biomarker development as well as new insights into the molecular mechanisms of the disease.

The organic cation transporter, OCT1, is a major hepatic transporter that mediates the uptake of many organic cations from the blood into the liver where the compounds may be metabolized or secreted into the bile. Because OCT1 interacts with a variety of structurally diverse organic cations, including clinically used drugs as well as toxic substances (e.g., N-methylpyridinium, MPP[sup +]), it is an important determinant of systemic exposure to many xenobiotics. To understand the genetic basis of extensive interindividual differences in xenobiotic disposition, we functionally characterized 15 protein-altering variants of the human liver organic cation transporter, OCT1, in Xenopus oocytes. All variants that reduced or eliminated function (OCT1R61C, OCT1-P341L, OCT1-G220V, OCT1-G401S, and OCT1-G465R) altered evolutionarily conserved amino acid residues. In general, variants with decreased function had amino acid substitutions that resulted in more radical chemical changes (higher Grantham values) and were less evolutionarily favorable (lower BLOSUM62 values) than variants that maintained function. A variant with increased function (OCT1-S14F) changed an amino acid residue such that the human protein matched the consensus of the OCT1 mammalian orthologs. Our results indicate that changes at evolutionarily conserved positions of OCT1 are strong predictors of decreased function and suggest that a combination of evolutionary conservation and chemical change might be a stronger predictor of function. [ABSTRACT FROM AUTHOR]

The gene encoding the human muscarinic receptor, type 1 ( CHRM1), was genotyped from 245 samples of the Coriell Collection (Coriell Institute for Medical Research, Camden, NJ). Fifteen single nucleotide polymorphisms (SNPs) were discovered, 9 of which are located in the coding region of the receptor. Of these, 8 represent synonymous SNPs, indicating that CHRM1 is highly conserved in humans. Only a single allele was found to contain a nonsynonymous SNP, which encodes an amino acid change of Cys to Arg at position 417. This may have functional consequences because a C417S point mutation in rat M1 was previously shown to affect receptor binding and coupling. Furthermore, 0 of 4 SNPs within CHRM1 previously deduced from sequencing of the human genome were found in this study despite a prediction that a majority of such inferred SNPs are accurate. The consensus sequence of CHRM1 obtained in our study differs from the deposited reference sequence (AC NM_000738) in 2 adjacent nucleotides, leading to a V173M change, suggesting a sequencing error in the reference sequence. The extraordinary sequence conservation of the CHRM1 gene-coding region was unexpected as M1-knockout mice show only minimal functional impairments. [ABSTRACT FROM AUTHOR]

Using DNA methylation profiles ( n = 15,456) from 348 mammalian species, we constructed phyloepigenetic trees that bear marked similarities to traditional phylogenetic ones. Using unsupervised clustering across all samples, we identified 55 distinct cytosine modules, of which 30 are related to traits such as maximum life span, adult weight, age, sex, and human mortality risk. Maximum life span is associated with methylation levels in HOXL subclass homeobox genes and developmental processes and is potentially regulated by pluripotency transcription factors. The methylation state of some modules responds to perturbations such as caloric restriction, ablation of growth hormone receptors, consumption of high-fat diets, and expression of Yamanaka factors. This study reveals an intertwined evolution of the genome and epigenome that mediates the biological characteristics and traits of different mammalian species.

The multidisciplinary UCSF Pharmacogenetics of Membrane Transporters project seeks to systematically identify sequence variants in transporters and to determine the functional significance of these variants through evaluation of relevant cellular and clinical phenotypes. The project is structured around four interacting cores: genomics, cellular phenotyping, clinical phenotyping, and bioinformatics. The bioinformatics core is responsible for collecting, storing, and analyzing the information obtained by the other cores and for presenting the results, in particular, for the genomic data. Most of this process is automated using locally developed software written in Python, an open source language well suited for rapid, modular development that meets requirements that are themselves constantly evolving. Here we present the details of transforming ABI trace file data into useful information for project investigators and a description of the types of data analysis and display that we have developed.