Genetic determinants of IgG antibody response to COVID-19 vaccination.

Human humoral immune responses to SARS-CoV-2 vaccines exhibit substantial inter-individual variability and have been linked to vaccine efficacy. To elucidate the underlying mechanism behind this variability, we conducted a genome-wide association study (GWAS) on the anti-spike IgG serostatus of UK Biobank participants who were previously uninfected by SARS-CoV-2 and had received either the first dose (n = 54,066) or the second dose (n = 46,232) of COVID-19 vaccines. Our analysis revealed significant genome-wide associations between the IgG antibody serostatus following the initial vaccine and human leukocyte antigen (HLA) class II alleles. Specifically, the HLA-DRB1∗13:02 allele (MAF = 4.0%, OR = 0.75, p = 2.34e-16) demonstrated the most statistically significant protective effect against IgG seronegativity. This protective effect was driven by an alteration from arginine (Arg) to glutamic acid (Glu) at position 71 on HLA-DRß1 (p = 1.88e-25), leading to a change in the electrostatic potential of pocket 4 of the peptide binding groove. Notably, the impact of HLA alleles on IgG responses was cell type specific, and we observed a shared genetic predisposition between IgG status and susceptibility/severity of COVID-19. These results were replicated within independent cohorts where IgG serostatus was assayed by two different antibody serology tests. Our findings provide insights into the biological mechanism underlying individual variation in responses to COVID-19 vaccines and highlight the need to consider the influence of constitutive genetics when designing vaccination strategies for optimizing protection and control of infectious disease across diverse populations.

A joint transcriptome-wide association study across multiple tissues identifies candidate breast cancer susceptibility genes.

Genome-wide association studies (GWASs) have identified more than 200 genomic loci for breast cancer risk, but specific causal genes in most of these loci have not been identified. In fact, transcriptome-wide association studies (TWASs) of breast cancer performed using gene expression prediction models trained in breast tissue have yet to clearly identify most target genes. To identify candidate genes, we performed a GWAS analysis in a breast cancer dataset from UK Biobank (UKB) and combined the results with the GWAS results of the Breast Cancer Association Consortium (BCAC) by a meta-analysis. Using the summary statistics from the meta-analysis, we performed a joint TWAS analysis that combined TWAS signals from multiple tissues. We used expression prediction models trained in 11 tissues that are potentially relevant to breast cancer from the Genotype-Tissue Expression (GTEx) data. In the GWAS analysis, we identified eight loci distinct from those reported previously. In the TWAS analysis, we identified 309 genes at 108 genomic loci to be significantly associated with breast cancer at the Bonferroni threshold. Of these, 17 genes were located in eight regions that were at least 1 Mb away from published GWAS hits. The remaining TWAS-significant genes were located in 100 known genomic loci from previous GWASs of breast cancer. We found that 21 genes located in known GWAS loci remained statistically significant after conditioning on previous GWAS index variants. Our study provides insights into breast cancer genetics through mapping candidate target genes in a large proportion of known GWAS loci and discovering multiple new loci.

Integrative splicing-quantitative-trait-locus analysis reveals risk loci for non-small-cell lung cancer.

Splicing quantitative trait loci (sQTLs) have been demonstrated to contribute to disease etiology by affecting alternative splicing. However, the role of sQTLs in the development of non-small-cell lung cancer (NSCLC) remains unknown. Thus, we performed a genome-wide sQTL study to identify genetic variants that affect alternative splicing in lung tissues from 116 individuals of Chinese ancestry, which resulted in the identification of 1,385 sQTL-harboring genes (sGenes) containing 378,210 significant variant-intron pairs. A comprehensive characterization of these sQTLs showed that they were enriched in actively transcribed regions, genetic regulatory elements, and splicing-factor-binding sites. Moreover, sQTLs were largely distinct from expression quantitative trait loci (eQTLs) and showed significant enrichment in potential risk loci of NSCLC. We also integrated sQTLs into NSCLC GWAS datasets (13,327 affected individuals and 13,328 control individuals) by using splice-transcriptome-wide association study (spTWAS) and identified alternative splicing events in 19 genes that were significantly associated with NSCLC risk. By using functional annotation and experiments, we confirmed an sQTL variant, rs35861926, that reduced the risk of lung adenocarcinoma (rs35861926-T, OR = 0.88, 95% confidence interval [CI]: 0.82-0.93, p = 1.87 × 10-5) by promoting FARP1 exon 20 skipping to downregulate the expression level of the long transcript FARP1-011. Transcript FARP1-011 promoted the migration and proliferation of lung adenocarcinoma cells. Overall, our study provided informative lung sQTL resources and insights into the molecular mechanisms linking sQTL variants to NSCLC risk.

