Genome-wide association analyses of breast cancer in women of African ancestry identify new susceptibility loci and improve risk prediction.

We performed genome-wide association studies of breast cancer including 18,034 cases and 22,104 controls of African ancestry. Genetic variants at 12 loci were associated with breast cancer risk (P < 5 × 10-8), including associations of a low-frequency missense variant rs61751053 in ARHGEF38 with overall breast cancer (odds ratio (OR) = 1.48) and a common variant rs76664032 at chromosome 2q14.2 with triple-negative breast cancer (TNBC) (OR = 1.30). Approximately 15.4% of cases with TNBC carried six risk alleles in three genome-wide association study-identified TNBC risk variants, with an OR of 4.21 (95% confidence interval = 2.66-7.03) compared with those carrying fewer than two risk alleles. A polygenic risk score (PRS) showed an area under the receiver operating characteristic curve of 0.60 for the prediction of breast cancer risk, which outperformed PRS derived using data from females of European ancestry. Our study markedly increases the population diversity in genetic studies for breast cancer and demonstrates the utility of PRS for risk prediction in females of African ancestry.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes.

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

Leveraging generative AI to prioritize drug repurposing candidates for Alzheimer's disease with real-world clinical validation.

Drug repurposing represents an attractive alternative to the costly and time-consuming process of new drug development, particularly for serious, widespread conditions with limited effective treatments, such as Alzheimer's disease (AD). Emerging generative artificial intelligence (GAI) technologies like ChatGPT offer the promise of expediting the review and summary of scientific knowledge. To examine the feasibility of using GAI for identifying drug repurposing candidates, we iteratively tasked ChatGPT with proposing the twenty most promising drugs for repurposing in AD, and tested the top ten for risk of incident AD in exposed and unexposed individuals over age 65 in two large clinical datasets: (1) Vanderbilt University Medical Center and (2) the All of Us Research Program. Among the candidates suggested by ChatGPT, metformin, simvastatin, and losartan were associated with lower AD risk in meta-analysis. These findings suggest GAI technologies can assimilate scientific insights from an extensive Internet-based search space, helping to prioritize drug repurposing candidates and facilitate the treatment of diseases.

Novel breast cancer susceptibility loci under linkage peaks identified in African ancestry consortia.

BACKGROUND: Expansion of genome-wide association studies across population groups is needed to improve our understanding of shared and unique genetic contributions to breast cancer. We performed association and replication studies guided by a priori linkage findings from African ancestry (AA) relative pairs. METHODS: We performed fixed-effect inverse-variance weighted meta-analysis under three significant AA breast cancer linkage peaks (3q26-27, 12q22-23, and 16q21-22) in 9241 AA cases and 10 193 AA controls. We examined associations with overall breast cancer as well as estrogen receptor (ER)-positive and negative subtypes (193,132 SNPs). We replicated associations in the African-ancestry Breast Cancer Genetic Consortium (AABCG). RESULTS: In AA women, we identified two associations on chr12q for overall breast cancer (rs1420647, OR = 1.15, p = 2.50×10-6; rs12322371, OR = 1.14, p = 3.15×10-6), and one for ER-negative breast cancer (rs77006600, OR = 1.67, p = 3.51×10-6). On chr3, we identified two associations with ER-negative disease (rs184090918, OR = 3.70, p = 1.23×10-5; rs76959804, OR = 3.57, p = 1.77×10-5) and on chr16q we identified an association with ER-negative disease (rs34147411, OR = 1.62, p = 8.82×10-6). In the replication study, the chr3 associations were significant and effect sizes were larger (rs184090918, OR: 6.66, 95% CI: 1.43, 31.01; rs76959804, OR: 5.24, 95% CI: 1.70, 16.16). CONCLUSION: The two chr3 SNPs are upstream to open chromatin ENSR00000710716, a regulatory feature that is actively regulated in mammary tissues, providing evidence that variants in this chr3 region may have a regulatory role in our target organ. Our study provides support for breast cancer variant discovery using prioritization based on linkage evidence.

