Applicability of spatial transcriptional profiling to cancer research.

Spatial transcriptional profiling provides gene expression information within the important anatomical context of tissue architecture. This approach is well suited to characterizing solid tumors, which develop within a complex landscape of malignant cells, immune cells, and stroma. In a single assay, spatial transcriptional profiling can interrogate the role of spatial relationships among these cell populations as well as reveal spatial patterns of relevant oncogenic genetic events. The broad utility of this approach is reflected in the array of strategies that have been developed for its implementation as well as in the recent commercial development of several profiling platforms. The flexibility to apply these technologies to both hypothesis-driven and discovery-driven studies allows widespread applicability in research settings. This review discusses available technologies for spatial transcriptional profiling and several applications for their use in cancer research.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Mentoring Minority Cancer Researchers of Tomorrow: Comparison of the Face-to-Face, Virtual, and Hybrid Training Methods of the CaRE2 Summer Cancer Research Education and Training Program.

In our previous publication, we reported a framework to develop an undergraduate cancer research training program at Florida A&M University (FAMU) under the umbrella of the Florida-California Cancer Research, Education, and Engagement (CaRE2) Health Equity Center activity by harnessing the resources available at FAMU, the University of Florida (UF), and the University of Southern California (USC) Cancer Centers. The implementation of the CaRE2 face-to-face training platform was dramatically affected by the COVID-19 pandemic during the summer of 2020 and 2021 training periods. However, a concerted effort was made to restructure the face-to-face training model into virtual and hybrid training methods to maintain the continuity of the program during the pandemic. This article compared the three methods to identify the best platform for training URM students in cancer disparity research. The program's effectiveness was measured through motivation, experiences, and knowledge gained by trainees during and one year after the completion of the program. The results showed that the participants were highly positive in their feedback about the professional and academic values of the program. Although the virtual and hybrid methods experienced significant challenges during the pandemic, the hybrid training module offered an "above average" effectiveness in performance, like the face-to-face mentoring platform in mentoring URM students in cancer disparity research.

Florida-California Cancer Research, Education and Engagement (CaRE2) Health Equity Center: Structure, Innovations, and Initial Outcomes.

INTRODUCTION: The Florida-California Cancer Research, Education, and Engagement (CaRE2) Health Equity Center is a triad partnership committed to increasing institutional capacity for cancer disparity research, the diversity of the cancer workforce, and community empowerment. This article provides an overview of the structure, process innovations, and initial outcomes from the first 4 years of the CaRE2 triad partnership. METHODS: CaRE2 serves diverse populations in Florida and California using a "molecule to the community and back" model. We prioritize research on the complex intersection of biological, environmental, and social determinants health, working together with scientific and health disparities communities, sharing expertise across institutions, bidirectional training, and community outreach. Partnership progress and outcomes were assessed using mixed methods and four Program Steering Committee meetings. RESULTS: Research capacity was increased through development of a Living Repository of 81 cancer model systems from minority patients for novel cancer drug development. CaRE2 funded 15 scientific projects resulting in 38 publications. Workforce diversity entailed supporting 94 cancer trainees (92 URM) and 34 ESIs (32 URM) who coauthored 313 CaRE2-related publications and received 48 grants. Community empowerment was promoted via outreaching to more than 3000 individuals, training 145 community cancer advocates (including 28 Community Scientist Advocates), and publishing 10 community reports. CaRE2 members and trainees together have published 639 articles, received 61 grants, and 57 awards. CONCLUSION: The CaRE2 partnership has achieved its initial aims. Infrastructure for translational cancer research was expanded at one partner institution, and cancer disparities research was expanded at the two cancer centers.

Results and lessons from dual extraction of DNA and RNA from formalin-fixed paraffin-embedded breast tumor tissues for a large Cancer epidemiologic study.