Hierarchical joint analysis of marginal summary statistics-Part II: High-dimensional instrumental analysis of omics data.

Instrumental variable (IV) analysis has been widely applied in epidemiology to infer causal relationships using observational data. Genetic variants can also be viewed as valid IVs in Mendelian randomization and transcriptome-wide association studies. However, most multivariate IV approaches cannot scale to high-throughput experimental data. Here, we leverage the flexibility of our previous work, a hierarchical model that jointly analyzes marginal summary statistics (hJAM), to a scalable framework (SHA-JAM) that can be applied to a large number of intermediates and a large number of correlated genetic variants-situations often encountered in modern experiments leveraging omic technologies. SHA-JAM aims to estimate the conditional effect for high-dimensional risk factors on an outcome by incorporating estimates from association analyses of single-nucleotide polymorphism (SNP)-intermediate or SNP-gene expression as prior information in a hierarchical model. Results from extensive simulation studies demonstrate that SHA-JAM yields a higher area under the receiver operating characteristics curve (AUC), a lower mean-squared error of the estimates, and a much faster computation speed, compared to an existing approach for similar analyses. In two applied examples for prostate cancer, we investigated metabolite and transcriptome associations, respectively, using summary statistics from a GWAS for prostate cancer with more than 140,000 men and high dimensional publicly available summary data for metabolites and transcriptomes.

Genome- and transcriptome-wide association studies of 386,000 Asian and European-ancestry women provide new insights into breast cancer genetics.

By combining data from 160,500 individuals with breast cancer and 226,196 controls of Asian and European ancestry, we conducted genome- and transcriptome-wide association studies of breast cancer. We identified 222 genetic risk loci and 137 genes that were associated with breast cancer risk at a p < 5.0 × 10-8 and a Bonferroni-corrected p < 4.6 × 10-6, respectively. Of them, 32 loci and 15 genes showed a significantly different association between ER-positive and ER-negative breast cancer after Bonferroni correction. Significant ancestral differences in risk variant allele frequencies and their association strengths with breast cancer risk were identified. Of the significant associations identified in this study, 17 loci and 14 genes are located 1Mb away from any of the previously reported breast cancer risk variants. Pathways analyses including 221 putative risk genes identified multiple signaling pathways that may play a significant role in the development of breast cancer. Our study provides a comprehensive understanding of and new biological insights into the genetics of this common malignancy.

Integrating transcriptomics, metabolomics, and GWAS helps reveal molecular mechanisms for metabolite levels and disease risk.

Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.

Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study.

BACKGROUND: Male-pattern baldness (MPB) is the most common cause of hair loss in men. It can be categorized into three types: type 2 (T2), type 3 (T3), and type 4 (T4), with type 1 (T1) being considered normal. Although various MPB-associated genetic variants have been suggested, a comprehensive study for linking these variants to gene expression regulation has not been performed to the best of our knowledge. RESULTS: In this study, we prioritized MPB-related tissue panels using tissue-specific enrichment analysis and utilized single-tissue panels from genotype-tissue expression version 8, as well as cross-tissue panels from context-specific genetics. Through a transcriptome-wide association study and colocalization analysis, we identified 52, 75, and 144 MPB associations for T2, T3, and T4, respectively. To assess the causality of MPB genes, we performed a conditional and joint analysis, which revealed 10, 11, and 54 putative causality genes for T2, T3, and T4, respectively. Finally, we conducted drug repositioning and identified potential drug candidates that are connected to MPB-associated genes. CONCLUSIONS: Overall, through an integrative analysis of gene expression and genotype data, we have identified robust MPB susceptibility genes that may help uncover the underlying molecular mechanisms and the novel drug candidates that may alleviate MPB.

