Joseph DeYoung

Evolutionary conservation predicts function of variants of the human organic cation transporter, OCT1

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+), 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 (OCT1-R61C, 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.

Natural variation in human membrane transporter genes reveals evolutionary and functional constraints

Membrane transporters maintain cellular and organismal homeostasis by importing nutrients and exporting toxic compounds. Transporters also play a crucial role in drug response, serving as drug targets and setting drug levels. As part of a pharmacogenetics project, we screened exons and flanking intronic regions for variation in a set of 24 membrane transporter genes (96 kb; 57% coding) in 247 DNA samples from ethnically diverse populations. We identified 680 single nucleotide polymorphisms (SNPs), of which 175 were synonymous and 155 caused amino acid changes, and 29 small insertions and deletions. Amino acid diversity (πNS) in transmembrane domains (TMDs) was significantly lower than in loop domains, suggesting that TMDs have special functional constraints. This difference was especially striking in the ATP-binding cassette superfamily and did not parallel evolutionary conservation: there was little variation in the TMDs, even in evolutionarily unconserved residues. We used allele frequency distribution to evaluate different scoring systems (Grantham, blosum62, SIFT, and evolutionarily conserved/evolutionarily unconserved) for their ability to predict which SNPs affect function. Our underlying assumption was that alleles that are functionally deleterious will be selected against and thus under represented at high frequencies and over represented at low frequencies. We found that evolutionary conservation of orthologous sequences, as assessed by evolutionarily conserved/evolutionarily unconserved and SIFT, was the best predictor of allele frequency distribution and hence of function. European Americans had an excess of high frequency alleles in comparison to African Americans, consistent with a historic bottleneck. In addition, African Americans exhibited a much higher frequency of population specific medium-frequency alleles than did European Americans.

Characterization of mutations in ATP8B1 associated with hereditary cholestasis

Progressive familial intrahepatic cholestasis (PFIC) and benign recurrent intrahepatic cholestasis (BRIC) are clinically distinct hereditary disorders. PFIC patients suffer from chronic cholestasis and develop liver fibrosis. BRIC patients experience intermittent attacks of cholestasis that resolve spontaneously. Mutations in ATP8B1 (previously FIC1) may result in PFIC or BRIC. We report the genomic organization of ATP8B1 and mutation analyses of 180 families with PFIC or BRIC that identified 54 distinct disease mutations, including 10 mutations predicted to disrupt splicing, 6 nonsense mutations, 11 small insertion or deletion mutations predicted to induce frameshifts, 1 large genomic deletion, 2 small inframe deletions, and 24 missense mutations. Most mutations are rare, occurring in 1–3 families, or are limited to specific populations. Many patients are compound heterozygous for 2 mutations. Mutation type or location correlates overall with clinical severity: missense mutations are more common in BRIC (58% vs. 38% in PFIC), while nonsense, frameshifting, and large deletion mutations are more common in PFIC (41% vs. 16% in BRIC). Some mutations, however, lead to a wide range of phenotypes, from PFIC to BRIC or even no clinical disease. ATP8B1 mutations were detected in 30% and 41%, respectively, of the PFIC and BRIC patients screened. Supplementary material for this article can be found on the HEPATOLOGY website (http://interscience.wiley.com/jpages/0270‐9139/suppmat/index.html) and at www.atp8b1‐primers.nl (HEPATOLOGY 2004;40:27–38.)

Polymorphisms in a human kidney xenobiotic transporter, OCT2, exhibit altered function

The completion of the Human Genome Project and the development of high-throughput polymorphism identification methods have allowed researchers to carry out full genetic analyses of many clinically relevant genes. However, few studies have combined genetic analysis with in vitro phenotyping to better understand the relationship between genetic variation and protein function. Many transporters in the kidney are thought to play key roles in defense against a variety of foreign substances. The goal of this study was to understand the relationship between variation in a gene encoding a major renal xenobiotic transporter, OCT2, and transporter function. We report a comprehensive genetic analysis and functional characterization of variants of OCT2. Twenty-eight variable sites in the OCT2 gene were identified in a collection of 247 ethnically diverse DNA samples. Eight caused non-synonymous amino acid changes, of which four were present at >/= 1% in an ethnic population. All four of these altered transporter function assayed in Xenopus laevis oocytes. Analysis of nucleotide diversity (pi) revealed a higher prevalence of synonymous (pi = 22.4 x 10-4) versus non-synonymous (pi = 2.1 x 10-4) changes in OCT2 than in other genes. In addition, the non-synonymous sites had a significant tendency to exhibit more skewed allele frequencies (more negative Tajimas D-values) compared to synonymous sites. The population-genetic analysis, together with the functional characterization, suggests that selection has acted against amino acid changes in OCT2. This selection may be due to a necessary role of OCT2 in the renal elimination of endogenous amines or xenobiotics, including environmental toxins, neurotoxic amines and therapeutic drugs.

Screening a large reference sample to identify very low frequency sequence variants: comparisons between two genes.

Most human sequence variation is in the form of single-nucleotide polymorphisms (SNPs). It has been proposed that coding-region SNPs (cSNPs) be used for direct association studies to determine the genetic basis of complex traits. The success of such studies depends on the frequency of disease-associated alleles, and their distribution in different ethnic populations. If disease-associated alleles are frequent in most populations, then direct genotyping of candidate variants could show robust associations in manageable study samples. This approach is less feasible if the genetic risk from a given candidate gene is due to many infrequent alleles. Previous studies of several genes demonstrated that most variants are relatively infrequent (<0.05). These surveys genotyped small samples (n<75) and thus had limited ability to identify rare alleles. Here we evaluate the prevalence and distribution of such rare alleles by genotyping an ethnically diverse reference sample that is more than six times larger than those used in previous studies (n=450). We screened for variants in the complete coding sequence and intron-exon junctions of two candidate genes for neuropsychiatric phenotypes: SLC6A4, encoding the serotonin transporter; and SLC18A2, encoding the vesicular monoamine transporter. Both genes have unique roles in neuronal transmission, and variants in either gene might be associated with neurobehavioral phenotypes.

Hereditary Vascular Retinopathy, Cerebroretinal Vasculopathy, and Hereditary Endotheliopathy with Retinopathy, Nephropathy, and Stroke Map to a Single Locus on Chromosome 3p21.1-p21.3

We performed a genomewide search for linkage in an extended Dutch family with hereditary vascular retinopathy associated with migraine and Raynaud phenomenon. Patients with vascular retinopathy are characterized by microangiopathy of the retina, accompanied by microaneurysms and telangiectatic capillaries. The genome search, using a high throughput capillary sequencer, revealed significant evidence of linkage to chromosome 3p21.1-p21.3 (maximum pairwise LOD score 5.25, with D3S1578). Testing of two additional families that had a similar phenotype, cerebroretinal vasculopathy, and hereditary endotheliopathy with retinopathy, nephropathy, and stroke, revealed linkage to the same chromosomal region (combined maximum LOD score 6.30, with D3S1588). Haplotype analysis of all three families defined a 3-cM candidate region between D3S1578 and D3S3564. Our study shows that three autosomal dominant vasculopathy syndromes with prominent cerebroretinal manifestations map to the same 3-cM interval on 3p21, suggesting a common locus.

Weighted gene co-expression network analysis of the peripheral blood from Amyotrophic Lateral Sclerosis patients

BackgroundAmyotrophic 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.ResultsOur 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.ConclusionThis 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.

Identification of brain transcriptional variation reproduced in peripheral blood: an approach for mapping brain expression traits

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.

Fine-resolution mapping by haplotype evaluation: the examples of PFIC1 and BRIC

Abstract Loci for two inherited liver diseases, benign recurrent intrahepatic cholestasis (BRIC) and progressive familial intrahepatic cholestasis type 1 (PFIC1), have previously been mapped to 18q21 by a search for shared haplotypes in patients in two isolated populations. This paper describes the use of further haplotype evaluation with a larger sample of patients for both disorders, drawn from several different populations. Our assessment places both loci in the same interval of less than 1 cM and has led to the discovery of the PFIC1/BRIC gene, FIC1; this discovery permits retrospective examination of the general utility of haplotype evaluation and highlights possible caveats regarding this method of genetic mapping.

Genome-wide association study of monoamine metabolite levels in human cerebrospinal fluid

Studying genetic determinants of intermediate phenotypes is a powerful tool to increase our understanding of genotype–phenotype correlations. Metabolic traits pertinent to the central nervous system (CNS) constitute a potentially informative target for genetic studies of intermediate phenotypes as their genetic underpinnings may elucidate etiological mechanisms. We therefore conducted a genome-wide association study (GWAS) of monoamine metabolite (MM) levels in cerebrospinal fluid (CSF) of 414 human subjects from the general population. In a linear model correcting for covariates, we identified one locus associated with MMs at a genome-wide significant level (standardized β=0.32, P=4.92 × 10−8), located 20 kb from SSTR1, a gene involved with brain signal transduction and glutamate receptor signaling. By subsequent whole-genome expression quantitative trait locus (eQTL) analysis, we provide evidence that this variant controls expression of PDE9A (β=0.21; Punadjusted=5.6 × 10−7; Pcorrected=0.014), a gene previously implicated in monoaminergic transmission, major depressive disorder and antidepressant response. A post hoc analysis of loci significantly associated with psychiatric disorders suggested that genetic variation at CSMD1, a schizophrenia susceptibility locus, plays a role in the ratio between dopamine and serotonin metabolites in CSF. The presented DNA and mRNA analyses yielded genome-wide and suggestive associations in biologically plausible genes, two of which encode proteins involved with glutamate receptor functionality. These findings will hopefully contribute to an exploration of the functional impact of the highlighted genes on monoaminergic transmission and neuropsychiatric phenotypes.