Leveraging base-pair mammalian constraint to understand genetic variation and human disease.

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease.

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Rare tandem repeat expansions associate with genes involved in synaptic and neuronal signaling functions in schizophrenia.

Tandem repeat expansions (TREs) are associated with over 60 monogenic disorders and have recently been implicated in complex disorders such as cancer and autism spectrum disorder. The role of TREs in schizophrenia is now emerging. In this study, we have performed a genome-wide investigation of TREs in schizophrenia. Using genome sequence data from 1154 Swedish schizophrenia cases and 934 ancestry-matched population controls, we have detected genome-wide rare (<0.1% population frequency) TREs that have motifs with a length of 2-20 base pairs. We find that the proportion of individuals carrying rare TREs is significantly higher in the schizophrenia group. There is a significantly higher burden of rare TREs in schizophrenia cases than in controls in genic regions, particularly in postsynaptic genes, in genes overlapping brain expression quantitative trait loci, and in brain-expressed genes that are differentially expressed between schizophrenia cases and controls. We demonstrate that TRE-associated genes are more constrained and primarily impact synaptic and neuronal signaling functions. These results have been replicated in an independent Canadian sample that consisted of 252 schizophrenia cases of European ancestry and 222 ancestry-matched controls. Our results support the involvement of rare TREs in schizophrenia etiology.

A New Method for Detecting Associations with Rare Copy-Number Variants.

Copy number variants (CNVs) play an important role in the etiology of many diseases such as cancers and psychiatric disorders. Due to a modest marginal effect size or the rarity of the CNVs, collapsing rare CNVs together and collectively evaluating their effect serves as a key approach to evaluating the collective effect of rare CNVs on disease risk. While a plethora of powerful collapsing methods are available for sequence variants (e.g., SNPs) in association analysis, these methods cannot be directly applied to rare CNVs due to the CNV-specific challenges, i.e., the multi-faceted nature of CNV polymorphisms (e.g., CNVs vary in size, type, dosage, and details of gene disruption), and etiological heterogeneity (e.g., heterogeneous effects of duplications and deletions that occur within a locus or in different loci). Existing CNV collapsing analysis methods (a.k.a. the burden test) tend to have suboptimal performance due to the fact that these methods often ignore heterogeneity and evaluate only the marginal effects of a CNV feature. We introduce CCRET, a random effects test for collapsing rare CNVs when searching for disease associations. CCRET is applicable to variants measured on a multi-categorical scale, collectively modeling the effects of multiple CNV features, and is robust to etiological heterogeneity. Multiple confounders can be simultaneously corrected. To evaluate the performance of CCRET, we conducted extensive simulations and analyzed large-scale schizophrenia datasets. We show that CCRET has powerful and robust performance under multiple types of etiological heterogeneity, and has performance comparable to or better than existing methods when there is no heterogeneity.

Endogenous retinoids in the pathogenesis of alopecia areata.

Alopecia areata (AA) is an autoimmune disease that attacks anagen hair follicles. Gene array in graft-induced C3H/HeJ mice revealed that genes involved in retinoic acid (RA) synthesis were increased, whereas RA degradation genes were decreased in AA compared with sham controls. This was confirmed by immunohistochemistry in biopsies from patients with AA and both mouse and rat AA models. RA levels were also increased in C3H/HeJ mice with AA. C3H/HeJ mice were fed a purified diet containing one of the four levels of dietary vitamin A or an unpurified diet 2 weeks before grafting and disease progression followed. High vitamin A accelerated AA, whereas mice that were not fed vitamin A had more severe disease by the end of the study. More hair follicles were in anagen in mice fed high vitamin A. Both the number and localization of granzyme B-positive cells were altered by vitamin A. IFNγ was also the lowest and IL13 highest in mice fed high vitamin A. Other cytokines were reduced and chemokines increased as the disease progressed, but no additional effects of vitamin A were seen. Combined, these results suggest that vitamin A regulates both the hair cycle and immune response to alter the progression of AA.

Chromosome neighborhood composition determines translocation outcomes after exposure to high-dose radiation in primary cells.

Radiation exposure is an occupational hazard for military personnel, some health care professionals, airport security screeners, and medical patients, with some individuals at risk for acute, high-dose exposures. Therefore, the biological effects of radiation, especially the potential for chromosome damage, are major occupational and health concerns. However, the biophysical mechanisms of chromosome instability subsequent to radiation-induced DNA damage are poorly understood. It is clear that interphase chromosomes occupy discrete structural and functional subnuclear domains, termed chromosome territories (CT), which may be organized into 'neighborhoods' comprising groups of specific CTs. We directly evaluated the relationship between chromosome positioning, neighborhood composition, and translocation partner choice in primary lymphocytes, using a cell-based system in which we could induce multiple, concentrated DNA breaks via high-dose irradiation. We critically evaluated mis-rejoining profiles and tested whether breaks occurring nearby were more likely to fuse than breaks occurring at a distance. We show that CT neighborhoods comprise heterologous chromosomes, within which inter-CT distances directly relate to translocation partner choice. These findings demonstrate that interphase chromosome arrangement is a principal factor in genomic instability outcomes in primary lymphocytes, providing a structural context for understanding the biological effects of radiation exposure, and the molecular etiology of tumor-specific translocation patterns.

Defective carbohydrate metabolism in mice homozygous for the tubby mutation.

Tub is a member of a small gene family, the tubby-like proteins (TULPs), with predominant expression in neurons. Mice carrying a mutation in Tub develop retinal and cochlear degeneration as well as late-onset obesity with insulin resistance. During behavioral and metabolic testing, we found that homozygous C57BL/6J-Tub(tub) mice have a lower respiratory quotient than C57BL/6J controls before the onset of obesity, indicating that tubby homozygotes fail to activate carbohydrate metabolism and instead rely on fat metabolism for energy needs. In concordance with this, tubby mice show higher excretion of ketone bodies and accumulation of glycogen in the liver. Quantitation of liver mRNA levels shows that, during the transition from light to dark period, tubby mice fail to induce glucose-6-phosphate dehydrogenase (G6pdh), the rate-limiting enzyme in the pentose phosphate pathway that normally supplies NADPH for de novo fatty acid synthesis and glutathione reduction. Reduced G6PDH protein levels and enzymatic activity in tubby mice lead accordingly to lower levels of NADPH and reduced glutathione (GSH), respectively. mRNA levels for the lipolytic enzymes acetyl-CoA synthetase and carnitine palmitoyltransferase are increased during the dark cycle and decreased during the light period, and several citric acid cycle genes are dysregulated in tubby mice. Examination of hypothalamic gene expression showed high levels of preproorexin mRNA leading to accumulation of orexin peptide in the lateral hypothalamus. We hypothesize that abnormal hypothalamic orexin expression leads to changes in liver carbohydrate metabolism and may contribute to the moderate obesity observed in tubby mice.