Susceptibility gene identification and risk evaluation model construction by transcriptome-wide association analysis for salt sensitivity of blood pressure.

BACKGROUND: Salt sensitivity of blood pressure (SSBP) is an intermediate phenotype of hypertension and is a predictor of long-term cardiovascular events and death. However, the genetic structures of SSBP are uncertain, and it is difficult to precisely diagnose SSBP in population. So, we aimed to identify genes related to susceptibility to the SSBP, construct a risk evaluation model, and explore the potential functions of these genes. METHODS AND RESULTS: A genome-wide association study of the systemic epidemiology of salt sensitivity (EpiSS) cohort was performed to obtain summary statistics for SSBP. Then, we conducted a transcriptome-wide association study (TWAS) of 12 tissues using FUSION software to predict the genes associated with SSBP and verified the genes with an mRNA microarray. The potential roles of the genes were explored. Risk evaluation models of SSBP were constructed based on the serial P value thresholds of polygenetic risk scores (PRSs), polygenic transcriptome risk scores (PTRSs) and their combinations of the identified genes and genetic variants from the TWAS. The TWAS revealed that 2605 genes were significantly associated with SSBP. Among these genes, 69 were differentially expressed according to the microarray analysis. The functional analysis showed that the genes identified in the TWAS were enriched in metabolic process pathways. The PRSs were correlated with PTRSs in the heart atrial appendage, adrenal gland, EBV-transformed lymphocytes, pituitary, artery coronary, artery tibial and whole blood. Multiple logistic regression models revealed that a PRS of P < 0.05 had the best predictive ability compared with other PRSs and PTRSs. The combinations of PRSs and PTRSs did not significantly increase the prediction accuracy of SSBP in the training and validation datasets. CONCLUSIONS: Several known and novel susceptibility genes for SSBP were identified via multitissue TWAS analysis. The risk evaluation model constructed with the PRS of susceptibility genes showed better diagnostic performance than the transcript levels, which could be applied to screen for SSBP high-risk individuals.

Brain-trait-associated variants impact cell-type-specific gene regulation during neurogenesis.

Interpretation of the function of non-coding risk loci for neuropsychiatric disorders and brain-relevant traits via gene expression and alternative splicing quantitative trait locus (e/sQTL) analyses is generally performed in bulk post-mortem adult tissue. However, genetic risk loci are enriched in regulatory elements active during neocortical differentiation, and regulatory effects of risk variants may be masked by heterogeneity in bulk tissue. Here, we map e/sQTLs, and allele-specific expression in cultured cells representing two major developmental stages, primary human neural progenitors (n = 85) and their sorted neuronal progeny (n = 74), identifying numerous loci not detected in either bulk developing cortical wall or adult cortex. Using colocalization and genetic imputation via transcriptome-wide association, we uncover cell-type-specific regulatory mechanisms underlying risk for brain-relevant traits that are active during neocortical differentiation. Specifically, we identified a progenitor-specific eQTL for CENPW co-localized with common variant associations for cortical surface area and educational attainment.

Multi-tissue transcriptome-wide association study reveals susceptibility genes and drug targets for insulin resistance-relevant phenotypes.