BACKGROUND: The use of archived formalin-fixed paraffin-embedded (FFPE) tumor tissues has become a common practice in clinical and epidemiologic genetic research. Simultaneous extraction of DNA and RNA from FFPE tissues is appealing but can be practically challenging. Here we report our results and lessons learned from processing FFPE breast tumor tissues for a large epidemiologic study. METHODS: Qiagen AllPrep DNA/RNA FFPE kit was adapted for dual extraction using tissue punches or sections from breast tumor tissues. The yield was quantified using Qubit and fragmentation analysis by Agilent Bioanalyzer. A subset of the DNA samples were used for genome-wide DNA methylation assays and RNA samples for sequencing. The QC metrices and performance of the assays were analyzed with pre-analytical variables. RESULTS: A total of 1859 FFPE breast tumor tissues were processed. We found it critical to adjust proteinase K digestion time based on tissue volume to achieve balanced yields of DNA and RNA. Tissue punches taken from tumor-enriched regions provided the most reliable output. A median of 1475 ng DNA and 1786 ng RNA per sample was generated. The median DNA integrity number (DIN) was 3.8 and median DV200 for RNA was 33.2. Of 1294 DNA samples used in DNA methylation assays, 97% passed quality check by qPCR and 92% generated data deemed high quality. Of the 130 RNA samples with DV200 ≥ 20% used in RNA-sequencing, all but 5 generated usable transcriptomic data with a mapping rate ≥ 60%. CONCLUSIONS: Dual DNA/RNA purification using Qiagen AllPrep FFPE extraction protocol is feasible for clinical and epidemiologic studies. We recommend tissue punches as a reliable source material and fine tuning of proteinase K digestion time based on tissue volume. IMPACT: Our protocol and recommendations may be adapted by future studies for successful extraction of archived tumor tissues.

Analysis of the genomic landscapes of Barbadian and Nigerian women with triple negative breast cancer.

PURPOSE: Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype that disproportionately affects women of African ancestry (WAA) and is often associated with poor survival. Although there is a high prevalence of TNBC across West Africa and in women of the African diaspora, there has been no comprehensive genomics study to investigate the mutational profile of ancestrally related women across the Caribbean and West Africa. METHODS: This multisite cross-sectional study used 31 formalin-fixed paraffin-embedded (FFPE) samples from Barbadian and Nigerian TNBC participants. High-resolution whole exome sequencing (WES) was performed on the Barbadian and Nigerian TNBC samples to identify their mutational profiles and comparisons were made to African American, European American and Asian American sequencing data obtained from The Cancer Genome Atlas (TCGA). Whole exome sequencing was conducted on tumors with an average of 382 × coverage and 4335 × coverage for pooled germline non-tumor samples. RESULTS: Variants detected at high frequency in our WAA cohorts were found in the following genes NBPF12, PLIN4, TP53 and BRCA1. In the TCGA TNBC cases, these genes had a lower mutation rate, except for TP53 (32% in our cohort; 63% in TCGA-African American; 67% in TCGA-European American; 63% in TCGA-Asian). For all altered genes, there were no differences in frequency of mutations between WAA TNBC groups including the TCGA-African American cohort. For copy number variants, high frequency alterations were observed in PIK3CA, TP53, FGFR2 and HIF1AN genes. CONCLUSION: This study provides novel insights into the underlying genomic alterations in WAA TNBC samples and shines light on the importance of inclusion of under-represented populations in cancer genomics and biomarker studies.

Translating RNA sequencing into clinical diagnostics: opportunities and challenges.

With the emergence of RNA sequencing (RNA-seq) technologies, RNA-based biomolecules hold expanded promise for their diagnostic, prognostic and therapeutic applicability in various diseases, including cancers and infectious diseases. Detection of gene fusions and differential expression of known disease-causing transcripts by RNA-seq represent some of the most immediate opportunities. However, it is the diversity of RNA species detected through RNA-seq that holds new promise for the multi-faceted clinical applicability of RNA-based measures, including the potential of extracellular RNAs as non-invasive diagnostic indicators of disease. Ongoing efforts towards the establishment of benchmark standards, assay optimization for clinical conditions and demonstration of assay reproducibility are required to expand the clinical utility of RNA-seq.