Prostein expression on circulating tumor cells as a prognostic marker in metastatic castration-resistant prostate cancer.

BACKGROUND: Identification of emerging molecular biomarkers on circulating tumor cells (CTCs) represents an attractive feature of liquid biopsy that facilitates precision and tailored medicine in the management of metastatic castration-resistant prostate cancer (mCRPC). Prostein is an androgen-regulated transmembrane protein with high prostate specificity. Prostein-positive circulating tumor cell (CTC) was recently suggested to have diagnostic potential; however, no study has been conducted to evaluate its prognostic value in mCRPC. METHODS: CTCs from mCRPC patients were enumerated using the CellSearch System. Prostein-positive CTCs were identified by immunostaining results. The relationships between prostein expression on CTCs and PSA response rate, PSA progression-free survival (PSA-PFS), radiographic progression-free survival (PFS), and overall survival (OS) were tested by Fisher's exact test or evaluated using Kaplan-Meier and multivariate Cox analyses. RESULTS: Prostein-positive CTCs were identified in 31 of 87 baseline samples from mCRPC patients and 16 of 51 samples collected at the first follow-up visit. PSA response rates were significantly lower in baseline prostein-positive patients (0%, 0/31) than in prostein-negative patients (19.6%, 11/56) (p = 0.007). The 31 prostein-positive patients had significantly shorter PSA-PFS (p < 0.001), radiographic PFS (p < 0.001), and OS (p = 0.018), compared to the 56 prostein-negative patients at baseline. The association with PSA-PFS maintained its significance (p = 0.028) in multivariate analyses. Analyzing prostein expression at the first follow-up as well as the conversion of prostein expression from baseline to follow-up samples not only confirmed the association with PSA-PFS, but also demonstrated prognostic significance with OS. CONCLUSION: Our study provides the first evidence to support the potential of prostein expression on CTCs to serve as a novel prognostic marker in mCRPC patients. Future large-scale prospective studies are needed to validate our findings.

Molecular signature incorporating the immune microenvironment enhances thyroid cancer outcome prediction.

Genomic and transcriptomic analysis has furthered our understanding of many tumors. Yet, thyroid cancer management is largely guided by staging and histology, with few molecular prognostic and treatment biomarkers. Here, we utilize a large cohort of 251 patients with 312 samples from two tertiary medical centers and perform DNA/RNA sequencing, spatial transcriptomics, and multiplex immunofluorescence to identify biomarkers of aggressive thyroid malignancy. We identify high-risk mutations and discover a unique molecular signature of aggressive disease, the Molecular Aggression and Prediction (MAP) score, which provides improved prognostication over high-risk mutations alone. The MAP score is enriched for genes involved in epithelial de-differentiation, cellular division, and the tumor microenvironment. The MAP score also identifies aggressive tumors with lymphocyte-rich stroma that may benefit from immunotherapy. Future clinical profiling of the stromal microenvironment of thyroid cancer could improve prognostication, inform immunotherapy, and support development of novel therapeutics for thyroid cancer and other stroma-rich tumors.

The impact of common variants on gene expression in the human brain: from RNA to protein to schizophrenia risk.

