DNA methylation signatures reveal that distinct combinations of transcription factors specify human immune cell epigenetic identity.

Epigenetic reprogramming underlies specification of immune cell lineages, but patterns that uniquely define immune cell types and the mechanisms by which they are established remain unclear. Here, we identified lineage-specific DNA methylation signatures of six immune cell types from human peripheral blood and determined their relationship to other epigenetic and transcriptomic patterns. Sites of lineage-specific hypomethylation were associated with distinct combinations of transcription factors in each cell type. By contrast, sites of lineage-specific hypermethylation were restricted mostly to adaptive immune cells. PU.1 binding sites were associated with lineage-specific hypo- and hypermethylation in different cell types, suggesting that it regulates DNA methylation in a context-dependent manner. These observations indicate that innate and adaptive immune lineages are specified by distinct epigenetic mechanisms via combinatorial and context-dependent use of key transcription factors. The cell-specific epigenomics and transcriptional patterns identified serve as a foundation for future studies on immune dysregulation in diseases and aging.

Canine hereditary ataxia in old english sheepdogs and gordon setters is associated with a defect in the autophagy gene encoding RAB24.

Old English Sheepdogs and Gordon Setters suffer from a juvenile onset, autosomal recessive form of canine hereditary ataxia primarily affecting the Purkinje neuron of the cerebellar cortex. The clinical and histological characteristics are analogous to hereditary ataxias in humans. Linkage and genome-wide association studies on a cohort of related Old English Sheepdogs identified a region on CFA4 strongly associated with the disease phenotype. Targeted sequence capture and next generation sequencing of the region identified an A to C single nucleotide polymorphism (SNP) located at position 113 in exon 1 of an autophagy gene, RAB24, that segregated with the phenotype. Genotyping of six additional breeds of dogs affected with hereditary ataxia identified the same polymorphism in affected Gordon Setters that segregated perfectly with phenotype. The other breeds tested did not have the polymorphism. Genome-wide SNP genotyping of Gordon Setters identified a 1.9 MB region with an identical haplotype to affected Old English Sheepdogs. Histopathology, immunohistochemistry and ultrastructural evaluation of the brains of affected dogs from both breeds identified dramatic Purkinje neuron loss with axonal spheroids, accumulation of autophagosomes, ubiquitin positive inclusions and a diffuse increase in cytoplasmic neuronal ubiquitin staining. These findings recapitulate the changes reported in mice with induced neuron-specific autophagy defects. Taken together, our results suggest that a defect in RAB24, a gene associated with autophagy, is highly associated with and may contribute to canine hereditary ataxia in Old English Sheepdogs and Gordon Setters. This finding suggests that detailed investigation of autophagy pathways should be undertaken in human hereditary ataxia.

SLC25A46 Mutations Associated with Autosomal Recessive Cerebellar Ataxia in North African Families.

BACKGROUND: Autosomal recessive cerebellar ataxias (ARCA) are a complex group of neurodegenerative disorders with high clinical and genetic heterogeneity. In most cases, the cerebellar ataxia is not pure, and complicating clinical features such as pyramidal signs or extraneurological features are found. OBJECTIVE: To identify the genetic origin of the cerebellar ataxia for 3 consanguineous North African families presenting with ARCA. METHODS: Genome-wide high-density SNP genotyping and whole-exome sequencing were performed followed by Sanger sequencing for mutation confirmation. RESULTS: Two variants were identified in SLC25A46. Mutations in this gene have been previously associated with Charcot-Marie-Tooth type 2 and optic atrophy. While the previously reported variant p.Arg340Cys seems to be consistently associated with the same clinical features such as childhood onset, optic atrophy, gait and speech difficulties, and wasting of the lower limbs, the patient with the novel mutation p.Trp160Ser did not present with optic atrophy and his ocular abnormalities were limited to nystagmus and saccadic pursuit. CONCLUSION: In this study, we report a novel variant (p.Trp160Ser) in SLC25A46 and we broaden the phenotypic spectrum associated with mutations in SLC25A46.