Multiethnic PDX models predict a possible immune signature associated with TNBC of African ancestry.

PURPOSE: Triple-negative breast cancer (TNBC) is an aggressive subtype most prevalent among women of Western Sub-Saharan African ancestry. It accounts for 15-25% of African American (AA) breast cancers (BC) and up to 80% of Ghanaian breast cancers, thus contributing to outcome disparities in BC for black women. The aggressive biology of TNBC has been shown to be regulated partially by breast cancer stem cells (BCSC) which mediate tumor recurrence and metastasis and are more abundant in African breast tumors. METHODS: We studied the biological differences between TNBC in women with African ancestry and those of Caucasian women by comparing the gene expression of the BCSC. From low-passage patient derived xenografts (PDX) from Ghanaian (GH), AA, and Caucasian American (CA) TNBCs, we sorted for and sequenced the stem cell populations and analyzed for differential gene enrichment. RESULTS: In our cohort of TNBC tumors, we observed that the ALDH expressing stem cells display distinct ethnic specific gene expression patterns, with the largest difference existing between the GH and AA ALDH+ cells. Furthermore, the tumors from the women of African ancestry [GH/AA] had ALDH stem cell (SC) enrichment for expression of immune related genes and processes. Among the significantly upregulated genes were CD274 (PD-L1), CXCR9, CXCR10 and IFI27, which could serve as potential drug targets. CONCLUSIONS: Further exploration of the role of immune regulated genes and biological processes in BCSC may offer insight into developing novel approaches to treating TNBC to help ameliorate survival disparities in women with African ancestry.

Developing a Novel Framework for an Undergraduate Cancer Research Education and Engagement Program for Underrepresented Minority Students: the Florida-California CaRE2 Research Education Core (REC) Training Program.

Lack of substantive research experiences and technical skills mentoring during undergraduate studies leaves many underrepresented minority (URM) students unprepared to apply to competitive graduate programs. As a part of our ongoing effort to increase the pipeline for the development and training of successful URM scientists in biomedical sciences with focus on reducing cancer health disparities, the Florida-California Cancer Research Education and Engagement (CaRE2) Health Equity Center was launched in 2018. Funded through an NIH/NCI U54 grant mechanism, the CaRE2 Center is a triad partnership among Florida Agricultural and Mechanical University (FAMU), a minority-serving institution, University of Florida (UF), and University of Southern California (USC) Cancer Center. One of the objectives of the triad partnership is to promote the coordination and implementation of the training of the next generation of Black and Latinx biomedical scientists in Florida and California. An important component of the CaRE2 program is the Research and Education Core (REC) designed to coordinate the training of URM students and researchers at different levels in their academic and professional developments. The undergraduate cancer research training program under FAMU-CaRE2 Center is a 3-year (2018-2021) project to identify, train, mentor, and provide the URM undergraduate students with the support network they need to flourish in the program and beyond. In its year-1 funding cycle, the program has made significant progress in developing a novel framework for an undergraduate cancer research education and engagement program at FAMU, one of the forefront minority institutions in the nation. The mentored research program is complemented with professional development and engagement activities, including cancer research seminars, workshops, and community outreach activities. The purpose of this paper is to discuss the strategies implemented for an effective partnership, the leadership and mentoring skills, and outcomes from the year-1 experiences. In addition, we present the progress made in advancing the pool of underrepresented minority students with scientific and academic career progression paths focused on cancer health disparities.

Cancer transcriptomic profiling from rapidly enriched circulating tumor cells.

