A large-scale microRNA transcriptome-wide association study identifies two susceptibility microRNAs, miR-1307-5p and miR-192-3p, for colorectal cancer risk.

Transcriptome-wide association studies (TWAS) have identified many putative susceptibility genes for colorectal cancer (CRC) risk. However, susceptibility miRNAs, critical dysregulators of gene expression, remain unexplored. We genotyped DNA samples from 313 CRC East Asian patients and performed small RNA sequencing in their normal colon tissues distant from tumors to build genetic models for predicting miRNA expression. We applied these models and data from genome-wide association studies (GWAS) including 23 942 cases and 217 267 controls of East Asian ancestry to investigate associations of predicted miRNA expression with CRC risk. Perturbation experiments separately by promoting and inhibiting miRNAs expressions and further in vitro assays in both SW480 and HCT116 cells were conducted. At a Bonferroni-corrected threshold of P < 4.5 × 10-4, we identified two putative susceptibility miRNAs, miR-1307-5p and miR-192-3p, located in regions more than 500 kb away from any GWAS-identified risk variants in CRC. We observed that a high predicted expression of miR-1307-5p was associated with increased CRC risk, while a low predicted expression of miR-192-3p was associated with increased CRC risk. Our experimental results further provide strong evidence of their susceptible roles by showing that miR-1307-5p and miR-192-3p play a regulatory role, respectively, in promoting and inhibiting CRC cell proliferation, migration, and invasion, which was consistently observed in both SW480 and HCT116 cells. Our study provides additional insights into the biological mechanisms underlying CRC development.

Regulation of the innate immune response and gut microbiome by p53.

p53 is a key tumor suppressor mutated in half of human cancers. In recent years, p53 was shown to regulate a wide variety of functions. From the transcriptome analysis of 24 tissues of irradiated mice, we identified 553 genes markedly induced by p53. Gene Ontology (GO) enrichment analysis found that the most associated biological process was innate immunity. 16S rRNA-seq analysis revealed that Akkermansia, which has anti-inflammatory properties and is involved in the regulation of intestinal barrier integrity, was decreased in p53-knockout (p53-/- ) mice after radiation. p53-/- mice were susceptible to radiation-induced GI toxicity and had a significantly shorter survival time than p53-wild-type (p53+/+ ) mice following radiation. However, administration of antibiotics resulted in a significant improvement in survival and protection against GI toxicity. Mbl2 and Lcn2, which have antimicrobial activity, were identified to be directly transactivated by p53 and secreted by liver into the circulatory system. We also found the expression of MBL2 and LCN2 was decreased in liver cancer tissues with p53 mutations compared with those without p53 mutations. These results indicate that p53 is involved in shaping the gut microbiome through its downstream targets related to the innate immune system, thus protecting the intestinal barrier.

East Asian-specific and cross-ancestry genome-wide meta-analyses provide mechanistic insights into peptic ulcer disease.

Peptic ulcer disease (PUD) refers to acid-induced injury of the digestive tract, occurring mainly in the stomach (gastric ulcer (GU)) or duodenum (duodenal ulcer (DU)). In the present study, we conducted a large-scale, cross-ancestry meta-analysis of PUD combining genome-wide association studies with Japanese and European studies (52,032 cases and 905,344 controls), and discovered 25 new loci highly concordant across ancestries. An examination of GU and DU genetic architecture demonstrated that GUs shared the same risk loci as DUs, although with smaller genetic effect sizes and higher polygenicity than DUs, indicating higher heterogeneity of GUs. Helicobacter pylori (HP)-stratified analysis found an HP-related host genetic locus. Integrative analyses using bulk and single-cell transcriptome profiles highlighted the genetic factors of PUD being enriched in the highly expressed genes in stomach tissues, especially in somatostatin-producing D cells. Our results provide genetic evidence that gastrointestinal cell differentiations and hormone regulations are critical in PUD etiology.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

Genotype imputation accuracy and the quality metrics of the minor ancestry in multi-ancestry reference panels.

Large-scale imputation reference panels are currently available and have contributed to efficient genome-wide association studies through genotype imputation. However, whether large-size multi-ancestry or small-size population-specific reference panels are the optimal choices for under-represented populations continues to be debated. We imputed genotypes of East Asian (180k Japanese) subjects using the Trans-Omics for Precision Medicine reference panel and found that the standard imputation quality metric (Rsq) overestimated dosage r2 (squared correlation between imputed dosage and true genotype) particularly in marginal-quality bins. Variance component analysis of Rsq revealed that the increased imputed-genotype certainty (dosages closer to 0, 1 or 2) caused upward bias, indicating some systemic bias in the imputation. Through systematic simulations using different template switching rates (θ value) in the hidden Markov model, we revealed that the lower θ value increased the imputed-genotype certainty and Rsq; however, dosage r2 was insensitive to the θ value, thereby causing a deviation. In simulated reference panels with different sizes and ancestral diversities, the θ value estimates from Minimac decreased with the size of a single ancestry and increased with the ancestral diversity. Thus, Rsq could be deviated from dosage r2 for a subpopulation in the multi-ancestry panel, and the deviation represents different imputed-dosage distributions. Finally, despite the impact of the θ value, distant ancestries in the reference panel contributed only a few additional variants passing a predefined Rsq threshold. We conclude that the θ value substantially impacts the imputed dosage and the imputation quality metric value.

Developing and validating polygenic risk scores for colorectal cancer risk prediction in East Asians.

