Understanding risk factors for endometrial cancer in young women.

BACKGROUND: The American Cancer Society recommends physicians inform average risk women about endometrial cancer (EC) risk on reaching menopause, but new diagnoses are rising fastest in women <50 years. Educating these women about EC risks requires knowledge of risk factors. However, EC in young women is rare and challenging to study in single study populations. METHODS: We included 13,846 incident EC patients (1,639 < 50 years) and 30,569 matched control individuals from the Epidemiology of Endometrial Cancer Consortium. We used generalized linear models to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for 6 risk factors and EC risk. We created a risk score to evaluate the combined associations and population attributable fractions of these factors. RESULTS: In younger and older women, we observed positive associations with BMI and diabetes, and inverse associations with age at menarche, oral contraceptive use, and parity. Current smoking was associated with reduced risk only in women ≥50 years (PHet<0.01). BMI was the strongest risk factor [OR≥35 vs <25 kg/m2=5.57 (95% CI:4.33-7.16) for <50 years; OR≥35 vs <25 kg/m2=4.68 (95% CI : 4.30-5.09) for ≥50 years; PHet=0.14]. Possessing ≥4 risk factors was associated with ∼9-fold increased risk in women <50 years and ∼4-fold increased risk in women ≥50 years (PHet<0.01). Together, 59.1% of ECs in women <50 and 55.6% in women ≥50 were attributable to these factors. CONCLUSIONS: Our data confirm younger and older women share common EC risk factors. Early educational efforts centered on these factors may help mitigate the rising EC burden in young women.

Integrating muti-omics data to identify tissue-specific DNA methylation biomarkers for cancer risk.

The relationship between tissue-specific DNA methylation and cancer risk remains inadequately elucidated. Leveraging resources from the Genotype-Tissue Expression consortium, here we develop genetic models to predict DNA methylation at CpG sites across the genome for seven tissues and apply these models to genome-wide association study data of corresponding cancers, namely breast, colorectal, renal cell, lung, ovarian, prostate, and testicular germ cell cancers. At Bonferroni-corrected P < 0.05, we identify 4248 CpGs that are significantly associated with cancer risk, of which 95.4% (4052) are specific to a particular cancer type. Notably, 92 CpGs within 55 putative novel loci retain significant associations with cancer risk after conditioning on proximal signals identified by genome-wide association studies. Integrative multi-omics analyses reveal 854 CpG-gene-cancer trios, suggesting that DNA methylation at 309 distinct CpGs might influence cancer risk through regulating the expression of 205 unique cis-genes. These findings substantially advance our understanding of the interplay between genetics, epigenetics, and gene expression in cancer etiology.

ERBB2 amplification in gastric cancer: a genomic insight into ethnic disparities.

Overall, gastric adenocarcinoma (GC) incidence rates have declined in recent years, but racial/ethnic disparities persist. Individuals who identify as Hispanic/Spanish/Latino are diagnosed with GC at younger ages and have poorer outcomes than non-Hispanic individuals. However, our understanding of GC biology across racial/ethnic groups remains limited. We assessed tumor genomic patterns by race/ethnicity among 1019 patients with primary GC in the AACR Project GENIE Consortium. Hispanic individuals presented with significantly higher rates of ERBB2/HER2 amplification vs other racial/ethnic groups (Hispanic: 13.9% vs 9.8% non-Hispanic White, 8.1% non-Hispanic Asian, and 11.0% non-Hispanic Black; p < .001, FDR adjusted q < 0.001). Hispanic patients also had higher odds of an ERBB2 amplification vs non-Hispanic whites in adjusted models (OR = 2.52, 95%CI = 1.20-5.33, p = .015). These findings underscore the important role of genomic factors in GC disparities. Ensuring equitable access to genomic profiling and targeted therapies, such as trastuzumab for HER2-overexpressing GC, is a promising avenue to mitigate GC disparities and improve outcomes.

Triterpenoid and trichothecenes from a rhizospheric soil-derived Paramyrothecium sp. KMU22107.

A new phthalide derivative named paramlyktone (1) and a new arborinane-type triterpenoid named paramyrpenoid (2), together with ten previously described trichothecenes derivatives (3-12) were isolated and identified from a rhizospheric soil-derived Paramyrothecium sp. KMU22107 associated with Delphinium yunnanense. Their structural elucidation was achieved by the comprehensive analysis of spectroscopic data and comparison with literature values. Notably, paramyrpenoid (2) was the first example of an arborinane-type triterpenoid with a double bond at Δ12(13) and an additional methyl motif at C-8. This was the first report of arborinane-type triterpenoids from a fungus belonging to Paramyrothecium genus. In pharmacological studies, paramyrpenoid (2) demonstrated significant cytotoxic activity against the HL-60, SW480, A-549, MDA-MB-231 and SMMC-7721 cell lines, with IC50 values from 2.0 to 16.1 µM. Compounds 1 and 2 were also evaluated for anti-inflammatory, anti-acetylcholinesterase (AChE), and protein tyrosine phosphatase 1B (PTP1B) inhibitory activities in vitro.

