Candidate Regulatory Variants Research Articles

Single-nucleus multi-omics of Parkinson’s disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks

The genetic architecture of Parkinson’s disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

FORGEdb: a tool for identifying candidate functional variants and uncovering target genes and mechanisms for complex diseases.

The majority of disease-associated variants identified through genome-wide association studies are located outside of protein-coding regions. Prioritizing candidate regulatory variants and gene targets to identify potential biological mechanisms for further functional experiments can be challenging. To address this challenge, we developed FORGEdb ( https://forgedb.cancer.gov/ ; https://forge2.altiusinstitute.org/files/forgedb.html ; and https://doi.org/10.5281/zenodo.10067458 ), a standalone and web-based tool that integrates multiple datasets, delivering information on associated regulatory elements, transcription factor binding sites, and target genes for over 37 million variants. FORGEdb scores provide researchers with a quantitative assessment of the relative importance of each variant for targeted functional experiments.

Sequenced-based GWAS for linear classification traits in Belgian Blue beef cattle reveals new coding variants in genes regulating body size in mammals.

Cohorts of individuals that have been genotyped and phenotyped for genomic selection programs offer the opportunity to better understand genetic variation associated with complex traits. Here, we performed an association study for traits related to body size and muscular development in intensively selected beef cattle. We leveraged multiple trait information to refine and interpret the significant associations. After a multiple-step genotype imputation to the sequence-level for 14,762 Belgian Blue beef (BBB) cows, we performed a genome-wide association study (GWAS) for 11 traits related to muscular development and body size. The 37 identified genome-wide significant quantitative trait loci (QTL) could be condensed in 11 unique QTL regions based on their position. Evidence for pleiotropic effects was found in most of these regions (e.g., correlated association signals, overlap between credible sets (CS) of candidate variants). Thus, we applied a multiple-trait approach to combine information from different traits to refine the CS. In several QTL regions, we identified strong candidate genes known to be related to growth and height in other species such as LCORL-NCAPG or CCND2. For some of these genes, relevant candidate variants were identified in the CS, including three new missense variants in EZH2, PAPPA2 and ADAM12, possibly two additional coding variants in LCORL, and candidate regulatory variants linked to CCND2 and ARMC12. Strikingly, four other QTL regions associated with dimension or muscular development traits were related to five (recessive) deleterious coding variants previously identified. Our study further supports that a set of common genes controls body size across mammalian species. In particular, we added new genes to the list of those associated with height in both humans and cattle. We also identified new strong candidate causal variants in some of these genes, strengthening the evidence of their causality. Several breed-specific recessive deleterious variants were identified in our QTL regions, probably as a result of the extreme selection for muscular development in BBB cattle.

Bioinformatics pipeline to identify candidate regulatory variants in late‐onset Alzheimer’s disease (LOAD) associated regions

AbstractBackgroundSmall insertions and deletions (indels) in the human genome are substantial contributors to genetic variation and impact human traits and diseases. The role of indels in late onset Alzheimer’s disease has been understudied. Few examples, such as the intronic poly‐T variant in the TOMM40 gene, suggest that systematic exploration of indels within LOAD risk regions will advance the understanding of the genetic architecture of LOAD. Previously we developed a bioinformatics pipeline that characterizes and prioritizes candidate regulatory SNPs in enhancers located in LOAD‐GWAS regions. Here we extend the pipeline to the analysis of indels. The proposed bioinformatics pipeline progresses from indels located in LOAD‐GWAS regions to a filtered set of regulatory variants that have a predicted strong effect on transcription factor (TFs) binding.MethodThe pipeline utilized publicly available functional genomics data sources. Primarily, candidate cis‐regulatory elements (cCREs) from ENCODE and single‐cell RNA‐seq data from LOAD patient samples (synapse: syn22079621). For TF binding analysis we employed motifs from MotifDb. In addition, we used various bioinformatics software including motifbreakR.ResultWe catalogued 1230 proximal CTCF‐bound candidate cis‐regulatory elements in LOAD‐GWAS regions, 912 showed epigenetic evidence in relevant brain tissue. We catalogued 426 indels in these cCREs. These indels disrupted 391 TFs, 362 of these had snRNA‐seq data from LOAD samples. Of note, TF motifs within the APOE‐TOMM40, SPI1 and MS4A2 regions were significantly disrupted by the candidate regulatory indels. Amongst these TFs are RUNX3, SPI1 and SMAD3. Interestingly, these significant findings with the APOE‐TOMM40, SPI1 and MS4A2 regions are consistent with our prior results for SNPs.ConclusionThis study provides an analytical framework to catalogue noncoding indel variation in cis‐regulatory elements located in LOAD‐GWAS loci and characterize their likelihood to perturb TF binding. The approach integrates multiple data types to prioritize genes and variants for validation experiments using disease models and gene editing technologies.

Genetic variants of the hypoxia-inducible factor 3 alpha subunit (Hif3a) gene in the Fat and Lean mouse selection lines.

Adipose tissue hypoxia and members of the hypoxia-inducible factor alpha (HIFA) are involved in development of obesity. However, the mechanism and functions of HIF3A, one of three HIFA paralogs, in fat deposition have not been sufficiently studied. In the present study, we investigated whether Hif3a sequence variants are associated with divergent fat deposition in mouse selection lines for fatness and leanness. Sequencing and RFLP were used to analyse sequence variants within Hif3a. To identify candidate regulatory variants, we performed literature screening and used databases and bioinformatics tools like Ensembl, MethPrimer, TargetScanMouse, miRDB, PolyAsite, RISE, LncRRIsearch, RNAfold, PredictProtein, CAIcal, and switches.ELM Resource. There are 90 sequence variants in Hif3a between the two mouse lines. While most Fat line variants locate within intronic regions, Lean line variants are mainly in 3' UTR. We constructed a map of Hif3a potential regulatory regions and identified 39 regulatory variants by integrating data on constrained and regulatory elements, CpGs, and miRNAs and lncRNAs binding sites. Moreover, 3' UTR and two exonic variants may influence mRNA stability, translation rate and protein functionality. We propose as priority candidates for further functional studies a missense (rs37398126) and synonymous (rs37739792) variants, and intronic (rs47471302) variant that overlap conserved element in promoter region and predicted lncRNAs binding site. The results indicate a potential involvement of Hif3a in fat deposition. Additionally, approach used in the present study may serve as a general guideline for constructing an integrative gene map for prioritizing candidate gene variants with phenotypic effects.

Abstract 829: Genetic risk factors for histiocytic sarcoma in a canine disease model

Abstract Histiocytic sarcoma is a rare, aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart in dogs affects two breeds at unusually high frequency, with 20-25% of Flatcoated Retrievers and Bernese Mountain Dogs developing histiocytic sarcoma, while the disease is extremely rare in other breeds. In this study, we leveraged the canine model system to investigate genetic risk factors for histiocytic sarcoma in the Flatcoated Retriever. Population bottlenecks and artificial selection during breed formation within last 200 years have created the unique breed structure in dogs, characterized by long linkage disequilibrium and low within-breed genetic diversity. This allows for the identification of disease-associated traits with relatively few markers and individuals. We conducted multiple genome-wide association studies in Flatcoated Retrievers in which we identified two loci conferring risk for histiocytic sarcoma (177 cases and 132 controls, total). The first locus, on canine chromosome (CFA) 5, also colocalizes with previously identified risk for hemangiosarcoma and B-cell lymphoma, two hematopoietic malignancies in the closely related golden retriever breed. The second locus is unique to the Flatcoated Retriever and increases risk for histiocytic sarcoma in combination with CFA5. There have been multiple reports of human histiocytic sarcoma diagnosis concurrent with other hematological cancers, including lymphoma, suggested to be the result of transdifferentiation or a common neoplastic progenitor. Here we link germline predisposition for histiocytic sarcoma and B-cell lymphoma to the same locus in retrievers. Subsequent whole genome, transcriptome, and ChIP sequencing, highlighted candidate genes at each locus with roles in cell migration, as well as candidate regulatory variants with predicted impact on transcription factor binding. This comprehensive, multi-omics approach has been previously largely unexplored in the canine model. Here we advance our knowledge of the biology of histiocytic sarcoma, which will have direct impact on canine health and elucidate the genetic mechanisms underlying a rare human cancer. Citation Format: Jacquelyn M. Evans, Heidi G. Parker, Jocelyn Plassais, Gerard R. Rutteman, Elaine A. Ostrander. Genetic risk factors for histiocytic sarcoma in a canine disease model [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 829.

1137-P: Liver Regulatory Variants and Genes at Metabolic Trait Loci

Genome-wide association studies (GWAS) have identified many loci that affect liver- and lipid-related metabolic traits, but the causal variants, target genes and mechanisms for most loci have not been characterized. GWAS loci are often enriched in tissue-specific transcriptional regulatory elements that drive physiological changes, and loci can be annotated with colocalized molecular quantitative traits for gene expression (eQTL) and accessible chromatin (caQTL) to identify candidate variants and genes. We previously profiled chromatin accessibility in 20 human liver samples, identified caQTL variants that colocalized with GWAS loci for metabolic traits, and predicted target genes based on promoter proximity, Hi-C chromatin contact with a promoter, and co-localization with liver eQTL. Here, we examined selected GWAS loci and tested variants for effects on gene regulation. We prioritized loci based on trait, gene function, and QTL characteristics and selected four variant-gene pairs to study. To evaluate effects of specific variants on gene expression, we performed transcriptional reporter assays in HepG2 cells. Three caQTL variants exhibited significant (P<0.05) allelic differences in transcriptional activity. rs13395911, near EFHD1 and associated with liver enzyme gamma glutamyl transferase, showed allelic differences in enhancer activity (P=0.008). In addition, rs9556404, near GPR180 and associated with triglycerides, and rs34003091, near ZNF329 and associated with LDL-cholesterol, also showed significant allelic effects on transcriptional activity (P<0.05). To study EFHD1 function in liver metabolism, we used gene-specific shRNA constructs in HepG2 cells and observed up to 90% reduction of EFHD1 mRNA expression compared to scramble shRNA. Subsequent metabolic assays will test role of EFHD1 on hepatic lipid and glucose metabolism. These data identified candidate regulatory variants that may be responsible for nearby GWAS loci and suggest genes as therapeutic targets for fatty liver disease. Disclosure G. Pandey: None. K. Currin: None. S. Vadlamudi: None. K. Broadaway: None. K. L. Mohlke: None.

Multi-omics approach identifies germline regulatory variants associated with hematopoietic malignancies in retriever dog breeds.

Histiocytic sarcoma is an aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart affects flat-coated retrievers (FCRs) at unusually high frequency, with 20% developing the lethal disease. The similar clinical presentation combined with the closed population structure of dogs, leading to high genetic homogeneity, makes dogs an excellent model for genetic studies of cancer susceptibility. To determine the genetic risk factors underlying histiocytic sarcoma in FCRs, we conducted multiple genome-wide association studies (GWASs), identifying two loci that confer significant risk on canine chromosomes (CFA) 5 (Pwald = 4.83x10-9) and 19 (Pwald = 2.25x10-7). We subsequently undertook a multi-omics approach that has been largely unexplored in the canine model to interrogate these regions, generating whole genome, transcriptome, and chromatin immunoprecipitation sequencing. These data highlight the PI3K pathway gene PIK3R6 on CFA5, and proximal candidate regulatory variants that are strongly associated with histiocytic sarcoma and predicted to impact transcription factor binding. The CFA5 association colocalizes with susceptibility loci for two hematopoietic malignancies, hemangiosarcoma and B-cell lymphoma, in the closely related golden retriever breed, revealing the risk contribution this single locus makes to multiple hematological cancers. By comparison, the CFA19 locus is unique to the FCR and harbors risk alleles associated with upregulation of TNFAIP6, which itself affects cell migration and metastasis. Together, these loci explain ~35% of disease risk, an exceptionally high value that demonstrates the advantages of domestic dogs for complex trait mapping and genetic studies of cancer susceptibility.

Genetic Control of Rod Bipolar Cell Number in the Mouse Retina.

Genetic variants modulate the numbers of various retinal cell types in mice. For instance, there is minimal variation in the number of rod bipolar cells (RBCs) in two inbred strains of mice (A/J and C57BL/6J), yet their F1 offspring contain significantly more RBCs than either parental strain. To investigate the genetic source of this variation, we mapped the variation in the number of RBCs across 24 genetically distinct recombinant inbred (RI) strains (the AXB/BXA strain-set), seeking to identify quantitative trait loci (QTL). We then sought to identify candidate genes and potential casual variants at those genomic loci. Variation in RBC number mapped to three genomic loci, each modulating cell number in excess of one-third of the range observed across the RI strains. At each of these loci, we identified candidate genes containing variants that might alter gene function or expression. The latter genes were also analyzed using a transcriptome database, revealing a subset for which expression correlated with variation in RBC number. Using an electroporation strategy, we demonstrate that early postnatal expression of one of them, Ggct (gamma-glutamyl cyclotransferase), modulates bipolar cell number. We identify candidate regulatory variants for this gene, finding a large structural variant (SV) in the putative promoter that reduces expression using a luciferase assay. This SV reducing Ggct expression in vitro is likely the causal variant within the gene associated with the variation in Ggct expression in vivo, implicating it as a quantitative trait variant (QTV) participating in the control of RBC number.

Abstract 4290: A study in locus heterogeneity: Targeted sequencing analysis of 6q reveals multiple significant loci as the source of a previous linkage peak in familial lung cancer

Abstract Lung cancer is the leading cancer-related cause of death in the United States. We had previously performed multipoint linkage analysis on families with a strong family history of lung cancer and found significant linkage to the 6q25 region. In order to find the source of that linkage signal, we performed targeted sequencing 6q23-6q27 on 75 individuals from the 9 most highly linked families. We performed two types of parametric linkage analysis, using an autosomal dominant mode of inheritance with 10% penetrance for carriers and a 1% phenocopy rate using a disease allele frequency of 1%. The first type was a two-point analysis using an Elston-Stewart algorithm. This approach did not lead to any genome-wide significant results but demonstrated significant heterogeneity throughout the nine families.In an effort to recover the power from the original multipoint analysis, we performed a regional linkage analysis using SEQLinkage and MERLIN. SEQLinkage built regional haplotypes that corresponded to a gene or an intergenic region - the regions were based on a customized map of our design. The regional markers were multiallelic which allowed for greater information content and were similar to the microsatellites that were used in the original multipoint analysis. The regional markers were then analyzed in a two-point linkage analysis via MERLIN. This allowed us to identify two genome-wide significant signals at PACRG-AS1 at 6q26 (HLOD = 3.4) and SAMD5 at 6q24 (3.3). The PACRG-AS1 is novel, though it is associated with the known lung cancer gene PARK2. SAMD5 may be involved in lung cancer cell proliferation. The heterogeneity of the two signals was particularly interesting. The PACRG-AS1 was driven primarily by two families. The SAMD5 marker was not being strongly driven by any of the families and appears to be a small cumulative effect across the nine families. In addition, several families have large, but non-significant LOD scores at other loci across the region. This further reinforces the locus heterogeneity within the region, and it is likely that our two significant signals here are not the only variants affecting the phenotype. We further attempted to localize the signals by running SEQLinkage using a custom map where genes are broken into exons and introns. Though this resulted in no significant markers, the highest signal was located in the intronic region of SAMD5 (HLOD = 3) and several other suggestive signals were localized to intronic regions of good candidate genes like SASH1 and ARID1B. Examination of predicted effects of the candidate regulatory variants using eQTL databases is ongoing. This is a region full of promising candidates, and it is likely that the two significant signals found here are just part of many that could be affecting lung cancer risk. We plan to do further analysis within the individual families to elucidate any more genes affecting this signal. Citation Format: Anthony M. Musolf, Claire L. Simpson, Bilal A. Moiz, Candace Middlebrooks, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, Yafang Li, Ming You, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Joan E. Bailey-Wilson. A study in locus heterogeneity: Targeted sequencing analysis of 6q reveals multiple significant loci as the source of a previous linkage peak in familial lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4290. doi:10.1158/1538-7445.AM2017-4290

A Functional Variant Associated with Atrial Fibrillation Regulates PITX2c Expression through TFAP2a

The most significantly associated genetic locus for atrial fibrillation (AF) is in chromosomal region 4q25, where four independent association signals have been identified. Although model-system studies suggest that altered PITX2c expression might underlie the association, the link between specific variants and the direction of effect on gene expression remains unknown for all four signals. In the present study, we analyzed the AF-associated region most proximal to PITX2 at 4q25. First, we identified candidate regulatory variants that might confer AF risk through a combination of mammalian conservation, DNase hypersensitivity, and histone modification from ENCODE and the Roadmap Epigenomics Project, as well as through invivo analysis of enhancer activity in embryonic zebrafish. Within candidate regions, we then identified a single associated SNP, rs2595104, which displayed dramatically reduced enhancer activity with the AF risk allele. CRISPR-Cas9-mediated deletion of the rs2595104 region and editing of the rs2595104 risk allele in human stem-cell-derived cardiomyocytes resulted in diminished PITX2c expression in comparison to that of the non-risk allele. This differential activity was mediated by activating enhancer binding protein 2 alpha (TFAP2a), which bound robustly to the non-risk allele at rs2595104, but not to the risk allele, in cardiomyocytes. In sum, we found that the AF-associated SNP rs2595104 altered PITX2c expression via interaction with TFAP2a. Such a pathway could ultimately contribute to AF susceptibility at the PITX2 locus associated with AF.

Integrative functional genomics identifies regulatory mechanisms at coronary artery disease loci.

Coronary artery disease (CAD) is the leading cause of mortality and morbidity, driven by both genetic and environmental risk factors. Meta-analyses of genome-wide association studies have identified >150 loci associated with CAD and myocardial infarction susceptibility in humans. A majority of these variants reside in non-coding regions and are co-inherited with hundreds of candidate regulatory variants, presenting a challenge to elucidate their functions. Herein, we use integrative genomic, epigenomic and transcriptomic profiling of perturbed human coronary artery smooth muscle cells and tissues to begin to identify causal regulatory variation and mechanisms responsible for CAD associations. Using these genome-wide maps, we prioritize 64 candidate variants and perform allele-specific binding and expression analyses at seven top candidate loci: 9p21.3, SMAD3, PDGFD, IL6R, BMP1, CCDC97/TGFB1 and LMOD1. We validate our findings in expression quantitative trait loci cohorts, which together reveal new links between CAD associations and regulatory function in the appropriate disease context.

Construction of an integrative regulatory element and variation map of the murine Tst locus.

BackgroundGiven the abundance of new genomic projects and gene annotations, researchers trying to pinpoint causal genetic variants are faced with a challenging task of how to efficiently integrate all current genomic information. The objective of the study was to develop an approach to integrate various genomic annotations for a recently positionally-cloned Tst gene (Thiosulfate Sulfur Transferase, synonym Rhodanese) responsible for the Fob3b2 QTL effect on leanness and improved metabolic parameters. The second aim was to identify and prioritize Tst genetic variants that may be causal for the phenotypic effects.ResultsA bioinformatics approach was developed to integrate existing knowledge of regulatory elements of the Tst gene. The entire Tst locus along with flanking segments was sequenced between our unique polygenic mouse Fat and Lean strains that were generated by divergent selection on adiposity for over 60 generations. The bioinformatics-generated regulatory element map of the Tst locus was then combined with genetic variants between the Fat and Lean mice and with comparative analyses of polymorphisms across 17 mouse strains in order to prioritise likely causal polymorphisms. Two candidate regulatory variants were identified, one overlapping an evolutionary constrained Tst intronic element and the other residing in the seed region of a predicted 3′UTR miRNA binding site.ConclusionsThis study developed a map of regulatory elements for the Tst locus in mice and identified candidate genetic variants with increased causal likelihood. This map provides a basis for experimental validation and functional analyses of this novel candidate leanness and antidiabetic gene. Our methodological approach is of general utility for analyzing regulation of loci that have limited annotations and experimental evidence and for identifying candidate causal regulatory genetic variants in post-GWAS or post-QTL- cloning studies.Electronic supplementary materialThe online version of this article (doi:10.1186/s12863-016-0381-6) contains supplementary material, which is available to authorized users.

Allele-specific transcription factor binding to common and rare variants associated with disease and gene expression.

Genome-wide association studies (GWAS) have identified a large number of disease-associated SNPs, but in few cases the functional variant and the gene it controls have been identified. To systematically identify candidate regulatory variants, we sequenced ENCODE cell lines and used public ChIP-seq data to look for transcription factors binding preferentially to one allele. We found 9962 candidate regulatory SNPs, of which 16 % were rare and showed evidence of larger functional effect than common ones. Functionally rare variants may explain divergent GWAS results between populations and are candidates for a partial explanation of the missing heritability. The majority of allele-specific variants (96 %) were specific to a cell type. Furthermore, by examining GWAS loci we found >400 allele-specific candidate SNPs, 141 of which were highly relevant in our cell types. Functionally validated SNPs support identification of an SNP in SYNGR1 which may expose to the risk of rheumatoid arthritis and primary biliary cirrhosis, as well as an SNP in the last intron of COG6 exposing to the risk of psoriasis. We propose that by repeating the ChIP-seq experiments of 20 selected transcription factors in three to ten people, the most common polymorphisms can be interrogated for allele-specific binding. Our strategy may help to remove the current bottleneck in functional annotation of the genome.Electronic supplementary materialThe online version of this article (doi:10.1007/s00439-016-1654-x) contains supplementary material, which is available to authorized users.

Maps of open chromatin highlight cell type–restricted patterns of regulatory sequence variation at hematological trait loci

Nearly three-quarters of the 143 genetic signals associated with platelet and erythrocyte phenotypes identified by meta-analyses of genome-wide association (GWA) studies are located at non-protein-coding regions. Here, we assessed the role of candidate regulatory variants associated with cell type–restricted, closely related hematological quantitative traits in biologically relevant hematopoietic cell types. We used formaldehyde-assisted isolation of regulatory elements followed by next-generation sequencing (FAIRE-seq) to map regions of open chromatin in three primary human blood cells of the myeloid lineage. In the precursors of platelets and erythrocytes, as well as in monocytes, we found that open chromatin signatures reflect the corresponding hematopoietic lineages of the studied cell types and associate with the cell type–specific gene expression patterns. Dependent on their signal strength, open chromatin regions showed correlation with promoter and enhancer histone marks, distance to the transcription start site, and ontology classes of nearby genes. Cell type–restricted regions of open chromatin were enriched in sequence variants associated with hematological indices. The majority (63.6%) of such candidate functional variants at platelet quantitative trait loci (QTLs) coincided with binding sites of five transcription factors key in regulating megakaryopoiesis. We experimentally tested 13 candidate regulatory variants at 10 platelet QTLs and found that 10 (76.9%) affected protein binding, suggesting that this is a frequent mechanism by which regulatory variants influence quantitative trait levels. Our findings demonstrate that combining large-scale GWA data with open chromatin profiles of relevant cell types can be a powerful means of dissecting the genetic architecture of closely related quantitative traits.

Modelling Human Regulatory Variation in Mouse: Finding the Function in Genome-Wide Association Studies and Whole-Genome Sequencing

An increasing body of literature from genome-wide association studies and human whole-genome sequencing highlights the identification of large numbers of candidate regulatory variants of potential therapeutic interest in numerous diseases. Our relatively poor understanding of the functions of non-coding genomic sequence, and the slow and laborious process of experimental validation of the functional significance of human regulatory variants, limits our ability to fully benefit from this information in our efforts to comprehend human disease. Humanized mouse models (HuMMs), in which human genes are introduced into the mouse, suggest an approach to this problem. In the past, HuMMs have been used successfully to study human disease variants; e.g., the complex genetic condition arising from Down syndrome, common monogenic disorders such as Huntington disease and β-thalassemia, and cancer susceptibility genes such as BRCA1. In this commentary, we highlight a novel method for high-throughput single-copy site-specific generation of HuMMs entitled High-throughput Human Genes on the X Chromosome (HuGX). This method can be applied to most human genes for which a bacterial artificial chromosome (BAC) construct can be derived and a mouse-null allele exists. This strategy comprises (1) the use of recombineering technology to create a human variant–harbouring BAC, (2) knock-in of this BAC into the mouse genome using Hprt docking technology, and (3) allele comparison by interspecies complementation. We demonstrate the throughput of the HuGX method by generating a series of seven different alleles for the human NR2E1 gene at Hprt. In future challenges, we consider the current limitations of experimental approaches and call for a concerted effort by the genetics community, for both human and mouse, to solve the challenge of the functional analysis of human regulatory variation.

Cattle gain stature

Identifying causal variants for complex traits and understanding their function remain arduous tasks. A new study combines the advantages of gene mapping in livestock with elegant genetic and functional analyses to address these challenges and identifies candidate regulatory variants affecting stature in cattle.

Research Highlights

Research Highlights