Mutations in GBA2 cause autosomal-recessive cerebellar ataxia with spasticity.

Autosomal-recessive cerebellar ataxia (ARCA) comprises a large and heterogeneous group of neurodegenerative disorders with more than 20 different forms currently recognized, many of which are also associated with increased tone and some of which have limb spasticity. Gaucher disease is a lysosomal storage disease resulting from a defect in the enzyme acid ß-glucosidase 1. ß-glucosidase 2 is an enzyme with similar glucosylceramidase activity but to date has not been associated with a monogenic disorder. We studied four unrelated consanguineous families of Tunisian decent diagnosed with cerebellar ataxia of unknown origin. We performed homozygosity mapping and whole-exome sequencing in an attempt to identify the genetic origin of their disorder. We were able to identify mutations responsible for autosomal-recessive ataxia in these families within the gene encoding ß-glucosidase 2, GBA2. Two nonsense mutations (c.363C>A [p.Tyr121(∗)] and c.1018C>T [p.Arg340(∗)]) and a substitution (c.2618G>A [p.Arg873His]) were identified, probably resulting in nonfunctional enzyme. This study suggests GBA2 mutations are a cause of recessive spastic ataxia and responsible for a form of glucosylceramide storage disease in humans.

A 7.5-Mb duplication at chromosome 11q21-11q22.3 is associated with a novel spastic ataxia syndrome.

BACKGROUND: The autosomal dominant spinocerebellar ataxias are most commonly caused by nucleotide repeat expansions followed by base-pair changes in functionally important genes. Structural variation has recently been shown to underlie spinocerebellar ataxia types 15 and 20. METHODS: We applied single-nucleotide polymorphism (SNP) genotyping to determine whether structural variation causes spinocerebellar ataxia in a family from France. RESULTS: We identified an approximately 7.5-megabasepair duplication on chromosome 11q21-11q22.3 that segregates with disease. This duplication contains an estimated 44 genes. Duplications at this locus were not found in control individuals. CONCLUSIONS: We have identified a new spastic ataxia syndrome caused by a genomic duplication, which we have denoted as spinocerebellar ataxia type 39. Finding additional families with this phenotype will be important to identify the genetic lesion underlying disease.

Resolving the polymorphism-in-probe problem is critical for correct interpretation of expression QTL studies.

Polymorphisms in the target mRNA sequence can greatly affect the binding affinity of microarray probe sequences, leading to false-positive and false-negative expression quantitative trait locus (QTL) signals with any other polymorphisms in linkage disequilibrium. We provide the most complete solution to this problem, by using the latest genome and exome sequence reference data to identify almost all common polymorphisms (frequency >1% in Europeans) in probe sequences for two commonly used microarray panels (the gene-based Illumina Human HT12 array, which uses 50-mer probes, and exon-based Affymetrix Human Exon 1.0 ST array, which uses 25-mer probes). We demonstrate the impact of this problem using cerebellum and frontal cortex tissues from 438 neuropathologically normal individuals. We find that although only a small proportion of the probes contain polymorphisms, they account for a large proportion of apparent expression QTL signals, and therefore result in many false signals being declared as real. We find that the polymorphism-in-probe problem is insufficiently controlled by previous protocols, and illustrate this using some notable false-positive and false-negative examples in MAPT and PRICKLE1 that can be found in many eQTL databases. We recommend that both new and existing eQTL data sets should be carefully checked in order to adequately address this issue.

MAPT expression and splicing is differentially regulated by brain region: relation to genotype and implication for tauopathies.

