Putative Causal Genes Research Articles

EQTL in diseased colon tissue identifies novel target genes associated with IBD.

Genome-wide association studies (GWAS) have identified over 300 loci associated with the inflammatory bowel diseases (IBD), but putative causal genes for most are unknown. We conducted the largest disease-focused expression quantitative trait loci (eQTL) analysis using colon tissue from 252 IBD patients to determine genetic effects on gene expression and potential contribution to IBD. Combined with two non-IBD colon eQTL studies, we identified 194 potential target genes for 108 GWAS loci. eQTL in IBD tissue were enriched for IBD GWAS loci colocalizations, provided novel evidence for IBD-associated genes such as ABO and TNFRSF14 , and identified additional target genes compared to non-IBD tissue eQTL. IBD-associated eQTL unique to diseased tissue had distinct regulatory and functional characteristics with increased effect sizes. Together, these highlight the importance of eQTL studies in diseased tissue for understanding functional consequences of genetic variants, and elucidating molecular mechanisms and regulation of key genes involved in IBD.

A functional variant rs912304 for late-onset T1D risk contributes to islet dysfunction by regulating proinflammatory cytokine-responsive gene STXBP6 expression

BackgroundOur previous genome‑wide association studies (GWAS) have suggested rs912304 in 14q12 as a suggestive risk variant for type 1 diabetes (T1D). However, the association between this risk region and T1D subgroups and related clinical risk features, the underlying causal functional variant(s), putative candidate gene(s), and related mechanisms are yet to be elucidated.MethodsWe assessed the association between variant rs912304 and T1D, as well as islet autoimmunity and islet function, stratified by the diagnosed age of 12. We used epigenome bioinformatics analyses, dual luciferase reporter assays, and expression quantitative trait loci (eQTL) analyses to prioritize the most likely functional variant and potential causal gene. We also performed functional experiments to evaluate the role of the causal gene on islet function and its related mechanisms.ResultsWe identified rs912304 as a risk variant for T1D subgroups with diagnosed age ≥ 12 but not < 12. This variant is associated with residual islet function but not islet-specific autoantibody positivity in T1D individuals. Bioinformatics analysis indicated that rs912304 is a functional variant exhibiting spatial overlaps with enhancer active histone marks (H3K27ac and H3K4me1) and open chromatin status (ATAC-seq) in the human pancreas and islet tissues. Luciferase reporter gene assays and eQTL analyses demonstrated that the biallelic sites of rs912304 had differential allele-specific enhancer activity in beta cell lines and regulated STXBP6 expression, which was defined as the most putative causal gene based on Open Targets Genetics, GTEx v8 and Tiger database. Moreover, Stxbp6 was upregulated by T1D-related proinflammatory cytokines but not high glucose/fat. Notably, Stxbp6 over-expressed INS-1E cells exhibited decreasing insulin secretion and increasing cell apoptosis through Glut1 and Gadd45β, respectively.ConclusionsThis study expanded the genomic landscape regarding late-onset T1D risk and supported islet function mechanistically connected to T1D pathogenesis.

A functional genomic framework to elucidate novel causal metabolic dysfunction-associated fatty liver disease genes.

Metabolic dysfunction-associated fatty liver disease (MASLD) is the most prevalent chronic liver pathology in western countries, with serious public health consequences. Efforts to identify causal genes for MASLD have been hampered by the relative paucity of human data from gold-standard magnetic resonance quantification of hepatic fat. To overcome insufficient sample size, genome-wide association studies using MASLD surrogate phenotypes have been used, but only a small number of loci have been identified to date. In this study, we combined GWAS of MASLD composite surrogate phenotypes with genetic colocalization studies followed by functional in vitro screens to identify bona fide causal genes for MASLD. We used the UK Biobank to explore the associations of our novel MASLD score, and genetic colocalization to prioritize putative causal genes for in vitro validation. We created a functional genomic framework to study MASLD genes in vitro using CRISPRi. Our data identify VKORC1, TNKS, LYPLAL1 and GPAM as regulators of lipid accumulation in hepatocytes and suggest the involvement of VKORC1 in the lipid storage related to the development of MASLD. Complementary genetic and genomic approaches are useful for the identification of MASLD genes. Our data supports VKORC1 as a bona fide MASLD gene. We have established a functional genomic framework to study at scale putative novel MASLD genes from human genetic association studies.