Background: The impact of genetic variants on gene expression has been intensely studied at the transcription level, yielding in valuable insights into the association between genes and the risk of complex disorders, such as schizophrenia (SCZ). However, the downstream impact of these variants and the molecular mechanisms connecting transcription variation to disease risk are not well understood. Results: We quantitated ribosome occupancy in prefrontal cortex samples of the BrainGVEX cohort. Together with transcriptomics and proteomics data from the same cohort, we performed cis-Quantitative Trait Locus (QTL) mapping and identified 3,253 expression QTLs (eQTLs), 1,344 ribosome occupancy QTLs (rQTLs), and 657 protein QTLs (pQTLs) out of 7,458 genes quantitated in all three omics types from 185 samples. Of the eQTLs identified, only 34% have their effects propagated to the protein level. Further analysis on the effect size of prefrontal cortex eQTLs identified from an independent dataset showed clear post-transcriptional attenuation of eQTL effects. To investigate the biological relevance of the attenuated eQTLs, we identified 70 expression-specific QTLs (esQTLs), 51 ribosome-occupancy-specific QTLs (rsQTLs), and 107 protein-specific QTLs (psQTLs). Five of these omics-specific QTLs showed strong colocalization with SCZ GWAS signals, three of them are esQTLs. The limited number of GWAS colocalization discoveries from omics-specific QTLs and the apparent prevalence of eQTL attenuation prompted us to take a complementary approach to investigate the functional relevance of attenuated eQTLs. Using S-PrediXcan we identified 74 SCZ risk genes, 34% of which were novel, and 67% of these risk genes were replicated in a MR-Egger test. Notably, 52 out of 74 risk genes were identified using eQTL data and 70% of these SCZ-risk-gene-driving eQTLs show little to no evidence of driving corresponding variations at the protein level. Conclusion: The effect of eQTLs on gene expression in the prefrontal cortex is commonly attenuated post-transcriptionally. Many of the attenuated eQTLs still correlate with SCZ GWAS signal. Further investigation is needed to elucidate a mechanistic link between attenuated eQTLs and SCZ disease risk.

Drug repurposing for Alzheimer's disease from 2012-2022-a 10-year literature review.

Background: Alzheimer's disease (AD) is a debilitating neurodegenerative condition with few treatment options available. Drug repurposing studies have sought to identify existing drugs that could be repositioned to treat AD; however, the effectiveness of drug repurposing for AD remains unclear. This review systematically analyzes the progress made in drug repurposing for AD throughout the last decade, summarizing the suggested drug candidates and analyzing changes in the repurposing strategies used over time. We also examine the different types of data that have been leveraged to validate suggested drug repurposing candidates for AD, which to our knowledge has not been previous investigated, although this information may be especially useful in appraising the potential of suggested drug repurposing candidates. We ultimately hope to gain insight into the suggested drugs representing the most promising repurposing candidates for AD. Methods: We queried the PubMed database for AD drug repurposing studies published between 2012 and 2022. 124 articles were reviewed. We used RxNorm to standardize drug names across the reviewed studies, map drugs to their constituent ingredients, and identify prescribable drugs. We used the Anatomical Therapeutic Chemical (ATC) Classification System to group drugs. Results: 573 unique drugs were proposed for repurposing in AD over the last 10 years. These suggested repurposing candidates included drugs acting on the nervous system (17%), antineoplastic and immunomodulating agents (16%), and drugs acting on the cardiovascular system (12%). Clozapine, a second-generation antipsychotic medication, was the most frequently suggested repurposing candidate (N = 6). 61% (76/124) of the reviewed studies performed a validation, yet only 4% (5/124) used real-world data for validation. Conclusion: A large number of potential drug repurposing candidates for AD has accumulated over the last decade. However, among these drugs, no single drug has emerged as the top candidate, making it difficult to establish research priorities. Validation of drug repurposing hypotheses is inconsistently performed, and real-world data has been critically underutilized for validation. Given the urgent need for new AD therapies, the utility of real-world data in accelerating identification of high-priority candidates for AD repurposing warrants further investigation.

Leveraging Generative AI to Prioritize Drug Repurposing Candidates: Validating Identified Candidates for Alzheimer's Disease in Real-World Clinical Datasets.