The MAPT (microtubule-associated protein tau) locus is one of the most remarkable in neurogenetics due not only to its involvement in multiple neurodegenerative disorders, including progressive supranuclear palsy, corticobasal degeneration, Parksinson's disease and possibly Alzheimer's disease, but also due its genetic evolution and complex alternative splicing features which are, to some extent, linked and so all the more intriguing. Therefore, obtaining robust information regarding the expression, splicing and genetic regulation of this gene within the human brain is of immense importance. In this study, we used 2011 brain samples originating from 439 individuals to provide the most reliable and coherent information on the regional expression, splicing and regulation of MAPT available to date. We found significant regional variation in mRNA expression and splicing of MAPT within the human brain. Furthermore, at the gene level, the regional distribution of mRNA expression and total tau protein expression levels were largely in agreement, appearing to be highly correlated. Finally and most importantly, we show that while the reported H1/H2 association with gene level expression is likely to be due to a technical artefact, this polymorphism is associated with the expression of exon 3-containing isoforms in human brain. These findings would suggest that contrary to the prevailing view, genetic risk factors for neurodegenerative diseases at the MAPT locus are likely to operate by changing mRNA splicing in different brain regions, as opposed to the overall expression of the MAPT gene.

Genome-wide association study of white blood cell count in 16,388 African Americans: the continental origins and genetic epidemiology network (COGENT).

Total white blood cell (WBC) and neutrophil counts are lower among individuals of African descent due to the common African-derived "null" variant of the Duffy Antigen Receptor for Chemokines (DARC) gene. Additional common genetic polymorphisms were recently associated with total WBC and WBC sub-type levels in European and Japanese populations. No additional loci that account for WBC variability have been identified in African Americans. In order to address this, we performed a large genome-wide association study (GWAS) of total WBC and cell subtype counts in 16,388 African-American participants from 7 population-based cohorts available in the Continental Origins and Genetic Epidemiology Network. In addition to the DARC locus on chromosome 1q23, we identified two other regions (chromosomes 4q13 and 16q22) associated with WBC in African Americans (P<2.5×10(-8)). The lead SNP (rs9131) on chromosome 4q13 is located in the CXCL2 gene, which encodes a chemotactic cytokine for polymorphonuclear leukocytes. Independent evidence of the novel CXCL2 association with WBC was present in 3,551 Hispanic Americans, 14,767 Japanese, and 19,509 European Americans. The index SNP (rs12149261) on chromosome 16q22 associated with WBC count is located in a large inter-chromosomal segmental duplication encompassing part of the hydrocephalus inducing homolog (HYDIN) gene. We demonstrate that the chromosome 16q22 association finding is most likely due to a genotyping artifact as a consequence of sequence similarity between duplicated regions on chromosomes 16q22 and 1q21. Among the WBC loci recently identified in European or Japanese populations, replication was observed in our African-American meta-analysis for rs445 of CDK6 on chromosome 7q21 and rs4065321 of PSMD3-CSF3 region on chromosome 17q21. In summary, the CXCL2, CDK6, and PSMD3-CSF3 regions are associated with WBC count in African American and other populations. We also demonstrate that large inter-chromosomal duplications can result in false positive associations in GWAS.

Multiple loci are associated with white blood cell phenotypes.

White blood cell (WBC) count is a common clinical measure from complete blood count assays, and it varies widely among healthy individuals. Total WBC count and its constituent subtypes have been shown to be moderately heritable, with the heritability estimates varying across cell types. We studied 19,509 subjects from seven cohorts in a discovery analysis, and 11,823 subjects from ten cohorts for replication analyses, to determine genetic factors influencing variability within the normal hematological range for total WBC count and five WBC subtype measures. Cohort specific data was supplied by the CHARGE, HeamGen, and INGI consortia, as well as independent collaborative studies. We identified and replicated ten associations with total WBC count and five WBC subtypes at seven different genomic loci (total WBC count-6p21 in the HLA region, 17q21 near ORMDL3, and CSF3; neutrophil count-17q21; basophil count- 3p21 near RPN1 and C3orf27; lymphocyte count-6p21, 19p13 at EPS15L1; monocyte count-2q31 at ITGA4, 3q21, 8q24 an intergenic region, 9q31 near EDG2), including three previously reported associations and seven novel associations. To investigate functional relationships among variants contributing to variability in the six WBC traits, we utilized gene expression- and pathways-based analyses. We implemented gene-clustering algorithms to evaluate functional connectivity among implicated loci and showed functional relationships across cell types. Gene expression data from whole blood was utilized to show that significant biological consequences can be extracted from our genome-wide analyses, with effect estimates for significant loci from the meta-analyses being highly corellated with the proximal gene expression. In addition, collaborative efforts between the groups contributing to this study and related studies conducted by the COGENT and RIKEN groups allowed for the examination of effect homogeneity for genome-wide significant associations across populations of diverse ancestral backgrounds.