Transcriptomic profiling of metastatic cancer can illuminate mechanisms of progression and lead to new therapies, but standard biopsy is invasive and reflects only a single metastatic site. In contrast, circulating tumor cell (CTC) profiling is noninvasive and repeatable, reflecting the dynamic and systemic nature of advanced disease. To date, transcriptomic profiling of CTCs has not delivered on its full potential, because white blood cells (WBCs) vastly outnumber CTCs. Current profiling strategies either lack cancer sensitivity and specificity or require specialized CTC capture protocols that are not readily scalable to large patient cohorts. Here, we describe a new strategy for rapid CTC enrichment and transcriptomic profiling using commercially available WBC depletion, microfluidic enrichment and RNA sequencing. When applied to blood samples from patients with advanced prostate cancer (PC), transcriptomes from enriched samples cluster with cancer positive controls and previously undetectable prostate-specific transcripts become readily measurable. Gene set enrichment analysis reveals multiple significantly enriched signaling pathways associated with PC, as well as novel pathways that merit further study. This accessible and scalable approach yields cancer-specific transcriptomic data and can be applied repeatedly and noninvasively in large cancer patient cohorts to discover new therapeutic targets in advanced disease.

Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma.

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. In contrast with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

Comprehensive molecular profiling of 718 Multiple Myelomas reveals significant differences in mutation frequencies between African and European descent cases.

Multiple Myeloma (MM) is a plasma cell malignancy with significantly greater incidence and mortality rates among African Americans (AA) compared to Caucasians (CA). The overall goal of this study is to elucidate differences in molecular alterations in MM as a function of self-reported race and genetic ancestry. Our study utilized somatic whole exome, RNA-sequencing, and correlated clinical data from 718 MM patients from the Multiple Myeloma Research Foundation CoMMpass study Interim Analysis 9. Somatic mutational analyses based upon self-reported race corrected for ancestry revealed significant differences in mutation frequency between groups. Of interest, BCL7A, BRWD3, and AUTS2 demonstrate significantly higher mutation frequencies among AA cases. These genes are all involved in translocations in B-cell malignancies. Moreover, we detected a significant difference in mutation frequency of TP53 and IRF4 with frequencies higher among CA cases. Our study provides rationale for interrogating diverse tumor cohorts to best understand tumor genomics across populations.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium.

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHRadjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHRadjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHRadjBMI and eight previously established loci at P < 5×10-8: seven for BMI, and one for WHRadjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHRadjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHRadjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHRadjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHRadjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHRadjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants.

The vast majority of coding variants are rare, and assessment of the contribution of rare variants to complex traits is hampered by low statistical power and limited functional data. Improved methods for predicting the pathogenicity of rare coding variants are needed to facilitate the discovery of disease variants from exome sequencing studies. We developed REVEL (rare exome variant ensemble learner), an ensemble method for predicting the pathogenicity of missense variants on the basis of individual tools: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons. REVEL was trained with recently discovered pathogenic and rare neutral missense variants, excluding those previously used to train its constituent tools. When applied to two independent test sets, REVEL had the best overall performance (p < 10-12) as compared to any individual tool and seven ensemble methods: MetaSVM, MetaLR, KGGSeq, Condel, CADD, DANN, and Eigen. Importantly, REVEL also had the best performance for distinguishing pathogenic from rare neutral variants with allele frequencies <0.5%. The area under the receiver operating characteristic curve (AUC) for REVEL was 0.046-0.182 higher in an independent test set of 935 recent SwissVar disease variants and 123,935 putatively neutral exome sequencing variants and 0.027-0.143 higher in an independent test set of 1,953 pathogenic and 2,406 benign variants recently reported in ClinVar than the AUCs for other ensemble methods. We provide pre-computed REVEL scores for all possible human missense variants to facilitate the identification of pathogenic variants in the sea of rare variants discovered as sequencing studies expand in scale.

Histological Validation of 11Carbon-Acetate Positron Emission Tomography/Computerized Tomography in Detecting Lymph Node Metastases in Prostate Cancer.

