Fine-mapping of genome-wide association study-identified risk loci for colorectal cancer in African Americans

Combining Asian and European genome-wide association studies of colorectal cancer improves risk prediction across racial and ethnic populations

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expand PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS are 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1681-3651 cases and 8696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They are significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values < 0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

HSPA12A was identified as a key driver in colorectal cancer GWAS loci 10q26.12 and modulated by an enhancer-promoter interaction

Although genome-wide association studies (GWASs) have identified over 100 colorectal cancer (CRC) risk loci, an understanding of causal genes or risk variants and their biological functions in these loci remain unclear. Recently, genomic loci 10q26.12 with lead SNP rs1665650 was identified as an essential CRC risk loci of Asian populations. However, the functional mechanism of this region has not been fully clarified. Here, we applied an RNA interfering-based on-chip approach to screen for the genes essential for cell proliferation in the CRC risk loci 10q26.12. Notably, HSPA12A had the most significant effect among the identified genes and functioned as a crucial oncogene facilitating cell proliferation. Moreover, we conducted an integrative fine-mapping analysis to identify putative casual variants and further explored their association with CRC risk in a large-scale Chinese population consisting of 4054 cases and 4054 controls and also independently validated in 5208 cases and 20,832 controls from the UK biobank cohort. We identified a risk SNP rs7093835 in the intron of HSPA12A that was significantly associated with an increased risk of CRC (OR 1.23, 95% CI 1.08-1.41, P = 1.92 × 10-3). Mechanistically, the risk variant could facilitate an enhancer-promoter interaction mediated by the transcriptional factor (TF) GRHL1 and ultimately upregulate HSPA12A expression, which provides functional evidence to support our population findings. Collectively, our study reveals the important role of HSPA12A in CRC development and illustrates a novel enhancer-promoter interaction module between HSPA12A and its regulatory elements rs7093835, providing new insights into the etiology of CRC.

Recurrent EWSR1::COLCA2 Fusions Define a Novel Sarcoma With Spindle/Round Cell Morphology and Strong Predilection for the Sinonasal Tract

The last 2 decades have attended a dynamic evolution in the nosology of poorly differentiated sinonasal tract malignancies, with several new molecularly defined entities having been described in addition to delineation of the genetic driver/s of some established older entities. These discoveries, however, mostly concerned epithelial-derived neoplasms (carcinomas). Adamantinoma-like Ewing sarcoma and biphenotypic sinonasal sarcoma are the major representatives of the newly defined mesenchymal categories. The colorectal cancer associated 2 (COLCA2) has been discovered recently as a colorectal cancer risk gene locus, but fusions involving this gene have not been well characterized. We, herein, describe clinicopathologic and molecular features of a novel sinonasal sarcoma characterized by undifferentiated spindle/round cell morphology and defined by recurrent EWSR1::COLCA2 fusions. All patients (n=5) were adults (3 female and 2 male) with a median age of 46 years (range, 23 to 60 y). The tumors originated in different subsites of the sinonasal tract with frequent multisite involvement. Original diagnoses were undifferentiated or unclassified round cell/spindle cell neoplasm/sarcoma (n=4) and neuroendocrine carcinoma (n=1). Surgery with or without adjuvant chemoradiation was the treatment in all cases. At the last follow-up, 1 patient developed multiple local recurrences over 21 years and another developed local recurrence and distant metastasis to bone 27 months after diagnosis. A third patient developed local recurrence 11 months later. Two patients were disease-free at 23, and 24 months. Histology showed nondescript highly cellular neoplasms with an admixture of spindled and round cells disposed into solid sheets and fascicles with brisk mitotic activity. Immunohistochemistry was negative for all lineage-specific markers with only limited focal membranous CD99 (4 of 5 cases) and weak pankeratin (1 of 5 cases) expression. Targeted RNA sequencing revealed an EWSR1::COLCA2 fusion, verified by EWSR1 fluorescence in situ hybridization, in all cases. This series identifies a novel member in the undifferentiated spindle/round cell sarcoma category with strong predilection for the sinonasal tract. None of >10,000 epithelial and mesenchymal neoplasms tested at the authors' centers during the same period showed this fusion, highlighting rarity of tumors carrying this gene fusion. Accordingly, molecular testing of unclassified sinonasal malignancies/sarcomas showing round and spindle cell morphology is recommended to enhance the identification and further characterization of this entity.

Genome-wide association study (GWAS) of circulating vitamin D outcomes among individuals of African ancestry

BACKGROUND

Vitamin D deficiency is more common among African-ancestry individuals and may be associated with adverse health outcomes. Vitamin D binding protein (VDBP) regulates concentrations of biologically active vitamin D.

OBJECTIVE