Distinct DNA methylation changes highly correlated with chronological age in the human brain.

Methylation at CpG sites is a critical epigenetic modification in mammals. Altered DNA methylation has been suggested to be a central mechanism in development, some disease processes and cellular senescence. Quantifying the extent and identity of epigenetic changes in the aging process is therefore potentially important for understanding longevity and age-related diseases. In the current study, we have examined DNA methylation at >27,000 CpG sites throughout the human genome, in frontal cortex, temporal cortex, pons and cerebellum from 387 human donors between the ages of 1 and 102 years. We identify CpG loci that show a highly significant, consistent correlation between DNA methylation and chronological age. The majority of these loci are within CpG islands and there is a positive correlation between age and DNA methylation level. Lastly, we show that the CpG sites where the DNA methylation level is significantly associated with age are physically close to genes involved in DNA binding and regulation of transcription. This suggests that specific age-related DNA methylation changes may have quite a broad impact on gene expression in the human brain.

Exome sequencing reveals riboflavin transporter mutations as a cause of motor neuron disease.

Brown-Vialetto-Van Laere syndrome was first described in 1894 as a rare neurodegenerative disorder characterized by progressive sensorineural deafness in combination with childhood amyotrophic lateral sclerosis. Mutations in the gene, SLC52A3 (formerly C20orf54), one of three known riboflavin transporter genes, have recently been shown to underlie a number of severe cases of Brown-Vialetto-Van Laere syndrome; however, cases and families with this disease exist that do not appear to be caused by SLC52A3 mutations. We used a combination of linkage and exome sequencing to identify the disease causing mutation in an extended Lebanese Brown-Vialetto-Van Laere kindred, whose affected members were negative for SLC52A3 mutations. We identified a novel mutation in a second member of the riboflavin transporter gene family (gene symbol: SLC52A2) as the cause of disease in this family. The same mutation was identified in one additional subject, from 44 screened. Within this group of 44 patients, we also identified two additional cases with SLC52A3 mutations, but none with mutations in the remaining member of this gene family, SLC52A1. We believe this strongly supports the notion that defective riboflavin transport plays an important role in Brown-Vialetto-Van Laere syndrome. Initial work has indicated that patients with SLC52A3 defects respond to riboflavin treatment clinically and biochemically. Clearly, this makes an excellent candidate therapy for the SLC52A2 mutation-positive patients identified here. Initial riboflavin treatment of one of these patients shows promising results.

Kinesin-associated protein 3 (KIFAP3) has no effect on survival in a population-based cohort of ALS patients.

It was recently reported that rs1541160 on chromosome 1q24.2 has a marked effect on survival of amyotrophic lateral sclerosis (ALS) patients by influencing KIFAP3 expression. The cohorts used in that study were collected from ALS specialty clinics. We attempted to replicate these findings in a population-based cohort of 504 Italian ALS patients. None of 140 SNPs genotyped within the KIFAP3 locus (including rs1541160) had an effect on survival (log-rank P value for rs1541160 = 0.47) or on gene expression in that region. These data illustrate the complexities associated with analyzing ALS phenotypes for association.

Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain.