Drug repurposing represents an attractive alternative to the costly and time-consuming process of new drug development, particularly for serious, widespread conditions with limited effective treatments, such as Alzheimer's disease (AD). Emerging generative artificial intelligence (GAI) technologies like ChatGPT offer the promise of expediting the review and summary of scientific knowledge. To examine the feasibility of using GAI for identifying drug repurposing candidates, we iteratively tasked ChatGPT with proposing the twenty most promising drugs for repurposing in AD, and tested the top ten for risk of incident AD in exposed and unexposed individuals over age 65 in two large clinical datasets: 1) Vanderbilt University Medical Center and 2) the All of Us Research Program. Among the candidates suggested by ChatGPT, metformin, simvastatin, and losartan were associated with lower AD risk in meta-analysis. These findings suggest GAI technologies can assimilate scientific insights from an extensive Internet-based search space, helping to prioritize drug repurposing candidates and facilitate the treatment of diseases.

Leveraging Generative AI to Prioritize Drug Repurposing Candidates: Validating Identified Candidates for Alzheimer's Disease in Real-World Clinical Datasets.

Drug repurposing represents an attractive alternative to the costly and time-consuming process of new drug development, particularly for serious, widespread conditions with limited effective treatments, such as Alzheimer's disease (AD). Emerging generative artificial intelligence (GAI) technologies like ChatGPT offer the promise of expediting the review and summary of scientific knowledge. To examine the feasibility of using GAI for identifying drug repurposing candidates, we iteratively tasked ChatGPT with proposing the twenty most promising drugs for repurposing in AD, and tested the top ten for risk of incident AD in exposed and unexposed individuals over age 65 in two large clinical datasets: 1) Vanderbilt University Medical Center and 2) the All of Us Research Program. Among the candidates suggested by ChatGPT, metformin, simvastatin, and losartan were associated with lower AD risk in meta-analysis. These findings suggest GAI technologies can assimilate scientific insights from an extensive Internet-based search space, helping to prioritize drug repurposing candidates and facilitate the treatment of diseases.

Multi-ancestry and multi-trait genome-wide association meta-analyses inform clinical risk prediction for systemic lupus erythematosus.

Systemic lupus erythematosus is a heritable autoimmune disease that predominantly affects young women. To improve our understanding of genetic etiology, we conduct multi-ancestry and multi-trait meta-analysis of genome-wide association studies, encompassing 12 systemic lupus erythematosus cohorts from 3 different ancestries and 10 genetically correlated autoimmune diseases, and identify 16 novel loci. We also perform transcriptome-wide association studies, computational drug repurposing analysis, and cell type enrichment analysis. We discover putative drug classes, including a histone deacetylase inhibitor that could be repurposed to treat lupus. We also identify multiple cell types enriched with putative target genes, such as non-classical monocytes and B cells, which may be targeted for future therapeutics. Using this newly assembled result, we further construct polygenic risk score models and demonstrate that integrating polygenic risk score with clinical lab biomarkers improves the diagnostic accuracy of systemic lupus erythematosus using the Vanderbilt BioVU and Michigan Genomics Initiative biobanks.

FAVOR: functional annotation of variants online resource and annotator for variation across the human genome.

Large biobank-scale whole genome sequencing (WGS) studies are rapidly identifying a multitude of coding and non-coding variants. They provide an unprecedented resource for illuminating the genetic basis of human diseases. Variant functional annotations play a critical role in WGS analysis, result interpretation, and prioritization of disease- or trait-associated causal variants. Existing functional annotation databases have limited scope to perform online queries and functionally annotate the genotype data of large biobank-scale WGS studies. We develop the Functional Annotation of Variants Online Resources (FAVOR) to meet these pressing needs. FAVOR provides a comprehensive multi-faceted variant functional annotation online portal that summarizes and visualizes findings of all possible nine billion single nucleotide variants (SNVs) across the genome. It allows for rapid variant-, gene- and region-level queries of variant functional annotations. FAVOR integrates variant functional information from multiple sources to describe the functional characteristics of variants and facilitates prioritizing plausible causal variants influencing human phenotypes. Furthermore, we provide a scalable annotation tool, FAVORannotator, to functionally annotate large-scale WGS studies and efficiently store the genotype and their variant functional annotation data in a single file using the annotated Genomic Data Structure (aGDS) format, making downstream analysis more convenient. FAVOR and FAVORannotator are available at https://favor.genohub.org.

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.