We conducted genome-wide association study (GWAS) of VDBP and 25-hydroxyvitamin D among African-ancestry individuals.

METHODS

Data were collected from 2,602 African American adults from the Southern Community Cohort Study (SCCS) and 6,934 African- or Caribbean-ancestry adults from the UK Biobank. Serum VDBP concentrations were available only in the SCCS and were measured by using the Polyclonal Human VDBP ELISA kit. Serum 25-hydroxyvitamin D concentrations for both study samples were measured by using Diasorin Liason, a chemiluminescent immunoassay. Participants were genotyped for single nucleotide polymorphisms (SNPs) with genome-wide coverage by using Illumina or Affymetrix platforms. Fine-mapping analysis was performed by using forward stepwise linear regression models including all variants with P value < 5 × 10-8 and within 250 kbps of a lead SNP.

RESULTS

We identified 4 loci notably associated with VDBP concentrations in the SCCS population: rs7041 (per allele β = 0.61 μg/mL, SE = 0.05, P = 1.4 × 10-48) and rs842998 (per allele β = 0.39 μg/mL, SE = 0.03, P = 4.0 × 10-31) in GC, rs8427873 (per allele β = 0.31 μg/mL, SE = 0.04, P = 3.0 × 10-14) near GC and rs11731496 (per allele β = 0.21 μg/mL, SE = 0.03, P = 3.6 × 10-11) in between GC and NPFFR2. In conditional analyses, which included the above-mentioned SNPs, only rs7041 remained notable (P = 4.1 × 10-21). SNP rs4588 in GC was the only GWAS-identified SNP associated with 25-hydroxyvitamin D concentration. Among UK Biobank participants: per allele β = -0.11 μg/mL, SE = 0.01, P = 1.5 × 10-13; in the SCCS: per allele β = -0.12 μg/mL, SE = 0.06, P = 2.8 × 10-02). rs7041 and rs4588 are functional SNPs that influence the binding affinity of VDBP to 25-hydroxyvitamin D.

CONCLUSIONS

Our results were in line with previous studies conducted in European-ancestry populations, showing that GC, the gene that directly encodes for VDBP, would be important for VDBP and 25-hydroxyvitamin D concentrations. The current study extends our knowledge of the genetics of vitamin D in diverse populations.

Combining Asian-European Genome-Wide Association Studies of Colorectal Cancer Improves Risk Prediction Across Race and Ethnicity

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expanded PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS were 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1,681-3,651 cases and 8,696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They were significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values<0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

Repurposing non-invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome

In our previous work, genomic data generated through non-invasive prenatal testing (NIPT) based on low-coverage massively parallel whole-genome sequencing of total plasma DNA of pregnant women in Slovakia was described as a valuable source of population specific data. In the present study, these data were used to determine the population allele frequency of common risk variants located in genes associated with colorectal cancer (CRC) and Lynch syndrome (LS). Allele frequencies of identified variants were compared with six world populations to detect significant differences between populations. Finally, variants were interpreted, functional consequences were searched for and clinical significance of variants was investigated using publicly available databases. Although the present study did not identify any pathogenic variants associated with CRC or LS in the Slovak population using NIPT data, significant differences were observed in the allelic frequency of risk CRC variants previously reported in genome-wide association studies and common variants located in genes associated with LS. As Slovakia is one of the leading countries with the highest incidence of CRC among male patients in the world, there is a need for studies dedicated to investigating the cause of such a high incidence of CRC in Slovakia. The present study also assumed that extensive cross-country data aggregation of NIPT results would represent an unprecedented source of information concerning human genome variation in cancer research.

Ability of known colorectal cancer susceptibility SNPs to predict colorectal cancer risk: A cohort study within the UK Biobank

Colorectal cancer risk stratification is crucial to improve screening and risk-reducing recommendations, and consequently do better than a one-size-fits-all screening regimen. Current screening guidelines in the UK, USA and Australia focus solely on family history and age for risk prediction, even though the vast majority of the population do not have any family history. We investigated adding a polygenic risk score based on 45 single-nucleotide polymorphisms to a family history model (combined model) to quantify how it improves the stratification and discriminatory performance of 10-year risk and full lifetime risk using a prospective population-based cohort within the UK Biobank. For both 10-year and full lifetime risk, the combined model had a wider risk distribution compared with family history alone, resulting in improved risk stratification of nearly 2-fold between the top and bottom risk quintiles of the full lifetime risk model. Importantly, the combined model can identify people (n = 72,019) who do not have family history of colorectal cancer but have a predicted risk that is equivalent to having at least one affected first-degree relative (n = 44,950). We also confirmed previous findings by showing that the combined full lifetime risk model significantly improves discriminatory accuracy compared with a simple family history model 0.673 (95% CI 0.664–0.682) versus 0.666 (95% CI 0.657–0.675), p = 0.0065. Therefore, a combined polygenic risk score and first-degree family history model could be used to improve risk stratified population screening programs.