Several polygenic risk scores (PRSs) have been developed to predict the risk of colorectal cancer (CRC) in European descendants. We used genome-wide association study (GWAS) data from 22 702 cases and 212 486 controls of Asian ancestry to develop PRSs and validated them in two case-control studies (1454 Korean and 1736 Chinese). Eleven PRSs were derived using three approaches: GWAS-identified CRC risk SNPs, CRC risk variants identified through fine-mapping of known risk loci and genome-wide risk prediction algorithms. Logistic regression was used to estimate odds ratios (ORs) and area under the curve (AUC). PRS115-EAS , a PRS with 115 GWAS-reported risk variants derived from East-Asian data, validated significantly better than PRS115-EUR derived from European descendants. In the Korea validation set, OR per SD increase of PRS115-EAS was 1.63 (95% CI = 1.46-1.82; AUC = 0.63), compared with OR of 1.44 (95% CI = 1.29-1.60, AUC = 0.60) for PRS115-EUR . PRS115-EAS/EUR derived using meta-analysis results of both populations slightly improved the AUC to 0.64. Similar but weaker associations were found in the China validation set. Individuals among the highest 5% of PRS115-EAS/EUR have a 2.52-fold elevated CRC risk compared with the medium (41-60th) risk group and have a 12% to 20% risk of developing CRC by age 85. PRSs constructed using results from fine-mapping and genome-wide algorithms did not perform as well as PRS115-EAS and PRS115-EAS/EUR in risk prediction, possibly due to a small sample size. Our results indicate that CRC PRSs are promising in predicting CRC risk in East Asians and highlights the importance of using population-specific data to build CRC risk prediction models.

Large-scale Integrated Analysis of Genetics and Metabolomic Data Reveals Potential Links Between Lipids and Colorectal Cancer Risk.

BACKGROUND: The etiology of colorectal cancer is not fully understood. METHODS: Using genetic variants and metabolomics data including 217 metabolites from the Framingham Heart Study (n = 1,357), we built genetic prediction models for circulating metabolites. Models with prediction R2 > 0.01 (Nmetabolite = 58) were applied to predict levels of metabolites in two large consortia with a combined sample size of approximately 46,300 cases and 59,200 controls of European and approximately 21,700 cases and 47,400 controls of East Asian (EA) descent. Genetically predicted levels of metabolites were evaluated for their associations with colorectal cancer risk in logistic regressions within each racial group, after which the results were combined by meta-analysis. RESULTS: Of the 58 metabolites tested, 24 metabolites were significantly associated with colorectal cancer risk [Benjamini-Hochberg FDR (BH-FDR) < 0.05] in the European population (ORs ranged from 0.91 to 1.06; P values ranged from 0.02 to 6.4 × 10-8). Twenty one of the 24 associations were replicated in the EA population (ORs ranged from 0.26 to 1.69, BH-FDR < 0.05). In addition, the genetically predicted levels of C16:0 cholesteryl ester was significantly associated with colorectal cancer risk in the EA population only (OREA: 1.94, 95% CI, 1.60-2.36, P = 2.6 × 10-11; OREUR: 1.01, 95% CI, 0.99-1.04, P = 0.3). Nineteen of the 25 metabolites were glycerophospholipids and triacylglycerols (TAG). Eighteen associations exhibited significant heterogeneity between the two racial groups (PEUR-EA-Het < 0.005), which were more strongly associated in the EA population. This integrative study suggested a potential role of lipids, especially certain glycerophospholipids and TAGs, in the etiology of colorectal cancer. CONCLUSIONS: This study identified potential novel risk biomarkers for colorectal cancer by integrating genetics and circulating metabolomics data. IMPACT: The identified metabolites could be developed into new tools for risk assessment of colorectal cancer in both European and EA populations.

Genetic susceptibility to hepatocellular carcinoma in chromosome 22q13.31, findings of a genome-wide association study.

BACKGROUND AND AIM: Chronic hepatitis C virus (HCV) infection, long-term alcohol use, cigarette smoking, and obesity are the major risk factors for hepatocellular carcinoma (HCC) in the United States, but the disease risk varies substantially among individuals with these factors, suggesting host susceptibility to and gene-environment interactions in HCC. To address genetic susceptibility to HCC, we conducted a genome-wide association study (GWAS). METHODS: Two case-control studies on HCC were conducted in the United States. DNA samples were genotyped using the Illumian microarray chip with over 710 000 single nucleotide polymorphisms (SNPs). We compared these SNPs between 705 HCC cases and 1455 population controls for their associations with HCC and verified our findings in additional studies. RESULTS: In this GWAS, we found that two SNPs were associated with HCC at P < 5E-8 and six SNPs at P < 5E-6 after adjusting for age, sex, and the top three principal components (PCs). Five of the SNPs in chromosome 22q13.31, three in PNPLA3 (rs2281135, rs2896019, and rs4823173) and two in SAMM50 (rs3761472, rs3827385), were replicated in a small US case-control study and a cohort study in Singapore. The associations remained significant after adjusting for body mass index and HCV infection. Meta-analysis of multiple datasets indicated that these SNPs were significantly associated with HCC. CONCLUSIONS: SNPs in PNPLA3 and SAMM50 are known risk loci for nonalcoholic fatty liver disease (NAFLD) and are suspected to be associated with HCC. Our GWAS demonstrated the associations of these SNPs with HCC in a US population. Biological mechanisms underlying the relationship remain to be elucidated.