Druggable genome-wide Mendelian randomization for identifying the role of integrated stress response in therapeutic targets of bipolar disorder

For bipolar disorder (BD), the inconsistency of treatment guidelines and the long phases of pharmacological adjustment remain major challenges. BD is known to be comorbid with many medical and psychiatric conditions and they may share inflammatory and stress-related aetiologies, which could give rise to this association. The integrated stress response (ISR) responds to various stress conditions that lead to alterations in cellular homeostasis. However, as a causative mechanism underlying cognitive deficits and neurodegeneration in a broad range of brain disorders, the impact of ISR on BD is understudied. Mendelian randomization has been widely used to repurpose licensed drugs and discover novel therapeutic targets. Thus, we aimed to identify novel therapeutic targets for BD and analyze their pathophysiological mechanisms, using the summary data-based Mendelian Randomization (SMR) and Bayesian colocalization (COLOC) methods to integrate the summary-level data of the GWAS on BD and the expression quantitative trait locus (eQTL) study in blood. We utilized the GWAS data including 41,917 BD cases and 371,549 controls from the Psychiatric Genomics Consortium and the eQTL data from 31,684 participants of predominantly European ancestry from the eQTLGen consortium. The SMR analysis identified the EIF2B5 gene that was associated with BD due to no linkage but pleiotropy or causality. The COLOC analysis strongly suggested that EIF2B5 and the trait of BD were affected by shared causal variants, and thus were colocalized. Utilizing data in EpiGraphDB we find other putative causal BD genes (EIF2AK4 and GSK3B) to prioritize potential alternative drug targets.

Genome-Wide Association Study Identifies Genes for Hair Growth and Patterning are Associated With Pilonidal Disease.

Pilonidal sinus disease is a highly morbid condition characterized by the formation of chronic sinus tracts throughout the sacrococcygeal region. Despite its commonality and strong association with family history, no prior investigation of genetic risk factors for pilonidal sinus disease exists. To identify genetic risk factors for pilonidal sinus disease. A genome-wide association study. The United Kingdom Biobank, FinnGen Biobank, and Penn Medicine BioBank. There were 772,072 participants. Genome-wide significant variants ( p < 5 × 10 -8 ) were mapped to genes using physical distance and gene expression in skin. Genetic correlation between pilonidal sinus disease and morphometric, androgen-driven, and hair phenotypes was estimated with linkage disequilibrium score regression. Finally, a genome-first approach to rare predicted deleterious variants in hair shaft genes TCHH , PADI3 , and TGM3 was conducted for association with pilonidal sinus disease via the Penn Medicine BioBank. A genome-wide association study comprising 2835 individuals with pilonidal sinus disease identified 5 genome-wide significant loci, prioritizing HDAC9, TBX15, WARS2, RP11-293M10.1 , PRKAR1B , TWIST1, GPATCH2L, NEK9 , and EIF2B2 , as putative causal genes; several of these genes have known roles in balding and hair patterning. There was a significant correlation between the genetic background of pilonidal sinus disease and the androgen-driven hair traits of male pattern baldness and young age at first facial hair. In a candidate analysis of genes associated with syndromic hair disorders, rare coding variants in TCHH , a monogenic cause of uncombable hair syndrome, were associated with increased prevalence of pilonidal sinus disease (OR 4.81 [95% CI, 2.06-11.2]). This study is limited to European ancestry. However, because there is a higher incidence of pilonidal sinus disease in men of European ancestry, this analysis is focused on the at-risk population. Genetic analysis of pilonidal sinus disease identified shared genetic architecture with hair biology and androgen-driven traits. As the first study investigating the genetic basis of pilonidal sinus disease, this provides biological insight into the long-appreciated connection between the disease state, male sex, and hair. See Video abstract. ANTECEDENTES:La enfermedad del seno pilonidal es una condición muy mórbida caracterizada por la formación de tractos sinusales crónicos en toda la región sacrococcígea. A pesar de su frecuencia y su fuerte asociación con los antecedentes familiares, no se han investigado previamente los factores de riesgo genéticos de la enfermedad sinusal pilonidal.OBJETIVO:Identificar factores genéticos de riesgo para la enfermedad del seno pilonidal.DISEÑO:Estudio de asociación de genoma completo.CONJUNTOS:Biobanco del Reino Unido, Biobanco FinnGen y Biobanco PennMedicine.PACIENTES:772.072 participantes.MEDIDA DE RESULTADO PRINCIPAL:Las variantes significativas en todo el genoma (p < 5x10-8) se asignaron a genes utilizando la distancia física y la expresión génica en la piel. La correlación genética entre la enfermedad del seno pilonidal y los fenotipos morfométricos, androgénicos y de cabello se estimó con regresión de puntuación LD. Por último, se realizó una aproximación genómica a variantes deletéreas raras predichas en los genes del tallo piloso TCHH, PADI3 y TGM3 para su asociación con la enfermedad del seno pilonidal a través del Biobanco PennMedicine.RESULTADOS:El estudio de asociación de todo el genoma, que incluyó a 2.835 individuos con enfermedad del seno pilonidal, identificó 5 loci significativos en todo el genoma, dando prioridad a HDAC9, TBX15, WARS2, RP11-293M10.1, PRKAR1B, TWIST1, GPATCH2L, NEK9 y EIF2B2, como genes causales putativos; varios de estos genes tienen funciones conocidas en la calvicie y el patrón del cabello. Se observó una correlación significativa entre los antecedentes genéticos de la enfermedad del seno pilonidal y los de los rasgos calvicie de patrón masculino y edad temprana del primer vello facial impulsados por andrógenos. En un análisis de genes candidatos asociados a trastornos capilares sindrómicos, las variantes raras de codificación en TCHH, una causa monogénica del síndrome capilar incombustible, se asociaron a una mayor prevalencia de la enfermedad del seno pilonidal (OR 4,81 [IC del 5%, 2,06-11,2]).LIMITACIONES:Este estudio se limita a la ascendencia europea. Sin embargo, debido a que hay una mayor incidencia de la enfermedad sinusal pilonidal en los hombres de ascendencia europea, este análisis se centra en la población de riesgo.CONCLUSIÓN:El análisis genético de la enfermedad del seno pilonidal identificó una arquitectura genética compartida con la biología del cabello y los rasgos impulsados por andrógenos. Siendo el primer estudio que investiga las bases genéticas de la enfermedad del seno pilonidal, esto proporciona una visión biológica de la conexión, apreciada desde hace tiempo, entre el estado de la enfermedad, el sexo masculino y el cabello. (Traducción-Dr. Aurian Garcia Gonzalez ).