Genetic risk factors for colorectal cancer in multiethnic Indonesians

Colorectal cancer is a common cancer in Indonesia, yet it has been understudied in this resource-constrained setting. We conducted a genome-wide association study focused on evaluation and preliminary discovery of colorectal cancer risk factors in Indonesians. We administered detailed questionnaires and collecting blood samples from 162 colorectal cancer cases throughout Makassar, Indonesia. We also established a control set of 193 healthy individuals frequency matched by age, sex, and ethnicity. A genome-wide association analysis was performed on 84 cases and 89 controls passing quality control. We evaluated known colorectal cancer genetic variants using logistic regression and established a genome-wide polygenic risk model using a Bayesian variable selection technique. We replicate associations for rs9497673, rs6936461 and rs7758229 on chromosome 6; rs11255841 on chromosome 10; and rs4779584, rs11632715, and rs73376930 on chromosome 15. Polygenic modeling identified 10 SNP associated with colorectal cancer risk. This work helps characterize the relationship between variants in the SCL22A3, SCG5, GREM1, and STXBP5-AS1 genes and colorectal cancer in a diverse Indonesian population. With further biobanking and international research collaborations, variants specific to colorectal cancer risk in Indonesians will be identified.

Variants in the 8q24 region associated with risk of breast cancer: Systematic research synopsis and meta-analysis

Breast cancer is a molecularly heterogeneous disorder associated with high lethal malignant tumors among women worldwide. Genetic factors play an important role in breast cancer development. Several single nucleotide polymorphisms in the 8q24 region associated with risk of breast cancer have been identified. Fifteen studies including 32,955 cases and 43,716 controls were collected to conduct a meta-analysis to evaluate the associations between variants in 8q24 region and risk of breast cancer. Our study showed that only rs13281615 is associated with breast cancer risk in this large-scale research synopsis and meta-analysis. Further studies are needed to explore the role of the 8q24 variants in the development of breast cancer.

Systematic Functional Interrogation of Genes in GWAS Loci Identified ATF1 as a Key Driver in Colorectal Cancer Modulated by a Promoter-Enhancer Interaction

Genome-wide association studies (GWASs) have identified approximately 100 colorectal cancer (CRC) risk loci. However, the causal genes in these loci have not been systematically interrogated. We conducted a high-throughput RNA-interference functional screen to identify the genes essential for proliferation in the CRC risk loci of Asian populations. We found that ATF1, located in the 12q13.12 region, functions as an oncogene that facilitates cell proliferation; ATF1 has the most significant effect of the identified genes and promotes CRC xenograft growth by affecting cell apoptosis. Next, by integrating a fine-mapping analysis, a two-stage affected-control study consisting of 6,213 affected individuals and 10,388 controls, and multipronged experiments, we elucidated that two risk variants, dbSNP: rs61926301 and dbSNP: rs7959129, that located in the ATF1 promoter and first intron, respectively, facilitate a promoter-enhancer interaction, mediated by the synergy of SP1 and GATA3, to upregulate ATF1 expression, thus synergistically predisposing to CRC risk (OR = 1.77, 95% CI = 1.42-2.21, p = 3.16 × 10^-7; P_{multiplicative-interaction} = 1.20 × 10^-22; P_{additive-interaction} = 6.50 × 10^-3). Finally, we performed RNA-seq and ChIP-seq assays in CRC cells treated with ATF1 overexpression in order to dissect the target programs of ATF1. Results showed that ATF1 activates a subset of genes, including BRAF, NRAS, MYC, BIRC2, DAAM1, MAML2, STAT1, ID1, and NKD2, related to apoptosis, Wnt, TGF-β, and MAPK pathways, and these effects could cooperatively increase the risk of CRC. These findings reveal the clinical potential of ATF1 in CRC development and illuminate a promoter-enhancer interaction module between the ATF1 regulatory elements dbSNP: rs61926301 and dbSNP: rs7959129, and they bring us closer to understanding the molecular drivers of cancer.

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.

A geographically matched control population efficiently limits the number of candidate disease-causing variants in an unbiased whole-genome analysis