AIM: Genome-wide association studies (GWAS) have identified multiple susceptibility loci associated with insulin resistance (IR)-relevant phenotypes. However, the genes responsible for these associations remain largely unknown. We aim to identify susceptibility genes for IR-relevant phenotypes via a transcriptome-wide association study. MATERIALS AND METHODS: We conducted a large-scale multi-tissue transcriptome-wide association study for IR (Insulin Sensitivity Index, homeostasis model assessment-IR, fasting insulin) and lipid-relevant traits (high-density lipoprotein cholesterol, triglycerides, low-density lipoprotein cholesterol and total cholesterol) using the largest GWAS summary statistics and precomputed gene expression weights of 49 human tissues. Conditional and joint analyses were implemented to identify significantly independent genes. Furthermore, we estimated the causal effects of independent genes by Mendelian randomization causal inference analysis. RESULTS: We identified 1190 susceptibility genes causally associated with IR-relevant phenotypes, including 58 genes that were not implicated in the original GWAS. Among them, 11 genes were further supported in differential expression analyses or a gene knockout mice database, such as KRIT1 showed both significantly differential expression and IR-related phenotypic effects in knockout mice. Meanwhile, seven proteins encoded by susceptibility genes were targeted by clinically approved drugs, and three of these genes (H6PD, CACNB2 and DRD2) have been served as drug targets for IR-related diseases/traits. Moreover, drug repurposing analysis identified four compounds with profiles opposing the expression of genes associated with IR risk. CONCLUSIONS: Our study provided new insights into IR aetiology and avenues for therapeutic development.

RNA profiling of human dorsal root ganglia reveals sex differences in mechanisms promoting neuropathic pain.

Neuropathic pain is a leading cause of high-impact pain, is often disabling and is poorly managed by current therapeutics. Here we focused on a unique group of neuropathic pain patients undergoing thoracic vertebrectomy where the dorsal root ganglia is removed as part of the surgery allowing for molecular characterization and identification of mechanistic drivers of neuropathic pain independently of preclinical models. Our goal was to quantify whole transcriptome RNA abundances using RNA-seq in pain-associated human dorsal root ganglia from these patients, allowing comprehensive identification of molecular changes in these samples by contrasting them with non-pain-associated dorsal root ganglia. We sequenced 70 human dorsal root ganglia, and among these 50 met inclusion criteria for sufficient neuronal mRNA signal for downstream analysis. Our expression analysis revealed profound sex differences in differentially expressed genes including increase of IL1B, TNF, CXCL14 and OSM in male and CCL1, CCL21, PENK and TLR3 in female dorsal root ganglia associated with neuropathic pain. Coexpression modules revealed enrichment in members of JUN-FOS signalling in males and centromere protein coding genes in females. Neuro-immune signalling pathways revealed distinct cytokine signalling pathways associated with neuropathic pain in males (OSM, LIF, SOCS1) and females (CCL1, CCL19, CCL21). We validated cellular expression profiles of a subset of these findings using RNAscope in situ hybridization. Our findings give direct support for sex differences in underlying mechanisms of neuropathic pain in patient populations.

Transcriptome-wide association study-derived genes as potential visceral adipose tissue-specific targets for type 2 diabetes.

AIMS/HYPOTHESIS: This study aimed to assess the causal relationship between visceral obesity and type 2 diabetes and subsequently to screen visceral adipose tissue (VAT)-specific targets for type 2 diabetes. METHODS: We examined the causal relationship between VAT and type 2 diabetes using bidirectional Mendelian randomisation (MR) followed by multivariable MR. We conducted a transcriptome-wide association study (TWAS) leveraging prediction models and a large-scale type 2 diabetes genome-wide association study (74,124 cases and 824,006 controls) to identify candidate genes in VAT and used summary-data-based MR (SMR) and co-localisation analysis to map causal genes. We performed enrichment and single-cell RNA-seq analyses to determine the cell-specific localisation of the TWAS-identified genes. We also conducted knockdown experiments in 3T3-L1 pre-adipocytes. RESULTS: MR analyses showed a causal relationship between genetically increased VAT mass and type 2 diabetes (inverse-variance weighted OR 2.48 [95% CI 2.21, 2.79]). Ten VAT-specific candidate genes were associated with type 2 diabetes after Bonferroni correction, including five causal genes supported by SMR and co-localisation: PABPC4 (1p34.3); CCNE2 (8q22.1); HAUS6 (9p22.1); CWF19L1 (10q24.31); and CCDC92 (12q24.31). Combined with enrichment analyses, clarifying cell-type specificity with single-cell RNA-seq data indicated that most TWAS-identified candidate genes appear more likely to be associated with adipocytes in VAT. Knockdown experiments suggested that Pabpc4 likely contributes to regulating differentiation and energy metabolism in 3T3-L1 adipocytes. CONCLUSIONS/INTERPRETATION: Our findings provide new insights into the genetic basis and biological processes of the association between VAT accumulation and type 2 diabetes and warrant investigation through further functional studies to validate these VAT-specific candidate genes.