A multi-ancestry cerebral cortex transcriptome-wide association study identifies genes associated with smoking behaviors.

Transcriptome-wide association studies (TWAS) have provided valuable insight in identifying genes that may impact cigarette smoking. Most of previous studies, however, mainly focused on European ancestry. Limited TWAS studies have been conducted across multiple ancestries to explore genes that may impact smoking behaviors. In this study, we used cis-eQTL data of cerebral cortex from multiple ancestries in MetaBrain, including European, East Asian, and African samples, as reference panels to perform multi-ancestry TWAS analyses on ancestry-matched GWASs of four smoking behaviors including smoking initiation, smoking cessation, age of smoking initiation, and number of cigarettes per day in GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN). Multiple-ancestry fine-mapping approach was conducted to identify credible gene sets associated with these four traits. Enrichment and module network analyses were further performed to explore the potential roles of these identified gene sets. A total of 719 unique genes were identified to be associated with at least one of the four smoking traits across ancestries. Among those, 249 genes were further prioritized as putative causal genes in multiple ancestry-based fine-mapping approach. Several well-known smoking-related genes, including PSMA4, IREB2, and CHRNA3, showed high confidence across ancestries. Some novel genes, e.g., TSPAN3 and ANK2, were also identified in the credible sets. The enrichment analysis identified a series of critical pathways related to smoking such as synaptic transmission and glutamate receptor activity. Leveraging the power of the latest multi-ancestry GWAS and eQTL data sources, this study revealed hundreds of genes and relevant biological processes related to smoking behaviors. These findings provide new insights for future functional studies on smoking behaviors.

The genetic landscape of substance use disorders.

