BMC Genomics

Single nucleus DNA sequencing of flow sorted archived frozen and formalin fixed paraffin embedded solid tumors

Single-cell and transcriptomic analyses reveal the influence of diabetes on ovarian cancer

Abstract Background There has been a significant surge in the global prevalence of diabetes mellitus (DM), which increases the susceptibility of individuals to ovarian cancer (OC). However, the relationship between DM and OC remains largely unexplored. The objective of this study is to provide preliminary insights into the shared molecular regulatory mechanisms and potential biomarkers between DM and OC. Methods Multiple datasets from the GEO database were utilized for bioinformatics analysis. Single cell datasets from the GEO database were analysed. Subsequently, immune cell infiltration analysis was performed on mRNA expression data. The intersection of these datasets yielded a set of common genes associated with both OC and DM. Using these overlapping genes and Cytoscape, a protein‒protein interaction (PPI) network was constructed, and 10 core targets were selected. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were then conducted on these core targets. Additionally, advanced bioinformatics analyses were conducted to construct a TF-mRNA-miRNA coregulatory network based on identified core targets. Furthermore, immunohistochemistry staining (IHC) and real-time quantitative PCR (RT-qPCR) were employed for the validation of the expression and biological functions of core proteins, including HSPAA1, HSPA8, SOD1, and transcription factors SREBF2 and GTAT2, in ovarian tumors. Results The immune cell infiltration analysis based on mRNA expression data for both DM and OC, as well as analysis using single-cell datasets, reveals significant differences in mononuclear cell levels. By intersecting the single-cell datasets, a total of 119 targets related to mononuclear cells in both OC and DM were identified. PPI network analysis further identified 10 hub genesincludingHSP90AA1, HSPA8, SNRPD2, UBA52, SOD1, RPL13A, RPSA, ITGAM, PPP1CC, and PSMA5, as potential targets of OC and DM. Enrichment analysis indicated that these genes are primarily associated with neutrophil degranulation, GDP-dissociation inhibitor activity, and the IL-17 signaling pathway, suggesting their involvement in the regulation of the tumor microenvironment. Furthermore, the TF-gene and miRNA-gene regulatory networks were validated using NetworkAnalyst. The identified TFs included SREBF2, GATA2, and SRF, while the miRNAs included miR-320a, miR-378a-3p, and miR-26a-5p. Simultaneously, IHC and RT-qPCR reveal differential expression of core targets in ovarian tumors after the onset of diabetes. RT-qPCR further revealed that SREBF2 and GATA2 may influence the expression of core proteins, including HSP90AA1, HSPA8, and SOD1. Conclusion This study revealed the shared gene interaction network between OC and DM and predicted the TFs and miRNAs associated with core genes in monocytes. Our research findings contribute to identifying potential biological mechanisms underlying the relationship between OC and DM.

Combined PD-1/PD-L1 and tumor-infiltrating immune cells redefined a unique molecular subtype of high-grade serous ovarian carcinoma

AbstractBackgroundHigh-grade serous ovarian carcinoma is highly heterogeneous, and although many studies have been conducted to identify high-grade serous ovarian carcinoma molecular subtypes that are sensitive to immunotherapy, no precise molecular subtype has been proposed to date. Immune cell infiltration and immune checkpoints are highly correlated with immunotherapy. Here, we investigated immune cell infiltration and immune checkpoint values for prognosis and precise immunotherapy for high-grade serous ovarian carcinoma based on molecular subtype classification.Results“High antigen-presenting cells infiltration molecular subtype of high-grade serous ovarian carcinoma” was identified in immune cell infiltration profiles. Each of the three immune cell infiltration clusters (A, B, and C) demonstrated distinct immune cell characterization, with immune cell infiltration cluster C exhibiting high antigen-presenting cell infiltration, improved prognosis, and higher sensitivity to immunotherapy. Programmed death-1/programmed death ligand 1 has a prognostic and predictive role that can help classify molecular subtypes.ConclusionsOur findings redefined a unique molecular subtype of high-grade serous ovarian carcinoma, suggesting that high-grade serous ovarian carcinoma patients with higher antigen-presenting cell infiltration and programmed death-1/programmed death ligand 1 expression can benefit from precise immunotherapy.