A splicing transcriptome-wide association study identifies novel altered splicing for Alzheimer's disease susceptibility.

Alzheimer's disease (AD) is a common neurodegenerative disease in aging individuals. Alternative splicing is reported to be relevant to AD development while their roles in etiology of AD remain largely elusive. We performed a comprehensive splicing transcriptome-wide association study (spTWAS) using intronic excision expression genetic prediction models of 12 brain tissues developed through three modelling strategies, to identify candidate susceptibility splicing introns for AD risk. A total of 111,326 (46,828 proxy) cases and 677,663 controls of European ancestry were studied. We identified 343 associations of 233 splicing introns (143 genes) with AD risk after Bonferroni correction (0.05/136,884 = 3.65 × 10-7). Fine-mapping analyses supported 155 likely causal associations corresponding to 83 splicing introns of 55 genes. Eighteen causal splicing introns of 15 novel genes (EIF2D, WDR33, SAP130, BYSL, EPHB6, MRPL43, VEGFB, PPP1R13B, TLN2, CLUHP3, LRRC37A4P, CRHR1, LINC02210, ZNF45-AS1, and XPNPEP3) were identified for the first time to be related to AD susceptibility. Our study identified novel genes and splicing introns associated with AD risk, which can improve our understanding of the etiology of AD.

Bayesian Genome-wide TWAS Method to Leverage both cis- and trans-eQTL Information through Summary Statistics.

Transcriptome-wide association studies (TWASs) have been widely used to integrate gene expression and genetic data for studying complex traits. Due to the computational burden, existing TWAS methods do not assess distant trans-expression quantitative trait loci (eQTL) that are known to explain important expression variation for most genes. We propose a Bayesian genome-wide TWAS (BGW-TWAS) method that leverages both cis- and trans-eQTL information for a TWAS. Our BGW-TWAS method is based on Bayesian variable selection regression, which not only accounts for cis- and trans-eQTL of the target gene but also enables efficient computation by using summary statistics from standard eQTL analyses. Our simulation studies illustrated that BGW-TWASs achieved higher power compared to existing TWAS methods that do not assess trans-eQTL information. We further applied BWG-TWAS to individual-level GWAS data (N = ∼3.3K), which identified significant associations between the genetically regulated gene expression (GReX) of ZC3H12B and Alzheimer dementia (AD) (p value = 5.42 × 10-13), neurofibrillary tangle density (p value = 1.89 × 10-6), and global measure of AD pathology (p value = 9.59 × 10-7). These associations for ZC3H12B were completely driven by trans-eQTL. Additionally, the GReX of KCTD12 was found to be significantly associated with ß-amyloid (p value = 3.44 × 10-8) which was driven by both cis- and trans-eQTL. Four of the top driven trans-eQTL of ZC3H12B are located within APOC1, a known major risk gene of AD and blood lipids. Additionally, by applying BGW-TWAS with summary-level GWAS data of AD (N = ∼54K), we identified 13 significant genes including known GWAS risk genes HLA-DRB1 and APOC1, as well as ZC3H12B.

Identifying Candidate Genes Associated with Sporadic Amyotrophic Lateral Sclerosis via Integrative Analysis of Transcriptome-Wide Association Study and Messenger RNA Expression Profile.