Large-scale integration of omics and electronic health records to identify potential risk protein biomarkers and therapeutic drugs for cancer prevention and intervention.

Identifying risk protein targets and their therapeutic drugs is crucial for effective cancer prevention. Here, we conduct integrative and fine-mapping analyses of large genome-wide association studies data for breast, colorectal, lung, ovarian, pancreatic, and prostate cancers, and characterize 710 lead variants independently associated with cancer risk. Through mapping protein quantitative trait loci (pQTL) for these variants using plasma proteomics data from over 75,000 participants, we identify 365 proteins associated with cancer risk. Subsequent colocalization analysis identifies 101 proteins, including 74 not reported in previous studies. We further characterize 36 potential druggable proteins for cancers or other disease indications. Analyzing >3.5 million electronic health records, we uncover five drugs (Haloperidol, Trazodone, Tranexamic Acid, Haloperidol, and Captopril) associated with increased cancer risk and two drugs (Caffeine and Acetazolamide) linked to reduced colorectal cancer risk. This study offers novel insights into therapeutic drugs targeting risk proteins for cancer prevention and intervention.

Case-Case Genome-Wide Analyses Identify Subtype-Informative Variants That Confer Risk for Breast Cancer.

Breast cancer includes several subtypes with distinct characteristic biological, pathologic, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate the development of improved prevention and treatment approaches. In this study, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case genome-wide association study (CC-GWAS) approach to summary statistics data of the Breast Cancer Association Consortium. The approach identified 13 statistically significant loci and eight suggestive loci, the majority of which were identified from comparisons between triple-negative breast cancer (TNBC) and luminal A breast cancer. Associations of lead variants in 12 loci remained statistically significant after accounting for previously reported breast cancer susceptibility variants, among which, two were genome-wide significant. Fine mapping implicated putative functional/causal variants and risk genes at several loci, e.g., 3q26.31/TNFSF10, 8q22.3/NACAP1/GRHL2, and 8q23.3/LINC00536/TRPS1, for TNBC as compared with luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates the enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to two-fold increased risk of being diagnosed with TNBC instead of luminal cancers. The CC-GWAS PRS remained statistically significant after adjusting for TNBC PRS derived from traditional case-control GWAS in The Cancer Genome Atlas and the African Ancestry Breast Cancer Genetic Consortium. The CC-GWAS PRS was also associated with overall survival and disease-specific survival among patients with breast cancer. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC. Significance: The discovery of subtype-informative genetic risk variants for breast cancer advances our understanding of the etiologic heterogeneity of breast cancer, which could accelerate the identification of targets and personalized strategies for prevention and treatment.

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

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

Large-Scale Alternative Polyadenylation-Wide Association Studies to Identify Putative Cancer Susceptibility Genes.

Alternative polyadenylation (APA) modulates mRNA processing in the 3'-untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. In this study, we conducted large APA-wide association studies to investigate associations between APA levels and cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA sequencing data from 1,337 samples from the Genotype-Tissue Expression project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European ancestry populations: breast, ovarian, prostate, colorectal, lung, and pancreatic cancers. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3'-UTR APA quantitative trait loci and colocalization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3'-UTR variants demonstrated that the risk alleles of 3'-UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the posttranscriptional activities of their target genes compared with reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers. Significance: Systematic evaluation of associations of alternative polyadenylation with cancer risk reveals 58 putative susceptibility genes, highlighting the contribution of genetically regulated alternative polyadenylation of 3'UTRs to genetic susceptibility to cancer.

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

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

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Circulating microRNAs in association with pancreatic cancer risk within 5 years.

