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

Breast cancer includes several subtypes with distinct characteristic biological, pathological, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate development of improved prevention and treatment approaches. Here, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case GWAS (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 to luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to 2-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 breast cancer patients. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC.

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.

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. Herein, we conducted large alternative polyadenylation-wide association studies (APA-WAS) to investigate associations of APA levels with 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, including breast, ovary, prostate, colorectum, lung, and pancreas. 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 co-localization 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 post-transcriptional activities of their target genes compared to 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.

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

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

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

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

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.

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.

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.

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.

Regulation of secondary metabolites in the endophytic fungus Penicillium sp. KMU18029 by the chemical epigenetic modifier 5-azacitidine.

The chemical epigenetic modifier 5-azacitidine (5-Aza C), a DNA methyltransferase inhibitor, was used to manipulate the endophytic fungus Penicillium sp. KMU18029. From its rice fermentation extract, a new polyketone compound (3S,4R)-3,4,8-trihydroxy-6-methyl-3,4-dihydronaphthalen-1(2H)-one (1), along with 13 known compounds, 3,4,8-trihydroxy-6-(hydroxymethyl)-3,4-dihydronaphthalen-1(2H)-one (2), decaturin B (3), 15-hydroxydecaturin A (4), oxalicine A (5), pileotin A (6), pyrandecarurin A (7), decaturenol A (8), decaturenoid (9), penisarins A (10), oxaline (11), (4E,8E)-N-D-2'-hydroxyocta-decanoyl-1-O-ß-D-glycopy-ranosyl-9-methyl-4,8-sphingadienine (12), ergosterol (13) and stigma-5-en-3-O-ß-glucoside (14), were separated. Among the known compounds, 2, 7, 12 and 14 were not found in our previous research on this strain. The structure of the new compound was identified by spectroscopic techniques such as HR-ESIMS, 1D NMR, 2D NMR and CD. Furthermore, all the isolated compounds were tested for their antimicrobial activities, and only compounds 1, 2 and 11 showed weak activities against S. aureus, with MICs of 128 µg/mL.

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.

High intake of cruciferous vegetables reduces the risk of gastrointestinal cancers: results from observational studies.

The relationship between cruciferous vegetables (CV) and the risk of gastrointestinal (GI) cancers has been extensively investigated. However, epidemiologic investigations have produced inconsistent results. This meta-analysis investigated the association between CV intake and the risk of GI cancers. Due to the heterogeneity, fixed- or random-effects models were used for the analyses. The final analysis included 81 articles covering 89 studies. In comparison to the lowest consumption categories, the highest consumption categories of CV were associated with a lower risk for all GI cancers [rate ratio (RR): 0.81, 95% confidence interval (95% CI) 0.76-0.87]. Compared to a CV intake of 75 g/day, subjects with CV intake <75 g/day experienced a 7% reduction in risk (RR: 0.93; 95% CI: 0.84-0.96) for each 50 g increase in consumption. A negative correlation was identified between CV intake and the risk of esophageal cancer, gastric cancer, liver cancer, pancreatic cancer, and colorectal cancer (CRC), but not gallbladder cancer (RR: 0.70; 95% CI: 0.38-1.27). High intake of broccoli and cabbage was associated with a decreased risk of gastric cancer (RR: 0.64; 95% CI: 0.47-0.87) and gallbladder cancer (RR: 0.46; 95% CI: 0.29-0.75). These results confirm the association between high intake of CV with a reduced risk of GI cancers.

cLD: Rare-variant linkage disequilibrium between genomic regions identifies novel genomic interactions.

Linkage disequilibrium (LD) is a fundamental concept in genetics; critical for studying genetic associations and molecular evolution. However, LD measurements are only reliable for common genetic variants, leaving low-frequency variants unanalyzed. In this work, we introduce cumulative LD (cLD), a stable statistic that captures the rare-variant LD between genetic regions, which reflects more biological interactions between variants, in addition to lack of recombination. We derived the theoretical variance of cLD using delta methods to demonstrate its higher stability than LD for rare variants. This property is also verified by bootstrapped simulations using real data. In application, we find cLD reveals an increased genetic association between genes in 3D chromatin interactions, a phenomenon recently reported negatively by calculating standard LD between common variants. Additionally, we show that cLD is higher between gene pairs reported in interaction databases, identifies unreported protein-protein interactions, and reveals interacting genes distinguishing case/control samples in association studies.

An exploratory study of knowledge, attitudes, and practices toward HPV associated anal cancer among Pakistani population.

Background: Anal cancer, mainly attributed to human papillomavirus (HPV) infection, is rising in prevalence among the general population in Pakistan. This study aimed to examine the knowledge, attitudes, and practices (KAP) towards anal cancer screening and HPV of the general population in Pakistan. Method: We surveyed anal cancer KAP using social media and snowball sampling from December 2022 to May 2023. The questionnaire had 16 knowledge, 12 attitudes, 6 practice questions, and socio-demographic variables. We applied validity criteria for inclusion and exclusion and used cutoffs ≥50% for each KAP category. We analyzed data in R with Guttman's λ2 for reliability, did univariate and bivariate analysis, and reported frequencies, percentages, p-values, coefficients, odds ratios, and 95% confidence intervals. Results: We surveyed 1620 people and discovered low awareness of HPV and anal cancer causes prevention, and screening (11%-24%), high stigma and embarrassment for screening (54%-70%), strong moral beliefs (89%), condom nonuse (91%), and low engagement in health services and programs (9.1%-14%). Knowledge (75.23%, OR = 1.0984, p = 0.05) was shaped by socio-demographic factors, attitude, and practice, with higher education enhancing knowledge (OR = 1.0984, p = 0.05). Attitude (78.45%, OR = 6.6052, p< 0.001) was influenced by socio-demographic factors, practice, and knowledge as well. Younger females, single, unemployed, students, living with more family members, earning more income, and residing in Islamabad had a more positive attitude (ORs from 1.0115 to 6.6052, p< 0.05), while religion did not affect attitude (p = 0.51). Practice (9.16%, OR = 0.1820, p< 0.001) was determined by socio-demographic factors, knowledge, and attitude. Older males, employed teachers, living with more family members, earning less income, and residing in Islamabad had better practice (ORs from 0.1323 to 3.8431, p< 0.05), but marital status and religion did not influence practice (p > 0.05). Conclusion: Pakistani young adults need more education, awareness, health services, and programs on HPV and anal cancer, as they have low awareness, high stigma, and socio-cultural challenges. In addition, it is recommended for more research and policy initiatives are needed to address socio-cultural factors and increase anal Pap to overcome anal cancer.