A fundamental challenge in the post-genome era is to understand and annotate the consequences of genetic variation, particularly within the context of human tissues. We present a set of integrated experiments that investigate the effects of common genetic variability on DNA methylation and mRNA expression in four human brain regions each from 150 individuals (600 samples total). We find an abundance of genetic cis regulation of mRNA expression and show for the first time abundant quantitative trait loci for DNA CpG methylation across the genome. We show peak enrichment for cis expression QTLs to be approximately 68,000 bp away from individual transcription start sites; however, the peak enrichment for cis CpG methylation QTLs is located much closer, only 45 bp from the CpG site in question. We observe that the largest magnitude quantitative trait loci occur across distinct brain tissues. Our analyses reveal that CpG methylation quantitative trait loci are more likely to occur for CpG sites outside of islands. Lastly, we show that while we can observe individual QTLs that appear to affect both the level of a transcript and a physically close CpG methylation site, these are quite rare. We believe these data, which we have made publicly available, will provide a critical step toward understanding the biological effects of genetic variation.

Epigenetic signature of human immune aging in the GESTALT study.

Age-associated DNA methylation in blood cells convey information on health status. However, the mechanisms that drive these changes in circulating cells and their relationships to gene regulation are unknown. We identified age-associated DNA methylation sites in six purified blood-borne immune cell types (naive B, naive CD4+ and CD8+ T cells, granulocytes, monocytes, and NK cells) collected from healthy individuals interspersed over a wide age range. Of the thousands of age-associated sites, only 350 sites were differentially methylated in the same direction in all cell types and validated in an independent longitudinal cohort. Genes close to age-associated hypomethylated sites were enriched for collagen biosynthesis and complement cascade pathways, while genes close to hypermethylated sites mapped to neuronal pathways. In silico analyses showed that in most cell types, the age-associated hypo- and hypermethylated sites were enriched for ARNT (HIF1ß) and REST transcription factor (TF) motifs, respectively, which are both master regulators of hypoxia response. To conclude, despite spatial heterogeneity, there is a commonality in the putative regulatory role with respect to TF motifs and histone modifications at and around these sites. These features suggest that DNA methylation changes in healthy aging may be adaptive responses to fluctuations of oxygen availability.

NeuroBooster Array: A Genome-Wide Genotyping Platform to Study Neurological Disorders Across Diverse Populations.

Genome-wide genotyping platforms have the capacity to capture genetic variation across different populations, but there have been disparities in the representation of population-dependent genetic diversity. The motivation for pursuing this endeavor was to create a comprehensive genome-wide array capable of encompassing a wide range of neuro-specific content for the Global Parkinson's Genetics Program (GP2) and the Center for Alzheimer's and Related Dementias (CARD). CARD aims to increase diversity in genetic studies, using this array as a tool to foster inclusivity. GP2 is the first supported resource project of the Aligning Science Across Parkinson's (ASAP) initiative that aims to support a collaborative global effort aimed at significantly accelerating the discovery of genetic factors contributing to Parkinson's disease and atypical parkinsonism by generating genome-wide data for over 200,000 individuals in a multi-ancestry context. Here, we present the Illumina NeuroBooster array (NBA), a novel, high-throughput and cost-effective custom-designed content platform to screen for genetic variation in neurological disorders across diverse populations. The NBA contains a backbone of 1,914,934 variants (Infinium Global Diversity Array) complemented with custom content of 95,273 variants implicated in over 70 neurological conditions or traits with potential neurological complications. Furthermore, the platform includes over 10,000 tagging variants to facilitate imputation and analyses of neurodegenerative disease-related GWAS loci across diverse populations. The NBA can identify low frequency variants and accurately impute over 15 million common variants from the latest release of the TOPMed Imputation Server as of August 2023 (reference of over 300 million variants and 90,000 participants). We envisage this valuable tool will standardize genetic studies in neurological disorders across different ancestral groups, allowing researchers to perform genetic research inclusively and at a global scale.

Integration of GWAS SNPs and tissue specific expression profiling reveal discrete eQTLs for human traits in blood and brain.