TVAR: assessing tissue-specific functional effects of non-coding variants with deep learning.

MOTIVATION: Analysis of whole-genome sequencing (WGS) for genetics is still a challenge due to the lack of accurate functional annotation of non-coding variants, especially the rare ones. As eQTLs have been extensively implicated in the genetics of human diseases, we hypothesize that rare non-coding variants discovered in WGS play a regulatory role in predisposing disease risk. RESULTS: With thousands of tissue- and cell-type-specific epigenomic features, we propose TVAR. This multi-label learning-based deep neural network predicts the functionality of non-coding variants in the genome based on eQTLs across 49 human tissues in the GTEx project. TVAR learns the relationships between high-dimensional epigenomics and eQTLs across tissues, taking the correlation among tissues into account to understand shared and tissue-specific eQTL effects. As a result, TVAR outputs tissue-specific annotations, with an average AUROC of 0.77 across these tissues. We evaluate TVAR's performance on four complex diseases (coronary artery disease, breast cancer, Type 2 diabetes and Schizophrenia), using TVAR's tissue-specific annotations, and observe its superior performance in predicting functional variants for both common and rare variants, compared with five existing state-of-the-art tools. We further evaluate TVAR's G-score, a scoring scheme across all tissues, on ClinVar, fine-mapped GWAS loci, Massive Parallel Reporter Assay (MPRA) validated variants and observe the consistently better performance of TVAR compared with other competing tools. AVAILABILITY AND IMPLEMENTATION: The TVAR source code and its scores on the ClinVar catalog, fine mapped GWAS Loci, high confidence eQTLs from GTEx dataset, and MPRA validated functional variants are available at https://github.com/haiyang1986/TVAR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Risperidone-induced changes in DNA methylation in peripheral blood from first-episode schizophrenia patients parallel changes in neuroimaging and cognitive phenotypes.

BACKGROUND: Second generation antipsychotics such as risperidone are first-line pharmacotherapy treatment choices for schizophrenia. However, our ability to reliably predict and monitor treatment reaction is impeded by the lack of relevant biomarkers. As a biomarker for the susceptibility of schizophrenia and clozapine treatment response, DNA methylation (DNAm) has been studied, but the impact of antipsychotics on DNAm has not been explored in drug-naïve patients. OBJECTIVE: The aim of the present study was to examine changes of DNAm after short-term antipsychotic therapy in first-episode drug-naïve schizophrenia (FES) to identify the beneficial and adverse effect of risperidone on DNAm and their relation to treatment outcome. METHODS: Thirty-eight never treated schizophrenia patients and 38 demographically matched individuals (healthy controls) were assessed at baseline and at 8-week follow-up with symptom ratings, and cognitive and functional imaging procedures, at which time a blood draw for DNAm studies was performed. During the 8-week period, patients received treatment with risperidone monotherapy. An independent data set was used as replication. RESULTS: We identified brain related pathways enriched in 4,888 FES-associated CpG sites relative to controls. Risperidone administration in patients altered DNAm in 5,979 CpG sites relative to baseline. Significant group differences in DNAm at follow-up were seen in FES patients at 6,760 CpG sites versus healthy controls. Through comparison of effect size, we found 87.54% out of the risperidone-associated changes in DNAm showed possible beneficial effect, while only 12.46% showed potential adverse effect. There were 580 DNAm sites in which changes shifted methylation levels to be indistinguishable from controls after risperidone treatment. The DNAm changes of some sites that normalized after treatment were correlated with treatment-related changes in symptom severity, spontaneous neurophysiological activity, and cognition. We replicated our results in an independent data set. CONCLUSION: The normalizing effect of risperidone monotherapy on gene DNAm, and its correlation with clinically relevant phenotypes, indicates that risperidone therapy is associated with DNAm changes that are related to changes in brain physiology, cognition and symptom severity.