Substance use disorders represent a significant public health concern with considerable socioeconomic implications worldwide. Twin and family-based studies have long established a heritable component underlying these disorders. In recent years, genome-wide association studies of large, broadly phenotyped samples have identified regions of the genome that harbour genetic risk variants associated with substance use disorders. These regions have enabled the discovery of putative causal genes and improved our understanding of genetic relationships among substance use disorders and other traits. Furthermore, the integration of these data with clinical information has yielded promising insights into how individuals respond to medications, allowing for the development of personalized treatment approaches based on an individual's genetic profile. This review article provides an overview of recent advances in the genetics of substance use disorders and demonstrates how genetic data may be used to reduce the burden of disease and improve public health outcomes.

Response eQTLs, chromatin accessibility, and 3D chromatin structure in chondrocytes provide mechanistic insight into osteoarthritis risk.

∘ Comprehensive analysis of sex- and age-related global gene expression in human chondrocytes revealed differences that correlate with osteoarthritis ∘ First response eQTLs in chondrocytes treated with an OA-related stimulus ∘ Deeply sequenced Hi-C in resting and activated chondrocytes helps connect OA risk variants to their putative causal genes ∘ Colocalization analysis reveals 13 (including 10 novel) putative OA risk genes.

GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson’s disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Rare copy-number variants as modulators of common disease susceptibility

BackgroundCopy-number variations (CNVs) have been associated with rare and debilitating genomic disorders (GDs) but their impact on health later in life in the general population remains poorly described.MethodsAssessing four modes of CNV action, we performed genome-wide association scans (GWASs) between the copy-number of CNV-proxy probes and 60 curated ICD-10 based clinical diagnoses in 331,522 unrelated white British UK Biobank (UKBB) participants with replication in the Estonian Biobank.ResultsWe identified 73 signals involving 40 diseases, all of which indicating that CNVs increased disease risk and caused earlier onset. We estimated that 16% of these associations are indirect, acting by increasing body mass index (BMI). Signals mapped to 45 unique, non-overlapping regions, nine of which being linked to known GDs. Number and identity of genes affected by CNVs modulated their pathogenicity, with many associations being supported by colocalization with both common and rare single-nucleotide variant association signals. Dissection of association signals provided insights into the epidemiology of known gene-disease pairs (e.g., deletions in BRCA1 and LDLR increased risk for ovarian cancer and ischemic heart disease, respectively), clarified dosage mechanisms of action (e.g., both increased and decreased dosage of 17q12 impacted renal health), and identified putative causal genes (e.g., ABCC6 for kidney stones). Characterization of the pleiotropic pathological consequences of recurrent CNVs at 15q13, 16p13.11, 16p12.2, and 22q11.2 in adulthood indicated variable expressivity of these regions and the involvement of multiple genes. Finally, we show that while the total burden of rare CNVs—and especially deletions—strongly associated with disease risk, it only accounted for ~ 0.02% of the UKBB disease burden. These associations are mainly driven by CNVs at known GD CNV regions, whose pleiotropic effect on common diseases was broader than anticipated by our CNV-GWAS.ConclusionsOur results shed light on the prominent role of rare CNVs in determining common disease susceptibility within the general population and provide actionable insights for anticipating later-onset comorbidities in carriers of recurrent CNVs.

Transcriptome-wide association analyses reveal the impact of regulatory variants on rice panicle architecture and causal gene regulatory networks

Panicle architecture is a key determinant of rice grain yield and is mainly determined at the 1-2 mm young panicle stage. Here, we investigated the transcriptome of the 1-2 mm young panicles from 275 rice varieties and identified thousands of genes whose expression levels were associated with panicle traits. Multimodel association studies suggested that many small-effect genetic loci determine spikelet per panicle (SPP) by regulating the expression of genes associated with panicle traits. We found that alleles at cis-expression quantitative trait loci of SPP-associated genes underwent positive selection, with a strong preference for alleles increasing SPP. We further developed a method that integrates the associations of cis- and trans-expression components of genes with traits to identify causal genes at even small-effect loci and construct regulatory networks. We identified 36 putative causal genes of SPP, including SDT (MIR156j) and OsMADS17, and inferred that OsMADS17 regulates SDT expression, which was experimentally validated. Our study reveals the impact of regulatory variants on rice panicle architecture and provides new insights into the gene regulatory networks of panicle traits.

Genome-wide association study of thyroid-stimulating hormone highlights new genes, pathways and associations with thyroid disease

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

