Multiple genes in a single GWAS risk locus synergistically mediate aberrant synaptic development and function in human neurons

The mechanistic tie between genome-wide association study (GWAS)-implicated risk variants and disease-relevant cellular phenotypes remains largely unknown. Here, using human induced pluripotent stem cell (hiPSC)-derived neurons as a neurodevelopmental model, we identify multiple schizophrenia (SZ) risk variants that display allele-specific open chromatin (ASoC) and are likely to be functional. Editing the strongest ASoC SNP, rs2027349, near vacuolar protein sorting 45 homolog (VPS45) alters the expression of VPS45, lncRNA AC244033.2, and a distal gene, C1orf54. Notably, the transcriptomic changes in neurons are associated with SZ and other neuropsychiatric disorders. Neurons carrying the risk allele exhibit increased dendritic complexity and hyperactivity. Interestingly, individual/combinatorial gene knockdown shows that these genes alter cellular phenotypes in a non-additive synergistic manner. Our study reveals that multiple genes at a single GWAS risk locus mediate a compound effect on neural function, providing a mechanistic link between a non-coding risk variant and disease-related cellular phenotypes.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Blood cell traits' GWAS loci colocalization with variation in PU.1 genomic occupancy prioritizes causal noncoding regulatory variants

Genome-wide association studies (GWASs) have uncovered numerous trait-associated loci across the human genome, most of which are located in noncoding regions, making interpretation difficult. Moreover, causal variants are hard to statistically fine-map at many loci because of widespread linkage disequilibrium. To address this challenge, we present a strategy utilizing transcription factor (TF) binding quantitative trait loci (bQTLs) for colocalization analysis to identify trait associations likely mediated by TF occupancy variation and to pinpoint likely causal variants using motif scores. We applied this approach to PU.1 bQTLs in lymphoblastoid cell lines and blood cell trait GWAS data. Colocalization analysis revealed 69 blood cell trait GWAS loci putatively driven by PU.1 occupancy variation. We nominate PU.1 motif-altering variants as the likely shared causal variants at 51 loci. Such integration of TF bQTL data with other GWAS data may reveal transcriptional regulatory mechanisms and causal noncoding variants underlying additional complex traits.

The role of noncoding genetic variants in cardiomyopathy

Cardiomyopathies remain one of the leading causes of morbidity and mortality worldwide. Environmental risk factors and genetic predisposition account for most cardiomyopathy cases. As with all complex diseases, there are significant challenges in the interpretation of the molecular mechanisms underlying cardiomyopathy-associated genetic variants. Given the technical improvements and reduced costs of DNA sequence technologies, an increasing number of patients are now undergoing genetic testing, resulting in a continuously expanding list of novel mutations. However, many patients carry noncoding genetic variants, and although emerging evidence supports their contribution to cardiac disease, their role in cardiomyopathies remains largely understudied. In this review, we summarize published studies reporting on the association of different types of noncoding variants with various types of cardiomyopathies. We focus on variants within transcriptional enhancers, promoters, intronic sites, and untranslated regions that are likely associated with cardiac disease. Given the broad nature of this topic, we provide an overview of studies that are relatively recent and have sufficient evidence to support a significant degree of causality. We believe that more research with additional validation of noncoding genetic variants will provide further mechanistic insights on the development of cardiac disease, and noncoding variants will be increasingly incorporated in future genetic screening tests.

High-yield identification of pathogenic NF1 variants by skin fibroblast transcriptome screening after apparently normal diagnostic DNA testing

Neurofibromatosis type 1 (NF1) is caused by inactivating mutations in NF1. Due to the size, complexity, and high mutation rate at the NF1 locus, the identification of causative variants can be challenging. To obtain a molecular diagnosis in 15 individuals meeting diagnostic criteria for NF1, we performed transcriptome analysis (RNA-seq) on RNA obtained from cultured skin fibroblasts. In each case, routine molecular DNA diagnostics had failed to identify a disease-causing variant in NF1. A pathogenic variant or abnormal mRNA splicing was identified in 13 cases: 6 deep intronic variants and 2 transposon insertions causing noncanonical splicing, 3 postzygotic changes, 1 branch point mutation and, in 1 case, abnormal splicing for which the responsible DNA change remains to be identified. These findings helped resolve the molecular findings for an additional 17 individuals in multiple families with NF1, demonstrating the utility of skin-fibroblast-based transcriptome analysis for molecular diagnostics. RNA-seq improves mutation detection in NF1 and provides a powerful complementary approach to DNA-based methods. Importantly, our approach is applicable to other genetic disorders, particularly those caused by a wide variety of variants in a limited number of genes and specifically for individuals in whom routine molecular DNA diagnostics did not identify the causative variant.