Genomic Aberrations in Circulating Tumor DNAs from Palbociclib-Treated Metastatic Breast Cancer Patients Reveal a Novel Resistance Mechanism.

Previously undescribed molecular mechanisms of resistance will emerge with the increased use of cyclin-dependent kinase 4/6 inhibitors in clinical settings. To identify genomic aberrations in circulating tumor DNA associated with treatment resistance in palbociclib-treated metastatic breast cancer (MBC) patients, we collected 35 pre- and post-treatment blood samples from 16 patients with estrogen receptor-positive (ER+) MBC, including 9 with inflammatory breast cancer (IBC). Circulating cell-free DNAs (cfDNAs) were isolated for sequencing using a targeted panel of 91 genes. Our data showed that FBXW7 and CDK6 were more frequently altered in IBC than in non-IBC, whereas conversely, PIK3CA was more frequently altered in non-IBC than in IBC. The cfDNA samples collected at follow-up harbored more mutations than baseline samples. By analyzing paired samples, we observed a higher percentage of patients with mutations in RB1, CCNE1, FBXW7, EZH2, and ARID1A, but a lower proportion of patients with mutated TSC2 at the post-treatment stage when they developed progression. Moreover, acquisition of CCNE1 mutations or loss of TSC2 mutations after treatment initiation conferred an unfavorable prognosis. These data provide insights into the relevance of novel genomic alterations in cfDNA to palbociclib resistance in MBC patients. Future large-scale prospective studies are warranted to confirm our findings.

A computational framework to unify orthogonal information in DNA methylation and copy number aberrations in cell-free DNA for early cancer detection.

Cell-free DNA (cfDNA) provides a convenient diagnosis avenue for noninvasive cancer detection. The current methods are focused on identifying circulating tumor DNA (ctDNA)s genomic aberrations, e.g. mutations, copy number aberrations (CNAs) or methylation changes. In this study, we report a new computational method that unifies two orthogonal pieces of information, namely methylation and CNAs, derived from whole-genome bisulfite sequencing (WGBS) data to quantify low tumor content in cfDNA. It implements a Bayes model to enrich ctDNA from WGBS data based on hypomethylation haplotypes, and subsequently, models CNAs for cancer detection. We generated WGBS data in a total of 262 samples, including high-depth (>20×, deduped high mapping quality reads) data in 76 samples with matched triplets (tumor, adjacent normal and cfDNA) and low-depth (~2.5×, deduped high mapping quality reads) data in 186 samples. We identified a total of 54 Mb regions of hypomethylation haplotypes for model building, a vast majority of which are not covered in the HumanMethylation450 arrays. We showed that our model is able to substantially enrich ctDNA reads (tens of folds), with clearly elevated CNAs that faithfully match the CNAs in the paired tumor samples. In the 19 hepatocellular carcinoma cfDNA samples, the estimated enrichment is as high as 16 fold, and in the simulation data, it can achieve over 30-fold enrichment for a ctDNA level of 0.5% with a sequencing depth of 600×. We also found that these hypomethylation regions are also shared among many cancer types, thus demonstrating the potential of our framework for pancancer early detection.

Integration of multidimensional splicing data and GWAS summary statistics for risk gene discovery.