Whole-genome sequencing is a promising approach for human autosomal dominant disease studies. However, the vast number of genetic variants observed by this method constitutes a challenge when trying to identify the causal variants. This is often handled by restricting disease studies to the most damaging variants, e.g. those found in coding regions, and overlooking the remaining genetic variation. Such a biased approach explains in part why the genetic causes of many families with dominantly inherited diseases, in spite of being included in whole-genome sequencing studies, are left unsolved today. Here we explore the use of a geographically matched control population to minimize the number of candidate disease-causing variants without excluding variants based on assumptions on genomic position or functional predictions. To exemplify the benefit of the geographically matched control population we apply a typical disease variant filtering strategy in a family with an autosomal dominant form of colorectal cancer. With the use of the geographically matched control population we end up with 26 candidate variants genome wide. This is in contrast to the tens of thousands of candidates left when only making use of available public variant datasets. The effect of the local control population is dual, it (1) reduces the total number of candidate variants shared between affected individuals, and more importantly (2) increases the rate by which the number of candidate variants are reduced as additional affected family members are included in the filtering strategy. We demonstrate that the application of a geographically matched control population effectively limits the number of candidate disease-causing variants and may provide the means by which variants suitable for functional studies are identified genome wide.

Novel Common Genetic Susceptibility Loci for Colorectal Cancer

BACKGROUND

Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10-8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk.

METHODS

We conducted a GWAS in European descent CRC cases and control subjects using a discovery-replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10-8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided.

RESULTS

The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10-8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0.

CONCLUSIONS

This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.

CDH1 rs9929218 variant at 16q22.1 contributes to colorectal cancer susceptibility

Colorectal cancer (CRC) is the third most common cancer. Large-scale genome-wide association studies (GWAS) have been performed and reported some novel CRC susceptibility variants in European ancestry including the CDH1 rs9929218. Following GWAS and candidate studies evaluated the association between the CDH1 rs9929218 polymorphism and CRC in European, Asian and American populations. However, these studies reported inconsistent associations. Evidence shows that rs9929218 may regulate different gene expressions in different human tissues. Here, we reevaluated this association using large-scale samples from 16 studies (n=131768) using a meta-analysis method. In heterogeneity test, we did not identify significant heterogeneity among these studies. Meta-analysis using fixed effect model showed significant association between rs9929218 and CRC (P=6.16E-21, odds ratio (OR) =0.92, 95% confidence interval (CI) 0.91-0.94). In order to validate the effect of rs9929218 variant on CDH1 expression, we further performed a functional analysis using two large-scale expression datasets. We identified significant regulation relation between rs9929218 variant and the expression of CDH1, ZFP90, RP11-354M1.2 and MCOLN2 by both cis-effect and trans-effect. In summary, our analysis highlights significant association between rs9929218 polymorphism and CRC susceptibility.

Colorectal Cancer Disparity in African Americans: Risk Factors and Carcinogenic Mechanisms

African Americans have the highest incidence and mortality rates of colorectal cancer (CRC) of any ethnic group in the United States. Although some of these disparities can be explained by differences in access to care, cancer screening, and other socioeconomic factors, disparities remain after adjustment for these factors. Consequently, an examination of recent advances in the understanding of ethnicity-specific factors, including genetic and environmental factors relating to risk of CRC, the biology of CRC progression, and the changes in screening and mortality, is important for evaluating our progress toward eliminating the disparities. An overarching limitation in this field is the number and sample size of studies performed to characterize the etiological bases of CRC incidence and mortality in African Americans. Despite this limitation, significant differences in etiology are manifest in many studies. These differences need validation, and their impacts on disparities need more detailed investigation. Perhaps most heartening, improvements in CRC screening can be attributed to the smallest difference in CRC incidence between African Americans and whites since the late 1980s. Cancer mortality, however, remains a persistent difference.

Novel colon cancer susceptibility variants identified from a genome-wide association study in African Americans

Genome-wide association studies (GWAS) in ethnic/racial minority populations can help to fine-map previously identified risk regions or discover new risk loci because of the genetic diversity in these populations. We conducted a GWAS of colorectal cancer (CRC) in 6,597 African Americans (1,894 cases and 4,703 controls) (Stage 1) and followed up the most promising markers in a replication set of 2,041 participants of African descent (891 cases and 1,150 controls) (Stage 2). We identified a novel variant, rs56848936 in the gene SYMPK at 19q13.3, associated with colon cancer risk (odds ratio 0.61 for the risk allele G, p = 2.4 × 10-8 ). The frequency of the G allele was 0.06 in African Americans, compared to <0.01 in Europeans, Asians and Amerindians in the 1000 Genomes project. In addition, a variant previously identified through fine-mapping in this GWAS in the region 19q13.1, rs7252505, was confirmed to be more strongly associated with CRC in the African American replication set than the variant originally reported in Europeans (rs10411210). The association between rs7252505 and CRC was of borderline significance (p = 0.05) in a Hispanic population GWAS with 1,611 CRC cases and 4,330 controls. With the three datasets combined, the odds ratio was 0.84 for the risk allele A (95% confidence interval 0.79-0.89, p = 3.7 × 10-8 ). This study further highlights the importance of conducting GWAS studies in diverse ancestry populations.