PURPOSE: Conventional imaging cannot definitively detect nodal metastases of prostate cancer. We histologically validated C-acetate positron emission tomography/computerized tomography to identify nodal metastases, examining prostate cancer factors that influence detection rates. MATERIALS AND METHODS: Patients with C-acetate avid positron emission tomography/computerized tomography imaged pelvic/retroperitoneal lymph nodes underwent high extended robotic lymphadenectomy. A standardized mapping template comprising 8 predetermined anatomical regions was dissected during lymphadenectomy, allowing for matched, region based analysis and comparison of imaging and histological data. RESULTS: In 25 patients a total of 2,149 lymph nodes were excised (mean 86 per patient, range 27 to 136) and 528 (22%) harbored metastases (mean 21 positive nodes per patient, range 0 to 109). A total of 174 anatomical regions had matching imaging histological data. C-acetate positron emission tomography/computerized tomography accurately identified 48 node-positive regions and accurately ruled out 88 regions as metastasis-free. C-acetate sensitivity, specificity, and positive and negative predictive values were 67%, 84%, 74% and 79%, respectively. An increasing, histologically measured metastatic lesion size in long axis diameter of 5 or less, 6 to 10, 11 to 15, 16 to 20 and 21 mm or greater correlated with improved C-acetate detection rates of 45%, 62%, 81%, 89% and 100%, respectively. Each standard uptake value unit increase correlated with a 1.9 mm increase in nodal long axis diameter (p <0.001) and a 1.2 mm increase in short axis diameter (p <0.001). Positive C-acetate positron emission tomography/computerized tomography findings correlated with histological lymph node size (long axis diameter 12 mm and short axis diameter 6 mm), metastatic lesion size (long axis diameter 11 mm and short axis diameter 6 mm) and extranodal extension (positive 88% vs false-negative 58%, p = 0.005). CONCLUSIONS: C-acetate positron emission tomography/computerized tomography can identify prostate cancer metastatic nodal disease. However, it underestimates the true cephalad extent of nodal involvement, performing better in the pelvis than in the retroperitoneum. Standard uptake value, histological nodal size, intranodal metastasis size and extranodal extension correlate with cancer bearing nodes.

E6201, an intravenous MEK1 inhibitor, achieves an exceptional response in BRAF V600E-mutated metastatic malignant melanoma with brain metastases.

Malignant melanoma (MM) exhibits a high propensity for central nervous system dissemination with ~50% of metastatic MM patients developing brain metastases (BM). Targeted therapies and immune checkpoint inhibitors have improved overall survival for MM patients with BM. However, responses are usually of short duration and new agents that effectively penetrate the blood brain barrier (BBB) are needed. Here, we report a MM patient with BM who experienced an exceptional response to E6201, an ATP-competitive MEK1 inhibitor, on a Phase 1 study, with ongoing near-complete response and overall survival extending beyond 8 years. Whole exome and transcriptome sequencing revealed a high mutational burden tumor (22 mutations/Megabase) with homozygous BRAF V600E mutation. Correlative preclinical studies demonstrated broad activity for E6201 across BRAF V600E mutant melanoma cell lines and effective BBB penetration in vivo. Together, these results suggest that E6201 may represent a potential new treatment option for BRAF-mutant MM patients with BM.

Mitochondrial biology and prostate cancer ethnic disparity.

Prostate cancer remains the second most prevalent cancer in men. Its incidence, progression and mortality profiles vary significantly by race and ethnicity, with African-American men having the highest incidence rate and mortality rate in the world. Although these disparities can be partially explained by socioeconomic factors, the underlying molecular causes are complex and require careful research. A considerable amount of literature exists, supporting the association between mitochondrial health and the incidence, aggression and risk of prostate cancer. Genetic alterations in mitochondrial DNA are frequent in prostate cancer; therefore, the resulting mitochondrial dysfunction and metabolic dysregulation may contribute to or indicate oncogenesis. Many of the prominent features of cancer cells are also closely related to mitochondrial functions, such as resistance to apoptosis, excess reactive oxygen species production and altered oxidative phosphorylation. In addition, prostate cancer ethnic disparity is influenced by environmental and lifestyle factors, which involves differences in mitochondrial metabolism and retrograde signaling events.