Genome sequencing reveals novel noncoding variants in PLA2G6 and LMNB1 causing progressive neurologic disease

Neurodegenerative disorders and leukodystrophies are progressive neurologic conditions that can occur following the disruption of intricately coordinated patterns of gene expression. Exome sequencing has been adopted as an effective diagnostic tool for determining the underlying genetic etiology of Mendelian neurologic disorders, however genome sequencing offer advantages in its ability to identify and characterize copy number, structural, and sequence variants in noncoding regions. Genome sequencing from peripheral leukocytes was performed on two patients with progressive neurologic disease of unknown etiology following negative genetic investigations including exome sequencing. RNA sequencing from peripheral blood was performed to determine gene expression patterns in one of the patients. Potential causative variants were matched to the patients' clinical presentation. The first proband was found to be heterozygous for a likely pathogenic missense variant in PLA2G6 (c.386T>C; p.Leu129Pro) and have an additional deep intronic variant in PLA2G6 (c.2035-926G>A). RNA sequencing indicated this latter variant created a splice acceptor site leading to the incorporation of a pseudo-exon introducing a premature termination codon. The second proband was heterozygous for a 261 kb deletion upstream of LMNB1 that included an enhancer region. Previous reports of copy number variants spanning this region of cis-acting regulatory elements corroborated its pathogenicity. When combined with clinical presentations, these findings led to a definitive diagnosis of autosomal recessive infantile neuroaxonal dystrophy and autosomal dominant adult-onset demyelinating leukodystrophy, respectively. In patients with progressive neurologic disease of unknown etiology, genome sequencing with the addition of RNA analysis where appropriate should be considered for the identification of causative noncoding pathogenic variants.

A single-cell atlas of chromatin accessibility in the human genome

Current catalogs of regulatory sequences in the human genome are still incomplete and lack cell type resolution. To profile the activity of gene regulatory elements in diverse cell types and tissues in the human body, we applied single-cell chromatin accessibility assays to 30 adult human tissue types from multiple donors. We integrated these datasets with previous single-cell chromatin accessibility data from 15 fetal tissue types to reveal the status of open chromatin for ∼1.2 million candidate cis-regulatory elements (cCREs) in 222 distinct cell types comprised of >1.3 million nuclei. We used these chromatin accessibility maps to delineate cell-type-specificity of fetal and adult human cCREs and to systematically interpret the noncoding variants associated with complex human traits and diseases. This rich resource provides a foundation for the analysis of gene regulatory programs in human cell types across tissues, life stages, and organ systems.

WEVar: a novel statistical learning framework for predicting noncoding regulatory variants

Understanding the functional consequence of noncoding variants is of great interest. Though genome-wide association studies or quantitative trait locus analyses have identified variants associated with traits or molecular phenotypes, most of them are located in the noncoding regions, making the identification of causal variants a particular challenge. Existing computational approaches developed for prioritizing noncoding variants produce inconsistent and even conflicting results. To address these challenges, we propose a novel statistical learning framework, which directly integrates the precomputed functional scores from representative scoring methods. It will maximize the usage of integrated methods by automatically learning the relative contribution of each method and produce an ensemble score as the final prediction. The framework consists of two modes. The first 'context-free' mode is trained using curated causal regulatory variants from a wide range of context and is applicable to predict regulatory variants of unknown and diverse context. The second 'context-dependent' mode further improves the prediction when the training and testing variants are from the same context. By evaluating the framework via both simulation and empirical studies, we demonstrate that it outperforms integrated scoring methods and the ensemble score successfully prioritizes experimentally validated regulatory variants in multiple risk loci.

A method for scoring the cell type-specific impacts of noncoding variants in personal genomes

Here we use the expression and accessibility data from a diverse set of cell types to learn a model for the dependence of the accessibility of a regulatory element on its DNA sequence and TF expression. Using GTEx samples with WGS data, we show that the noncoding variants predicted to affect accessibility are more strongly associated with the expression of nearby genes. To interpret a personal genome, we combine the sequence information with context-specific TF expression to prioritize variants and regulatory elements in any genomic region of interest. This approach should be helpful in the study of risk loci previously identified by GWAS. Results from analysis of height and WGS data from the GTEx project support this hypothesis.

