Systematic localization of common disease-associated variation in regulatory DNA

TFTG: A comprehensive database for human transcription factors and their targets

Transcription factors (TFs) are major contributors to gene transcription, especially in controlling cell-specific gene expression and disease occurrence and development. Uncovering the relationship between TFs and their target genes is critical to understanding the mechanism of action of TFs. With the development of high-throughput sequencing techniques, a large amount of TF-related data has accumulated, which can be used to identify their target genes. In this study, we developed TFTG (Transcription Factor and Target Genes) database (http://tf.liclab.net/TFTG), which aimed to provide a large number of available human TF-target gene resources by multiple strategies, besides performing a comprehensive functional and epigenetic annotations and regulatory analyses of TFs. We identified extensive available TF-target genes by collecting and processing TF-associated ChIP-seq datasets, perturbation RNA-seq datasets and motifs. We also obtained experimentally confirmed relationships between TF and target genes from available resources. Overall, the target genes of TFs were obtained through integrating the relevant data of various TFs as well as fourteen identification strategies. Meanwhile, TFTG was embedded with user-friendly search, analysis, browsing, downloading and visualization functions. TFTG is designed to be a convenient resource for exploring human TF-target gene regulations, which will be useful for most users in the TF and gene expression regulation research.

Identification of candidate genes associated with milk production and mastitis based on transcriptome-wide association study

Genetic research for the assessment of mastitis and milk production traits simultaneously has a long history. The main issue that arises in this context is the known existence of a positive correlation between the risk of mastitis and lactation performance due to selection. The transcriptome-wide association study (TWAS) approach endeavors to combine the expression quantitative trait loci and genome-wide association study summary statistics to decode complex traits or diseases. Accordingly, we used the farmgtex project results as a complete bovine database for mastitis and milk production. The results of colocalization and TWAS approaches were used for the detection of functional associated candidate genes with milk production and mastitis traits on multiple tissue-based transcriptome records. Also, we used the david database for gene ontology to identify significant terms and associated genes. For the identification of interaction networks, the genemania and string databases were used. Also, the available z-scores in TWAS results were used for the calculation of the correlation between tissues. Therefore, the present results confirm that LYNX1, DGAT1, C14H8orf33, and LY6E were identified as significant genes associated with milk production in eight, six, five, and five tissues, respectively. Also, FBXL6 was detected as a significant gene associated with mastitis trait. CLN3 and ZNF34 genes emerged via both the colocalization and TWAS approaches as significant genes for milk production trait. It is expected that TWAS and colocalization can improve our perception of the potential health status control mechanism in high-yielding dairy cows.

Integrative analysis of genetics, epigenetics and RNA expression data reveal three susceptibility loci for smoking behavior in Chinese Han population

Despite numerous studies demonstrate that genetics and epigenetics factors play important roles on smoking behavior, our understanding of their functional relevance and coordinated regulation remains largely unknown. Here we present a multiomics study on smoking behavior for Chinese smoker population with the goal of not only identifying smoking-associated functional variants but also deciphering the pathogenesis and mechanism underlying smoking behavior in this under-studied ethnic population. After whole-genome sequencing analysis of 1329 Chinese Han male samples in discovery phase and OpenArray analysis of 3744 samples in replication phase, we discovered that three novel variants located near FOXP1 (rs7635815), and between DGCR6 and PRODH (rs796774020), and in ARVCF (rs148582811) were significantly associated with smoking behavior. Subsequently cis-mQTL and cis-eQTL analysis indicated that these variants correlated significantly with the differential methylation regions (DMRs) or differential expressed genes (DEGs) located in the regions where these variants present. Finally, our in silico multiomics analysis revealed several hub genes, like DRD2, PTPRD, FOXP1, COMT, CTNNAP2, to be synergistic regulated each other in the etiology of smoking.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Genetic regulation of cell type-specific chromatin accessibility shapes brain disease etiology

Nucleotide variants in cell type-specific gene regulatory elements in the human brain are risk factors for human disease. We measured chromatin accessibility in 1932 aliquots of sorted neurons and non-neurons from 616 human postmortem brains and identified 34,539 open chromatin regions with chromatin accessibility quantitative trait loci (caQTLs). Only 10.4% of caQTLs are shared between neurons and non-neurons, which supports cell type-specific genetic regulation of the brain regulome. Incorporating allele-specific chromatin accessibility improves statistical fine-mapping and refines molecular mechanisms that underlie disease risk. Using massively parallel reporter assays in induced excitatory neurons, we screened 19,893 brain QTLs and identified the functional impact of 476 regulatory variants. Combined, this comprehensive resource captures variation in the human brain regulome and provides insights into disease etiology.

Genetics of Hypertension: From Monogenic Analysis to GETomics

Arterial hypertension is the most frequent cardiovascular risk factor all over the world, and it is one of the leading drivers of the risk of cardiovascular events and death. It is a complex trait influenced by heritable and environmental factors. To date, the World Health Organization estimates that 1.28 billion adults aged 30-79 years worldwide have arterial hypertension (defined by European guidelines as office systolic blood pressure ≥ 140 mmHg or office diastolic blood pressure ≥ 90 mmHg), and 7.1 million die from this disease. The molecular genetic basis of primary arterial hypertension is the subject of intense research and has recently yielded remarkable progress. In this review, we will discuss the genetics of arterial hypertension. Recent studies have identified over 900 independent loci associated with blood pressure regulation across the genome. Comprehending these mechanisms not only could shed light on the pathogenesis of the disease but also hold the potential for assessing the risk of developing arterial hypertension in the future. In addition, these findings may pave the way for novel drug development and personalized therapeutic strategies.

Robust differential expression testing for single-cell CRISPR screens at low multiplicity of infection

Single-cell CRISPR screens (perturb-seq) link genetic perturbations to phenotypic changes in individual cells. The most fundamental task in perturb-seq analysis is to test for association between a perturbation and a count outcome, such as gene expression. We conduct the first-ever comprehensive benchmarking study of association testing methods for low multiplicity-of-infection (MOI) perturb-seq data, finding that existing methods produce excess false positives. We conduct an extensive empirical investigation of the data, identifying three core analysis challenges: sparsity, confounding, and model misspecification. Finally, we develop an association testing method - SCEPTRE low-MOI - that resolves these analysis challenges and demonstrates improved calibration and power.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Genetic regulation of m6A RNA methylation and its contribution in human complex diseases

N6-methyladenosine (m6A) has been established as the most prevalent chemical modification in message RNA (mRNA), playing an essential role in determining the fate of RNA molecules. Dysregulation of m6A has been revealed to lead to abnormal physiological conditions and cause various types of human diseases. Recent studies have delineated the genetic regulatory maps for m6A methylation by mapping the quantitative trait loci of m6A (m6A-QTLs), thereby building up the regulatory circuits linking genetic variants, m6A, and human complex traits. Here, we review the recent discoveries concerning the genetic regulatory maps of m6A, describing the methodological and technical details of m6A-QTL identification, and introducing the key findings of the cis- and trans-acting drivers of m6A. We further delve into the tissue- and ethnicity-specificity of m6A-QTL, the association with other molecular phenotypes in light of genetic regulation, the regulators underlying m6A genetics, and importantly, the functional roles of m6A in mediating human complex diseases. Lastly, we discuss potential research avenues that can accelerate the translation of m6A genetics studies toward the development of therapies for human genetic diseases.

Fundamentals for predicting transcriptional regulations from DNA sequence patterns

Cell-type-specific regulatory elements, cataloged through extensive experiments and bioinformatics in large-scale consortiums, have enabled enrichment analyses of genetic associations that primarily utilize positional information of the regulatory elements. These analyses have identified cell types and pathways genetically associated with human complex traits. However, our understanding of detailed allelic effects on these elements' activities and on-off states remains incomplete, hampering the interpretation of human genetic study results. This review introduces machine learning methods to learn sequence-dependent transcriptional regulation mechanisms from DNA sequences for predicting such allelic effects (not associations). We provide a concise history of machine-learning-based approaches, the requirements, and the key computational processes, focusing on primers in machine learning. Convolution and self-attention, pivotal in modern deep-learning models, are explained through geometrical interpretations using dot products. This facilitates understanding of the concept and why these have been used for machine learning for DNA sequences. These will inspire further research in this genetics and genomics field.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Advances in computational and experimental approaches for deciphering transcriptional regulatory networks: Understanding the roles of cis-regulatory elements is essential, and recent research utilizing MPRAs, STARR-seq, CRISPR-Cas9, and machine learning has yielded valuable insights

Understanding the influence of cis-regulatory elements on gene regulation poses numerous challenges given complexities stemming from variations in transcription factor (TF) binding, chromatin accessibility, structural constraints, and cell-type differences. This review discusses the role of gene regulatory networks in enhancing understanding of transcriptional regulation and covers construction methods ranging from expression-based approaches to supervised machine learning. Additionally, key experimental methods, including MPRAs and CRISPR-Cas9-based screening, which have significantly contributed to understanding TF binding preferences and cis-regulatory element functions, are explored. Lastly, the potential of machine learning and artificial intelligence to unravel cis-regulatory logic is analyzed. These computational advances have far-reaching implications for precision medicine, therapeutic target discovery, and the study of genetic variations in health and disease.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Regulatory landscape enrichment analysis (RLEA): a computational toolkit for non-coding variant enrichment and cell type prioritization

BACKGROUND

As genomic studies continue to implicate non-coding sequences in disease, testing the roles of these variants requires insights into the cell type(s) in which they are likely to be mediating their effects. Prior methods for associating non-coding variants with cell types have involved approaches using linkage disequilibrium or ontological associations, incurring significant processing requirements. GaiaAssociation is a freely available, open-source software that enables thousands of genomic loci implicated in a phenotype to be tested for enrichment at regulatory loci of multiple cell types in minutes, permitting insights into the cell type(s) mediating the studied phenotype.

RESULTS

In this work, we present Regulatory Landscape Enrichment Analysis (RLEA) by GaiaAssociation and demonstrate its capability to test the enrichment of 12,133 variants across the cis-regulatory regions of 44 cell types. This analysis was completed in 134.0 ± 2.3 s, highlighting the efficient processing provided by GaiaAssociation. The intuitive interface requires only four inputs, offers a collection of customizable functions, and visualizes variant enrichment in cell-type regulatory regions through a heatmap matrix. GaiaAssociation is available on PyPi for download as a command line tool or Python package and the source code can also be installed from GitHub at https://github.com/GreallyLab/gaiaAssociation .

CONCLUSIONS

GaiaAssociation is a novel package that provides an intuitive and efficient resource to understand the enrichment of non-coding variants across the cis-regulatory regions of different cells, empowering studies seeking to identify disease-mediating cell types.

A functional schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome-wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes, yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated four isogenic clones homozygous for the risk allele and five clones homozygous for the non-risk allele into neural progenitor cells. Employing RNA sequencing, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at a false discovery rate (FDR) of <0.05 and 259 genes at a FDR of <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

Chromatin accessibility during human first-trimester neurodevelopment

The human brain develops through a tightly organized cascade of patterning events, induced by transcription factor expression and changes in chromatin accessibility. Although gene expression across the developing brain has been described at single-cell resolution1, similar atlases of chromatin accessibility have been primarily focused on the forebrain2-4. Here we describe chromatin accessibility and paired gene expression across the entire developing human brain during the first trimester (6-13 weeks after conception). We defined 135 clusters and used multiomic measurements to link candidate cis-regulatory elements to gene expression. The number of accessible regions increased both with age and along neuronal differentiation. Using a convolutional neural network, we identified putative functional transcription factor-binding sites in enhancers characterizing neuronal subtypes. We applied this model to cis-regulatory elements linked to ESRRB to elucidate its activation mechanism in the Purkinje cell lineage. Finally, by linking disease-associated single nucleotide polymorphisms to cis-regulatory elements, we validated putative pathogenic mechanisms in several diseases and identified midbrain-derived GABAergic neurons as being the most vulnerable to major depressive disorder-related mutations. Our findings provide a more detailed view of key gene regulatory mechanisms underlying the emergence of brain cell types during the first trimester and a comprehensive reference for future studies related to human neurodevelopment.

Genetic mapping across autoimmune diseases reveals shared associations and mechanisms

Autoimmune and inflammatory diseases are polygenic disorders of the immune system. Many genomic loci harbor risk alleles for several diseases, but the limited resolution of genetic mapping prevents determining whether the same allele is responsible, indicating a shared underlying mechanism. Here, using a collection of 129,058 cases and controls across 6 diseases, we show that ~40% of overlapping associations are due to the same allele. We improve fine-mapping resolution for shared alleles twofold by combining cases and controls across diseases, allowing us to identify more expression quantitative trait loci driven by the shared alleles. The patterns indicate widespread sharing of pathogenic mechanisms but not a single global autoimmune mechanism. Our approach can be applied to any set of traits and is particularly valuable as sample collections become depleted.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Identification of bone mineral density associated genes with shared genetic architectures across multiple tissues: Functional insights for EPDR1, PKDCC, and SPTBN1

Recent studies suggest a shared genetic architecture between muscle and bone, yet the underlying molecular mechanisms remain elusive. This study aims to identify the functionally annotated genes with shared genetic architecture between muscle and bone using the most up-to-date genome-wide association study (GWAS) summary statistics from bone mineral density (BMD) and fracture-related genetic variants. We employed an advanced statistical functional mapping method to investigate shared genetic architecture between muscle and bone, focusing on genes highly expressed in muscle tissue. Our analysis identified three genes, EPDR1, PKDCC, and SPTBN1, which are highly expressed in muscle tissue and previously unlinked to bone metabolism. About 90% and 85% of filtered Single-Nucleotide Polymorphisms were in the intronic and intergenic regions for the threshold at P≤5×10-8 and P≤5×10-100, respectively. EPDR1 was highly expressed in multiple tissues, including muscles, adrenal glands, blood vessels, and the thyroid. SPTBN1 was highly expressed in all 30 tissue types except blood, while PKDCC was highly expressed in all 30 tissue types except the brain, pancreas, and skin. Our study provides a framework for using GWAS findings to highlight functional evidence of crosstalk between multiple tissues based on shared genetic architecture between muscle and bone. Further research should focus on functional validation, multi-omics data integration, gene-environment interactions, and clinical relevance in musculoskeletal disorders.

A systems biology-based identification and in vivo functional screening of Alzheimer's disease risk genes reveal modulators of memory function

Genome-wide association studies (GWASs) have uncovered over 75 genomic loci associated with risk for late-onset Alzheimer's disease (LOAD), but identification of the underlying causal genes remains challenging. Studies of induced pluripotent stem cell (iPSC)-derived neurons from LOAD patients have demonstrated the existence of neuronal cell-intrinsic functional defects. Here, we searched for genetic contributions to neuronal dysfunction in LOAD using an integrative systems approach that incorporated multi-evidence-based gene mapping and network-analysis-based prioritization. A systematic perturbation screening of candidate risk genes in Caenorhabditis elegans (C. elegans) revealed that neuronal knockdown of the LOAD risk gene orthologs vha-10 (ATP6V1G2), cmd-1 (CALM3), amph-1 (BIN1), ephx-1 (NGEF), and pho-5 (ACP2) alters short-/intermediate-term memory function, the cognitive domain affected earliest during LOAD progression. These results highlight the impact of LOAD risk genes on evolutionarily conserved memory function, as mediated through neuronal endosomal dysfunction, and identify new targets for further mechanistic interrogation.

Correlation analysis of serum TLR4 protein levels and TLR4 gene polymorphisms in gouty arthritis patients

OBJECTIVE

The Toll-like receptor (TLR) 4-mediated nuclear factor kappa B (NF-κB) signaling pathway regulates the production of inflammatory factors and plays a key role in the pathogenesis of gouty arthritis. The aim of the present study was to investigate the link among TLR4 gene polymorphisms at various loci, protein expression, and gouty arthritis susceptibility.

METHODS

Between 2016 and 2021, a case-control study was used to collect a total of 1207 study subjects, including 317 male patients with gouty arthritis (gout group) and 890 healthy males (control group). The association between gout susceptibility and different genetic models was analyzed by typing three loci of the TLR4 gene (rs2149356, rs2737191, and rs10759932) using a multiplex point mutation rapid assay, and the association between protein expression and gout was confirmed by measuring TLR4 protein concentrations using enzyme-linked immunosorbent assays (ELISAs).

RESULTS

In a codominant models AA and AG, the rs2737191 polymorphism in the gout group increased the risk of gout compared to the AA genotype (OR = 2.249, 95%CI 1.010~5.008), and the risk of gout was higher for those carrying the G allele compared to the A allele (OR = 2.227, 95%CI 1.006~4.932). TLR4 protein expression was different between the two groups with different locus genotypes. The differences in TLR4 protein expression between the gout group and control group were statistically significant between the following genotypes: the GG and GT genotypes of the rs2149356 polymorphism; the AA and AG genotypes of the rs2737191 polymorphism; and the TT and TC genotypes of the rs10759932 polymorphism(P<0.05). The TLR4 protein level in the gout group (19.19±3.09 ng/ml) was significantly higher than that in the control group (15.85±4.75 ng/ml).

CONCLUSION

The AG genotype of the TLR4 gene rs2737191 polymorphism may be correlated with the development of gouty arthritis. The level of TLR4 protein expression is significantly higher in patients with gouty arthritis than in controls, and there is a correlation between high TLR4 protein expression and the development of gouty arthritis.

Systematic druggable genome-wide Mendelian randomization identifies therapeutic targets for sarcopenia

BACKGROUND

There are no effective pharmacological treatments for sarcopenia. We aim to identify potential therapeutic targets for sarcopenia by integrating various publicly available datasets.

METHODS

We integrated druggable genome data, cis-eQTL/cis-pQTL from human blood and skeletal muscle tissue, and GWAS summary data of sarcopenia-related traits to analyse the potential causal relationships between drug target genes and sarcopenia using the Mendelian Randomization (MR) method. Sensitivity analyses and Bayesian colocalization were employed to validate the causal relationships. We also assessed the side effects or additional indications of the identified drug targets using a phenome-wide MR (Phe-MR) approach and investigated actionable drugs for target genes using available databases.

RESULTS

MR analysis identified 17 druggable genes with potential causation to sarcopenia in human blood or skeletal muscle tissue. Six of them (HP, HLA-DRA, MAP 3K3, MFGE8, COL15A1, and AURKA) were further confirmed by Bayesian colocalization (PPH4 > 90%). The up-regulation of HP [higher ALM (beta: 0.012, 95% CI: 0.007-0.018, P = 1.2*10-5) and higher grip strength (OR: 0.96, 95% CI: 0.94-0.98, P = 4.2*10-5)], MAP 3K3 [higher ALM (beta: 0.24, 95% CI: 0.21-0.26, P = 1.8*10-94), higher grip strength (OR: 0.82, 95% CI: 0.75-0.90, P = 2.1*10-5), and faster walking pace (beta: 0.03, 95% CI: 0.02-0.05, P = 8.5*10-6)], and MFGE8 [higher ALM (muscle eQTL, beta: 0.09, 95% CI: 0.06-0.11, P = 6.1*10-13; blood pQTL, beta: 0.05, 95% CI: 0.03-0.07, P = 3.8*10-09)], as well as the down-regulation of HLA-DRA [lower ALM (beta: -0.09, 95% CI: -0.11 to -0.08, P = 5.4*10-36) and lower grip strength (OR: 1.13, 95% CI: 1.07-1.20, P = 1.8*10-5)] and COL15A1 [higher ALM (muscle eQTL, beta: -0.07, 95% CI: -0.10 to -0.04, P = 3.4*10-07; blood pQTL, beta: -0.05, 95% CI: -0.06 to -0.03, P = 1.6*10-07)], decreased the risk of sarcopenia. AURKA in blood (beta: -0.16, 95% CI: -0.22 to -0.09, P = 2.1*10-06) and skeletal muscle (beta: 0.03, 95% CI: 0.02 to 0.05, P = 5.3*10-05) tissues showed an inverse relationship with sarcopenia risk. The Phe-MR indicated that the six potential therapeutic targets for sarcopenia had no significant adverse effects. Drug repurposing analysis supported zinc supplementation and collagenase clostridium histolyticum might be potential therapeutics for sarcopenia by activating HP and inhibiting COL15A1, respectively.

CONCLUSIONS

Our research indicated MAP 3K3, MFGE8, COL15A1, HP, and HLA-DRA may serve as promising targets for sarcopenia, while the effectiveness of zinc supplementation and collagenase clostridium histolyticum for sarcopenia requires further validation.

Robust differential expression testing for single-cell CRISPR screens at low multiplicity of infection

Single-cell CRISPR screens (perturb-seq) link genetic perturbations to phenotypic changes in individual cells. The most fundamental task in perturb-seq analysis is to test for association between a perturbation and a count outcome, such as gene expression. We conduct the first-ever comprehensive benchmarking study of association testing methods for low multiplicity-of-infection (MOI) perturb-seq data, finding that existing methods produce excess false positives. We conduct an extensive empirical investigation of the data, identifying three core analysis challenges: sparsity, confounding, and model misspecification. Finally, we develop an association testing method - SCEPTRE low-MOI - that resolves these analysis challenges and demonstrates improved calibration and power.

Deletion of a conserved genomic region associated with adolescent idiopathic scoliosis leads to vertebral rotation in mice

Adolescent idiopathic scoliosis (AIS) is the most common form of scoliosis, in which spinal curvature develops in adolescence, and 90% of patients are female. Scoliosis is a debilitating disease that often requires bracing or surgery in severe cases. AIS affects 2%-5.2% of the population; however, the biological origin of the disease remains poorly understood. In this study, we aimed to determine the function of a highly conserved genomic region previously linked to AIS using a mouse model generated by CRISPR-CAS9 gene editing to knockout this area of the genome to understand better its contribution to AIS, which we named AIS_CRMΔ. We also investigated the upstream factors that regulate the activity of this enhancer in vivo, whether the spatial expression of the LBX1 protein would change with the loss of AIS-CRM function, and whether any phenotype would arise after deletion of this region. We found a significant increase in mRNA expression in the developing neural tube at E10.5, and E12.5, for not only Lbx1 but also other neighboring genes. Adult knockout mice showed vertebral rotation and proprioceptive deficits, also observed in human AIS patients. In conclusion, our study sheds light on the elusive biological origins of AIS, by targeting and investigating a highly conserved genomic region linked to AIS in humans. These findings provide valuable insights into the function of the investigated region and contribute to our understanding of the underlying causes of this debilitating disease.

OCA-B promotes autoimmune demyelination through control of stem-like CD4 + T cells

Stem-like T cell populations can selectively contribute to autoimmunity, but the activities that promote and sustain these populations are incompletely understood. Here, we show that T cell-intrinsic loss of the transcription cofactor OCA-B protects mice from experimental autoimmune encephalomyelitis (EAE) while preserving responses to CNS infection. In adoptive transfer EAE models driven by multiple antigen encounters, OCA-B deletion nearly eliminates CNS infiltration, proinflammatory cytokine production and clinical disease. OCA-B-expressing CD4 + T cells within the CNS of mice with EAE comprise a minority of the population but display a memory phenotype and preferentially confer disease. In a relapsing-remitting EAE model, OCA-B T cell deficiency specifically protects mice from relapse. During remission, OCA-B promotes the expression of Tcf7 , Slamf6 , and Sell in proliferating T cell populations. At relapse, OCA-B loss results in both the accumulation of an immunomodulatory CD4 + T cell population expressing Ccr9 and Bach2 , and the loss of pro-inflammatory gene expression from Th17 cells. These results identify OCA-B as a driver of pathogenic stem-like T cells.

LINC00998 Modulating M2 Macrophage Activation in Allergic Rhinitis by Stabilizing BOB.1 mRNA

Background

Allergic rhinitis (AR) is globally recognized as a considerable threat to human health with a rising prevalence and a substantial medical and socioeconomic burden. Numerous studies have emphasized the significance of long noncoding RNAs (lncRNAs) in allergic responses. Hence, this research dealt with exploring the involvement of the lncRNA LINC00998 in the mechanism of AR.

Methods

LINC00998 expression was assessed by qRT-PCR in peripheral blood mononuclear cells acquired from individuals with AR. Additionally, the potential relationship between LINC00998 and macrophage polarization was observed in vitro. Then we constructed AR mice model and macrophage polarization models using THP-1 cells as well as primary human macrophages to verify the M2 shift in AR and the low expression level of LINC00998 in M2 macrophages. We used gain- and loss-of-function experiments to explore the modification of LINC00998 in macrophage polarization. Furthermore, we explored the underlying mechanism of LINC00998 mediates through qRT-PCR, flow cytometry, and Western blot.

Results

The analysis revealed a significant decrease in LINC00998 expression in the samples obtained from patients with AR. LINC00998 is markedly increased in M1 macrophages whereas decreased in M2 macrophages in vitro. Furthermore, suppression of LINC00998 caused a remarkable enhancement in M2 polarization, whereas its overexpression led to its attenuation. Knockdown of LINC00998 led to a remarkable downregulation of BOB.1 mRNA and protein, while overexpression of LINC00998 upregulated their expression. Moreover, it was found that BOB.1 modulated macrophage polarization through the PU.1/IL-1β axis. Meanwhile, the modulation of LINC00098 overexpression on macrophage polarization and PU.1/ IL-1β can be reversed by BOB.1 siRNA.

Conclusion

This research revealed the lncRNA LINC00998 altered M2 macrophage polarization by regulating the BOB.1/PU.1/IL-1β axis, which open up new avenues for studying the pathogenesis of AR.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

CWAS-Plus: Estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human brain pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS) employs both whole-genome sequencing and user-provided functional data to enhance noncoding variant analysis, with a faster and more efficient execution of the CWAS workflow. Here, we used single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type specific enhancers and promoters. Examining autism spectrum disorder whole-genome sequencing data (n = 7,280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease whole-genome sequencing data (n = 1,087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale whole-genome sequencing data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/ , respectively.

KEY POINTS

CWAS-Plus efficiently identifies noncoding associations in WGS data, supporting user-friendly categorization and burden enrichment tests.CWAS-Plus integrates various functional datasets, emphasizing cell-type-specific noncoding associations.CWAS-Plus provides a novel approach for multiple testing correction, enhancing the reliability of the results.Autism spectrum disorder risk noncoding variants are identified as enriched with transcription factors, suggesting their role in the pathology.Rare variant analysis with Alzheimer's disease samples reveals strong association with microglia, supporting the reliability of the results produced by CWAS-Plus.

Chromatin accessibility variation provides insights into missing regulation underlying immune-mediated diseases

Most genetic loci associated with complex traits and diseases through genome-wide association studies (GWAS) are noncoding, suggesting that the causal variants likely have gene regulatory effects. However, only a small number of loci have been linked to expression quantitative trait loci (eQTLs) detected currently. To better understand the potential reasons for many trait-associated loci lacking eQTL colocalization, we investigated whether chromatin accessibility QTLs (caQTLs) in lymphoblastoid cell lines (LCLs) explain immune-mediated disease associations that eQTLs in LCLs did not. The power to detect caQTLs was greater than that of eQTLs and was less affected by the distance from the transcription start site of the associated gene. Meta-analyzing LCL eQTL data to increase the sample size to over a thousand led to additional loci with eQTL colocalization, demonstrating that insufficient statistical power is still likely to be a factor. Moreover, further eQTL colocalization loci were uncovered by surveying eQTLs of other immune cell types. Altogether, insufficient power and context-specificity of eQTLs both contribute to the 'missing regulation.'

Liver regulatory mechanisms of noncoding variants at lipid and metabolic trait loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for liver disease and lipid-related metabolic traits, although identifying their target genes and molecular mechanisms remains challenging. We predicted target genes at GWAS signals by integrating them with molecular quantitative trait loci for liver gene expression (eQTL) and liver chromatin accessibility QTL (caQTL). We predicted specific regulatory caQTL variants at four GWAS signals located near EFHD1, LITAF, ZNF329, and GPR180. Using transcriptional reporter assays, we determined that caQTL variants rs13395911, rs11644920, rs34003091, and rs9556404 exhibit allelic differences in regulatory activity. We also performed a protein binding assay for rs13395911 and found that FOXA2 differentially interacts with the alleles of rs13395911. For variants rs13395911 and rs11644920 in putative enhancer regulatory elements, we used CRISPRi to demonstrate that repression of the enhancers altered the expression of the predicted target and/or nearby genes. Repression of the element at rs13395911 reduced the expression of EFHD1, and repression of the element at rs11644920 reduced the expression of LITAF, SNN, and TXNDC11. Finally, we showed that EFHD1 is a metabolically active gene in HepG2 cells. Together, these results provide key steps to connect genetic variants with cellular mechanisms and help elucidate the causes of liver disease.

Improved estimation of functional enrichment in SNP heritability using feasible generalized least squares

Functional enrichment results typically implicate tissue or cell-type-specific biological pathways in disease pathogenesis and as therapeutic targets. We propose generalized linkage disequilibrium score regression (g-LDSC) that requires only genome-wide association studies (GWASs) summary-level data to estimate functional enrichment. The method adopts the same assumptions and regression model formulation as stratified linkage disequilibrium score regression (s-LDSC). Although s-LDSC only partially uses LD information, our method uses the whole LD matrix, which accounts for possible correlated error structure via a feasible generalized least-squares estimation. We demonstrate through simulation studies under various scenarios that g-LDSC provides more precise estimates of functional enrichment than s-LDSC, regardless of model misspecification. In an application to GWAS summary statistics of 15 traits from the UK Biobank, estimates of functional enrichment using g-LDSC were lower and more realistic than those obtained from s-LDSC. In addition, g-LDSC detected more significantly enriched functional annotations among 24 functional annotations for the 15 traits than s-LDSC (118 vs. 51).

Subset-based method for cross-tissue transcriptome-wide association studies improves power and interpretability

Integrating results from genome-wide association studies (GWASs) and studies of molecular phenotypes such as gene expressions can improve our understanding of the biological functions of trait-associated variants and can help prioritize candidate genes for downstream analysis. Using reference expression quantitative trait locus (eQTL) studies, several methods have been proposed to identify gene-trait associations, primarily based on gene expression imputation. To increase the statistical power by leveraging substantial eQTL sharing across tissues, meta-analysis methods aggregating such gene-based test results across multiple tissues or contexts have been developed as well. However, most existing meta-analysis methods have limited power to identify associations when the gene has weaker associations in only a few tissues and cannot identify the subset of tissues in which the gene is "activated." For this, we developed a cross-tissue subset-based transcriptome-wide association study (CSTWAS) meta-analysis method that improves power under such scenarios and can extract the set of potentially associated tissues. To improve applicability, CSTWAS uses only GWAS summary statistics and pre-computed correlation matrices to identify a subset of tissues that have the maximal evidence of gene-trait association. Through numerical simulations, we found that CSTWAS can maintain a well-calibrated type-I error rate, improves power especially when there is a small number of associated tissues for a gene-trait association, and identifies an accurate associated tissue set. By analyzing GWAS summary statistics of three complex traits and diseases, we demonstrate that CSTWAS could identify biological meaningful signals while providing an interpretation of disease etiology by extracting a set of potentially associated tissues.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

Trans-eQTL mapping in gene sets identifies network effects of genetic variants

Nearly all trait-associated variants identified in genome-wide association studies (GWASs) are noncoding. The cis regulatory effects of these variants have been extensively characterized, but how they affect gene regulation in trans has been the subject of fewer studies because of the difficulty in detecting trans-expression quantitative loci (eQTLs). We developed trans-PCO for detecting trans effects of genetic variants on gene networks. Our simulations demonstrate that trans-PCO substantially outperforms existing trans-eQTL mapping methods. We applied trans-PCO to two gene expression datasets from whole blood, DGN (N = 913) and eQTLGen (N = 31,684), and identified 14,985 high-quality trans-eSNP-module pairs associated with 197 co-expression gene modules and biological processes. We performed colocalization analyses between GWAS loci of 46 complex traits and the trans-eQTLs. We demonstrated that the identified trans effects can help us understand how trait-associated variants affect gene regulatory networks and biological pathways.

Genetic determinants of plasma protein levels in the Estonian population

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

Alteration of DNA Methylation and Epigenetic Scores Associated With Features of Schizophrenia and Common Variant Genetic Risk

BACKGROUND

Unpacking molecular perturbations associated with features of schizophrenia is a critical step toward understanding phenotypic heterogeneity in this disorder. Recent epigenome-wide association studies have uncovered pervasive dysregulation of DNA methylation in schizophrenia; however, clinical features of the disorder that account for a large proportion of phenotypic variability are relatively underexplored.

METHODS

We comprehensively analyzed patterns of DNA methylation in a cohort of 381 individuals with schizophrenia from the deeply phenotyped Australian Schizophrenia Research Bank. Epigenetic changes were investigated in association with cognitive status, age of onset, treatment resistance, Global Assessment of Functioning scores, and common variant polygenic risk scores for schizophrenia. We subsequently explored alterations within genes previously associated with psychiatric illness, phenome-wide epigenetic covariance, and epigenetic scores.

RESULTS

Epigenome-wide association studies of the 5 primary traits identified 662 suggestively significant (p < 6.72 × 10-5) differentially methylated probes, with a further 432 revealed after controlling for schizophrenia polygenic risk on the remaining 4 traits. Interestingly, we uncovered many probes within genes associated with a variety of psychiatric conditions as well as significant epigenetic covariance with phenotypes and exposures including acute myocardial infarction, C-reactive protein, and lung cancer. Epigenetic scores for treatment-resistant schizophrenia strikingly exhibited association with clozapine administration, while epigenetic proxies of plasma protein expression, such as CCL17, MMP10, and PRG2, were associated with several features of schizophrenia.

CONCLUSIONS

Our findings collectively provide novel evidence suggesting that several features of schizophrenia are associated with alteration of DNA methylation, which may contribute to interindividual phenotypic variation in affected individuals.

Multicenter integrated analysis of noncoding CRISPRi screens

The ENCODE Consortium's efforts to annotate noncoding cis-regulatory elements (CREs) have advanced our understanding of gene regulatory landscapes. Pooled, noncoding CRISPR screens offer a systematic approach to investigate cis-regulatory mechanisms. The ENCODE4 Functional Characterization Centers conducted 108 screens in human cell lines, comprising >540,000 perturbations across 24.85 megabases of the genome. Using 332 functionally confirmed CRE-gene links in K562 cells, we established guidelines for screening endogenous noncoding elements with CRISPR interference (CRISPRi), including accurate detection of CREs that exhibit variable, often low, transcriptional effects. Benchmarking five screen analysis tools, we find that CASA produces the most conservative CRE calls and is robust to artifacts of low-specificity single guide RNAs. We uncover a subtle DNA strand bias for CRISPRi in transcribed regions with implications for screen design and analysis. Together, we provide an accessible data resource, predesigned single guide RNAs for targeting 3,275,697 ENCODE SCREEN candidate CREs with CRISPRi and screening guidelines to accelerate functional characterization of the noncoding genome.

Comprehensive functional interrogation of susceptibility loci in GWASs identified KIAA0391 as a novel oncogenic driver via regulating pyroptosis in NSCLC

Approximately 51 non-small-cell lung cancer (NSCLC) risk loci have been identified by genome-wide association studies (GWASs). We conducted a high throughput RNA-interference (RNAi) screening to identify the candidate causal genes in NSCLC risk loci. KIAA0391 at 14q13.1 had the highest score and could promote proliferation and metastasis of NSCLC in vitro and in vivo. We next prioritized rs3783313 as a causal variant at 14q13.1, by integrating a large-scale population study consisting of 27,120 lung cancer cases and 27,355 controls, functional annotation, and expression quantitative trait locus (eQTL) analysis. Then we found that rs3783313 could facilitate a promoter-enhancer interaction to upregulate KIAA0391 expression by affecting the affinity of transcription factor NFYA. Mechanistically, KIAA0391 knockdown dramatically influenced pyroptosis-related pathways and increased the expression of CASP1. And KIAA0391 transcriptionally repressed CASP1 by binding to SMAD2 and induced an anti-pyroptosis phenotype, promoting tumorigenesis of NSCLC, which provides new insights and potential target for NSCLC.

Genetic control of DNA methylation is largely shared across European and East Asian populations

DNA methylation is an ideal trait to study the extent of the shared genetic control across ancestries, effectively providing hundreds of thousands of model molecular traits with large QTL effect sizes. We investigate cis DNAm QTLs in three European (n = 3701) and two East Asian (n = 2099) cohorts to quantify the similarities and differences in the genetic architecture across populations. We observe 80,394 associated mQTLs (62.2% of DNAm probes with significant mQTL) to be significant in both ancestries, while 28,925 mQTLs (22.4%) are identified in only a single ancestry. mQTL effect sizes are highly conserved across populations, with differences in mQTL discovery likely due to differences in allele frequency of associated variants and differing linkage disequilibrium between causal variants and assayed SNPs. This study highlights the overall similarity of genetic control across ancestries and the value of ancestral diversity in increasing the power to detect associations and enhancing fine mapping resolution.

Identification of candidate causal cis -regulatory variants underlying electrocardiographic QT interval GWAS loci

Identifying causal variants among tens or hundreds of associated variants at each locus mapped by genome-wide association studies (GWAS) of complex traits is a challenge. As vast majority of GWAS variants are noncoding, sequence variation at cis -regulatory elements affecting transcriptional expression of specific genes is a widely accepted molecular hypothesis. Following this cis -regulatory hypothesis and combining it with the observation that nucleosome-free open chromatin is a universal hallmark of all types of cis -regulatory elements, we aimed to identify candidate causal regulatory variants underlying electrocardiographic QT interval GWAS loci. At a dozen loci, selected for higher effect sizes and a better understanding of the likely causal gene, we identified and included all common variants in high linkage disequilibrium with the GWAS variants as candidate variants. Using ENCODE DNase-seq and ATAC-seq from multiple human adult cardiac left ventricle tissue samples, we generated genome-wide maps of open chromatin regions marking putative regulatory elements. QT interval associated candidate variants were filtered for overlap with cardiac left ventricle open chromatin regions to identify candidate causal cis -regulatory variants, which were further assessed for colocalizing with a known cardiac GTEx expression quantitative trait locus variant as additional evidence for their causal role. Together, these efforts have generated a comprehensive set of candidate causal variants that are expected to be enriched for cis -regulatory potential and thereby, explaining the observed genetic associations.

Population-enriched innate immune variants may identify candidate gene targets at the intersection of cancer and cardio-metabolic disease

Both cancer and cardio-metabolic disease disparities exist among specific populations in the US. For example, African Americans experience the highest rates of breast and prostate cancer mortality and the highest incidence of obesity. Native and Hispanic Americans experience the highest rates of liver cancer mortality. At the same time, Pacific Islanders have the highest death rate attributed to type 2 diabetes (T2D), and Asian Americans experience the highest incidence of non-alcoholic fatty liver disease (NAFLD) and cancers induced by infectious agents. Notably, the pathologic progression of both cancer and cardio-metabolic diseases involves innate immunity and mechanisms of inflammation. Innate immunity in individuals is established through genetic inheritance and external stimuli to respond to environmental threats and stresses such as pathogen exposure. Further, individual genomes contain characteristic genetic markers associated with one or more geographic ancestries (ethnic groups), including protective innate immune genetic programming optimized for survival in their corresponding ancestral environment(s). This perspective explores evidence related to our working hypothesis that genetic variations in innate immune genes, particularly those that are commonly found but unevenly distributed between populations, are associated with disparities between populations in both cancer and cardio-metabolic diseases. Identifying conventional and unconventional innate immune genes that fit this profile may provide critical insights into the underlying mechanisms that connect these two families of complex diseases and offer novel targets for precision-based treatment of cancer and/or cardio-metabolic disease.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Integrating functional scoring and regulatory data to predict the effect of non-coding SNPs in a complex neurological disease

Most SNPs associated with complex diseases seem to lie in non-coding regions of the genome; however, their contribution to gene expression and disease phenotype remains poorly understood. Here, we established a workflow to provide assistance in prioritising the functional relevance of non-coding SNPs of candidate genes as susceptibility loci in polygenic neurological disorders. To illustrate the applicability of our workflow, we considered the multifactorial disorder migraine as a model to follow our step-by-step approach. We annotated the overlap of selected SNPs with regulatory elements and assessed their potential impact on gene expression based on publicly available prediction algorithms and functional genomics information. Some migraine risk loci have been hypothesised to reside in non-coding regions and to be implicated in the neurotransmission pathway. In this study, we used a set of 22 non-coding SNPs from neurotransmission and synaptic machinery-related genes previously suggested to be involved in migraine susceptibility based on our candidate gene association studies. After prioritising these SNPs, we focused on non-reported ones that demonstrated high regulatory potential: (1) VAMP2_rs1150 (3' UTR) was predicted as a target of hsa-mir-5010-3p miRNA, possibly disrupting its own gene expression; (2) STX1A_rs6951030 (proximal enhancer) may affect the binding affinity of zinc-finger transcription factors (namely ZNF423) and disturb TBL2 gene expression; and (3) SNAP25_rs2327264 (distal enhancer) expected to be in a binding site of ONECUT2 transcription factor. This study demonstrated the applicability of our practical workflow to facilitate the prioritisation of potentially relevant non-coding SNPs and predict their functional impact in multifactorial neurological diseases.

SUMMIT-FA: a new resource for improved transcriptome imputation using functional annotations

Transcriptome-wide association studies (TWAS) integrate gene expression prediction models and genome-wide association studies (GWAS) to identify gene-trait associations. The power of TWAS is determined by the sample size of GWAS and the accuracy of the expression prediction model. Here, we present a new method, the Summary-level Unified Method for Modeling Integrated Transcriptome using Functional Annotations (SUMMIT-FA), which improves gene expression prediction accuracy by leveraging functional annotation resources and a large expression quantitative trait loci (eQTL) summary-level dataset. We build gene expression prediction models in whole blood using SUMMIT-FA with the comprehensive functional database MACIE and eQTL summary-level data from the eQTLGen consortium. We apply these models to GWAS for 24 complex traits and show that SUMMIT-FA identifies significantly more gene-trait associations and improves predictive power for identifying "silver standard" genes compared to several benchmark methods. We further conduct a simulation study to demonstrate the effectiveness of SUMMIT-FA.

Two redundant transcription factor binding sites in a single enhancer are essential for mammalian sex determination

Male development in mammals depends on the activity of the two SOX gene: Sry and Sox9, in the embryonic testis. As deletion of Enhancer 13 (Enh13) of the Sox9 gene results in XY male-to-female sex reversal, we explored the critical elements necessary for its function and hence, for testis and male development. Here, we demonstrate that while microdeletions of individual transcription factor binding sites (TFBS) in Enh13 lead to normal testicular development, combined microdeletions of just two SRY/SOX binding motifs can alone fully abolish Enh13 activity leading to XY male-to-female sex reversal. This suggests that for proper male development to occur, these few nucleotides of non-coding DNA must be intact. Interestingly, we show that depending on the nature of these TFBS mutations, dramatically different phenotypic outcomes can occur, providing a molecular explanation for the distinct clinical outcomes observed in patients harboring different variants in the same enhancer.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD)

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.