A common strategy for the functional interpretation of genome-wide association study (GWAS) findings has been the integrative analysis of GWAS and expression data. Using this strategy, many association methods (e.g., PrediXcan and FUSION) have been successful in identifying trait-associated genes via mediating effects on RNA expression. However, these approaches often ignore the effects of splicing, which can carry as much disease risk as expression. Compared to expression data, one challenge to detect associations using splicing data is the large multiple testing burden due to multidimensional splicing events within genes. Here, we introduce a multidimensional splicing gene (MSG) approach, which consists of two stages: 1) we use sparse canonical correlation analysis (sCCA) to construct latent canonical vectors (CVs) by identifying sparse linear combinations of genetic variants and splicing events that are maximally correlated with each other; and 2) we test for the association between the genetically regulated splicing CVs and the trait of interest using GWAS summary statistics. Simulations show that MSG has proper type I error control and substantial power gains over existing multidimensional expression analysis methods (i.e., S-MultiXcan, UTMOST, and sCCA+ACAT) under diverse scenarios. When applied to the Genotype-Tissue Expression Project data and GWAS summary statistics of 14 complex human traits, MSG identified on average 83%, 115%, and 223% more significant genes than sCCA+ACAT, S-MultiXcan, and UTMOST, respectively. We highlight MSG's applications to Alzheimer's disease, low-density lipoprotein cholesterol, and schizophrenia, and found that the majority of MSG-identified genes would have been missed from expression-based analyses. Our results demonstrate that aggregating splicing data through MSG can improve power in identifying gene-trait associations and help better understand the genetic risk of complex traits.

Evaluating breast cancer predisposition genes in women of African ancestry.

PURPOSE: Studies conducted primarily among European ancestry women reported 12 breast cancer predisposition genes. However, etiologic roles of these genes in breast cancer among African ancestry women have been less well-investigated. METHODS: We conducted a case-control study in African American women, which included 1117 breast cancer cases and 2169 cancer-free controls, and a pooled analysis, which included 7096 cases and 8040 controls of African descent. Odds ratios of associations with breast cancer risk were estimated. RESULTS: Using sequence data, we identified 61 pathogenic variants in 12 breast cancer predisposition genes, including 11 pathogenic variants not yet reported in previous studies. Pooled analysis showed statistically significant associations of breast cancer risk with pathogenic variants in BRCA1, BRCA2, PALB2, ATM, CHEK2, TP53, NF1, RAD51C, and RAD51D (all P < .05). The associations with BRCA1, PALB2, and RAD51D were stronger for estrogen receptor (ER)-negative than for ER-positive breast cancer (P heterogeneity < .05), whereas the association with CHEK2 was stronger for ER-positive than for ER-negative breast cancer. CONCLUSION: Our study confirmed previously identified associations of breast cancer risk with BRCA1, BRCA2, PALB2, ATM, TP53, NF1, and CHEK2 and provided new evidence to extend the associations of breast cancer risk with RAD51C and RAD51D, which was identified previously in European ancestry populations, to African ancestry women.

A Bayesian framework to integrate multi-level genome-scale data for Autism risk gene prioritization.

BACKGROUND: Autism spectrum disorder (ASD) is a group of complex neurodevelopment disorders with a strong genetic basis. Large scale sequencing studies have identified over one hundred ASD risk genes. Nevertheless, the vast majority of ASD risk genes remain to be discovered, as it is estimated that more than 1000 genes are likely to be involved in ASD risk. Prioritization of risk genes is an effective strategy to increase the power of identifying novel risk genes in genetics studies of ASD. As ASD risk genes are likely to exhibit distinct properties from multiple angles, we reason that integrating multiple levels of genomic data is a powerful approach to pinpoint genuine ASD risk genes. RESULTS: We present BNScore, a Bayesian model selection framework to probabilistically prioritize ASD risk genes through explicitly integrating evidence from sequencing-identified ASD genes, biological annotations, and gene functional network. We demonstrate the validity of our approach and its improved performance over existing methods by examining the resulting top candidate ASD risk genes against sets of high-confidence benchmark genes and large-scale ASD genome-wide association studies. We assess the tissue-, cell type- and development stage-specific expression properties of top prioritized genes, and find strong expression specificity in brain tissues, striatal medium spiny neurons, and fetal developmental stages. CONCLUSIONS: In summary, we show that by integrating sequencing findings, functional annotation profiles, and gene-gene functional network, our proposed BNScore provides competitive performance compared to current state-of-the-art methods in prioritizing ASD genes. Our method offers a general and flexible strategy to risk gene prioritization that can potentially be applied to other complex traits as well.