Genome-wide association studies of cancer: current insights and future perspectives

CD147 TagSNP is associated with the vulnerability to lung cancer in the Chinese population: a case-control study

BACKGROUND

Lung cancer, with its high morbidity and mortality, presents a major significant public health challenge. CD147, linked to cancer progression and metastasis, is a promising therapeutic target, including for lung cancer. The genetic variation may influence the expression of the gene and consequently the risk of lung cancer. This study aims to investigate single nucleotide polymorphisms (SNPs) in CD147 to understand their association with the risk of developing lung cancer in the Han Chinese population.

METHODS

A hospital-based case-control investigation was conducted, enrolling 700 lung cancer patients and 700 cancer-free controls. TagSNPs were selected using Haploview v4.2, and genotype data from the 1000 Genomes Project database were utilized. The selected SNPs (rs28992491, rs67945626, and rs79361899) within the CD147 gene were evaluated using the improved multiple ligation detection reaction method. Statistical analysis included chi-square tests, logistic regression models, and interaction analyses.

RESULTS

Baseline characteristics of the study population showed no significant differences in gender distribution between cases and controls, but there was a notable difference in smoking rates. No significant associations were found between the three TagSNPs and lung cancer susceptibility in the codominant model. However, stratification analyses revealed interesting findings. Among females, the rs79361899 AA/AG genotype was associated with an increased risk of lung cancer. In individuals aged ≥ 65 years old, the rs28992491 GG and rs79361899 AA genotypes were linked to a higher susceptibility. Furthermore, an interaction analysis demonstrated significant genotype × gender interactions in the rs79361899 recessive model, indicating an increased lung cancer risk in female carriers of the heterozygous or homozygous polymorphic genotype.

CONCLUSIONS

CD147 polymorphisms play an important role in lung cancer development, particularly in specific subgroup of age and gender. These findings highlight the significance of incorporating genetic variations and their interactions with demographic factors in comprehending the intricate etiology of lung cancer.

Co-expression in tissue-specific gene networks links genes in cancer-susceptibility loci to known somatic driver genes

BACKGROUND

The genetic background of cancer remains complex and challenging to integrate. Many somatic mutations within genes are known to cause and drive cancer, while genome-wide association studies (GWAS) of cancer have revealed many germline risk factors associated with cancer. However, the overlap between known somatic driver genes and positional candidate genes from GWAS loci is surprisingly small. We hypothesised that genes from multiple independent cancer GWAS loci should show tissue-specific co-regulation patterns that converge on cancer-specific driver genes.

RESULTS

We studied recent well-powered GWAS of breast, prostate, colorectal and skin cancer by estimating co-expression between genes and subsequently prioritising genes that show significant co-expression with genes mapping within susceptibility loci from cancer GWAS. We observed that the prioritised genes were strongly enriched for cancer drivers defined by COSMIC, IntOGen and Dietlein et al. The enrichment of known cancer driver genes was most significant when using co-expression networks derived from non-cancer samples of the relevant tissue of origin.

CONCLUSION

We show how genes within risk loci identified by cancer GWAS can be linked to known cancer driver genes through tissue-specific co-expression networks. This provides an important explanation for why seemingly unrelated sets of genes that harbour either germline risk factors or somatic mutations can eventually cause the same type of disease.

Machine Learning of Three-Dimensional Protein Structures to Predict the Functional Impacts of Genome Variation

Research in the human genome sciences generates a substantial amount of genetic data for hundreds of thousands of individuals, which concomitantly increases the number of variants of unknown significance (VUS). Bioinformatic analyses can successfully reveal rare variants and variants with clear associations with disease-related phenotypes. These studies have had a significant impact on how clinical genetic screens are interpreted and how patients are stratified for treatment. There are few, if any, computational methods for variants comparable to biological activity predictions. To address this gap, we developed a machine learning method that uses protein three-dimensional structures from AlphaFold to predict how a variant will influence changes to a gene's downstream biological pathways. We trained state-of-the-art machine learning classifiers to predict which protein regions will most likely impact transcriptional activities of two proto-oncogenes, nuclear factor erythroid 2 (NFE2L2)-related factor 2 (NRF2) and c-Myc. We have identified classifiers that attain accuracies higher than 80%, which have allowed us to identify a set of key protein regions that lead to significant perturbations in c-Myc or NRF2 transcriptional pathway activities.

A Genome-Wide Association Study of Endoxifen Serum Concentrations and Adjuvant Tamoxifen Efficacy in Early-Stage Breast Cancer Patients

Tamoxifen is part of the standard of care of endocrine therapy for adjuvant treatment of breast cancer. However, survival outcomes with tamoxifen are highly variable. The concentration of endoxifen, the 30-100 times more potent metabolite of tamoxifen and bioactivated by the CYP2D6 enzyme, has been described as the most relevant metabolite of tamoxifen metabolism. A genome-wide association study (GWAS) was performed with the objective to identify genetic polymorphisms associated with endoxifen serum concentration levels and clinical outcome in early-stage breast cancer patients receiving tamoxifen. A GWAS was conducted in 608 women of the CYPTAM study (NTR1509/PMID: 30120701). Germline DNA and clinical and survival characteristics were readily available. Genotyping was performed on Infinium Global Screening Array (686,082 markers) and single nucleotide polymorphism (SNP) imputation by using 1000 Genomes. Relapse-free survival during tamoxifen (RFSt) was defined the primary clinical outcome. Endoxifen serum concentration was analyzed as a continuous variable. Several genetic variants reached genome-wide significance (P value: ≤5 × 10-8). Endoxifen concentrations analysis identified 430 variants, located in TCF20 and WBP2NL genes (chromosome 22), which are in strong linkage disequilibrium with CYP2D6 variants. In the RFSt analysis, several SNP were identified (LPP gene: rs77693286, HR 18.3, 95% CI: 15.2-21.1; rs6790761, OR 18.2, 95% CI: 15.5-21.1). Endoxifen concentrations have a strong association with the chromosome 22, which contains the CYP2D6 gene.

Genome-wide enhancer RNA profiling adds molecular links between genetic variation and human cancers

BACKGROUND

Dysregulation of enhancer transcription occurs in multiple cancers. Enhancer RNAs (eRNAs) are transcribed products from enhancers that play critical roles in transcriptional control. Characterizing the genetic basis of eRNA expression may elucidate the molecular mechanisms underlying cancers.

METHODS

Initially, a comprehensive analysis of eRNA quantitative trait loci (eRNAQTLs) was performed in The Cancer Genome Atlas (TCGA), and functional features were characterized using multi-omics data. To establish the first eRNAQTL profiles for colorectal cancer (CRC) in China, epigenomic data were used to define active enhancers, which were subsequently integrated with transcription and genotyping data from 154 paired CRC samples. Finally, large-scale case-control studies (34,585 cases and 69,544 controls) were conducted along with multipronged experiments to investigate the potential mechanisms by which candidate eRNAQTLs affect CRC risk.

RESULTS

