Repeat polymorphisms underlie top genetic risk loci for glaucoma and colorectal cancer.

Many regions in the human genome vary in length among individuals due to variable numbers of tandem repeats (VNTRs). To assess the phenotypic impact of VNTRs genome-wide, we applied a statistical imputation approach to estimate the lengths of 9,561 autosomal VNTR loci in 418,136 unrelated UK Biobank participants and 838 GTEx participants. Association and statistical fine-mapping analyses identified 58 VNTRs that appeared to influence a complex trait in UK Biobank, 18 of which also appeared to modulate expression or splicing of a nearby gene. Non-coding VNTRs at TMCO1 and EIF3H appeared to generate the largest known contributions of common human genetic variation to risk of glaucoma and colorectal cancer, respectively. Each of these two VNTRs associated with a >2-fold range of risk across individuals. These results reveal a substantial and previously unappreciated role of non-coding VNTRs in human health and gene regulation.

Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases.

Genetic studies have revealed many variant loci that are associated with immune-mediated diseases. To elucidate the disease pathogenesis, it is essential to understand the function of these variants, especially under disease-associated conditions. Here, we performed a large-scale immune cell gene-expression analysis, together with whole-genome sequence analysis. Our dataset consists of 28 distinct immune cell subsets from 337 patients diagnosed with 10 categories of immune-mediated diseases and 79 healthy volunteers. Our dataset captured distinctive gene-expression profiles across immune cell types and diseases. Expression quantitative trait loci (eQTL) analysis revealed dynamic variations of eQTL effects in the context of immunological conditions, as well as cell types. These cell-type-specific and context-dependent eQTLs showed significant enrichment in immune disease-associated genetic variants, and they implicated the disease-relevant cell types, genes, and environment. This atlas deepens our understanding of the immunogenetic functions of disease-associated variants under in vivo disease conditions.

A Thalamic Orphan Receptor Drives Variability in Short-Term Memory.

Working memory is a form of short-term memory that involves maintaining and updating task-relevant information toward goal-directed pursuits. Classical models posit persistent activity in prefrontal cortex (PFC) as a primary neural correlate, but emerging views suggest additional mechanisms may exist. We screened ∼200 genetically diverse mice on a working memory task and identified a genetic locus on chromosome 5 that contributes to a substantial proportion (17%) of the phenotypic variance. Within the locus, we identified a gene encoding an orphan G-protein-coupled receptor, Gpr12, which is sufficient to drive substantial and bidirectional changes in working memory. Molecular, cellular, and imaging studies revealed that Gpr12 enables high thalamus-PFC synchrony to support memory maintenance and choice accuracy. These findings identify an orphan receptor as a potent modifier of short-term memory and supplement classical PFC-based models with an emerging thalamus-centric framework for the mechanistic understanding of working memory.

Histone Acetylome-wide Association Study of Autism Spectrum Disorder.

The association of histone modification changes with autism spectrum disorder (ASD) has not been systematically examined. We conducted a histone acetylome-wide association study (HAWAS) by performing H3K27ac chromatin immunoprecipitation sequencing (ChIP-seq) on 257 postmortem samples from ASD and matched control brains. Despite etiological heterogeneity, ≥68% of syndromic and idiopathic ASD cases shared a common acetylome signature at >5,000 cis-regulatory elements in prefrontal and temporal cortex. Similarly, multiple genes associated with rare genetic mutations in ASD showed common "epimutations." Acetylome aberrations in ASD were not attributable to genetic differentiation at cis-SNPs and highlighted genes involved in synaptic transmission, ion transport, epilepsy, behavioral abnormality, chemokinesis, histone deacetylation, and immunity. By correlating histone acetylation with genotype, we discovered >2,000 histone acetylation quantitative trait loci (haQTLs) in human brain regions, including four candidate causal variants for psychiatric diseases. Due to the relative stability of histone modifications postmortem, we anticipate that the HAWAS approach will be applicable to multiple diseases.

Shaping immunity: The influence of natural selection on population immune diversity.

Humans exhibit considerable variability in their immune responses to the same immune challenges. Such variation is widespread and affects individual and population-level susceptibility to infectious diseases and immune disorders. Although the factors influencing immune response diversity are partially understood, what mechanisms lead to the wide range of immune traits in healthy individuals remain largely unexplained. Here, we discuss the role that natural selection has played in driving phenotypic differences in immune responses across populations and present-day susceptibility to immune-related disorders. Further, we touch on future directions in the field of immunogenomics, highlighting the value of expanding this work to human populations globally, the utility of modeling the immune response as a dynamic process, and the importance of considering the potential polygenic nature of natural selection. Identifying loci acted upon by evolution may further pinpoint variants critically involved in disease etiology, and designing studies to capture these effects will enrich our understanding of the genetic contributions to immunity and immune dysregulation.