Lack of APC somatic mutation is associated with early-onset colorectal cancer in African Americans.

African Americans (AAs) have higher incidence and mortality rates of colorectal cancer (CRC) compared with other US populations. They present with more right-sided, microsatellite stable disease and are diagnosed at earlier ages compared with non-Hispanic Whites (NHWs). To gain insight into these trends, we conducted exome sequencing (n = 45), copy number (n = 33) and methylation analysis (n = 11) of microsatellite stable AA CRCs. Results were compared with data from The Cancer Genome Atlas (TCGA). Two of the 45 tumors contained POLE mutations. In the remaining 43 tumors, only 27 (63%) contained loss-of-function mutations in APC compared with 80% of TCGA NHW CRCs. APC-mutation-negative CRCs were associated with an earlier onset of CRC (P = 0.01). They were also associated with lower overall mutation burden, fewer copy number variants and a DNA methylation signature that was distinct from the CpG island methylator phenotype characterized in microsatellite unstable disease. Three of the APC-mutation-negative CRCs had loss-of-function mutations in BCL9L. Mutations in driver genes identified by TCGA exome analysis were less frequent in AA CRC cases than TCGA NHWs. Genes that regulate the WNT signaling pathway, including SOX9, GATA6, TET1, GLIS1 and FAT1, were differentially hypermethylated in APC-mutation-negative CRCs, suggesting a novel mechanism for cancer development in these tumors. In summary, we have identified a subtype of CRC that is associated with younger age of diagnosis, lack of APC mutation, microsatellite and chromosome stability, lower mutation burden and distinctive methylation changes.

gsSKAT: Rapid gene set analysis and multiple testing correction for rare-variant association studies using weighted linear kernels.

Next-generation sequencing technologies have afforded unprecedented characterization of low-frequency and rare genetic variation. Due to low power for single-variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel-machine regression and adaptive testing methods for aggregative rare-variant association testing have been demonstrated to be powerful approaches for pathway-level analysis, although these methods tend to be computationally intensive at high-variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare-variant analysis using component gene-level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family-wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case-control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open-source R code for public use to facilitate easy application of our methods to existing rare-variant analysis results.

Post hoc Analysis for Detecting Individual Rare Variant Risk Associations Using Probit Regression Bayesian Variable Selection Methods in Case-Control Sequencing Studies.

Rare variants (RVs) have been shown to be significant contributors to complex disease risk. By definition, these variants have very low minor allele frequencies and traditional single-marker methods for statistical analysis are underpowered for typical sequencing study sample sizes. Multimarker burden-type approaches attempt to identify aggregation of RVs across case-control status by analyzing relatively small partitions of the genome, such as genes. However, it is generally the case that the aggregative measure would be a mixture of causal and neutral variants, and these omnibus tests do not directly provide any indication of which RVs may be driving a given association. Recently, Bayesian variable selection approaches have been proposed to identify RV associations from a large set of RVs under consideration. Although these approaches have been shown to be powerful at detecting associations at the RV level, there are often computational limitations on the total quantity of RVs under consideration and compromises are necessary for large-scale application. Here, we propose a computationally efficient alternative formulation of this method using a probit regression approach specifically capable of simultaneously analyzing hundreds to thousands of RVs. We evaluate our approach to detect causal variation on simulated data and examine sensitivity and specificity in instances of high RV dimensionality as well as apply it to pathway-level RV analysis results from a prostate cancer (PC) risk case-control sequencing study. Finally, we discuss potential extensions and future directions of this work.