Early detection is critical for improving pancreatic cancer prognosis. Our study aims to identify circulating microRNAs (miRNAs) associated with pancreatic cancer risk. The two-stage study used plasma samples collected ≤5 years prior to cancer diagnosis, from case-control studies nested in five prospective cohort studies. The discovery stage included 185 case-control pairs from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Replication stage samples comprised 277 pairs from Shanghai Women's Health Study/Shanghai Men's Health Study, Southern Community Cohort Study, and Multiethnic Cohort Study. Seven hundred and ninety-eight miRNAs were measured using the NanoString nCounter Analysis System. Odds ratios (OR) and 95% confidence intervals (CI) for per 10% change in miRNAs in association with pancreatic cancer risk were derived from conditional logistic regression analysis in discovery and replication studies, separately, and then meta-analyzed. Stratified analysis was conducted by age at diagnosis (<65/≥65 years) and time interval between sample collection and diagnosis (≤2/>2 years). In the discovery stage, 120 risk associated miRNAs were identified at p < .05. Three were validated in the replication stage: hsa-miR-199a-3p/hsa-miR-199b-3p, hsa-miR-767-5p, and hsa-miR-191-5p, with respective ORs (95% CI) being 0.89 (0.84-0.95), 1.08 (1.02-1.13), and 0.90 (0.85-0.95). Five additional miRNAs, hsa-miR-640, hsa-miR-874-5p, hsa-miR-1299, hsa-miR-22-3p, and hsa-miR-449b-5p, were validated among patients diagnosed at ≥65 years, with OR (95% CI) of 1.23 (1.09-1.39), 1.33 (1.16-1.52), 1.25 (1.09-1.43), 1.28 (1.12-1.46), 0.76 (0.65-0.89), and 1.22 (1.07-1.39), respectively. The miRNA targets were enriched in pancreatic carcinogenesis/progression-related pathways. Our study suggests that circulating miRNAs may identify individuals at high risk for pancreatic cancer ≤5 years prior to diagnosis, indicating its potential utility in cancer screening and surveillance.

Drug-target Mendelian randomization revealed a significant association of genetically proxied metformin effects with increased prostate cancer risk.

The association between metformin use and risk of prostate cancer remains controversial, while data from randomized trials is lacking. We aim to evaluate the association of genetically proxied metformin effects with prostate cancer risk using a drug-target Mendelian randomization (MR) approach. Summary statistics for prostate cancer were obtained from the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome Consortium (79,148 cases and 61,106 controls). Cis-expression quantitative trait loci (cis-eQTL) variants in the gene targets of metformin were identified in the GTEx project and eQTLGen consortium. We also obtained male-specific genome-wide association study data for type 2 diabetes, body mass index (BMI), total testosterone, bioavailable testosterone, estradiol, and sex hormone binding globulin for mediation analysis. Inverse-variance weighted (IVW) regression, weighted median, MR-Egger regression, and MR-PRESSO were performed in the main MR analysis. Multivariable MR was used to identify potential mediators and genetic colocalization analysis was performed to assess any shared genetic basis between two traits of interest. We found that genetically proxied metformin effects (1-SD HbA1c reduction, equivalent to 6.75 mmol/mol) were associated with higher risk of prostate cancer (odds ratioIVW [ORIVW]: 1.55, 95% confidence interval, CI: 1.23-1.96, p = 3.0 × 10-3). Two metformin targets, mitochondrial complex I (ORIVW: 1.48, 95% CI: 1.07-2.03, p = 0.016) and gamma-secretase complex (ORIVW: 2.58, 95%CI :1.47-4.55, p = 0.001), showed robust associations with prostate cancer risk, and their effects were partly mediated through BMI (16.4%) and total testosterone levels (34.3%), respectively. These results were further supported by colocalization analysis that expressions of NDUFA13 and BMI, APH1A, and total testosterone may be influenced by shared genetic factors, respectively. In summary, our study indicated that genetically proxied metformin effects may be associated with an increased risk of prostate cancer. Repurposing metformin for prostate cancer prevention in general populations is not supported by our findings.

One-year efficacy of myopia control by the defocus distributed multipoint lens: a multicentric randomised controlled trial.

AIMS: To report the 1-year results of the efficacy of a defocus distributed multipoint (DDM) lens in controlling myopia progression in a multicentre, randomised controlled trial. METHODS: Overall, 168 children aged 6-13 years were recruited and randomly assigned to wear a DDM lens (n=84) or single-vision (SV) lens (n=84) in three centres. Cycloplegic autorefraction (spherical equivalent refraction (SER)) and axial length (AL) were measured. Linear mixed model analysis was performed to compare between-group SER and AL changes. Logistic regression analysis was used to analyse the between-group difference in rapid myopia progression (SER increase≥0.75 D per year or AL growth≥0.40 mm per year). RESULTS: After 1 year, mean changes in SER were significantly lower in the DDM group (-0.47±0.37 D) than in the SV group (-0.71±0.42 D) (p<0.001). Similarly, mean changes in AL were significantly lower in the DDM group (0.21±0.17 mm) than in the SV group (0.34±0.16 mm) (p<0.001). After adjusting for age, sex, daily wearing time and parental myopia, rapid myopia progression risk was higher in the SV group than in the DDM group (OR=3.51, 95% CI: 1.77 to 6.99), especially for children who wore a lens for >12 hours per day, boys and younger children (6-9 years) with ORs (95% CIs) of 10.82 (3.22 to 36.37), 5.34 (1.93 to 14.78) and 8.73 (2.6 to 29.33), respectively. CONCLUSIONS: After 1 year, DDM lenses effectively retarded myopia progression in children. Longer daily wearing time of DDM lens improved the efficacy of myopia control. Future long-term studies are needed for validation. TRIAL REGISTRATION NUMBER: NCT05340699.