Fine-Mapping of Common Genetic Variants Associated with Colorectal Tumor Risk Identified Potential Functional Variants

Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) associated with colorectal cancer risk. These SNPs may tag correlated variants with biological importance. Fine-mapping around GWAS loci can facilitate detection of functional candidates and additional independent risk variants. We analyzed 11,900 cases and 14,311 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. To fine-map genomic regions containing all known common risk variants, we imputed high-density genetic data from the 1000 Genomes Project. We tested single-variant associations with colorectal tumor risk for all variants spanning genomic regions 250-kb upstream or downstream of 31 GWAS-identified SNPs (index SNPs). We queried the University of California, Santa Cruz Genome Browser to examine evidence for biological function. Index SNPs did not show the strongest association signals with colorectal tumor risk in their respective genomic regions. Bioinformatics analysis of SNPs showing smaller P-values in each region revealed 21 functional candidates in 12 loci (5q31.1, 8q24, 11q13.4, 11q23, 12p13.32, 12q24.21, 14q22.2, 15q13, 18q21, 19q13.1, 20p12.3, and 20q13.33). We did not observe evidence of additional independent association signals in GWAS-identified regions. Our results support the utility of integrating data from comprehensive fine-mapping with expanding publicly available genomic databases to help clarify GWAS associations and identify functional candidates that warrant more onerous laboratory follow-up. Such efforts may aid the eventual discovery of disease-causing variant(s).

Genome-wide association study of colorectal cancer in Hispanics

This manuscript describes the first large-scale genome-wide association study of colorectal cancer in Hispanics and Latinos. Our results demonstrate the broad replication of known susceptibility regions and the importance of fine-mapping in ethnic minority populations.

Genome-wide association studies (GWAS) have identified 58 susceptibility alleles across 37 regions associated with the risk of colorectal cancer (CRC) with P < 5×10⁻⁸. Most studies have been conducted in non-Hispanic whites and East Asians; however, the generalizability of these findings and the potential for ethnic-specific risk variation in Hispanic and Latino (HL) individuals have been largely understudied. We describe the first GWAS of common genetic variation contributing to CRC risk in HL (1611 CRC cases and 4330 controls). We also examine known susceptibility alleles and implement imputation-based fine-mapping to identify potential ethnicity-specific association signals in known risk regions. We discovered 17 variants across 4 independent regions that merit further investigation due to suggestive CRC associations (P < 1×10⁻⁶) at 1p34.3 (rs7528276; Odds Ratio (OR) = 1.86 [95% confidence interval (CI): 1.47–2.36); P = 2.5×10⁻⁷], 2q23.3 (rs1367374; OR = 1.37 (95% CI: 1.21–1.55); P = 4.0×10⁻⁷), 14q24.2 (rs143046984; OR = 1.65 (95% CI: 1.36–2.01); P = 4.1×10⁻⁷) and 16q12.2 [rs142319636; OR = 1.69 (95% CI: 1.37–2.08); P=7.8×10⁻⁷]. Among the 57 previously published CRC susceptibility alleles with minor allele frequency ≥1%, 76.5% of SNPs had a consistent direction of effect and 19 (33.3%) were nominally statistically significant (P < 0.05). Further, rs185423955 and rs60892987 were identified as novel secondary susceptibility variants at 3q26.2 (P = 5.3×10^–5) and 11q12.2 (P = 6.8×10⁻⁵), respectively. Our findings demonstrate the importance of fine mapping in HL. These results are informative for variant prioritization in functional studies and future risk prediction modeling in minority populations.

Interaction between Red Meat Intake and NAT2 Genotype in Increasing the Risk of Colorectal Cancer in Japanese and African Americans

Heterocyclic aromatic amines formed in cooked meat may be an underlying mechanism for the red meat-colorectal cancer (CRC) association. These compounds require bioactivaction by N-acetyltransferase 2 (NAT2). An interaction effect between red meat consumption and NAT2 in increasing CRC risk has been inconsistently reported in whites. We investigated this interaction in two populations in which the high-activity rapid NAT2 phenotype is 10- and 2-fold more common than in whites. We meta-analyzed four studies of Japanese (2,217 cases, 3,788 controls) and three studies of African Americans (527 cases, 4,527 controls). NAT2 phenotype was inferred from an optimized seven-SNP genotyping panel. Processed and total red meat intakes were associated with an increased CRC risk in Japanese and in both ethnic groups combined (P’s ≤ 0.002). We observed an interaction between processed meat intake and NAT2 in Japanese (P = 0.04), African Americans (P = 0.02), and in both groups combined (P = 0.006). The association of processed meat with CRC was strongest among individuals with the rapid NAT2 phenotype (combined analysis, OR for highest vs. lowest quartile: 1.62, 95% CI: 1.28–2.05; P_trend = 8.0×10⁻⁵), intermediate among those with the intermediate NAT2 phenotype (1.29, 95% CI: 1.05–1.59; P_trend = 0.05) and null among those with the slow phenotype (P_trend = 0.45). A similar interaction was found for NAT2 and total red meat (P_interaction = 0.03). Our findings support a role for NAT2 in modifying the association between red meat consumption and CRC in Japanese and African Americans.