Amyotrophic lateral sclerosis, a fatal neurodegeneration disease affecting motor neurons in the brain and spinal cord, is difficult to diagnose and treat. The objective of this study is to identify novel candidate genes related to ALS. Transcriptome-wide association study of ALS was conducted by integrating the genome-wide association study summary data (including 1234 ALS patients and 2850 controls) and pre-computed gene expression weights of different tissues. The ALS-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the mRNA expression profiles of the sporadic ALS. Functional enrichment and annotation analysis of identified genes were performed by an R package and the functional mapping and annotation software. TWAS identified 761 significant genes (PTWAS < 0.05), 627 Gene ontology terms, and 8 Kyoto Encyclopedia of Genes and Genomes pathways for ALS, such as C9orf72, with three expression quantitative trait loci were found significantly: rs2453554 (PTWAS CBRS = 4.68 × 10-10, PTWAS CBRS = 2.54 × 10-9), rs10967976 (PTWAS CBRS = 7.85 × 10-10, PTWAS CBRS = 8.91 × 10-9, PTWAS CBRS = 1.49 × 10-7, PTWAS CBRS = 5.59 × 10-7), rs3849946 (PTWAS CBRS = 7.69 × 10-4, PTWAS YBL = 4.02 × 10-2), Mitochondrion (Padj = 4.22 × 10-16), and Cell cycle (Padj = 2.04 × 10-3). Moreover, 107 common genes, 4 KEGG pathways and 41 GO terms were detected by integrating mRNA expression profiles of sALS, such as CPVL (FC = 2.06, PmRNA = 6.99 × 10-6, PTWAS CBR = 2.88 × 10-2, PTWAS CBR = 4.37 × 10-2), Pyrimidine Metabolism (Padj = 2.43 × 10-2), and Cell Activation (Padj = 5.54 × 10-3). Multiple candidate genes and pathways were detected for ALS. Our findings may provide novel clues for understanding the genetic mechanism of ALS.

A Transcriptome-Wide Analysis of Psoriasis: Identifying the Potential Causal Genes and Drug Candidates.

Psoriasis is a chronic inflammatory skin disease characterized by cutaneous eruptions and pruritus. Because the genetic backgrounds of psoriasis are only partially revealed, an integrative and rigorous study is necessary. We conducted a transcriptome-wide association study (TWAS) with the new Genotype-Tissue Expression version 8 reference panels, including some tissue and multi-tissue panels that were not used previously. We performed tissue-specific heritability analyses on genome-wide association study data to prioritize the tissue panels for TWAS analysis. TWAS and colocalization (COLOC) analyses were performed with eight tissues from the single-tissue panels and the multi-tissue panels of context-specific genetics (CONTENT) to increase tissue specificity and statistical power. From TWAS, we identified the significant associations of 101 genes in the single-tissue panels and 64 genes in the multi-tissue panels, of which 26 genes were replicated in the COLOC. Functional annotation and network analyses identified that the genes were associated with psoriasis and/or immune responses. We also suggested drug candidates that interact with jointly significant genes through a conditional and joint analysis. Together, our findings may contribute to revealing the underlying genetic mechanisms and provide new insights into treatments for psoriasis.

Identifying causal genes for migraine by integrating the proteome and transcriptome.

BACKGROUND: While previous genome-wide association studies (GWAS) have identified multiple risk variants for migraine, there is a lack of evidence about how these variants contribute to the development of migraine. We employed an integrative pipeline to efficiently transform genetic associations to identify causal genes for migraine. METHODS: We conducted a proteome-wide association study (PWAS) by combining data from the migraine GWAS data with proteomic data from the human brain and plasma to identify proteins that may play a role in the risk of developing migraine. We also combined data from GWAS of migraine with a novel joint-tissue imputation (JTI) prediction model of 17 migraine-related human tissues to conduct transcriptome-wide association studies (TWAS) together with the fine mapping method FOCUS to identify disease-associated genes. RESULTS: We identified 13 genes in the human brain and plasma proteome that modulate migraine risk by regulating protein abundance. In addition, 62 associated genes not reported in previous migraine TWAS studies were identified by our analysis of migraine using TWAS and fine mapping. Five genes including ICA1L, TREX1, STAT6, UFL1, and B3GNT8 showed significant associations with migraine at both the proteome and transcriptome, these genes are mainly expressed in ependymal cells, neurons, and glial cells, and are potential target genes for prevention of neuronal signaling and inflammatory responses in the pathogenesis of migraine. CONCLUSIONS: Our proteomic and transcriptome findings have identified disease-associated genes that may give new insights into the pathogenesis and potential therapeutic targets for migraine.