A total of 300,112 eRNAQTLs were identified across 30 different cancer types, which exert their influence on eRNA transcription by modulating chromatin status, binding affinity to transcription factors and RNA-binding proteins. These eRNAQTLs were found to be significantly enriched in cancer risk loci, explaining a substantial proportion of cancer heritability. Additionally, tumor-specific eRNAQTLs exhibited high responsiveness to the development of cancer. Moreover, the target genes of these eRNAs were associated with dysregulated signaling pathways and immune cell infiltration in cancer, highlighting their potential as therapeutic targets. Furthermore, multiple ethnic population studies have confirmed that an eRNAQTL rs3094296-T variant decreases the risk of CRC in populations from China (OR = 0.91, 95%CI 0.88-0.95, P = 2.92 × 10-7) and Europe (OR = 0.92, 95%CI 0.88-0.95, P = 4.61 × 10-6). Mechanistically, rs3094296 had an allele-specific effect on the transcription of the eRNA ENSR00000155786, which functioned as a transcriptional activator promoting the expression of its target gene SENP7. These two genes synergistically suppressed tumor cell proliferation. Our curated list of variants, genes, and drugs has been made available in CancereRNAQTL ( http://canernaqtl.whu.edu.cn/#/ ) to serve as an informative resource for advancing this field.

CONCLUSION

Our findings underscore the significance of eRNAQTLs in transcriptional regulation and disease heritability, pinpointing the potential of eRNA-based therapeutic strategies in cancers.

Metabolites Associated with Polygenic Risk of Breast Cancer

While hundreds of germline genetic variants have been associated with breast cancer risk, the mechanisms underlying the impacts of most of these variants on breast cancer remain uncertain. Metabolomics may offer valuable insights into the mechanisms underlying genetic risks of breast cancer. Among 143 cancer-free female participants, we used linear regression analyses to explore associations between the genetic risk of breast cancer, as determined by a previously developed polygenic risk score (PRS) that included 266 single-nucleotide polymorphisms (SNPs), and 223 measures of metabolites obtained from blood samples using nuclear magnetic resonance (NMR). A false discovery rate of 10% was applied to account for multiple comparisons. PRS was statistically significantly associated with 45 metabolite measures. These were primarily measures of very low-density lipoproteins (VLDLs) and high-density lipoproteins (HDLs), including triglycerides, cholesterol, and phospholipids. For example, the strongest effect was observed with the percent ratio of medium VLDL triglycerides to total lipids (0.53 unit increase in mean-standardized ln-transformed percent ratio per unit increase in PRS; q = 0.1). While larger-scale studies are needed to confirm these results, this exploratory study presents biologically plausible findings that are consistent with previously reported associations between lipids and breast cancer risk. If confirmed, these lipids could be targeted for lifestyle and pharmaceutical interventions among women at increased genetic risk of breast cancer.

Association Between Polymorphisms in DNA Damage Repair Pathway Genes and Female Breast Cancer Risk

Breast cancer risk have been discussed to be associated with polymorphisms in genes as well as abnormal DNA damage repair function. This study aims to assess the relationship between genes single nucleotide polymorphisms (SNPs) related to DNA damage repair and female breast cancer risk in Chinese population. A case-control study containing 400 patients and 400 healthy controls was conducted. Genotype was identified using the sequence MassARRAY method and expression of estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor-2 (HER-2) in tumor tissues was analyzed by immunohistochemistry assay. The results revealed that ATR rs13091637 decreased breast cancer risk influenced by ER, PR (CT/TT vs. CC: adjusted odds ratio [OR] = 1.54, 95% confidence interval [CI]: 1.04-2.27, p = 0.032; CT/TT vs. CC: adjusted OR = 1.63, 95%CI: 1.14-2.35, p = 0.008) expression. Stratified analysis revealed that PALB2 rs16940342 increased breast cancer risk in response to menstrual status (AG/GG vs. AA: adjusted OR = 1.72, 95%CI: 1.13-2.62, p = 0.011) and age of menarche (AG/GG vs. AA: adjusted OR = 1.54, 95%CI: 1.03-2.31, p = 0.037), whereas ATM rs611646 and Ku70 rs132793 were associated with reduced breast cancer risk influenced by menarche (GA/AA vs. GG: adjusted OR = 0.50, 95%CI: 0.30-0.95, p = 0.033). In a summary, PALB2 rs16940342, ATR rs13091637, ATM rs611646, and Ku70 rs132793 were associated with breast cancer risk.

Genetic Risk, Health-Associated Lifestyle, and Risk of Early-onset Total Cancer and Breast Cancer

Importance

Early-onset cancer (diagnosed under 50 years of age) is associated with aggressive disease characteristics and its rising incidence is a global concern. The association between healthy lifestyle and early-onset cancer and whether it varies by common genetic variants is unknown.

Objective

To examine the associations between genetic risk, lifestyle, and risk of early-onset cancers.

Design Setting and Participants

We analyzed a prospective cohort of 66,308 white British participants who were under age 50 and free of cancer at baseline in the UK Biobank.

Exposures

Sex-specific composite total cancer polygenic risk scores (PRSs), a breast cancer-specific PRS, and sex-specific health-associated lifestyle scores (HLSs, which summarize smoking status, body mass index [males only], physical activity, alcohol consumption, and diet).

Main Outcomes and Measures

Hazard ratios (HRs) and 95% confidence intervals (CIs) for early-onset total and breast cancer.

Results

A total of 1,247 incident invasive early-onset cancer cases (female: 820, male: 427, breast: 386) were documented. In multivariable-adjusted analyses with 2-year latency, higher genetic risk (highest vs. lowest tertile of PRS) was associated with significantly increased risks of early-onset total cancer in females (HR, 95% CI: 1.85, 1.50-2.29) and males (1.94, 1.45-2.59) as well as early-onset breast cancer in females (3.06, 2.20-4.25). An unfavorable lifestyle (highest vs. lowest category of HLS) was associated with higher risk of total cancer and breast cancer in females across genetic risk categories; the association with total cancer was stronger in the highest genetic risk category than the lowest: HRs in females and men were 1.85 (1.02, 3.36), 3.27 (0.78, 13.72) in the highest genetic risk category and 1.15 (0.44, 2.98), 1.16 (0.39, 3.40) in the lowest.

Conclusions and Relevance

Both genetic and lifestyle factors were independently associated with early-onset total and breast cancer risk. Compared to those with low genetic risk, individuals with a high genetic risk may benefit more from adopting a healthy lifestyle in preventing early-onset cancer.

Combined SNPs sequencing and allele specific proteomics capture reveal functional causality underpinning the 2p25 prostate cancer susceptibility locus

Genome wide association studies (GWASs) have identified numerous risk loci associated with prostate cancer, yet unraveling their functional significance remains elusive. Leveraging our high-throughput SNPs-seq method, we pinpointed rs4519489 within the multi-ancestry GWAS-discovered 2p25 locus as a potential functional SNP due to its significant allelic differences in protein binding. Here, we conduct a comprehensive analysis of rs4519489 and its associated gene, NOL10, employing diverse cohort data and experimental models. Clinical findings reveal a synergistic effect between rs4519489 genotype and NOL10 expression on prostate cancer prognosis and severity. Through unbiased proteomics screening, we reveal that the risk allele A of rs4519489 exhibits enhanced binding to USF1, a novel oncogenic transcription factor (TF) implicated in prostate cancer progression and prognosis, resulting in elevated NOL10 expression. Furthermore, we elucidate that NOL10 regulates cell cycle pathways, fostering prostate cancer progression. The concurrent expression of NOL10 and USF1 correlates with aggressive prostate cancer characteristics and poorer prognosis. Collectively, our study offers a robust strategy for functional SNP screening and TF identification through high-throughput SNPs-seq and unbiased proteomics, highlighting the rs4519489-USF1-NOL10 regulatory axis as a promising biomarker or therapeutic target for clinical diagnosis and treatment of prostate cancer.

The molecular map of CLL and Richter's syndrome

Clonal expansion of B-cells, from the early stages of monoclonal B-cell lymphocytosis through to chronic lymphocytic leukemia (CLL), and then in some cases to Richter's syndrome (RS) provides a comprehensive model of cancer evolution, notable for the marked morphological transformation and distinct clinical phenotypes. High-throughput sequencing of large cohorts of patients and single-cell studies have generated a molecular map of CLL and more recently, of RS, yielding fundamental insights into these diseases and of clonal evolution. A selection of CLL driver genes have been functionally interrogated to yield novel insights into the biology of CLL. Such findings have the potential to impact patient care through risk stratification, treatment selection and drug discovery. However, this molecular map remains incomplete, with extant questions concerning the origin of the B-cell clone, the role of the TME, inter- and intra-compartmental heterogeneity and of therapeutic resistance mechanisms. Through the application of multi-modal single-cell technologies across tissues, disease states and clinical contexts, these questions can now be addressed with the answers holding great promise of generating translatable knowledge to improve patient care.

Functional TET2 gene polymorphisms increase the risk of neuroblastoma in Chinese children

The 5-methylcytosine (m5C) is the key chemical modification in RNAs. As one of the demethylases in m5C, TET2 has been shown as a tumor suppressor. However, the impact of TET2 gene polymorphisms on neuroblastoma has not been elucidated. 402 neuroblastoma patients and 473 controls were genotyped for TET2 gene polymorphisms using the TaqMan method. The impact of TET2 gene polymorphisms on neuroblastoma susceptibility was determined using multivariate logistic regression analysis. We also adopted genotype-tissue expression database to explore the impact of TET2 gene polymorphisms on the expression of host and nearby genes. We used the R2 platform and Sangerbox tool to analyze the relationship between gene expression and neuroblastoma risk and prognosis through non-parametric testing and Kaplan-Meier analysis, respectively. We found the TET2 gene polymorphisms (rs10007915 G > C and rs7670522 A > C) and the combined 2-5 risk genotypes can significantly increase neuroblastoma risk. Stratification analysis showed that these significant associations were more prominent in certain subgroups. TET2 rs10007915 G > C and rs7670522 A > C are significantly associated with reduced expression of TET2 mRNA. Moreover, lower expression of TET2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastoma. The rs10007915 G > C and rs7670522 A > C are significantly related to the increased expression of inorganic pyrophosphatase 2 mRNA, and higher expression of PPA2 gene is associated with high risk, MYCN amplification, and poor prognosis of neuroblastomas. In summary, TET2 rs10007915 G > C and rs7670522 A > C significantly confer neuroblastoma susceptibility, and further research is needed to investigate the underlying mechanisms.

Prioritization of genes for translation: a computational approach

INTRODUCTION

Gene identification for genetic diseases is critical for the development of new diagnostic approaches and personalized treatment options. Prioritization of gene translation is an important consideration in the molecular biology field, allowing researchers to focus on the most promising candidates for further investigation.

AREAS COVERED

In this paper, we discussed different approaches to prioritize genes for translation, including the use of computational tools and machine learning algorithms, as well as experimental techniques such as knockdown and overexpression studies. We also explored the potential biases and limitations of these approaches and proposed strategies to improve the accuracy and reliability of gene prioritization methods. Although numerous computational methods have been developed for this purpose, there is a need for computational methods that incorporate tissue-specific information to enable more accurate prioritization of candidate genes. Such methods should provide tissue-specific predictions, insights into underlying disease mechanisms, and more accurate prioritization of genes.

EXPERT OPINION

Using advanced computational tools and machine learning algorithms to prioritize genes, we can identify potential targets for therapeutic intervention of complex diseases. This represents an up-and-coming method for drug development and personalized medicine.

Germline determinants of aberrant signaling pathways in cancer

Cancer is a complex disease influenced by a heterogeneous landscape of both germline genetic variants and somatic aberrations. While there is growing evidence suggesting an interplay between germline and somatic variants, and a substantial number of somatic aberrations in specific pathways are now recognized as hallmarks in many well-known forms of cancer, the interaction landscape between germline variants and the aberration of those pathways in cancer remains largely unexplored. Utilizing over 8500 human samples across 33 cancer types characterized by TCGA and considering binary traits defined using a large collection of somatic aberration profiles across ten well-known oncogenic signaling pathways, we conducted a series of GWAS and identified genome-wide and suggestive associations involving 276 SNPs. Among these, 94 SNPs revealed cis-eQTL links with cancer-related genes or with genes functionally correlated with the corresponding traits' oncogenic pathways. GWAS summary statistics for all tested traits were then used to construct a set of polygenic scores employing a customized computational strategy. Polygenic scores for 24 traits demonstrated significant performance and were validated using data from PCAWG and CCLE datasets. These scores showed prognostic value for clinical variables and exhibited significant effectiveness in classifying patients into specific cancer subtypes or stratifying patients with cancer-specific aggressive phenotypes. Overall, we demonstrate that germline genetics can describe patients' genetic liability to develop specific cancer molecular and clinical profiles.

The evolutionary impact of childhood cancer on the human gene pool

Germline pathogenic variants associated with increased childhood mortality must be subject to natural selection. Here, we analyze publicly available germline genetic metadata from 4,574 children with cancer [11 studies; 1,083 whole exome sequences (WES), 1,950 whole genome sequences (WGS), and 1,541 gene panel] and 141,456 adults [125,748 WES and 15,708 WGS]. We find that pediatric cancer predisposition syndrome (pCPS) genes [n = 85] are highly constrained, harboring only a quarter of the loss-of-function variants that would be expected. This strong indication of selective pressure on pCPS genes is found across multiple lines of germline genomics data from both pediatric and adult cohorts. For six genes [ELP1, GPR161, VHL and SDHA/B/C], a clear lack of mutational constraint calls the pediatric penetrance and/or severity of associated cancers into question. Conversely, out of 23 known pCPS genes associated with biallelic risk, two [9%, DIS3L2 and MSH2] show significant constraint, indicating that they may monoallelically increase childhood cancer risk. In summary, we show that population genetic data provide empirical evidence that heritable childhood cancer leads to natural selection powerful enough to have significantly impacted the present-day gene pool.

A distinct class of pan-cancer susceptibility genes revealed by an alternative polyadenylation transcriptome-wide association study

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3' untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3' untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

Anti-cancer bioactivity of sweet basil leaf derived extracellular vesicles on pancreatic cancer cells

Most living organisms secrete tiny lipid bilayer particles encapsulating various biomolecular entities, including nucleic acids and proteins. These secreted extracellular vesicles (EVs) are shown to aid in communication between cells and their environment. EVs are mainly involved in the signalling and manipulation of physiological processes. Plant EVs display similar functional activity as seen in mammalian EVs. Medicinal plants have many bioactive constituents with potential applications in cancer treatment. Particularly, Basil (Ocimum basilicum), has wide therapeutic properties including anti-inflammatory, anti-cancer, and anti-infection, among others. In this study, we focused on using EVs purified from Apoplast Washing Fluid (AWF) of Basil plant leaves as a biological therapeutic agent against cancer. Characterization of Basil EVs revealed a size range of 100-250 nm, which were later assessed for their cell uptake and apoptosis inducing abilities in pancreatic cancer cells. Basil plant EVs (BasEVs) showed a significant cytotoxic effect on pancreatic cancer cell line MIA PaCa-2 at a concentration of 80 and 160 μg/mL in cell viability, as well as clonogenic assays. Similarly, RT-PCR and western blot analysis has shown up regulation in apoptotic gene and protein expression of Bax, respectively, in BasEV treatment groups compared to untreated controls of MIA PaCa-2. Overall, our results suggest that EVs from basil plants have potent anti-cancer effects in pancreatic cancer cells and can serve as a drug delivery system, demanding an investigation into the therapeutic potential of other medicinal plant EVs.

A Multiethnic Germline-Somatic Association Database Deciphers Multilayered and Interconnected Genetic Mutations in Cancer

Inherited germline and acquired somatic alterations can both promote human tumor development. Elucidating the cooperation between somatic and germline genetic alterations that drive tumorigenesis could help inform precision cancer prevention and treatment strategies. Here, leveraging genomic genotyping and sequencing data from 9,029 patients with cancer with European, East Asian, and African ancestry, we performed a pan-cancer analysis to evaluate the associations between germline SNPs and somatic alterations, including single-nucleotide variant and small insertion/deletion mutations, copy-number variation, tumor mutational burden, and mutational signatures. Genome-wide significant germline-somatic pairs were abundant, and most of the associations were observed in one cancer type and one ancestry group. A user-friendly interactive Multiethnic Germline-Somatic Association (MGSA) database (http://wanglab-hkust.cn:3838/MGSA/) was developed, which can be used to query, browse, and download the results of the association analyses. Moreover, the MGSA database offers additional survival analysis and functional annotation. Together, this work provides a resource for uncovering the clinical and biological roles of associations between germline variants and somatic alterations in human cancer.

SIGNIFICANCE

Comprehensive analysis of connections between germline variants and somatic events in cancer offers a resource for investigating the functional significance of genetic mutations and exploring genetic factors contributing to racial disparities.

Differentiated genomic footprints suggest isolation and long-distance migration of Hmong-Mien populations

BACKGROUND

The underrepresentation of Hmong-Mien (HM) people in Asian genomic studies has hindered our comprehensive understanding of the full landscape of their evolutionary history and complex trait architecture. South China is a multi-ethnic region and indigenously settled by ethnolinguistically diverse HM, Austroasiatic (AA), Tai-Kadai (TK), Austronesian (AN), and Sino-Tibetan (ST) people, which is regarded as East Asia's initial cradle of biodiversity. However, previous fragmented genetic studies have only presented a fraction of the landscape of genetic diversity in this region, especially the lack of haplotype-based genomic resources. The deep characterization of demographic history and natural-selection-relevant genetic architecture of HM people was necessary.

RESULTS

We reported one HM-specific genomic resource and comprehensively explored the fine-scale genetic structure and adaptative features inferred from the genome-wide SNP data of 440 HM individuals from 33 ethnolinguistic populations, including previously unreported She. We identified solid genetic differentiation between HM people and Han Chinese at 7.64‒15.86 years ago (kya) and split events between southern Chinese inland (Miao/Yao) and coastal (She) HM people in the middle Bronze Age period and the latter obtained more gene flow from Ancient Northern East Asians. Multiple admixture models further confirmed that extensive gene flow from surrounding ST, TK, and AN people entangled in forming the gene pool of Chinese coastal HM people. Genetic findings of isolated shared unique ancestral components based on the sharing alleles and haplotypes deconstructed that HM people from the Yungui Plateau carried the breadth of previously unknown genomic diversity. We identified a direct and recent genetic connection between Chinese inland and Southeast Asian HM people as they shared the most extended identity-by-descent fragments, supporting the long-distance migration hypothesis. Uniparental phylogenetic topology and network-based phylogenetic relationship reconstruction found ancient uniparental founding lineages in southwestern HM people. Finally, the population-specific biological adaptation study identified the shared and differentiated natural selection signatures among inland and coastal HM people associated with physical features and immune functions. The allele frequency spectrum of cancer susceptibility alleles and pharmacogenomic genes showed significant differences between HM and northern Chinese people.

CONCLUSIONS

Our extensive genetic evidence combined with the historical documents supported the view that ancient HM people originated from the Yungui regions associated with ancient "Three-Miao tribes" descended from the ancient Daxi-Qujialing-Shijiahe people. Then, some have recently migrated rapidly to Southeast Asia, and some have migrated eastward and mixed respectively with Southeast Asian indigenes, Liangzhu-related coastal ancient populations, and incoming southward ST people. Generally, complex population migration, admixture, and adaptation history contributed to the complicated patterns of population structure of geographically diverse HM people.

A robust and inducible precise genome editing via an all-in-one prime editor in human pluripotent stem cells

Prime editing (PE) allows for precise genome editing in human pluripotent stem cells (hPSCs), such as introducing single nucleotide modifications, small deletions, or insertions at a specific genomic locus, a strategy that shows great promise for creating "Disease in a dish" models. To improve the effectiveness of prime editing in hPSCs, we systematically compared and combined the "inhibition of mismatch repair pathway and p53" on top of the "PEmax" to generate an all-in-one "PE-Plus" prime editor. We show that PE-Plus conducts the most efficient editing among the current PE tools in hPSCs. We further established an inducible prime editing platform in hPSCs by incorporating the all-in-one PE vector into a safe-harbor locus and demonstrated temporal control of precise editing in both hPSCs and differentiated cells. By evaluating disease-associated mutations, we show that this platform allows efficient creation of both monoallelic and biallelic disease-relevant mutations in hPSCs. In addition, this platform enables the efficient introduction of single or multiple edits in one step, demonstrating potential for multiplex editing. Therefore, our method presents an efficient and controllable multiplex prime editing tool in hPSCs and their differentiated progeny.

Leveraging GWAS: Path to Prevention?

Developing novel cancer prevention medication strategies is important for reducing mortality. Identification of common genetic variants associated with cancer risk suggests the potential to leverage these discoveries to define causal targets for cancer interception. Although each risk variant confers small increases in risk, researchers propose that blocking those that produce causal carcinogenic effects might have large impacts on cancer prevention. While a promising concept, we describe potential hurdles that may need to be scaled to reach this goal, including: (i) understanding the complexity of risk; (ii) achieving statistical power in studies with binary outcomes (cancer development: yes or no); (iii) characterization of cancer precursors; (iv) heterogeneity of cancer subtypes and the populations in which these diseases occur; (v) impact of static genetic markers across complex events of the life course; (vi) defining gene-gene and gene-environment interactions and (vii) demonstrating functional effects of markers in human populations. We assess short-term prospects for this research against the backdrop of these challenges and the potential to prevent cancer through other means. See related commentary by Peters and Tomlinson, p. 7.

Current trends and future prospects of drug repositioning in gastrointestinal oncology

Gastrointestinal (GI) cancers comprise a significant number of cancer cases worldwide and contribute to a high percentage of cancer-related deaths. To improve survival rates of GI cancer patients, it is important to find and implement more effective therapeutic strategies with better prognoses and fewer side effects. The development of new drugs can be a lengthy and expensive process, often involving clinical trials that may fail in the early stages. One strategy to address these challenges is drug repurposing (DR). Drug repurposing is a developmental strategy that involves using existing drugs approved for other diseases and leveraging their safety and pharmacological data to explore their potential use in treating different diseases. In this paper, we outline the existing therapeutic strategies and challenges associated with GI cancers and explore DR as a promising alternative approach. We have presented an extensive review of different DR methodologies, research efforts and examples of repurposed drugs within various GI cancer types, such as colorectal, pancreatic and liver cancers. Our aim is to provide a comprehensive overview of employing the DR approach in GI cancers to inform future research endeavors and clinical trials in this field.

Increased hsa-miR-100-5p Expression Improves Hepatocellular Carcinoma Prognosis in the Asian Population with PLK1 Variant rs27770A>G

Hepatocellular carcinoma (HCC) has the highest incidence and mortality in the Asian population, and race is an independent risk factor affecting survival time in liver cancer. Micro RNAs (miRNAs) are remarkably dysregulated in HCC and closely associated with HCC prognosis. Recent studies show that genetic variability between ethnic groups may result in differences in the specificity of HCC miRNA biomarkers. Here, we reveal a high expression level of hsa-miR-100-5p, an HCC prognosis-related miRNA, which improves HCC prognosis in the Asian Population with Polo-like kinase 1 (PLK1) variant rs27770A>G. In this study, we discovered that hsa-miR-100-5p was downregulated in various HCC cell lines. While mimics transient transfection and mouse liver cancer model confirmed the interaction between hsa-miR-100-5p and PLK1, a stratified analysis based on the Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) data suggest both low hsa-miR-100-5p expression level and high PLK1 expression level associated with poor HCC prognosis, especially in the Asian population. According to the 1000 Genomes Project database, the SNP rs27770 located in 3'UTR of PLK1 had a significantly higher G allele frequency in the East Asian population. Bioinformatics analysis suggested that rs27770 A>G affects PLK1 mRNA secondary structure and alters the hsa-miR-100-5p/PLK1 interaction by forming an additional seedless binding site. This racial variation caused PLK1 to be more vulnerable to hsa-miR-100-5p inhibition, resulting in hsa-miR-100-5p being more favorable for HCC prognosis in the Asian population.

Genetic modifiers of repeat expansion disorders

Repeat expansion disorders (REDs) are monogenic diseases caused by a sequence of repetitive DNA expanding above a pathogenic threshold. A common feature of the REDs is a strong genotype-phenotype correlation in which a major determinant of age at onset (AAO) and disease progression is the length of the inherited repeat tract. Over a disease-gene carrier's life, the length of the repeat can expand in somatic cells, through the process of somatic expansion which is hypothesised to drive disease progression. Despite being monogenic, individual REDs are phenotypically variable, and exploring what genetic modifying factors drive this phenotypic variability has illuminated key pathogenic mechanisms that are common to this group of diseases. Disease phenotypes are affected by the cognate gene in which the expansion is found, the location of the repeat sequence in coding or non-coding regions and by the presence of repeat sequence interruptions. Human genetic data, mouse models and in vitro models have implicated the disease-modifying effect of DNA repair pathways via the mechanisms of somatic mutation of the repeat tract. As such, developing an understanding of these pathways in the context of expanded repeats could lead to future disease-modifying therapies for REDs.

Causal association of gut microbiota and esophageal cancer: a Mendelian randomization study

Introduction

Despite the growing body of evidence, the link between the gut microbiota and different types of tumors, such as colorectal, gastric, and liver cancer, is becoming more apparent. The gut microbiota can be used as a reference for evaluating various diseases, including cancer, and can also act as risk factors or preventive factors. However, the specific connection between the gut microbiota and the advancement of esophageal cancer has yet to be investigated. Therefore, the aim of this research is to clarify the possible causal influence of intestinal microorganisms on the vulnerability to esophageal cancer through the utilization of Mendelian randomization (MR) studies.

Methods

In this study, we employed a two-sample Mendelian randomization approach to evaluate the unbiased causal association between 150 different gut microbiota types and the occurrence of esophageal cancer. Following the selection from the IEU GWAS database and SNP filtration, we utilized various MR statistical techniques on the suitable instrumental variables. These included IVW methods, employing inverse variance weighting. Additionally, we performed a range of sensitivity analyses to confirm the heterogeneity and pleiotropy of the instrumental variables, thus ensuring the reliability of the outcomes.

Results

The increased likelihood of developing esophageal cancer is linked to the genetically predicted high levels of Gordonibacter, Oxalobacter, Coprobacter, Veillonella, Ruminiclostridium 5, Ruminococcus 1, and Senegalimasilia genera. Conversely, a decreased risk of esophageal cancer is associated with the high abundance of Turicibacter, Eubacterium oxidoreducens group, Romboutsia, and Prevotella 9 genera. No heterogeneity and pleiotropy were detected in the sensitivity analysis.

Discussion

We found that 11 types of gut microbial communities are associated with esophageal cancer, thereby confirming that the gut microbiota plays a significant role in the path.

Causal link between gut microbiota and four types of pancreatitis: a genetic association and bidirectional Mendelian randomization study

BACKGROUND

A number of recent observational studies have indicated a correlation between the constitution of gut microbiota and the incidence of pancreatitis. Notwithstanding, observational studies are unreliable for inferring causality because of their susceptibility to confounding, bias, and reverse causality, the causal relationship between specific gut microbiota and pancreatitis is still unclear. Therefore, our study aimed to investigate the causal relationship between gut microbiota and four types of pancreatitis.

METHODS

An investigative undertaking encompassing a genome-wide association study (GWAS) comprising 18,340 participants was undertaken with the aim of discerning genetic instrumental variables that exhibit associations with gut microbiota, The aggregated statistical data pertaining to acute pancreatitis (AP), alcohol-induced AP (AAP), chronic pancreatitis (CP), and alcohol-induced CP (ACP) were acquired from the FinnGen Consortium. The two-sample bidirectional Mendelian randomization (MR) approach was utilized. Utilizing the Inverse-Variance Weighted (IVW) technique as the cornerstone of our primary analysis. The Bonferroni analysis was used to correct for multiple testing, In addition, a number of sensitivity analysis methodologies, comprising the MR-Egger intercept test, the Cochran's Q test, MR polymorphism residual and outlier (MR-PRESSO) test, and the leave-one-out test, were performed to evaluate the robustness of our findings.

RESULTS

A total of 28 intestinal microflora were ascertained to exhibit significant associations with diverse outcomes of pancreatitis. Among them, Class Melainabacteria (OR = 1.801, 95% CI: 1.288-2.519, p = 0.008) has a strong causality with ACP after the Bonferroni-corrected test, in order to assess potential reverse causation effects, we used four types of pancreatitis as the exposure variable and scrutinized its impact on gut microbiota as the outcome variable, this analysis revealed associations between pancreatitis and 30 distinct types of gut microflora. The implementation of Cochran's Q test revealed a lack of substantial heterogeneity among the various single nucleotide polymorphisms (SNP).

CONCLUSION

Our first systematic Mendelian randomization analysis provides evidence that multiple gut microbiota taxa may be causally associated with four types of pancreatitis disease. This discovery may contribute significant biomarkers conducive to the preliminary, non-invasive identification of Pancreatitis. Additionally, it could present viable targets for potential therapeutic interventions in the disease's treatment.

Clustering of colon, lung, and other cancer susceptibility genes with protein tyrosine phosphatases and protein kinases in multiple short genomic regions

Interactions of large gene families are poorly understood. We found that human, mouse, and rat colon and lung cancer susceptibility genes, presently considered as separate gene families, were frequently pairwise linked. The orthologous mouse map positions of 142 of 159 early discovered colon and lung cancer susceptibility genes formed 41 genomic clusters conserved >70 million years. These linked gene pairs concordantly affected both tumors and their majority was linked with two other gene families - protein tyrosine phosphatases and cancer driver protein kinases. 25% of both protein tyrosine phosphatases and protein kinases mapped <1 cM from a colon or lung cancer susceptibility gene, and 50% in <3 cM. Similar linkage was detected with most other human susceptibility genes that controlled 29 different cancer types. This concentration of tumor susceptibility genes with protein tyrosine phosphatases and driver protein kinases in multiple relatively short genomic regions suggests their possible functional diversity.

Assessing a Polygenic Risk Score for Lung Cancer Susceptibility in Non-Hispanic White and Black Populations

BACKGROUND

Polygenic risk scores (PRS) have become an increasingly popular approach to evaluate cancer susceptibility, but have not adequately represented Black populations in model development.

METHODS

We used a previously published lung cancer PRS on the basis of 80 SNPs associated with lung cancer risk in the OncoArray cohort and validated in UK Biobank. The PRS was evaluated for association with lung cancer risk adjusting for age, sex, total pack-years, family history of lung cancer, history of chronic obstructive pulmonary disease, and the top five principal components for genetic ancestry.

RESULTS

Among the 80 PRS SNPs included in the score, 14 were significantly associated with lung cancer risk (P < 0.05) in INHALE White participants, while there were no significant SNPs among INHALE Black participants. After adjusting for covariates, the PRS was significantly associated with risk in Whites (continuous score P = 0.007), but not in Blacks (continuous score P = 0.88). The PRS remained a statistically significant predictor of lung cancer risk in Whites ineligible for lung cancer screening under current U.S. Preventive Services Task Force guidelines (P = 0.02).

CONCLUSIONS

Using a previously validated PRS, we did find some predictive ability for lung cancer in INHALE White participants beyond traditional risk factors. However, this effect was not observed in Black participants, indicating the need to develop and validate ancestry-specific lung cancer risk models.

IMPACT

While a previously published lung cancer PRS was able to stratify White participants into different levels of risk, the model was not predictive in Blacks. Our findings highlight the need to develop and validate ancestry-specific lung cancer risk models.

Phenotypic variance partitioning by transcriptomic gene expression levels and environmental variables for anthropometric traits using GTEx data

Phenotypic variation in human is the results of genetic variation and environmental influences. Understanding the contribution of genetic and environmental components to phenotypic variation is of great interest. The variance explained by genome-wide single nucleotide polymorphisms (SNPs) typically represents a small proportion of the phenotypic variance for complex traits, which may be because the genome is only a part of the whole biological process to shape the phenotypes. In this study, we propose to partition the phenotypic variance of three anthropometric traits, using gene expression levels and environmental variables from GTEx data. We use the gene expression of four tissues that are deemed relevant for the anthropometric traits (two adipose tissues, skeletal muscle tissue and blood tissue). Additionally, we estimate the transcriptome-environment correlation that partly underlies the phenotypes of the anthropometric traits. We found that genetic factors play a significant role in determining body mass index (BMI), with the proportion of phenotypic variance explained by gene expression levels of visceral adipose tissue being 0.68 (SE = 0.06). However, we also observed that environmental factors such as age, sex, ancestry, smoking status, and drinking alcohol status have a small but significant impact (0.005, SE = 0.001). Interestingly, we found a significant negative correlation between the transcriptomic and environmental effects on BMI (transcriptome-environment correlation = -0.54, SE = 0.14), suggesting an antagonistic relationship. This implies that individuals with lower genetic profiles may be more susceptible to the effects of environmental factors on BMI, while those with higher genetic profiles may be less susceptible. We also show that the estimated transcriptomic variance varies across tissues, e.g., the gene expression levels of whole blood tissue and environmental variables explain a lower proportion of BMI phenotypic variance (0.16, SE = 0.05 and 0.04, SE = 0.004 respectively). We observed a significant positive correlation between transcriptomic and environmental effects (1.21, SE = 0.23) for this tissue. In conclusion, phenotypic variance partitioning can be done using gene expression and environmental data even with a small sample size (n = 838 from GTEx data), which can provide insights into how the transcriptomic and environmental effects contribute to the phenotypes of the anthropometric traits.

Association of Vitamin D Genetic Risk Score with Noncommunicable Diseases: A Systematic Review

Background and Aims: The genetic risk score (GRS) is an important tool for estimating the total genetic contribution or susceptibility to a certain outcome of interest in an individual, taking into account their genetic risk alleles. This study aims to systematically review the association between the GRS of low vitamin D with different noncommunicable diseases/markers. Methods: The article was first registered in PROSPERO CRD42023406929. PubMed and Embase were searched from the time of inception until March 2023 to capture all the literature related to the vitamin D genetic risk score (vD-GRS) in association with noncommunicable diseases. This was performed using comprehensive search terms including "Genetic Risk Score" OR "Genetics risk assessment" OR "Genome-wide risk score" AND "Vitamin D" OR 25(HO)D OR "25-hydroxyvitamin D". Results: Eleven eligible studies were included in this study. Three studies reported a significant association between vD-GRS and metabolic parameters, including body fat percentage, body mass index, glycated hemoglobin, and fasting blood glucose. Moreover, colorectal cancer overall mortality and the risk of developing arterial fibrillation were also found to be associated with genetically deprived vitamin D levels. Conclusions: This systematic review highlights the genetic contribution of low-vitamin-D-risk single nucleotides polymorphisms (SNPs) as an accumulative factor associated with different non-communicable diseases/markers, including cancer mortality and the risk of developing obesity, type 2 diabetes, and cardiovascular diseases such as arterial fibrillation.

Gene-based burden scores identify rare variant associations for 28 blood biomarkers

BACKGROUND

A relevant part of the genetic architecture of complex traits is still unknown; despite the discovery of many disease-associated common variants. Polygenic risk score (PRS) models are based on the evaluation of the additive effects attributable to common variants and have been successfully implemented to assess the genetic susceptibility for many phenotypes. In contrast, burden tests are often used to identify an enrichment of rare deleterious variants in specific genes. Both kinds of genetic contributions are typically analyzed independently. Many studies suggest that complex phenotypes are influenced by both low effect common variants and high effect rare deleterious variants. The aim of this paper is to integrate the effect of both common and rare functional variants for a more comprehensive genetic risk modeling.

METHODS

We developed a framework combining gene-based scores based on the enrichment of rare functionally relevant variants with genome-wide PRS based on common variants for association analysis and prediction models. We applied our framework on UK Biobank dataset with genotyping and exome data and considered 28 blood biomarkers levels as target phenotypes. For each biomarker, an association analysis was performed on full cohort using gene-based scores (GBS). The cohort was then split into 3 subsets for PRS construction and feature selection, predictive model training, and independent evaluation, respectively. Prediction models were generated including either PRS, GBS or both (combined).

RESULTS

Association analyses of the cohort were able to detect significant genes that were previously known to be associated with different biomarkers. Interestingly, the analyses also revealed heterogeneous effect sizes and directionality highlighting the complexity of the blood biomarkers regulation. However, the combined models for many biomarkers show little or no improvement in prediction accuracy compared to the PRS models.

CONCLUSION

This study shows that rare variants play an important role in the genetic architecture of complex multifactorial traits such as blood biomarkers. However, while rare deleterious variants play a strong role at an individual level, our results indicate that classical common variant based PRS might be more informative to predict the genetic susceptibility at the population level.

GATA2 co-opts TGFβ1/SMAD4 oncogenic signaling and inherited variants at 6q22 to modulate prostate cancer progression

BACKGROUND

Aberrant somatic genomic alteration including copy number amplification is a hallmark of cancer genomes. We previously profiled genomic landscapes of prostate cancer (PCa), yet the underlying causal genes with prognostic potential has not been defined. It remains unclear how a somatic genomic event cooperates with inherited germline variants contribute to cancer predisposition and progression.

METHODS

We applied integrated genomic and clinical data, experimental models and bioinformatic analysis to identify GATA2 as a highly prevalent metastasis-associated genomic amplification in PCa. Biological roles of GATA2 in PCa metastasis was determined in vitro and in vivo. Global chromatin co-occupancy and co-regulation of GATA2 and SMAD4 was investigated by coimmunoprecipitation, ChIP-seq and RNA-seq assays. Tumor cellular assays, qRT-PCR, western blot, ChIP, luciferase assays and CRISPR-Cas9 editing methods were performed to mechanistically understand the cooperation of GATA2 with SMAD4 in promoting TGFβ1 and AR signaling and mediating inherited PCa risk and progression.

RESULTS

In this study, by integrated genomics and experimental analysis, we identified GATA2 as a prevalent metastasis-associated genomic amplification to transcriptionally augment its own expression in PCa. Functional experiments demonstrated that GATA2 physically interacted and cooperated with SMAD4 for genome-wide chromatin co-occupancy and co-regulation of PCa genes and metastasis pathways like TGFβ signaling. Mechanistically, GATA2 was cooperative with SMAD4 to enhance TGFβ and AR signaling pathways, and activated the expression of TGFβ1 via directly binding to a distal enhancer of TGFβ1. Strinkingly, GATA2 and SMAD4 globally mediated inherited PCa risk and formed a transcriptional complex with HOXB13 at the PCa risk-associated rs339331/6q22 enhancer, leading to increased expression of the PCa susceptibility gene RFX6.

CONCLUSIONS

Our study prioritizes causal genomic amplification genes with prognostic values in PCa and reveals the pivotal roles of GATA2 in transcriptionally activating the expression of its own and TGFβ1, thereby co-opting to TGFβ1/SMAD4 signaling and RFX6 at 6q22 to modulate PCa predisposition and progression.

Pan-cancer and cross-population genome-wide association studies dissect shared genetic backgrounds underlying carcinogenesis

Integrating genomic data of multiple cancers allows de novo cancer grouping and elucidating the shared genetic basis across cancers. Here, we conduct the pan-cancer and cross-population genome-wide association study (GWAS) meta-analysis and replication studies on 13 cancers including 250,015 East Asians (Biobank Japan) and 377,441 Europeans (UK Biobank). We identify ten cancer risk variants including five pleiotropic associations (e.g., rs2076295 at DSP on 6p24 associated with lung cancer and rs2525548 at TRIM4 on 7q22 nominally associated with six cancers). Quantifying shared heritability among the cancers detects positive genetic correlations between breast and prostate cancer across populations. Common genetic components increase the statistical power, and the large-scale meta-analysis of 277,896 breast/prostate cancer cases and 901,858 controls identifies 91 newly genome-wide significant loci. Enrichment analysis of pathways and cell types reveals shared genetic backgrounds across said cancers. Focusing on genetically correlated cancers can contribute to enhancing our insights into carcinogenesis.

Germline variants of the genes involved in NF-kB activation are associated with the risk of COPD and lung cancer development

Chronic obstructive pulmonary disease (COPD) and lung cancer (LC) are closely related diseases associated with smoking history and dysregulated immune response. However, not all smokers develop the disease, indicating that genetic susceptibility could be important. Therefore, the aim of this study was to search for the potential overlapping genetic biomarkers, with a focus on single nucleotide polymorphisms (SNPs) located in the regulatory regions of immune-related genes. Additionally, the aim was to see if an identified SNP has potentially an effect on proinflamma-tory cytokine concentration in the serum of COPD patients. We extracted summary data of variants in 1511 immune-related genes from COPD and LC genome-wide association studies (GWAS) from the UK Biobank. The LC data had 203 cases, patients diagnosed with LC, and 360 938 controls, while COPD data had 1 897 cases and 359 297 controls. Assuming 1 association/gene, SNPs with a p-value < 3.3 × 10^-5 were considered statistically significantly associated with the disease. We identified seven SNPs located in different genes (BAG6, BTNL2, TNF, HCP5, MICB, NCR3, ABCF1, TCF7L1) to be associated with the COPD risk and two with the LC risk (HLA-C, HLA-B), with statistical significance. We also identified two SNPs located in the IL2RA gene associated with LC (rs2386841; p = 1.86 × 10^-4) and COPD (rs11256442; p = 9.79 × 10^-3) but with lower significance. Functional studies conducted on COPD patients showed that RNA expression of IL2RA, IFNγ and related proinflammatory cytokines in blood serum did not correlate with a specific genotype. Although results presented in this study do not fully support our hypothesis, it is worth to mention that the identified genes/SNPs that were associated with either COPD or LC risk, all were involved in the activation of the NF-κB transcription factor which is closely related to the regulation of the inflammatory response, a condition associated with both pathologies.

Development of a Somatic Variant Registry in a National Cancer Center: towards Molecular Real World Data preparedness

The Biomedical Research field is currently advancing to develop Clinical Trials and translational projects based on Real World Evidence. To make this transition feasible, clinical centers need to work toward Data Accessibility and Interoperability. This task is particularly challenging when applied to Genomics, that entered in routinary screening in the last years via mostly amplicon-based Next-Generation Sequencing panels. Said experiments produce up to hundreds of features per patient, and their summarized results are often stored in static clinical reports, making critical information inaccessible to automated access and Federated Search consortia. In this study, we present a reanalysis of 4620 solid tumor sequencing samples in five different histology settings. Furthermore, we describe all the Bioinformatics and Data Engineering processes that were put in place in order to create a Somatic Variant Registry able to deal with the large biotechnological variability of routinary Genomics Profiling.

Inherited mutations affecting the SRCAP complex are central in moderate-penetrance predisposition to uterine leiomyomas

Uterine leiomyomas (ULs) are benign smooth muscle tumors that are common in premenopausal women. Somatic alterations in MED12, HMGA2, FH, genes encoding subunits of the SRCAP complex, and genes involved in Cullin 3-RING E3 ligase neddylation are mutually exclusive UL drivers. Established predisposition genes explain only partially the estimated heritability of leiomyomas. Here, we examined loss-of-function variants across 18,899 genes in a cohort of 233,614 White European women, revealing variants in four genes encoding SRCAP complex subunits (YEATS4, ZNHIT1, DMAP1, and ACTL6A) with a significant association to ULs, and YEATS4 and ZNHIT1 strikingly rank first and second, respectively. Positive mutation status was also associated with younger age at diagnosis and hysterectomy. Moderate-penetrance UL risk was largely attributed to rare non-synonymous mutations affecting the SRCAP complex. To examine this disease phenotype more closely, we set out to identify inherited mutations affecting the SRCAP complex in our in-house sample collection of Finnish individuals with ULs (n = 860). We detected one individual with an ACTL6A splice-site mutation, two individuals with a YEATS4 missense mutation, and four individuals with DMAP1 mutations: one splice-site, one nonsense, and two missense variants. These individuals had large and/or multiple ULs, were often diagnosed at an early age, and many had family history of ULs. When a somatic second hit was found, ACTL6A and DMAP1 were silenced in tumors by somatic mutation and YEATS4 by promoter hypermethylation. Decreased H2A.Z staining was observed in the tumors, providing further evidence for the pathogenic nature of the germline mutations. Our results establish inactivation of genes encoding SRCAP complex subunits as a central contributor to moderate-penetrance UL predisposition.

Exome sequencing in a Swedish family with PMS2 mutation with varying penetrance of colorectal cancer: investigating the presence of genetic risk modifiers in colorectal cancer risk

OBJECTIVE

Lynch syndrome is caused by germline mutations in the mismatch repair (MMR) genes, such as the PMS2 gene, and is characterised by a familial accumulation of colorectal cancer. The penetrance of cancer in PMS2 carriers is still not fully elucidated as a colorectal cancer risk has been shown to vary between PMS2 carriers, suggesting the presence of risk modifiers.

METHODS

Whole exome sequencing was performed in a Swedish family carrying a PMS2 missense mutation [c.2113G>A, p.(Glu705Lys)]. Thirteen genetic sequence variants were further selected and analysed in a case-control study (724 cases and 711 controls).

RESULTS

The most interesting variant was an 18 bp deletion in gene BAG1. BAG1 has been linked to colorectal tumour progression with poor prognosis and is thought to promote colorectal tumour cell survival through increased NF-κB activity.

CONCLUSIONS

We conclude the genetic architecture behind the incomplete penetrance of PMS2 is complicated and must be assessed in a genome wide manner using large families and multifactorial analysis.

Prognostic evaluation of polygenic risk score underlying pan-cancer analysis: evidence from two large-scale cohorts

BACKGROUND

Polygenic risk score (PRS) has been demonstrated to be effective in identifying individuals at high risk of developing cancer, but its prognostic value remains unclear.

METHODS

We constructed site-specific PRSs by aggregating the risk effect of independent variants derived from previous genome-wide association studies (GWASs) across 17 cancer types. The Cox proportional hazards model was used to evaluate the association of each PRS with cancer survival, leveraging data from two prospective European cohorts, namely the UK Biobank involving 19,628 incident cases and The Cancer Genome Atlas involving 7079 prevalent cases. The combined PRS (CPRS), determined by merging site-specific PRSs, was further used to assess the prognostic effect of PRS on overall cancer in a sex-specific manner.

FINDINGS

We discovered that the cancer risk-related PRS was associated with neither overall survival (OS) nor cancer-specific survival (CSS) of each site-specific cancer with an underlying false discovery rate (FDR) P > 0.05, as evidenced by consistent findings from the two cohorts. Furthermore, the fixed-effect meta-analysis of the two cohorts provided no evidence to support for an association between CPRS and overall cancer survival in both males [OS: hazard ratio (HR)_meta = 1.00, P_meta = 0.760; CSS: HR_meta = 1.01, P_meta = 0.447] and females (OS: HR_meta = 0.97, P_meta = 0.067; CSS: HR_meta = 0.96, P_meta = 0.054). Similar results were observed across multiple sensitivity analyses.

INTERPRETATION

Our findings indicate that the risk-specific PRS might not be a clinically useful tool in cancer prognosis prediction and further studies focusing on the development of polygenic prognostic score are warranted.

FUNDING

This project was funded by the National Natural Science Foundation of China (82173601 and 82073631), and Priority Academic Program Development of Jiangsu Higher Education Institutions (Public Health and Preventive Medicine).

Causal relationship between gut microbiota and cancers: a two-sample Mendelian randomisation study

BACKGROUND

Evidence from observational studies and clinical trials suggests that the gut microbiota is associated with cancer. However, the causal association between gut microbiota and cancer remains to be determined.

METHODS

We first identified two sets of gut microbiota based on phylum, class, order, family, and genus level information, and cancer data were obtained from the IEU Open GWAS project. We then performed two-sample Mendelian randomisation (MR) to determine whether the gut microbiota is causally associated with eight cancer types. Furthermore, we performed a bi-directional MR analysis to examine the direction of the causal relations.

RESULTS

We identified 11 causal relationships between genetic liability in the gut microbiome and cancer, including those involving the genus Bifidobacterium. We found 17 strong associations between genetic liability in the gut microbiome and cancer. Moreover, we found 24 associations between genetic liability in the gut microbiome and cancer using multiple datasets.

CONCLUSIONS

Our MR analysis revealed that the gut microbiota was causally associated with cancers and may be useful in providing new insights for further mechanistic and clinical studies of microbiota-mediated cancer.

Learning from the nexus of autoimmunity and cancer

The immune system plays critical roles in both autoimmunity and cancer, diseases at opposite ends of the immune spectrum. Autoimmunity arises from loss of T cell tolerance against self, while in cancer, poor immunity against transformed self fails to control tumor growth. Blockade of pathways that preserve self-tolerance is being leveraged to unleash immunity against many tumors; however, widespread success is hindered by the autoimmune-like toxicities that arise in treated patients. Knowledge gained from the treatment of autoimmunity can be leveraged to treat these toxicities in patients. Further, the understanding of how T cell dysfunction arises in cancer can be leveraged to induce a similar state in autoreactive T cells. Here, we review what is known about the T cell response in autoimmunity and cancer and highlight ways in which we can learn from the nexus of these two diseases to improve the application, efficacy, and management of immunotherapies.

A genome-wide cross-cancer meta-analysis highlights the shared genetic links of five solid cancers

Breast, ovarian, prostate, lung, and head/neck cancers are five solid cancers with complex interrelationships. However, the shared genetic factors of the five cancers were often revealed either by the combination of individual genome-wide association study (GWAS) approach or by the fixed-effect model-based meta-analysis approach with practically impossible assumptions. Here, we presented a random-effect model-based cross-cancer meta-analysis framework for identifying the genetic variants jointly influencing the five solid cancers. A comprehensive genetic correlation analysis (genome-wide, partitioned, and local) approach was performed by using GWAS summary statistics of the five cancers, and we observed three cancer pairs with significant genetic correlation: breast-ovarian cancer (r _g = 0.221, p = 0.0003), breast-lung cancer (r _g = 0.234, p = 7.6 × 10^-6), and lung-head/neck cancer (r _g = 0.652, p = 0.010). Furthermore, a random-effect model-based cross-trait meta-analysis was conducted for each significant cancer pair, and we found 27 shared genetic loci between breast and ovarian cancers, 18 loci between breast and lung cancers, and three loci between lung and head/neck cancers. Functional analysis indicates that the shared genes are enriched in human T-cell leukemia virus 1 infection (HTLV-1) and antigen processing and presentation (APP) pathways. Our study investigates the shared genetic links across five solid cancers and will help to reveal their potential molecular mechanisms.

Deciphering the chromatin spatial organization landscapes during BMMSC differentiation

The differentiation imbalance in bone marrow mesenchymal stem cells (BMMSCs) is critical for the development of bone density diseases as the population ages. BMMSCs are precursor cells for osteoblasts and adipocytes; however, the chromatin organization landscapes during BMMSC differentiation remain elusive. In this study, we systematically delineate the four-dimensional (4D) genome and dynamic epigenetic atlas of BMMSCs by RNA sequencing (RNA-seq), assay for transposase-accessible chromatin sequencing (ATAC-seq), and high-throughput chromosome conformation capture (Hi-C). The structure analyses reveal 17.5% common and 28.5%-30% specific loops among BMMSCs, osteoblasts, and adipocytes. The subsequent correlation of genome-wide association studies (GWAS) and expression quantitative trait locus (eQTL) data with multi-omics analysis reveal 274 genes and 3634 single nucleotide polymorphisms (SNPs) associated with bone degeneration and osteoporosis (OP). We hypothesize that SNP mutations affect transcription factor (TF) binding sites, thereby affecting changes in gene expression. Furthermore, 26 motifs, 260 TFs, and 291 SNPs are identified to affect the eQTL. Among these genes, DAAM2, TIMP2, and TMEM241 were found to be essential for diseases such as bone degeneration and OP and may serve as potential drug targets.

Risk single-nucleotide polymorphism-mediated enhancer-promoter interaction drives keloids through long noncoding RNA down expressed in keloids

BACKGROUND

Keloids represent one extreme of aberrant dermal wound healing and are characterized by fibroblast hyperproliferation and excessive deposition of extracellular matrix. Genetics is a major factor for predisposition to keloids and genome-wide association study has identified a single-nucleotide polymorphism (SNP) rs873549 at 1q41 as a susceptibility locus. The SNP rs873549, and the SNPs in strong linkage disequilibrium (LD) with rs873549, may be involved in keloid development. However, the functional significance of these SNPs in keloid pathogenesis remains elusive.

OBJECTIVES

To investigate the function and mechanism of SNP rs873549 and the SNPs in strong LD with rs873549 in keloids.

METHODS

SNPs in strong LD with rs873549 were analysed using Haploview. The expression levels of the genes near the susceptibility locus were analysed using quantitative real-time polymerase chain reaction. The interaction between rs1348270-containing enhancer and the long noncoding RNA down expressed in keloids (DEIK) (formerly RP11-400N13.1) promoter in fibroblasts was investigated using chromosome conformation capture. The enhancer activity of the rs1348270 locus was evaluated using luciferase reporter assay. Knockdown experiments were used to explore the function of DEIK in keloids. RNA-Seq was performed to investigate the mechanism by which DEIK regulates the expression of collagens POSTN and COMP.

RESULTS

rs1348270, an enhancer-located SNP in strong LD with rs873549, mediated looping with the promoter of DEIK. The risk variant was associated with decreased enhancer-promoter interaction and DEIK down-expression in keloids. Mechanistically, downregulation of DEIK increased the expression of collagens POSTN and COMP through upregulating BMP2. Furthermore, correlation analysis revealed that DEIK expression was inversely correlated with BMP2, POSTN and COMP expression in both keloid and normal fibroblasts.

CONCLUSIONS

Our findings suggest that the risk variant rs1348270 is located in an enhancer and is associated with the downregulation of DEIK in keloids, and that downregulation of DEIK increases the expression of collagens POSTN and COMP through BMP2 in keloid fibroblasts. These findings will help to provide a more thorough understanding of the role played by genetic factors in keloid development and may lead to new strategies for screening and therapy in keloid-susceptible populations.

A genome-wide association study of germline variation and melanoma prognosis

Background

The high mortality of cutaneous melanoma (CM) is partly due to unpredictable patterns of disease progression in patients with early-stage lesions. The reliable prediction of advanced disease risk from early-stage CM, is an urgent clinical need, especially given the recent expansion of immune checkpoint inhibitor therapy to the adjuvant setting. In our study, we comprehensively investigated the role of germline variants as CM prognostic markers.

Methods

We performed a genome-wide association analysis in two independent cohorts of N=551 (discovery), and N=550 (validation) early-stage immunotherapy-naïve melanoma patients. A multivariable Cox proportional hazard regression model was used to identify associations with overall survival in the discovery group, followed by a validation analysis. Transcriptomic profiling and survival analysis were used to elucidate the biological relevance of candidate genes associated with CM progression.

Results

We found two independent associations of germline variants with melanoma prognosis. The alternate alleles of these two SNPs were both associated with an increased risk of death [rs60970102 in MELK: HR=3.14 (2.05-4.81), p=1.48×10^-7; and rs77480547 in SH3BP4: HR=3.02 (2.02-4.52), p=7.58×10^-8, both in the pooled cohort]. The addition of the combined risk alleles (CRA) of the identified variants into the prognostic model improved the predictive power, as opposed to a model of clinical covariates alone.

Conclusions

Our study provides suggestive evidence of novel melanoma germline prognostic markers, implicating two candidate genes: an oncogene MELK and a tumor suppressor SH3BP4, both previously suggested to affect CM progression. Pending further validation, these findings suggest that the genetic factors may improve the prognostic stratification of high-risk early-stage CM patients, and propose putative biological insights for potential therapeutic investigation of these targets to prevent aggressive outcome from early-stage melanoma.

Association of several loci of SMAD7 with colorectal cancer: A meta-analysis based on case-control studies

BACKGROUND

Sma-and mad-related protein 7 (SMAD7) can affect tumor progression by closing transforming growth factor-beta intracellular signaling channels. Despite the extensive research on the correlation between SMAD7 polymorphisms and colorectal cancer (CRC), the conclusions of studies are still contradictory. We conducted a study focusing on the association of SMAD7 polymorphisms rs4939827, rs4464148, and rs12953717 with CRC.

METHODS

We searched through 5 databases for articles and used odd ratios (ORs) and 95% confidence intervals (CIs) to discuss the correlation of SMAD7 polymorphisms with CRC risk. The heterogeneity will be appraised by subgroup analysis and meta-regression. Contour-enhanced funnel plot, Begg test and Egger test were utilized to estimate publication bias, and the sensitivity analysis illustrates the reliability of the outcomes. We performed False-positive report probability and trial sequential analysis methods to verify results. We also used public databases for bioinformatics analysis.

RESULTS

We conclusively included 34 studies totaling 173251 subjects in this study. The minor allele (C) of rs4939827 is a protective factor of CRC (dominant, OR/[95% CI] = 0.89/[0.83-0.97]; recessive, OR/[95% CI] = 0.89/[0.83-0.96]; homozygous, OR/[95% CI] = 0.84/[0.76-0.93]; heterozygous, OR/[95% CI] = 0.91/[0.85-0.97]; additive, OR/[95% CI] = 0.91/[0.87-0.96]). the T allele of rs12953717 (recessive, OR/[95% CI] = 1.22/[1.15-1.28]; homozygous, OR/[95% CI] = 1.25/[1.13-1.38]; additive, OR/[95% CI] = 1.11/[1.05-1.17]) and the C allele of rs4464148 (heterozygous, OR/[95% CI] = 1.13/[1.04-1.24]) can enhance the risk of CRC.

CONCLUSION

Rs4939827 (T > C) can decrease the susceptibility to CRC. However, the rs4464148 (T > C) and rs12953717 (C > T) variants were connected with an enhanced risk of CRC.

SUMMER: a Mendelian randomization interactive server to systematically evaluate the causal effects of risk factors and circulating biomarkers on pan-cancer survival

Genome-wide association studies (GWASs) underlying case-control design have uncovered hundreds of genetic loci involved in tumorigenesis and provided rich resources for identifying risk factors and biomarkers associated with cancer susceptibility. However, the application of GWAS in determining the genetic architecture of cancer survival remains unestablished. Here, we systematically evaluated genetic effects at the genome-wide level on cancer survival that included overall survival (OS) and cancer-specific survival (CSS), leveraging data deposited in the UK Biobank cohort of a total of 19 628 incident patients across 17 cancer types. Furthermore, we assessed the causal effects of risk factors and circulating biomarkers on cancer prognosis via a Mendelian randomization (MR) analytic framework, which integrated cancer survival GWAS dataset, along with phenome-wide association study (PheWAS) and blood genome-wide gene expression/DNA methylation quantitative trait loci (eQTL/meQTL) datasets. On average, more than 10 traits, 700 genes, and 4,500 CpG sites were prone to cancer prognosis. Finally, we developed a user-friendly online database, SUrvival related cancer Multi-omics database via MEndelian Randomization (SUMMER; http://njmu-edu.cn:3838/SUMMER/), to help users query, browse, and download cancer survival results. In conclusion, SUMMER provides an important resource to assist the research community in understanding the genetic mechanisms of cancer survival.

Extensive germline-somatic interplay contributes to prostate cancer progression through HNF1B co-option of TMPRSS2-ERG

Genome-wide association studies have identified 270 loci conferring risk for prostate cancer (PCa), yet the underlying biology and clinical impact remain to be investigated. Here we observe an enrichment of transcription factor genes including HNF1B within PCa risk-associated regions. While focused on the 17q12/HNF1B locus, we find a strong eQTL for HNF1B and multiple potential causal variants involved in the regulation of HNF1B expression in PCa. An unbiased genome-wide co-expression analysis reveals PCa-specific somatic TMPRSS2-ERG fusion as a transcriptional mediator of this locus and the HNF1B eQTL signal is ERG fusion status dependent. We investigate the role of HNF1B and find its involvement in several pathways related to cell cycle progression and PCa severity. Furthermore, HNF1B interacts with TMPRSS2-ERG to co-occupy large proportion of genomic regions with a remarkable enrichment of additional PCa risk alleles. We finally show that HNF1B co-opts ERG fusion to mediate mechanistic and biological effects of the PCa risk-associated locus 17p13.3/VPS53/FAM57A/GEMIN4. Taken together, we report an extensive germline-somatic interaction between TMPRSS2-ERG fusion and genetic variations underpinning PCa risk association and progression.

SNP-Target Genes Interaction Perturbing the Cancer Risk in the Post-GWAS

Cancer ranks as the second leading cause of death worldwide, and, being a genetic disease, it is highly heritable. Over the past few decades, genome-wide association studies (GWAS) have identified many risk-associated loci harboring hundreds of single nucleotide polymorphisms (SNPs). Some of these cancer-associated SNPs have been revealed as causal, and the functional characterization of the mechanisms underlying the cancer risk association has been illuminated in some instances. In this review, based on the different positions of SNPs and their modes of action, we discuss the mechanisms underlying how SNPs regulate the expression of target genes to consequently affect tumorigenesis and the development of cancer.

The current state of genetic risk models for the development of kidney cancer: a review and validation

OBJECTIVE

To review the current state of genetic risk models for predicting the development of kidney cancer, by identifying and comparing the performance of published models.

METHODS

Risk models were identified from a recent systematic review and the Cancer-PRS web directory. A narrative synthesis of the models, previous validation studies and related genome-wide association studies (GWAS) was carried out. The discrimination and calibration of the identified models was then assessed and compared in the UK Biobank (UKB) cohort (cases, 452; controls, 487 925).

RESULTS

A total of 39 genetic models predicting the development of kidney cancer were identified and 31 were validated in the UKB. Several of the genetic-only models (seven of 25) and most of the mixed genetic-phenotypic models (five of six) had some discriminatory ability (area under the receiver operating characteristic curve >0.5) in this cohort. In general, models containing a larger number of genetic variants identified in GWAS performed better than models containing a small number of variants associated with known causal pathways. However, the performance of the included models was consistently poorer than genetic risk models for other cancers.

CONCLUSIONS

Although there is potential for genetic models to identify those at highest risk of developing kidney cancer, their performance is poorer than the best genetic risk models for other cancers. This may be due to the comparatively small number of genetic variants associated with kidney cancer identified in GWAS to date. The development of improved genetic risk models for kidney cancer is dependent on the identification of more variants associated with this disease. Whether these will have utility within future kidney cancer screening pathways is yet to determined.

Research on the effects of rs1800566 C/T polymorphism of NAD(P)H quinone oxidoreductase 1 gene on cancer risk involves analysis of 43,736 cancer cases and 56,173 controls

Objective

A two-electron reductase known as NQO1 [NAD(P)H quinone oxidoreductase 1] is regarded as an excellent anticancer target. Studies have found that rs1800566 polymorphism of NQO1 is linked to different cancers, but their associations remain controversial.

Methods

In the present work, we selected to do a comprehensive meta-analysis to analyze their correlation. We performed searches on PubMed, Embase, Google Scholar, Chinese database, and Web of Science. The results we obtained covered all publications before April 3, 2022.

Results

There were 176 case-control studies among them, with 56,173 corresponding controls and 43,736 cancer cases. We determined that the NQO1 rs1800566 polymorphism was not related to the cancer risk by calculating 95% confidence intervals and odds ratios. However, stratified genotyping showed that this polymorphism was protective against hepatocellular carcinoma, renal cell carcinoma, and gastric cancer. In addition, on dividing cancer into six systems, the association with gastrointestinal cancer decreased. In the race-based subgroup, a decreasing trend was observed in Asians, while an increasing trend was found among Caucasians, Africans, and mixed populations. The decreased correlation in the hospital-based subgroup was also detected.

Conclusion

Current study shows that rs1800566 polymorphism of NQO1 was linked to cancer susceptibility and maybe as a tumor marker in their development.