Association of 8q23-24 region (8q23.3 loci and 8q24.21 loci) with susceptibility to colorectal cancer: a systematic and updated meta-analysis

BACKGROUND

Several single nucleotide polymorphisms (SNPs), rs16892766 in the 8q23.3 region and rs6983267, rs10505477, rs7014346 and rs7837328 in the 8q24.21 region, have been identified by genome-wide association studies (GWAS) and a number of case-control studies to be closely associated with risk of colorectal cancer (CRC). In the present study, a meta-analysis was performed to confirm if these loci are risk factors for susceptibility to CRC, taking heterogeneity of population into consideration.

METHODS

The whole literature search was conducted via database of MEDLINE and Embase, through which 33 articles with 49 studies (141,899 cases and 157,536 controls) were finally included in this meta-analysis to evaluate the association between the 5 polymorphisms and risk of CRC under allelic model.

RESULTS

A meta-analysis of the pooled data showed that the G allele of rs6983267, the A allele of rs7014346, the T allele of rs10505477, the C allele of rs16892766 and the A allele of rs7837328 were associated with significantly increased risk of CRC under allelic model. Additionally, subgroup analyses of four SNPs (rs7837328 excluded) by ethnicity witnessed a notable association between the G allele of rs6983267 and increased risk of CRC among Caucasians, Asians and Africans. Furthermore, the C allele of rs16892766 was strongly linked with elevated risk of CRC among Caucasians and Africans. However, the A allele of rs7014346 and T allele of rs10505477 only heightened risk for CRC among Caucasians and showed no effects among Asians.

CONCLUSION

In summary, rs6983267 is a risk factor for CRC among Caucasians, Asians and Africans; rs7014346 and rs10505477 are risky genetic polymorphisms only among Caucasians; rs16892766 is a hazardous element among populations with Caucasian and African ancestry; and rs7837328 could elevate the susceptibility to CRC in a multinational group. However, more potential factors related with CRC risk should be investigated in further studies.

Genetic architecture of colorectal cancer

Colorectal cancer (CRC) is a complex disease that develops as a consequence of both genetic and environmental risk factors. A small proportion (3-5%) of cases arise from hereditary syndromes predisposing to early onset CRC as a result of mutations in over a dozen well defined genes. In contrast, CRC is predominantly a late onset 'sporadic' disease, developing in individuals with no obvious hereditary syndrome. In recent years, genome wide association studies have discovered that over 40 genetic regions are associated with weak effects on sporadic CRC, and it has been estimated that increasingly large genome wide scans will identify many additional novel genetic regions. Subsequent experimental validations have identified the causally related variant(s) in a limited number of these genetic regions. Further biological insight could be obtained through ethnically diverse study populations, larger genetic sequencing studies and development of higher throughput functional experiments. Along with inherited variation, integration of the tumour genome may shed light on the carcinogenic processes in CRC. In addition to summarising the genetic architecture of CRC, this review discusses genetic factors that modify environmental predictors of CRC, as well as examples of how genetic insight has improved clinical surveillance, prevention and treatment strategies. In summary, substantial progress has been made in uncovering the genetic architecture of CRC, and continued research efforts are expected to identify additional genetic risk factors that further our biological understanding of this disease. Subsequently these new insights will lead to improved treatment and prevention of colorectal cancer.

Enrichment of inflammatory bowel disease and colorectal cancer risk variants in colon expression quantitative trait loci

Background

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with diseases of the colon including inflammatory bowel diseases (IBD) and colorectal cancer (CRC). However, the functional role of many of these SNPs is largely unknown and tissue-specific resources are lacking. Expression quantitative trait loci (eQTL) mapping identifies target genes of disease-associated SNPs. This study provides a comprehensive eQTL map of distal colonic samples obtained from 40 healthy African Americans and demonstrates their relevance for GWAS of colonic diseases.

Results