Large-Scale Cancer Genomic Analysis Reveals Significant Disparities between Microsatellite Instability and Tumor Mutational Burden.

BACKGROUND: Microsatellite instability (MSI) and tumor mutational burden (TMB) are predictive biomarkers for pan-cancer immunotherapy. The interrelationship between MSI-high (MSI-H) and TMB-high (TMB-H) in human cancers and their predictive value for immunotherapy in lung cancer remain unclear. METHODS: We analyzed somatic mutation data from the Genomics Evidence Neoplasia Information Exchange (n = 46,320) to determine the relationship between MSI-H and TMB-H in human cancers using adjusted multivariate regression models. Patient survival was examined using the Cox proportional hazards model. The association between MSI and genetic mutations was assessed. RESULTS: Patients (31-89%) with MSI-H had TMB-low phenotypes across 22 cancer types. Colorectal and stomach cancers showed the strongest association between TMB and MSI. TMB-H patients with lung cancer who received immunotherapy exhibited significantly higher overall survival [HR, 0.61; 95% confidence interval (CI), 0.44-0.86] and progression-free survival (HR, 0.65; 95% CI, 0.47-0.91) compared to the TMB-low group; no significant benefit was observed in the MSI-H group. Patients with TMB and MSI phenotypes showed further improvement in overall survival and PFS. We identified several mutated genes associated with MSI-H phenotypes, including known mismatch repair genes and novel mutated genes, such as ARID1A and ARID1B. CONCLUSIONS: Our results demonstrate that TMB-H and/or a combination of MSI-H can serve as biomarkers for immunotherapies in lung cancer. IMPACT: These findings suggest that distinct or combined biomarkers should be considered for immunotherapy in human cancers because notable discrepancies exist between MSI-H and TMB-H across different cancer types.

An expression-directed linear mixed model discovering low-effect genetic variants.

Detecting genetic variants with low-effect sizes using a moderate sample size is difficult, hindering downstream efforts to learn pathology and estimating heritability. In this work, by utilizing informative weights learned from training genetically predicted gene expression models, we formed an alternative approach to estimate the polygenic term in a linear mixed model. Our linear mixed model estimates the genetic background by incorporating their relevance to gene expression. Our protocol, expression-directed linear mixed model, enables the discovery of subtle signals of low-effect variants using moderate sample size. By applying expression-directed linear mixed model to cohorts of around 5,000 individuals with either binary (WTCCC) or quantitative (NFBC1966) traits, we demonstrated its power gain at the low-effect end of the genetic etiology spectrum. In aggregate, the additional low-effect variants detected by expression-directed linear mixed model substantially improved estimation of missing heritability. Expression-directed linear mixed model moves precision medicine forward by accurately detecting the contribution of low-effect genetic variants to human diseases.

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

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

Novel insights into genetic susceptibility for colorectal cancer from transcriptome-wide association and functional investigation.

BACKGROUND: Transcriptome-wide association studies have been successful in identifying candidate susceptibility genes for colorectal cancer (CRC). To strengthen susceptibility gene discovery, we conducted a large transcriptome-wide association study and an alternative splicing transcriptome-wide association study in CRC using improved genetic prediction models and performed in-depth functional investigations. METHODS: We analyzed RNA-sequencing data from normal colon tissues and genotype data from 423 European descendants to build genetic prediction models of gene expression and alternative splicing and evaluated model performance using independent RNA-sequencing data from normal colon tissues of the Genotype-Tissue Expression Project. We applied the verified models to genome-wide association studies (GWAS) summary statistics among 58 131 CRC cases and 67 347 controls of European ancestry to evaluate associations of genetically predicted gene expression and alternative splicing with CRC risk. We performed in vitro functional assays for 3 selected genes in multiple CRC cell lines. RESULTS: We identified 57 putative CRC susceptibility genes, which included the 48 genes from transcriptome-wide association studies and 15 genes from splicing transcriptome-wide association studies, at a Bonferroni-corrected P value less than .05. Of these, 16 genes were not previously implicated in CRC susceptibility, including a gene PDE7B (6q23.3) at locus previously not reported by CRC GWAS. Gene knockdown experiments confirmed the oncogenic roles for 2 unreported genes, TRPS1 and METRNL, and a recently reported gene, C14orf166. CONCLUSION: This study discovered new putative susceptibility genes of CRC and provided novel insights into the biological mechanisms underlying CRC development.

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.