Genome-wide association studies have nominated many genetic variants for common human traits, including diseases, but in many cases the underlying biological reason for a trait association is unknown. Subsets of genetic polymorphisms show a statistical association with transcript expression levels, and have therefore been nominated as expression quantitative trait loci (eQTL). However, many tissue and cell types have specific gene expression patterns and so it is not clear how frequently eQTLs found in one tissue type will be replicated in others. In the present study we used two appropriately powered sample series to examine the genetic control of gene expression in blood and brain. We find that while many eQTLs associated with human traits are shared between these two tissues, there are also examples where blood and brain differ, either by restricted gene expression patterns in one tissue or because of differences in how genetic variants are associated with transcript levels. These observations suggest that design of eQTL mapping experiments should consider tissue of interest for the disease or other traits studied.

A two-stage genome-wide association study of sporadic amyotrophic lateral sclerosis.

The cause of sporadic amyotrophic lateral sclerosis (ALS) is largely unknown, but genetic factors are thought to play a significant role in determining susceptibility to motor neuron degeneration. To identify genetic variants altering risk of ALS, we undertook a two-stage genome-wide association study (GWAS): we followed our initial GWAS of 545 066 SNPs in 553 individuals with ALS and 2338 controls by testing the 7600 most associated SNPs from the first stage in three independent cohorts consisting of 2160 cases and 3008 controls. None of the SNPs selected for replication exceeded the Bonferroni threshold for significance. The two most significantly associated SNPs, rs2708909 and rs2708851 [odds ratio (OR) = 1.17 and 1.18, and P-values = 6.98 x 10(-7) and 1.16 x 10(-6)], were located on chromosome 7p13.3 within a 175 kb linkage disequilibrium block containing the SUNC1, HUS1 and C7orf57 genes. These associations did not achieve genome-wide significance in the original cohort and failed to replicate in an additional independent cohort of 989 US cases and 327 controls (OR = 1.18 and 1.19, P-values = 0.08 and 0.06, respectively). Thus, we chose to cautiously interpret our data as hypothesis-generating requiring additional confirmation, especially as all previously reported loci for ALS have failed to replicate successfully. Indeed, the three loci (FGGY, ITPR2 and DPP6) identified in previous GWAS of sporadic ALS were not significantly associated with disease in our study. Our findings suggest that ALS is more genetically and clinically heterogeneous than previously recognized. Genotype data from our study have been made available online to facilitate such future endeavors.

A genome-wide association study identifies protein quantitative trait loci (pQTLs).

There is considerable evidence that human genetic variation influences gene expression. Genome-wide studies have revealed that mRNA levels are associated with genetic variation in or close to the gene coding for those mRNA transcripts - cis effects, and elsewhere in the genome - trans effects. The role of genetic variation in determining protein levels has not been systematically assessed. Using a genome-wide association approach we show that common genetic variation influences levels of clinically relevant proteins in human serum and plasma. We evaluated the role of 496,032 polymorphisms on levels of 42 proteins measured in 1200 fasting individuals from the population based InCHIANTI study. Proteins included insulin, several interleukins, adipokines, chemokines, and liver function markers that are implicated in many common diseases including metabolic, inflammatory, and infectious conditions. We identified eight Cis effects, including variants in or near the IL6R (p = 1.8x10(-57)), CCL4L1 (p = 3.9x10(-21)), IL18 (p = 6.8x10(-13)), LPA (p = 4.4x10(-10)), GGT1 (p = 1.5x10(-7)), SHBG (p = 3.1x10(-7)), CRP (p = 6.4x10(-6)) and IL1RN (p = 7.3x10(-6)) genes, all associated with their respective protein products with effect sizes ranging from 0.19 to 0.69 standard deviations per allele. Mechanisms implicated include altered rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA), variation in gene copy number (CCL4L1) and altered transcription (GGT1). We identified one novel trans effect that was an association between ABO blood group and tumour necrosis factor alpha (TNF-alpha) levels (p = 6.8x10(-40)), but this finding was not present when TNF-alpha was measured using a different assay , or in a second study, suggesting an assay-specific association. Our results show that protein levels share some of the features of the genetics of gene expression. These include the presence of strong genetic effects in cis locations. The identification of protein quantitative trait loci (pQTLs) may be a powerful complementary method of improving our understanding of disease pathways.