Finding genes and pathways that underlie coral adaptation.

Mass coral bleaching is one of the clearest threats of climate change to the persistence of marine biodiversity. Despite the negative impacts of bleaching on coral health and survival, some corals may be able to rapidly adapt to warming ocean temperatures. Thus, a significant focus in coral research is identifying the genes and pathways underlying coral heat adaptation. Here, we review state-of-the-art methods that may enable the discovery of heat-adaptive loci in corals and identify four main knowledge gaps. To fill these gaps, we describe an experimental approach combining seascape genomics with CRISPR/Cas9 gene editing to discover and validate heat-adaptive loci. Finally, we discuss how information on adaptive genotypes could be used in coral reef conservation and management strategies.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis.

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease.

A major fraction of loci identified by genome-wide association studies (GWASs) mediate alternative splicing, but mechanistic interpretation is hindered by the technical limitations of short-read RNA sequencing (RNA-seq), which cannot directly link splicing events to full-length protein isoforms. Long-read RNA-seq represents a powerful tool to characterize transcript isoforms, and recently, infer protein isoform existence. Here, we present an approach that integrates information from GWASs, splicing quantitative trait loci (sQTLs), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes that colocalized with BMD associations (H4PP ≥ 0.75). We generated PacBio Iso-Seq data (N = ∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were unannotated. By casting the sQTLs onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense-mediated decay and 190 that potentially resulted in the expression of unannotated protein isoforms. Finally, we functionally validated colocalizing sQTLs in TPM2, in which siRNA-mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization but exhibited no effect upon knockdown of the entire gene. Our approach should be to generalize across diverse clinical traits and to provide insights into protein isoform activities modulated by GWAS loci.

Polygenic sex determination in vertebrates - is there any such thing?

Genetic sex determination (SD) in most vertebrates is controlled by a single master sex gene, which ensures a 1:1 sex ratio. However, more complex systems abound, and several have been ascribed to polygenic SD (PSD), in which many genes at different loci interact to produce the sexual phenotype. Here we examine claims for PSD in vertebrates, finding that most constitute transient states during sex chromosome turnover, or aberrant systems in species hybrids. To avoid confusion about terminology, we propose a consistent nomenclature for genetic SD systems.

Chromatin modules and their implication in genomic organization and gene regulation.

Regulation of gene expression is a complex but highly guided process. While genomic technologies and computational approaches have allowed high-throughput mapping of cis-regulatory elements (CREs) and their interactions in 3D, their precise role in regulating gene expression remains obscure. Recent complementary observations revealed that interactions between CREs frequently result in the formation of small-scale functional modules within topologically associating domains. Such chromatin modules likely emerge from a complex interplay between regulatory machineries assembled at CREs, including site-specific binding of transcription factors. Here, we review the methods that allow identifying chromatin modules, summarize possible mechanisms that steer CRE interactions within these modules, and discuss outstanding challenges to uncover how chromatin modules fit in our current understanding of the functional 3D genome.

Genome-wide aggregated trans-effects on risk of type 1 diabetes: A test of the "omnigenic" sparse effector hypothesis of complex trait genetics.

The "omnigenic" hypothesis postulates that the polygenic effects of common SNPs on a typical complex trait are mediated through trans-effects on expression of a relatively sparse set of effector ("core") genes. We tested this hypothesis in a study of 4,964 cases of type 1 diabetes (T1D) and 7,497 controls by using summary statistics to calculate aggregated (excluding the HLA region) trans-scores for gene expression in blood. From associations of T1D with aggregated trans-scores, nine putative core genes were identified, of which three-STAT1, CTLA4 and FOXP3-are genes in which variants cause monogenic forms of autoimmune diabetes. Seven of these genes affect the activity of regulatory T cells, and two are involved in immune responses to microbial lipids. Four T1D-associated genomic regions could be identified as master regulators via trans-effects on gene expression. These results support the sparse effector hypothesis and reshape our understanding of the genetic architecture of T1D.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension.

BACKGROUND: Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms. METHODS: We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension. RESULTS: We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus. CONCLUSIONS: Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

Yield-related quantitative trait loci identification and lint percentage hereditary dissection under salt stress in upland cotton.

Salinity is frequently mentioned as a major constraint in worldwide agricultural production. Lint percentage (LP) is a crucial yield-component in cotton lint production. While the genetic factors affect cotton yield in saline soils are still unclear. Here, we employed a recombinant inbred line population in upland cotton (Gossypium hirsutum L.) and investigated the effects of salt stress on five yield and yield component traits, including seed cotton yield per plant, lint yield per plant, boll number per plant, boll weight, and LP. Between three datasets of salt stress (E1), normal growth (E2), and the difference values dataset of salt stress and normal conditions (D-value), 87, 82, and 55 quantitative trait loci (QTL) were detectable, respectively. In total, five QTL (qLY-Chr6-2, qBNP-Chr4-1, qBNP-Chr12-1, qBNP-Chr15-5, qLP-Chr19-2) detected in both in E1 and D-value were salt related QTL, and three stable QTL (qLP-Chr5-3, qLP-Chr13-1, qBW-Chr5-5) were detected both in E1 and E2 across 3 years. Silencing of nine genes within a stable QTL (qLP-Chr5-3) highly expressed in fiber developmental stages increased LP and decreased fiber length (FL), indicating that multiple minor-effect genes clustered on Chromosome 5 regulate LP and FL. Additionally, the difference in LP caused by Gh_A05G3226 is mainly in transcription level rather than in the sequence difference. Moreover, silencing of salt related gene (GhDAAT) within qBNP-Chr4-1 decreased salt tolerance in cotton. Our findings shed light on the regulatory mechanisms underlining cotton salt tolerance and fiber initiation.

The Genetics of Plant Metabolism.

Plant metabolic studies have traditionally focused on the role and regulation of the enzymes catalyzing key reactions within specific pathways. Within the past 20 years, reverse genetic approaches have allowed direct determination of the effects of the deficiency, or surplus, of a given protein on the biochemistry of a plant. In parallel, top-down approaches have also been taken, which rely on screening broad, natural genetic diversity for metabolic diversity. Here, we compare and contrast the various strategies that have been adopted to enhance our understanding of the natural diversity of metabolism. We also detail how these approaches have enhanced our understanding of both specific and global aspects of the genetic regulation of metabolism. Finally, we discuss how such approaches are providing important insights into the evolution of plant secondary metabolism.

Genetic control of RNA editing in neurodegenerative disease.

A-to-I RNA editing diversifies human transcriptome to confer its functional effects on the downstream genes or regulations, potentially involving in neurodegenerative pathogenesis. Its variabilities are attributed to multiple regulators, including the key factor of genetic variants. To comprehensively investigate the potentials of neurodegenerative disease-susceptibility variants from the view of A-to-I RNA editing, we analyzed matched genetic and transcriptomic data of 1596 samples across nine brain tissues and whole blood from two large consortiums, Accelerating Medicines Partnership-Alzheimer's Disease and Parkinson's Progression Markers Initiative. The large-scale and genome-wide identification of 95 198 RNA editing quantitative trait loci revealed the preferred genetic effects on adjacent editing events. Furthermore, to explore the underlying mechanisms of the genetic controls of A-to-I RNA editing, several top RNA-binding proteins were pointed out, such as EIF4A3, U2AF2, NOP58, FBL, NOP56 and DHX9, since their regulations on multiple RNA-editing events were probably interfered by these genetic variants. Moreover, these variants may also contribute to the variability of other molecular phenotypes associated with RNA editing, including the functions of 3 proteins, expressions of 277 genes and splicing of 449 events. All the analyses results shown in NeuroEdQTL (https://relab.xidian.edu.cn/NeuroEdQTL/) constituted a unique resource for the understanding of neurodegenerative pathogenesis from genotypes to phenotypes related to A-to-I RNA editing.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes.

BACKGROUND: Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes. METHODS: We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome. RESULTS: We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs. CONCLUSIONS: Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Connecting dementia risk loci to the CSF proteome identifies pathophysiological leads for dementia.

Genome-wide association studies have successfully identified many genetic risk loci for dementia, but exact biological mechanisms through which genetic risk factors contribute to dementia remains unclear. Integrating CSF proteomic data with dementia risk loci could reveal intermediate molecular pathways connecting genetic variance to the development of dementia. We tested to what extent effects of known dementia risk loci can be observed in CSF levels of 665 proteins [proximity extension-based (PEA) immunoassays] in a deeply-phenotyped mixed memory clinic cohort [n = 502, mean age (standard deviation, SD) = 64.1 (8.7) years, 181 female (35.4%)], including patients with Alzheimer's disease (AD, n = 213), dementia with Lewy bodies (DLB, n = 50) and frontotemporal dementia (FTD, n = 93), and controls (n = 146). Validation was assessed in independent cohorts (n = 99 PEA platform, n = 198, mass reaction monitoring-targeted mass spectroscopy and multiplex assay). We performed additional analyses stratified according to diagnostic status (AD, DLB, FTD and controls separately), to explore whether associations between CSF proteins and genetic variants were specific to disease or not. We identified four AD risk loci as protein quantitative trait loci (pQTL): CR1-CR2 (rs3818361, P = 1.65 × 10-8), ZCWPW1-PILRB (rs1476679, P = 2.73 × 10-32), CTSH-CTSH (rs3784539, P = 2.88 × 10-24) and HESX1-RETN (rs186108507, P = 8.39 × 10-8), of which the first three pQTLs showed direct replication in the independent cohorts. We identified one AD-specific association between a rare genetic variant of TREM2 and CSF IL6 levels (rs75932628, P = 3.90 × 10-7). DLB risk locus GBA showed positive trans effects on seven inter-related CSF levels in DLB patients only. No pQTLs were identified for FTD loci, either for the total sample as for analyses performed within FTD only. Protein QTL variants were involved in the immune system, highlighting the importance of this system in the pathophysiology of dementia. We further identified pQTLs in stratified analyses for AD and DLB, hinting at disease-specific pQTLs in dementia. Dissecting the contribution of risk loci to neurobiological processes aids in understanding disease mechanisms underlying dementia.

Deep learning predicts DNA methylation regulatory variants in the human brain and elucidates the genetics of psychiatric disorders.

There is growing evidence for the role of DNA methylation (DNAm) quantitative trait loci (mQTLs) in the genetics of complex traits, including psychiatric disorders. However, due to extensive linkage disequilibrium (LD) of the genome, it is challenging to identify causal genetic variations that drive DNAm levels by population-based genetic association studies. This limits the utility of mQTLs for fine-mapping risk loci underlying psychiatric disorders identified by genome-wide association studies (GWAS). Here we present INTERACT, a deep learning model that integrates convolutional neural networks with transformer, to predict effects of genetic variations on DNAm levels at CpG sites in the human brain. We show that INTERACT-derived DNAm regulatory variants are not confounded by LD, are concentrated in regulatory genomic regions in the human brain, and are convergent with mQTL evidence from genetic association analysis. We further demonstrate that predicted DNAm regulatory variants are enriched for heritability of brain-related traits and improve polygenic risk prediction for schizophrenia across diverse ancestry samples. Finally, we applied predicted DNAm regulatory variants for fine-mapping schizophrenia GWAS risk loci to identify potential novel risk genes. Our study shows the power of a deep learning approach to identify functional regulatory variants that may elucidate the genetic basis of complex traits.

The genetic regulation of the gastric transcriptome is associated with metabolic and obesity-related traits and diseases.

Tissue-specific gene expression and gene regulation lead to a better understanding of tissue-specific physiology and pathophysiology. We analyzed the transcriptome and genetic regulatory profiles of two distinct gastric sites, corpus and antrum, to identify tissue-specific gene expression and its regulation. Gastric corpus and antrum mucosa biopsies were collected during routine gastroscopies from up to 431 healthy individuals. We obtained genotype and transcriptome data and performed transcriptome profiling and expression quantitative trait locus (eQTL) studies. We further used data from genome-wide association studies (GWAS) of various diseases and traits to partition their heritability and to perform transcriptome-wide association studies (TWAS). The transcriptome data from corpus and antral mucosa highlights the heterogeneity of gene expression in the stomach. We identified enriched pathways revealing distinct and common physiological processes in gastric corpus and antrum. Furthermore, we found an enrichment of the single nucleotide polymorphism (SNP)-based heritability of metabolic, obesity-related, and cardiovascular traits and diseases by considering corpus- and antrum-specifically expressed genes. Particularly, we could prioritize gastric-specific candidate genes for multiple metabolic traits, like NQO1 which is involved in glucose metabolism, MUC1 which contributes to purine and protein metabolism or RAB27B being a regulator of weight and body composition. Our findings show that gastric corpus and antrum vary in their transcriptome and genetic regulatory profiles indicating physiological differences which are mostly related to digestion and epithelial protection. Moreover, our findings demonstrate that the genetic regulation of the gastric transcriptome is linked to biological mechanisms associated with metabolic, obesity-related, and cardiovascular traits and diseases. NEW & NOTEWORTHY We analyzed the transcriptomes and genetic regulatory profiles of gastric corpus and for the first time also of antrum mucosa in 431 healthy individuals. Through tissue-specific gene expression and eQTL analyses, we uncovered unique and common physiological processes across both primary gastric sites. Notably, our findings reveal that stomach-specific eQTLs are enriched in loci associated with metabolic traits and diseases, highlighting the pivotal role of gene expression regulation in gastric physiology and potential pathophysiology.