8.4 million imputed SNPs were tested for their associations with 16,252 expression probes representing 12,363 unique genes. 1,941 significant cis-eQTL, corresponding to 122 independent signals, were identified at a false discovery rate (FDR) of 0.01. Overall, among colon cis-eQTL, there was significant enrichment for GWAS variants for IBD (Crohn’s disease [CD] and ulcerative colitis [UC]) and CRC as well as type 2 diabetes and body mass index. ERAP2, ADCY3, INPP5E, UBA7, SFMBT1, NXPE1 and REXO2 were identified as target genes for IBD-associated variants. The CRC-associated eQTL rs3802842 was associated with the expression of C11orf93 (COLCA2). Enrichment of colon eQTL near transcription start sites and for active histone marks was demonstrated, and eQTL with high population differentiation were identified.

Conclusions

Through the comprehensive study of eQTL in the human colon, this study identified novel target genes for IBD- and CRC-associated genetic variants. Moreover, bioinformatic characterization of colon eQTL provides a tissue-specific tool to improve understanding of biological differences in diseases between different ethnic groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1292-z) contains supplementary material, which is available to authorized users.

Common genetic variant on BMP4 contributes to colorectal adenoma and cancer: A meta-analysis based on 15 studies

Epidemiological studies indicate a genetic contribution to colorectal cancer (CRC), but specific genetic variants remain unknown. Genome-wide association studies have identified rs4444235 at BMP4 as a new colorectal cancer (CRC) and colorectal adenoma (CRA) susceptibility locus in populations of European descent. After that, several validation studies have been conducted among various ethnic populations to investigate if the SNP was associated with CRC/CRA, but the results have been inconsistent. To investigate this inconsistency and derive a more precise estimation of the relationship, a meta-analysis involving 54,631 CRC cases, 3995 CRA cases and 88,098 controls from 15 studies was performed. Potential sources of heterogeneity including ethnicity, sample size, study design and endpoint were also assessed. Overall, the summary OR of CRC was 1.06 (95% CI: 1.04-1.08, P<10(-5)). In the subgroup analysis by ethnicity, significantly increased risks were found in East Asians (OR=1.07, 95% CI: 1.01-1.12, P=0.01) and Caucasians (OR=1.07, 95% CI: 1.05-1.10, P<10(-5)); while no significant associations were found among African Americans and other ethnic populations in all genetic models. In addition, significant associations were also detected for CRA with per-allele OR of 1.09 (95% CI: 1.03-1.14, P=0.001). Our findings demonstrated that BMP4-rs4444235 is a risk factor associated with increased CRC and CRA susceptibility, but these associations vary in different ethnic populations.

Analyzing large-scale samples confirms the association between rs16892766 polymorphism and colorectal cancer susceptibility

Colorectal cancer (CRC) is a common complex disease caused by the combination of genetic variants and environmental factors. Genome-wide association studies (GWAS) have been performed and reported some novel CRC susceptibility variants. The rs16892766 (8q23.3) polymorphism was first identified to be significantly associated with CRC in European ancestry. The following studies investigated this association in Chinese, Japanese, Romanian, Swedish, African American, European American, and Croatian populations. These studies reported consistent and inconsistent results. Here, we reevaluated this association using the relatively large-scale samples from 13 studies (N = 59737, 26237 cases and 33500 controls) using a meta-analysis by searching the PubMed, Google Scholar and CRCgene databases. We observed no significant heterogeneity among the included studies. Our results showed significant association between rs16892766 polymorphism and CRC (P = 1.33E-35, OR = 1.23, 95% CI 1.20-1.27). Collectively, our analysis further supports previous findings that the rs16892766 polymorphism is significantly associated with CRC susceptibility. We believe that our findings will be very useful for future genetic studies on CRC.

Red meat intake, NAT2, and risk of colorectal cancer: a pooled analysis of 11 studies

BACKGROUND

Red meat intake has been associated with risk of colorectal cancer, potentially mediated through heterocyclic amines. The metabolic efficiency of N-acetyltransferase 2 (NAT2) required for the metabolic activation of such amines is influenced by genetic variation. The interaction between red meat intake, NAT2 genotype, and colorectal cancer has been inconsistently reported.

METHODS

We used pooled individual-level data from the Colon Cancer Family Registry and the Genetics and Epidemiology of Colorectal Cancer Consortium. Red meat intake was collected by each study. We inferred NAT2 phenotype based on polymorphism at rs1495741, highly predictive of enzyme activity. Interaction was assessed using multiplicative interaction terms in multivariate-adjusted models.

RESULTS

From 11 studies, 8,290 colorectal cancer cases and 9,115 controls were included. The highest quartile of red meat intake was associated with increased risk of colorectal cancer compared with the lowest quartile [OR, 1.41; 95% confidence interval (CI), 1.29-1.55]. However, a significant association was observed only for studies with retrospective diet data, not for studies with diet prospectively assessed before cancer diagnosis. Combining all studies, high red meat intake was similarly associated with colorectal cancer in those with a rapid/intermediate NAT2 genotype (OR, 1.38; 95% CI, 1.20-1.59) as with a slow genotype (OR, 1.43; 95% CI, 1.28-1.61; P interaction = 0.9).