Genome-wide and transcriptome-wide association studies of mammographic density phenotypes reveal novel loci.

BACKGROUND: Mammographic density (MD) phenotypes, including percent density (PMD), area of dense tissue (DA), and area of non-dense tissue (NDA), are associated with breast cancer risk. Twin studies suggest that MD phenotypes are highly heritable. However, only a small proportion of their variance is explained by identified genetic variants. METHODS: We conducted a genome-wide association study, as well as a transcriptome-wide association study (TWAS), of age- and BMI-adjusted DA, NDA, and PMD in up to 27,900 European-ancestry women from the MODE/BCAC consortia. RESULTS: We identified 28 genome-wide significant loci for MD phenotypes, including nine novel signals (5q11.2, 5q14.1, 5q31.1, 5q33.3, 5q35.1, 7p11.2, 8q24.13, 12p11.2, 16q12.2). Further, 45% of all known breast cancer SNPs were associated with at least one MD phenotype at p < 0.05. TWAS further identified two novel genes (SHOX2 and CRISPLD2) whose genetically predicted expression was significantly associated with MD phenotypes. CONCLUSIONS: Our findings provided novel insight into the genetic background of MD phenotypes, and further demonstrated their shared genetic basis with breast cancer.

A transcriptome-wide association study identifies novel candidate susceptibility genes for prostate cancer risk.

A large proportion of heritability for prostate cancer risk remains unknown. Transcriptome-wide association study combined with validation comparing overall levels will help to identify candidate genes potentially playing a role in prostate cancer development. Using data from the Genotype-Tissue Expression Project, we built genetic models to predict normal prostate tissue gene expression using the statistical framework PrediXcan, a modified version of the unified test for molecular signatures and Joint-Tissue Imputation. We applied these prediction models to the genetic data of 79 194 prostate cancer cases and 61 112 controls to investigate the associations of genetically determined gene expression with prostate cancer risk. Focusing on associated genes, we compared their expression in prostate tumor vs normal prostate tissue, compared methylation of CpG sites located at these loci in prostate tumor vs normal tissue, and assessed the correlations between the differentiated genes' expression and the methylation of corresponding CpG sites, by analyzing The Cancer Genome Atlas (TCGA) data. We identified 573 genes showing an association with prostate cancer risk at a false discovery rate (FDR) ≤ 0.05, including 451 novel genes and 122 previously reported genes. Of the 573 genes, 152 showed differential expression in prostate tumor vs normal tissue samples. At loci of 57 genes, 151 CpG sites showed differential methylation in prostate tumor vs normal tissue samples. Of these, 20 CpG sites were correlated with expression of 11 corresponding genes. In this TWAS, we identified novel candidate susceptibility genes for prostate cancer risk, providing new insights into prostate cancer genetics and biology.

eQTLs play critical roles in regulating gene expression and identifying key regulators in rice.

The regulation of gene expression plays an essential role in both the phenotype and adaptation of plants. Transcriptome sequencing enables simultaneous identification of exonic variants and quantification of gene expression. Here, we sequenced the leaf transcriptomes of 287 rice accessions from around the world and obtained a total of 177 853 high-quality single nucleotide polymorphisms after filtering. Genome-wide association study identified 44 354 expression quantitative trait loci (eQTLs), which regulate the expression of 13 201 genes, as well as 17 local eQTL hotspots and 96 distant eQTL hotspots. Furthermore, a transcriptome-wide association study screened 21 candidate genes for starch content in the flag leaves at the heading stage. HS002 was identified as a significant distant eQTL hotspot with five downstream genes enriched for diterpene antitoxin synthesis. Co-expression analysis, eQTL analysis, and linkage mapping together demonstrated that bHLH026 acts as a key regulator to activate the expression of downstream genes. The transgenic assay revealed that bHLH026 is an important regulator of diterpenoid antitoxin synthesis and enhances the disease resistance of rice. These findings improve our knowledge of the regulatory mechanisms of gene expression variation and complex regulatory networks of the rice genome and will facilitate genetic improvement of cultivated rice varieties.