Using Portuguese BRCA pathogenic variation as a model to study the impact of human admixture on human health

Abstract Background Admixture occurs between different ethnic human populations. The global colonization in recent centuries by Europeans led to the most significant admixture in human history. While admixture may enhance genetic diversity for better fitness, it may also impact on human health by transmitting genetic variants for disease susceptibility in the admixture population. The admixture by Portuguese global exploration initiated in the 15th century has reached over 20 million of Portuguese-heritage population worldwide. It provides a valuable model to study the impact of admixture on human health. BRCA1 and BRCA2 (BRCA) are two of the important tumor suppressor genes. The pathogenic variation (PV) in BRCA is well determined to cause high risk of hereditary breast and ovarian cancer. Tracing the distribution of Portuguese BRCA PV in Portuguese-heritage population will help to understand the impact of admixture on cancer susceptibility in modern humans. In this study, we analyzed the distribution of the Portuguese-originated BRCA variation in Brazilian population, which has high degree Portuguese-heritage. Methods By comprehensive data mining, standardization and annotation, we generated a Portuguese-derived BRCA variation dataset and a Brazilian-derived BRCA variation dataset. We compared the two BRCA variation datasets to identify the BRCA variants shared between the two populations. Results The Portuguese-derived BRCA variation dataset consists of 220 BRCA variants including 78 PVs from 11,482 Portuguese cancer patients, 93 (42.2%) in BRCA1 and 127 (57.7%) in BRCA2. Of the 556 Portuguese BRCA PV carriers carrying the 78 PVs, 331 (59.5%) carried the three Portuguese-BRCA founder PVs of BRCA1 c.2037delinsCC, BRCA1 c.3331_3334del and BRCA2 c.156_157insAlu. The Brazilian-derived BRCA variation dataset consists of 255 BRCA PVs from 7,711 cancer patients, 136 (53.3%) in BRCA1 and 119 (46.6%) in BRCA2. We developed an open database named dbBRCA-Portuguese (https://genemutation.fhs.um.edu.mo/dbbrca-portuguese/) and an open database named dbBRCA-Brazilian (https://genemutation.fhs.um.edu.mo/dbbrca-brazilian) to host the BRCA variation data from Portuguese and Brazilian populations. We compared the BRCA PV datasets between Portuguese and Brazilian populations, and identified 29 Portuguese-specific BRCA PVs shared between Portuguese and Brazilian populations, 14 in BRCA1 including the Portuguese founder BRCA1 c.3331_3334del and BRCA1 c.2037delinsCC, and 15 in BRCA2 including the Portuguese founder BRCA2 c.156_157insAlu. Searching the 78 Portuguese BRCA PVs in over 5,000 ancient human genomes identified evolution origin for only 8 PVs in Europeans dated between 37,470 and 3,818 years before present, confirming the Portuguese-specificity of Portuguese BRCA PVs; comparing the 78 Portuguese BRCA PVs Portuguese, 255 Brazilian BRCA PVs, and 134 African BRCA PVs showed little overlapping, ruling out the possibility that the BRCA PVs shared between Portuguese and Brazilian may also be contributed by African. Conclusion Our study provides evidence that the admixture in recent human history contributed to cancer susceptibility in modern humans.

Upregulation of CALD1 predicted a poor prognosis for platinum-treated ovarian cancer and revealed it as a potential therapeutic resistance target

Abstract Background Ovarian cancer (OC) has the worst prognosis among gynecological malignancies, most of which are found to be in advanced stage. Cell reduction surgery based on platinum-based chemotherapy is the current standard of treatment for OC, but patients are prone to relapse and develop drug resistance. The objective of this study was to identify a specific molecular target responsible for platinum chemotherapy resistance in OC. Results We screened the protein-coding gene Caldesmon (CALD1), expressed in cisplatin-resistant OC cells in vitro. The prognostic value of CALD1 was evaluated using survival curve analysis in OC patients treated with platinum therapy. The diagnostic value of CALD1 was verified by drawing a Receiver Operating Characteristic (ROC) curve using clinical samples from OC patients. This study analyzed data from various databases including Gene Expression Omnibus (GEO), Human Protein Atlas (HPA), The Cancer Cell Line Encyclopedia (CCLE), The Cancer Genome Atlas (TCGA), GEPIA 2, UALCAN, Kaplan–Meier (KM) plotter, LinkedOmics database, and String. Different expression genes (DEGs) between cisplatin-sensitive and cisplatin-resistant cells were acquired respectively from 5 different datasets of GEO. CALD1 was selected as a common gene from 5 groups DEGs. Online data analysis of HPA and CCLE showed that CALD1 was highly expressed in both normal ovarian tissue and OC. In TCGA database, high expression of CALD1 was associated with disease stage and venous invasion in OC. Patients with high CALD1 expression levels had a worse prognosis under platinum drug intervention, according to Kaplan–Meier (KM) plotter analysis. Analysis of clinical sample data from GEO showed that CALD1 had superior diagnostic value in distinguishing patients with platinum "resistant" and platinum "sensitive" (AUC = 0.816), as well as patients with worse progression-free survival (AUC = 0.741), and those with primary and omental metastases (AUC = 0.811) in ovarian tumor. At last, CYR61 was identified as a potential predictive molecule that may play an important role alongside CALD1 in the development of platinum resistance in OC. Conclusions CALD1, as a member of cytoskeletal protein, was associated with poor prognosis of platinum resistance in OC, and could be used as a target protein for mechanism study of platinum resistance in OC.

Long-read sequencing reveals the structural complexity of genomic integration of HPV DNA in cervical cancer cell lines

Abstract Background Cervical cancer (CC) causes more than 311,000 deaths annually worldwide. The integration of human papillomavirus (HPV) is a crucial genetic event that contributes to cervical carcinogenesis. Despite HPV DNA integration is known to disrupt the genomic architecture of both the host and viral genomes in CC, the complexity of this process remains largely unexplored. Results In this study, we conducted whole-genome sequencing (WGS) at 55-65X coverage utilizing the PacBio long-read sequencing platform in SiHa and HeLa cells, followed by comprehensive analyses of the sequence data to elucidate the complexity of HPV integration. Firstly, our results demonstrated that PacBio long-read sequencing effectively identifies HPV integration breakpoints with comparable accuracy to targeted-capture Next-generation sequencing (NGS) methods. Secondly, we constructed detailed models of complex integrated genome structures that included both the HPV genome and nearby regions of the human genome by utilizing PacBio long-read WGS. Thirdly, our sequencing results revealed the occurrence of a wide variety of genome-wide structural variations (SVs) in SiHa and HeLa cells. Additionally, our analysis further revealed a potential correlation between changes in gene expression levels and SVs on chromosome 13 in the genome of SiHa cells. Conclusions Using PacBio long-read sequencing, we have successfully constructed complex models illustrating HPV integrated genome structures in SiHa and HeLa cells. This accomplishment serves as a compelling demonstration of the valuable capabilities of long-read sequencing in detecting and characterizing HPV genomic integration structures within human cells. Furthermore, these findings offer critical insights into the complex process of HPV16 and HPV18 integration and their potential contribution to the development of cervical cancer.

A general index for linear and nonlinear correlations for high dimensional genomic data

Abstract Background With the advance of high throughput sequencing, high-dimensional data are generated. Detecting dependence/correlation between these datasets is becoming one of most important issues in multi-dimensional data integration and co-expression network construction. RNA-sequencing data is widely used to construct gene regulatory networks. Such networks could be more accurate when methylation data, copy number aberration data and other types of data are introduced. Consequently, a general index for detecting relationships between high-dimensional data is indispensable. Results We proposed a Kernel-Based RV-coefficient, named KBRV, for testing both linear and nonlinear correlation between two matrices by introducing kernel functions into RV2 (the modified RV-coefficient). Permutation test and other validation methods were used on simulated data to test the significance and rationality of KBRV. In order to demonstrate the advantages of KBRV in constructing gene regulatory networks, we applied this index on real datasets (ovarian cancer datasets and exon-level RNA-Seq data in human myeloid differentiation) to illustrate its superiority over vector correlation. Conclusions We concluded that KBRV is an efficient index for detecting both linear and nonlinear relationships in high dimensional data. The correlation method for high dimensional data has possible applications in the construction of gene regulatory network.

Causal association between 637 human metabolites and ovarian cancer: a mendelian randomization study

Abstract Background Current evidence suggests a significant association between metabolites and ovarian cancer (OC); however, the causal relationship between the two remains unclear. This study employs Mendelian randomization (MR) to investigate the causal effects between different metabolites and OC. Methods In this study, a total of 637 metabolites were selected as the exposure variables from the Genome-wide Association Study (GWAS) database (http://gwas.mrcieu.ac.uk/datasets/). The OC related GWAS dataset (ieu-b-4963) was chosen as the outcome variable. R software and the TwoSampleMR package were utilized for the analysis in this study. MR analysis employed the inverse variance-weighted method (IVW), MR-Egger and weighted median (WM) for regression fitting, taking into consideration potential biases caused by linkage disequilibrium and weak instrument variables. Metabolites that did not pass the tests for heterogeneity and horizontal pleiotropy were considered to have no significant causal effect on the outcome. Steiger’s upstream test was used to determine the causal direction between the exposure and outcome variables. Results The results from IVW analysis revealed that a total of 31 human metabolites showed a significant causal effect on OC (P < 0.05). Among them, 9 metabolites exhibited consistent and stable causal effects, which were confirmed by Steiger’s upstream test (P < 0.05). Among these 9 metabolites, Androsterone sulfate, Propionylcarnitine, 5alpha-androstan-3beta,17beta-diol disulfate, Total lipids in medium VLDL and Concentration of medium VLDL particles demonstrated a significant positive causal effect on OC, indicating that these metabolites promote the occurrence of OC. On the other hand, X-12,093, Octanoylcarnitine, N2,N2-dimethylguanosine, and Cis-4-decenoyl carnitine showed a significant negative causal association with OC, suggesting that these metabolites can inhibit the occurrence of OC. Conclusions The study revealed the complex effect of metabolites on OC through Mendelian randomization. As promising biomarkers, these metabolites are worthy of further clinical validation.

Fine mapping of RNA isoform diversity using an innovative targeted long-read RNA sequencing protocol with novel dedicated bioinformatics pipeline

Solving the structure of mRNA transcripts is a major challenge for both research and molecular diagnostic purposes. Current approaches based on short-read RNA sequencing and RT-PCR techniques cannot fully explore the complexity of transcript structure. The emergence of third-generation long-read sequencing addresses this problem by solving this sequence directly. However, genes with low expression levels are difficult to study with the whole transcriptome sequencing approach. To fix this technical limitation, we propose a novel method to capture transcripts of a gene panel using a targeted enrichment approach suitable for Pacific Biosciences and Oxford Nanopore Technologies platforms. We designed a set of probes to capture transcripts of a panel of genes involved in hereditary breast and ovarian cancer syndrome. We present SOSTAR (iSofOrmS annoTAtoR), a versatile pipeline to assemble, quantify and annotate isoforms from long read sequencing using a new tool specially designed for this application. The significant enrichment of transcripts by our capture protocol, together with the SOSTAR annotation, allowed the identification of 1,231 unique transcripts within the gene panel from the eight patients sequenced. The structure of these transcripts was annotated with a resolution of one base relative to a reference transcript. All major alternative splicing events of the BRCA1 and BRCA2 genes described in the literature were found. Complex splicing events such as pseudoexons were correctly annotated. SOSTAR enabled the identification of abnormal transcripts in the positive controls. In addition, a case of unexplained inheritance in a family with a history of breast and ovarian cancer was solved by identifying an SVA retrotransposon in intron 13 of the BRCA1 gene. We have validated a new protocol for the enrichment of transcripts of interest using probes adapted to the ONT and PacBio platforms. This protocol allows a complete description of the alternative structures of transcripts, the estimation of their expression and the identification of aberrant transcripts in a single experiment. This proof-of-concept opens new possibilities for RNA structure exploration in both research and molecular diagnostics.