CONCLUSION

We found that high red meat intake was associated with increased risk of colorectal cancer only from retrospective case-control studies and not modified by NAT2 enzyme activity.

IMPACT

Our results suggest no interaction between NAT2 genotype and red meat intake in mediating risk of colorectal cancer.

Common genetic variants (rs4779584 and rs10318) at 15q13.3 contributes to colorectal adenoma and colorectal cancer susceptibility: evidence based on 22 studies

Several genome-wide association studies on colorectal cancer (CRC) have reported similar findings of a new susceptibility locus, 15q13.3. After that, a number of studies have reported that the rs4779584 and rs10318 polymorphisms at chromosome 15q13.3 have been implicated in CRC and colorectal adenoma (CRA) risk; however, these studies have yielded inconsistent results. To investigate this inconsistency, we performed a meta-analysis of 22 studies involving a total of 48,468 CRC cases, 4,189 CRA cases, and 85,105 controls for the two polymorphisms to evaluate its effect on genetic susceptibility for CRC/CRA. Potential sources of heterogeneity and publication bias were also systematically explored. Overall, the summary odds ratio (OR) of rs4779584-T variant for CRC was 1.13 (95 % CI 1.09-1.16, P < 10(-5)) and 1.15 (95 % CI 1.04-1.28, P = 0.006) for CRA. After stratified by ethnicity, significantly increased CRC risks were found for rs4779584 polymorphism among East Asians and Caucasians, while no significant associations were detected among African American and other ethnic populations. A meta-analysis of studies on the rs10318 polymorphism also showed significant overall association with CRC, yielding a per-allele OR of 1.13 (95 % CI 1.02-1.24, P = 0.02). In the subgroup analysis by ethnicity, significantly increased CRC risks were found in Caucasians; whereas no significant associations were found among East Asians and African Americans. This meta-analysis demonstrated that the rs4779584 and rs10318 polymorphism at 15q13.3 is a risk factor associated with increased CRC/CRA susceptibility, but these associations vary in different ethnic populations.

Trans-ethnic genome-wide association study of colorectal cancer identifies a new susceptibility locus in VTI1A

The genetic basis of sporadic colorectal cancer (CRC) is not well explained by known risk polymorphisms. Here we perform a meta-analysis of two genome-wide association studies in 2,627 cases and 3,797 controls of Japanese ancestry and 1,894 cases and 4,703 controls of African ancestry, to identify genetic variants that contribute to CRC susceptibility. We replicate genome-wide statistically significant associations (P<5 × 10(-8)) in 16,823 cases and 18,211 controls of European ancestry. This study reveals a new pan-ethnic CRC risk locus at 10q25 (rs12241008, intronic to VTI1A; P=1.4 × 10(-9)), providing additional insight into the aetiology of CRC and highlighting the value of association mapping in diverse populations.

Shared and independent colorectal cancer risk alleles in TGFβ-related genes in African and European Americans

Genome-wide association studies (GWAS) in colorectal cancer (CRC) identified five regions near transforming growth factor β-related genes BMP4, GREM1, CDH1, SMAD7 and RPHN2. The true risk alleles remain to be identified in these regions, and their role in CRC risk in non-European populations has been understudied. Our previous work noted significant genetic heterogeneity between African Americans (AAs) and European Americans (EAs) for single nucleotide polymorphisms (SNPs) identified in GWAS. We hypothesized that associations may not have been replicated in AAs due to differential or independent genetic structures. In order to test this hypothesis, we genotyped 195 tagging SNPs across these five gene regions in 1194 CRC cases (795 AAs and 399 EAs) and 1352 controls (985 AAs and 367 EAs). Imputation was performed, and association testing of genotyped and imputed SNPs included ancestry, age and sex as covariates. In two of the five genes originally associated with CRC, we found evidence for association in AAs including rs1862748 in CDH1 (OR(Add) = 0.82, P = 0.02) and in GREM1 the SNPs rs10318 (OR(Rec) = 60.1, P = 0.01), rs11632715 (OR(Rec) = 2.36; P = 0.004) and rs12902616 (OR(Rec) = 1.28, P = 0.005), the latter which is in linkage disequilibrium with the previously identified SNP rs4779584. Testing more broadly for associations in these gene regions in AAs, we noted three statistically significant association peaks in GREM1 and RHPN2 that were not identified in EAs. We conclude that some CRC risk alleles are shared between EAs and AAs and others are population specific.