A person’s genome typically contains millions of variants which represent the differences between this personal genome and the reference human genome. The interpretation of these variants, i.e., the assessment of their potential impact on a person’s phenotype, is currently of great interest in human genetics and medicine. We have developed a prioritization tool called OpenCausal which takes as inputs 1) a personal genome and 2) a reference context-specific TF expression profile and returns a list of noncoding variants prioritized according to their impact on chromatin accessibility for any given genomic region of interest. We applied OpenCausal to 6,430 samples across 18 tissues derived from the GTEx project and found that the variants prioritized by OpenCausal are highly enriched for eQTLs and caQTLs. We further propose a strategy to integrate the predicted open scores with genome-wide association studies (GWAS) data to prioritize putative causal variants and regulatory elements for a given risk locus (i.e., fine-mapping analysis). As an initial example, we applied this method to a GWAS dataset of human height and found that the prioritized putative variants and elements are correlated with the phenotype (i.e., heights of individuals) better than others.

The Role of Noncoding Variants in Heritable Disease

The genetic basis of disease has largely focused on coding regions. However, it has become clear that a large proportion of the noncoding genome is functional and harbors genetic variants that contribute to disease etiology. Here, we review recent examples of inherited noncoding alterations that are responsible for Mendelian disorders or act to influence complex traits. We explore both rare and common genetic variants and discuss the wide range of mechanisms by which they affect gene regulation to promote disease. We also debate the challenges and progress associated with identifying and interpreting the functional and clinical significance of genetic variation in the context of the noncoding regulatory landscape.

Limitations of exome sequencing in detecting rare and undiagnosed diseases

While exome sequencing (ES) is commonly the final diagnostic step in clinical genetics, it may miss diagnoses. To clarify the limitations of ES, we investigated the diagnostic yield of genetic tests beyond ES in our Undiagnosed Diseases Network (UDN) participants. We reviewed the yield of additional genetic testing including genome sequencing (GS), copy number variant (CNV), noncoding variant (NCV), repeat expansion (RE), or methylation testing in UDN cases with nondiagnostic ES results. Overall, 36/54 (67%) of total diagnoses were based on clinical findings and coding variants found by ES and 3/54 (6%) were based on clinical findings only. The remaining 15/54 (28%) required testing beyond ES. Of these, 7/15 (47%) had NCV, 6/15 (40%) CNV, and 2/15 (13%) had a RE or a DNA methylation disorder. Thus 18/54 (33%) of diagnoses were not solved exclusively by ES. Several methods were needed to detect and/or confirm the functional effects of the variants missed by ES, and in some cases by GS. These results indicate that tests to detect elusive variants should be considered after nondiagnostic preliminary steps. Further studies are needed to determine the cost-effectiveness of tests beyond ES that provide diagnoses and insights to possible treatment.

Variants regulating ZBTB4 are associated with age-at-onset of Alzheimer's disease

The identification of novel genetic modifiers of age-at-onset (AAO) of Alzheimer's disease (AD) could advance our understanding of AD and provide novel therapeutic targets. A previous genome scan for modifiers of AAO among families affected by early-onset AD caused by the PSEN2 N141I variant identified 2 loci with significant evidence for linkage: 1q23.3 and 17p13.2. Here, we describe the fine-mapping of these 2 linkage regions, and test for replication in 6 independent datasets. By fine-mapping these linkage signals in a single large family, we reduced the linkage regions to 11% their original size and nominated 54 candidate variants. Among the 11 variants associated with AAO of AD in a larger sample of Germans from Russia, the strongest evidence implicated promoter variants influencing NCSTN on 1q23.3 and ZBTB4 on 17p13.2. The association between ZBTB4 and AAO of AD was replicated by multiple variants in independent, trans-ethnic datasets. Our results show association between AAO of AD and both ZBTB4 and NCSTN. ZBTB4 is a transcriptional repressor that regulates the cell cycle, including the apoptotic response to amyloid beta, while NCSTN is part of the gamma secretase complex, known to influence amyloid beta production. These genes therefore suggest important roles for amyloid beta and cell cycle pathways in AAO of AD.