Pervasive Sharing of Causal Genetic Risk Factors Contributes to Clinical and Molecular Overlap between Sjögren's Disease and Systemic Lupus Erythematosus.

SjD (Sjögren's Disease) and SLE (Systemic Lupus Erythematosus) are similar diseases. There is extensive overlap between the two in terms of both clinical features and pathobiologic mechanisms. Shared genetic risk is a potential explanation of this overlap. In this study, we evaluated whether these diseases share causal genetic risk factors. We compared the causal genetic risk for SLE and SjD using three complementary approaches. First, we examined the published GWAS results for these two diseases by analyzing the predicted causal gene protein-protein interaction networks of both diseases. Since this method does not account for overlapping risk intervals, we examined whether such intervals also overlap. Third, we used two-sample Mendelian randomization (two sample MR) using GWAS summary statistics to determine whether risk variants for SLE are causal for SjD and vice versa. We found that both the putative causal genes and the genomic risk intervals for SLE and SjD overlap 28- and 130-times more than expected by chance (p < 1.1 × 10-24 and p < 1.1 × 10-41, respectively). Further, two sample MR analysis confirmed that alone or in aggregate, SLE is likely causal for SjD and vice versa. [SjD variants predicting SLE: OR = 2.56; 95% CI (1.98-3.30); p < 1.4 × 10-13, inverse-variance weighted; SLE variants predicting SjD: OR = 1.36; 95% CI (1.26-1.47); p < 1.6 × 10-11, inverse-variance weighted]. Notably, some variants have disparate impact in terms of effect size across disease states. Overlapping causal genetic risk factors were found for both diseases using complementary approaches. These observations support the hypothesis that shared genetic factors drive the clinical and pathobiologic overlap between these diseases. Our study has implications for both differential diagnosis and future genetic studies of these two conditions.

Genome-wide meta-analysis implicates variation affecting mast cell biology in urticaria

Urticaria is characterized by inappropriate mast cell degranulation leading to the development of wheals and/or angioedema. Twin and family studies indicate that there is a substantial heritable component to urticaria risk. Our aim was to identify genomic loci at which common genetic variation influences urticaria susceptibility. Genome-wide association studies of urticaria (including all subtypes) from 3 European cohorts (UK Biobank, FinnGen, and the Trøndelag Health Study [HUNT]) were combined through statistical meta-analysis (14,306 urticaria cases and 650,664 controls). Cases were identified via electronic health care records from primary and/or secondary care. To identify putative causal variants and genes, statistical fine-mapping, colocalization, and interrogation of publicly available single-cell transcriptome sequencing resources were performed. Genome-wide significant associations (P< 5×10-8) were identified at 6 independent loci. These included 2 previously reported association signals at 1q44 and the human leucocyte antigen region on chromosome 6. Genes with expected or established roles in mast cell biology were associated with the 4 other genome-wide association signals (GCSAML, FCER1A, TPSAB1, and CBLB). Colocalization of association signals consistent with the presence of shared causal variants was observed between urticaria susceptibility and increased expression of GCSAML (posterior probability of colocalization [PPcoloc]= 0.89) and FCER1A (PPcoloc= 0.91) in skin. Common genetic variation influencing the risk of developing urticaria was identified at 6 genomic loci. The relationship between genes with roles in mast cell biology and several association signals implicates genetic variability of specific components of mast cell function in the development of urticaria.

Genome-wide association study identifies functional genomic variants associated with young stock survival in Nordic Red Dairy Cattle

Identifying quantitative trait loci (QTL) associated with calf survival is essential for both reducing economic loss in cattle industry and understanding the genetic basis of the trait. To identify mutations and genes underlying young stock survival (YSS), we performed GWAS using de-regressed estimated breeding values of a YSS index and its component traits defined by sex and age in 3,077 Nordic Red Dairy Cattle (RDC) bulls and 2 stillbirth traits (first lactation and later lactations) in 5,141 RDC bulls. Two associated QTL regions on Bos taurus autosome (BTA) 4 and 6 were identified for the YSS index. The results of 4 YSS component traits indicate that same QTL regions were associated with bull and heifer calf mortality, but the effects were different over the growing period and suggested an additional QTL on BTA23. The GWAS on stillbirth identified 3 additional QTL regions on BTA5, 14, and 24 compared with YSS and its component traits. The conditional test of BTA6 showed at least 2 closely located QTL segregating for YSS component traits and stillbirth. We found 2 independent QTL for stillbirth on BTA23. The post-GWAS revealed LCORL, PPM1K, SSP1, MED28, and LAP3 are putative causal genes on BTA6, and a frame shift variant within LCORL, BTA6:37401770 (rs384548488) could be the putative causal variant. On BTA4, the GRB10 gene is the putative causal gene and BTA4:5296018 is the putative causal variant. In addition, NDUFA9 and FGF23 on BTA5, LYN on BTA14, and KCNK5 on BTA23 are putative causal genes for QTL for stillbirth. The gene analysis also proposed several candidate genes. Our findings shed new light on the candidate genes affecting calf survival, and the knowledge could be utilized to reduce calf mortality and thereby enhance welfare of dairy cattle.

Cross-trait genome-wide association analysis of C-reactive protein level and psychiatric disorders

C-reactive protein (CRP) tends to be elevated in individuals with psychiatric disorders. Recent findings have suggested a protective effect of the genetic liability to elevated CRP on schizophrenia risk and a causative effect on depression despite weak genetic correlations, while causal relationships with bipolar disorder were inconclusive. We investigated the shared genetic underpinnings of psychiatric disorders and variation in CRP levels. Genome-wide association studies for CRP (n = 575,531), bipolar disorder (n = 413,466), depression (n = 480,359), and schizophrenia (n = 130,644) were used in causal mixture models to compare CRP with psychiatric disorders based on polygenicity, discoverability, and genome-wide genetic overlap. The conjunctional false discovery rate method was used to identify specific shared genetic loci. Shared variants were mapped to putative causal genes, which were tested for overrepresentation among gene ontology gene-sets. CRP was six to ten times less polygenic (n = 1400 vs 8600–14,500 variants) and had a discoverability one to two orders of magnitude higher than psychiatric disorders. Most CRP-associated variants were overlapping with psychiatric disorders. We identified 401 genetic loci jointly associated with CRP and psychiatric disorders with mixed effect directions. Gene-set enrichment analyses identified predominantly CNS-related gene sets for CRP and each of depression and schizophrenia, and basic cellular processes for CRP and bipolar disorder. In conclusion, CRP has a markedly different genetic architecture to psychiatric disorders, but the majority of CRP associated variants are also implicated in psychiatric disorders. Shared genetic loci implicated CNS-related processes to a greater extent than immune processes, which may have implications for how we conceptualise causal relationships between CRP and psychiatric disorders.

A genome-wide association study of coat color in Chinese Rex rabbits.

Coat color is an important phenotypic characteristic of the domestic rabbit (Oryctolagus cuniculus) and has specific economic importance in the Rex rabbit industry. Coat color varies considerably among different populations of rabbits, and several causal genes for this variation have been thoroughly studied. Nevertheless, the candidate genes affecting coat color variation in Chinese Rex rabbits remained to be investigated. In this study, we collected blood samples from 250 Chinese Rex rabbits with six different coat colors. We performed genome sequencing using a restriction site-associated DNA sequencing approach. A total of 91,546 single nucleotide polymorphisms (SNPs), evenly distributed among 21 autosomes, were identified. Genome-wide association studies (GWAS) were performed using a mixed linear model, in which the individual polygenic effect was fitted as a random effect. We detected a total of 24 significant SNPs that were located within a genomic region on chromosome 4 (OCU4). After re-fitting the most significant SNP (OCU4:13,434,448, p = 1.31e-12) as a covariate, another near-significant SNP (OCU4:11,344,946, p = 7.03e-07) was still present. Hence, we conclude that the 2.1-Mb genomic region located between these two significant SNPs is significantly associated with coat color in Chinese Rex rabbits. The well-studied coat-color-associated agouti signaling protein (ASIP) gene is located within this region. Furthermore, low genetic differentiation was also observed among the six coat color varieties. In conclusion, our results confirmed that ASIP is a putative causal gene affecting coat color variation in Chinese Rex rabbits.

Cotton ethylene response factor GhERF91 is involved in the defense against Verticillium dahliae

Verticillium dahliae causes significant losses in cotton production. To reveal the mechanism of the defense response to V. dahliae in cotton, transcriptomic analyses were performed using cotton cultivars M138 (V. dahliae-resistant) and P2 (V. dahliae-susceptible). The results revealed 11,076 and 6,640 differentially expressed genes (DEGs) in response to V. dahliae, respectively. The weighted gene co-expression network analysis of 4,633 transcription factors (TFs) indicated a “MEblue” module containing 654 TFs that strongly correlate with resistance to V. dahliae. Among these TFs, the ethylene response factor Ghi_A05G10166 (GhERF91) was identified as a putative hub gene with a defense response against V. dahliae. A virus-induced gene silencing assay and exogenous application of ethephon showed that GhERF91 is activated by ethylene and positively regulates the response to V. dahliae exposure in cotton. This study provides fundamental transcriptome data and a putative causal gene (GhERF91) associated with resistance to V. dahliae, as well as genetic resources for breeding V. dahliae-resistant cotton.

Identification of drug repurposing candidates for the treatment of anxiety: A genetic approach

Anxiety disorders are a group of prevalent and heritable neuropsychiatric diseases. We previously conducted a genome-wide association study (GWAS) which identified genomic loci associated with anxiety; however, the biological consequences underlying the genetic associations are largely unknown. Integrating GWAS and functional genomic data may improve our understanding of the genetic effects on intermediate molecular phenotypes such as gene expression. This can provide an opportunity for the discovery of drug targets for anxiety via drug repurposing. We used the GWAS summary statistics to determine putative causal genes for anxiety using MAGMA and colocalization analyses. A transcriptome-wide association study was conducted to identify genes with differential genetically regulated levels of gene expression in human brain tissue. The genes were integrated with a large drug-gene expression database (Connectivity Map), discovering compounds that are predicted to “normalise” anxiety-associated expression changes. The study identified 64 putative causal genes associated with anxiety (35 genes upregulated; 29 genes downregulated). Drug mechanisms adrenergic receptor agonists, sigma receptor agonists, and glutamate receptor agonists gene targets were enriched in anxiety-associated genetic signal and exhibited an opposing effect on the anxiety-associated gene expression signature. The significance of the project demonstrated genetic links for novel drug candidates to potentially advance anxiety therapeutics.

Identifying potential risk genes and pathways for neuropsychiatric and substance use disorders using intermediate molecular mediator information.

Neuropsychiatric and substance use disorders (NPSUDs) have a complex etiology that includes environmental and polygenic risk factors with significant cross-trait genetic correlations. Genome-wide association studies (GWAS) of NPSUDs yield numerous association signals. However, for most of these regions, we do not yet have a firm understanding of either the specific risk variants or the effects of these variants. Post-GWAS methods allow researchers to use GWAS summary statistics and molecular mediators (transcript, protein, and methylation abundances) infer the effect of these mediators on risk for disorders. One group of post-GWAS approaches is commonly referred to as transcriptome/proteome/methylome-wide association studies, which are abbreviated as T/P/MWAS (or collectively as XWAS). Since these approaches use biological mediators, the multiple testing burden is reduced to the number of genes (∼20,000) instead of millions of GWAS SNPs, which leads to increased signal detection. In this work, our aim is to uncover likely risk genes for NPSUDs by performing XWAS analyses in two tissues-blood and brain. First, to identify putative causal risk genes, we performed an XWAS using the Summary-data-based Mendelian randomization, which uses GWAS summary statistics, reference xQTL data, and a reference LD panel. Second, given the large comorbidities among NPSUDs and the shared cis-xQTLs between blood and the brain, we improved XWAS signal detection for underpowered analyses by performing joint concordance analyses between XWAS results i) across the two tissues and ii) across NPSUDs. All XWAS signals i) were adjusted for heterogeneity in dependent instruments (HEIDI) (non-causality) p-values and ii) used to test for pathway enrichment. The results suggest that there were widely shared gene/protein signals within the major histocompatibility complex region on chromosome 6 (BTN3A2 and C4A) and elsewhere in the genome (FURIN, NEK4, RERE, and ZDHHC5). The identification of putative molecular genes and pathways underlying risk may offer new targets for therapeutic development. Our study revealed an enrichment of XWAS signals in vitamin D and omega-3 gene sets. So, including vitamin D and omega-3 in treatment plans may have a modest but beneficial effect on patients with bipolar disorder.