Identification of CELSR2 as a novel prognostic biomarker for hepatocellular carcinoma

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

A prognostic mathematical model based on tumor microenvironment-related genes expression for breast cancer patients

Background

Tumor microenvironment (TME) status is closely related to breast cancer (BC) prognosis and systemic therapeutic effects. However, to date studies have not considered the interactions of immune and stromal cells at the gene expression level in BC as a whole. Herein, we constructed a predictive model, for adjuvant decision-making, by mining TME molecular expression information related to BC patient prognosis and drug treatment sensitivity.

Methods

Clinical information and gene expression profiles were extracted from The Cancer Genome Atlas (TCGA), with patients divided into high- and low-score groups according to immune/stromal scores. TME-related prognostic genes were identified using Kaplan-Meier analysis, functional enrichment analysis, and protein-protein interaction (PPI) networks, and validated in the Gene Expression Omnibus (GEO) database. Least absolute shrinkage and selection operator (LASSO) Cox regression analysis was used to construct and verify a prognostic model based on TME-related genes. In addition, the patients' response to chemotherapy and immunotherapy was assessed by survival outcome and immunohistochemistry (IPS). Immunohistochemistry (IHC) staining laid a solid foundation for exploring the value of novel therapeutic target genes.

Results

By dividing patients into low- and high-risk groups, a significant distinction in overall survival was found (p < 0.05). The risk model was independent of multiple clinicopathological parameters and accurately predicted prognosis in BC patients (p < 0.05). The nomogram-integrated risk score had high prediction accuracy and applicability, when compared with simple clinicopathological features. As predicted by the risk model, regardless of the chemotherapy regimen, the survival advantage of the low-risk group was evident in those patients receiving chemotherapy (p < 0.05). However, in patients receiving anthracycline (A) therapy, outcomes were not significantly different when compared with those receiving no-A therapy (p = 0.24), suggesting these patients may omit from A-containing adjuvant chemotherapy. Our risk model also effectively predicted tumor mutation burden (TMB) and immunotherapy efficacy in BC patients (p < 0.05).

Conclusion

The prognostic score model based on TME-related genes effectively predicted prognosis and chemotherapy effects in BC patients. The model provides a theoretical basis for novel driver-gene discover in BC and guides the decision-making for the adjuvant treatment of early breast cancer (eBC).

METTL3-dependent m6A methylation facilitates uterine receptivity and female fertility via balancing estrogen and progesterone signaling

Infertility is a worldwide reproductive health problem and there are still many unknown etiologies of infertility. In recent years, increasing evidence emerged and confirmed that epigenetic regulation played a leading role in reproduction. However, the function of m⁶A modification in infertility remains unknown. Here we report that METTL3-dependent m⁶A methylation plays an essential role in female fertility via balancing the estrogen and progesterone signaling. Analysis of GEO datasets reveal a significant downregulation of METTL3 expression in the uterus of infertile women with endometriosis or recurrent implantation failure. Conditional deletion of Mettl3 in female reproductive tract by using a Pgr-Cre driver results in infertility due to compromised uterine endometrium receptivity and decidualization. m⁶A-seq analysis of the uterus identifies the 3'UTR of several estrogen-responsive genes with METTL3-dependent m⁶A modification, like Elf3 and Celsr2, whose mRNAs become more stable upon Mettl3 depletion. However, the decreased expression levels of PR and its target genes, including Myc, in the endometrium of Mettl3 cKO mice indicate a deficiency in progesterone responsiveness. In vitro, Myc overexpression could partially compensate for uterine decidualization failure caused by Mettl3 deficiency. Collectively, this study reveals the role of METTL3-dependent m⁶A modification in female fertility and provides insight into the pathology of infertility and pregnancy management.

Association between Genetic Variants of CELSR2-PSRC1-SORT1 and Cardiovascular Diseases: A Systematic Review and Meta-Analysis

A cluster of three genes CELSR2, PSRC1, and SORT1 has been associated with cardiovascular diseases. Thus, the aim of this study was (i) to perform a systematic review and updated meta-analysis of the association of three polymorphisms (rs646776, rs599839, and rs464218) of this cluster with cardiovascular diseases, and (ii) to explore by PheWAS signals of the three SNPs in cardiovascular diseases and to evaluate the effect of rs599839 with tissue expression by in silico tools. Three electronic databases were searched to identify eligible studies. The meta-analysis showed that the rs599839 (allelic OR 1.19, 95% CI 1.13–1.26, dominant OR 1.22, 95% CI 1.06–1.39, recessive OR 1.23, 95% CI 1.15–1.32), rs646776 (allelic OR 1.46, 95% CI 1.17–1.82) polymorphisms showed an increased risk for cardiovascular diseases. PheWas analysis showed associations with coronary artery disease and total cholesterol. Our results suggest a possible involvement of the CELSR2-PSRC1-SORT1 cluster variants in the risk association of cardiovascular diseases, particularly coronary artery disease.

Genome-wide characterization leading to simple sequence repeat (SSR) markers development in Shorea robusta

Tropical rainforests in Southeast Asia are enriched by multifarious biota dominated by Dipterocarpaceae. In this family, Shorea robusta is an ecologically sensitive and economically important timber species whose genomic diversity and phylogeny remain understudied due to lack of datasets on genetic resources. Smattering availability of molecular markers impedes population genetic studies indicating a necessity to develop genomic databases and species-specific markers in S. robusta. Accordingly, the present study focused on fostering de novo low-depth genome sequencing, identification of reliable microsatellites markers, and their validation in various populations of S. robusta in Uttarakhand Himalayas. With 69.88 million raw reads assembled into 1,97,489 contigs (read mapped to 93.2%) and a genome size of 357.11 Mb (29 × coverage), Illumina paired-end sequencing technology arranged a library of sequence data of ~ 10 gigabases (Gb). From 57,702 microsatellite repeats, a total of 35,049 simple sequence repeat (SSR) primer pairs were developed. Afterward, among randomly selected 60 primer pairs, 50 showed successful amplification and 24 were found as polymorphic. Out of which, nine polymorphic loci were further used for genetic analysis in 16 genotypes each from three different geographical locations of Uttarakhand (India). Prominently, the average number of alleles per locus (Na), observed heterozygosity (Ho), expected heterozygosity (He), and the polymorphism information content (PIC) were recorded as 2.44, 0.324, 0.277 and 0.252, respectively. The accessibility of sequence information and novel SSR markers potentially enriches the current knowledge of the genomic background for S. robusta and to be utilized in various genetic studies in species under tribe Shoreae.

Mouse-INtraDuctal (MIND): an in vivo model for studying the underlying mechanisms of DCIS malignancy

Due to widespread adoption of screening mammography, there has been a significant increase in new diagnoses of ductal carcinoma in situ (DCIS). However, DCIS prognosis remains unclear. To address this gap, we developed an in vivo model, Mouse-INtraDuctal (MIND), in which patient-derived DCIS epithelial cells are injected intraductally and allowed to progress naturally in mice. Similar to human DCIS, the cancer cells formed in situ lesions inside the mouse mammary ducts and mimicked all histologic subtypes including micropapillary, papillary, cribriform, solid, and comedo. Among 37 patient samples injected into 202 xenografts, at median duration of 9 months, 20 samples (54%) injected into 95 xenografts showed in vivo invasive progression, while 17 (46%) samples injected into 107 xenografts remained non-invasive. Among the 20 samples that showed invasive progression, nine samples injected into 54 xenografts exhibited a mixed pattern in which some xenografts showed invasive progression while others remained non-invasive. Among the clinically relevant biomarkers, only elevated progesterone receptor expression in patient DCIS and the extent of in vivo growth in xenografts predicted an invasive outcome. The Tempus XT assay was used on 16 patient DCIS formalin-fixed, paraffin-embedded sections including eight DCISs that showed invasive progression, five DCISs that remained non-invasive, and three DCISs that showed a mixed pattern in the xenografts. Analysis of the frequency of cancer-related pathogenic mutations among the groups showed no significant differences (KW: p > 0.05). There were also no differences in the frequency of high, moderate, or low severity mutations (KW; p > 0.05). These results suggest that genetic changes in the DCIS are not the primary driver for the development of invasive disease. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Insights into the Association Between QSER1 and M2 Macrophages and Remarkable Malignancy Characteristics in Hepatocellular Carcinoma

Purpose

Glutamine and serine rich 1 (QSER1), as a DNA methylation modulator, play a crucial role in transforming tumor cells. Previous studies have shown that QSER1 plays a role in regulating the progression of various malignancies and that QSER1 dysfunction is connected with precancerous lesions of hepatocellular carcinoma (HCC) as well as HCC prognosis. However, little is known about the detailed contribution of QSER1 in HCC.

Patients and Methods

Various statistical methods such as Kaplan–Meier method, AUC analysis, GSEA, and immune-infiltration analysis were used to evaluate the relationship between QSER1 expression and clinical features, prognostic factors, and potential functional mechanisms of QSER1.

Results

QSER1 expression was negatively correlated with clinicopathological features (clinical stage, pathological grade, TP53 mutation, lymph node metastasis) and clinical outcome (overall survival versus recurrence). Functional enrichment analysis further suggested that QSER1 is involved in multiple pathways related to DNA replication and tumor immunity. TIMER analysis indicated that high QSER1 expression was significantly associated with higher macrophage infiltration and poorer macrophage-related outcomes. In particular, QSER1 was significantly more associated with M2 macrophages than M1 macrophages.

Conclusion

Overall, elevated QSER1 is a potential prognostic marker for HCC and is associated with immune infiltration in HCC.

CELSR2 deficiency suppresses lipid accumulation in hepatocyte by impairing the UPR and elevating ROS level

Cadherin EGF LAG seven-pass G-type receptor 2 (CELSR2), a mammalian orthologue of drosophila flamingo, belongs to the cadherin subfamily. CELSR2 mainly function in neural development and cilium polarity. Recent studies showed that the CELSR2 gene is related to many human diseases, including coronary artery disease, idiopathic scoliosis, and cancer. Genome-Wide Association Studies data showed that SNP in the CELSR2-PSRC1-SORT1 gene loci has a strong association with circulating lipid levels and coronary artery disease. However, the function and underlying mechanism of CELSR2 in hepatic lipid metabolism remain unknown. Here, we found that CELSR2 expression is decreased in the liver of NAFLD/NASH patients and db/db mice. Depletion of CELSR2 significantly decreased the lipid accumulation in hepatocytes by suppressing the expression of lipid synthesis enzymes. Moreover, CELSR2 deficiency impaired the physiological unfolded protein response (UPR), which damages the ER homeostasis, and elevates the reactive oxygen species (ROS) level by decreasing the antioxidant expression. Scavenging of ROS by N-acetylcysteine treatment could restore the decreased lipid accumulation of CELSR2 knockdown cells. Furthermore, CELSR2 loss impaired cell survival by suppressing cell proliferation and promoting apoptosis. Our results uncovered a new role of CELSR2 in regulating lipid homeostasis and UPR, suggesting CELSR2 may be a new therapeutic target for non-alcoholic fatty liver disease.

Genome skimming-based simple sequence repeat (SSR) marker discovery and characterization in Grevillea robusta

Proteaceae, a largely southern hemisphere family consisting of 80 genera distributed in Australia and southern Africa as its centres of greatest diversity, also extends well in northern and southern America. Under this family, Grevillea robusta is a fast-growing species got popularity in farm and avenue plantations. Despite the ecological and economic importance, the species has not yet been investigated for its genetic improvement and genome-based studies. Only a few molecular markers are available for the species or its close relatives, which hinders  genomic and population genetics studies. Genetic markers have been intensively applied for the main strategies in breeding programs, especially for the economically important traits. Hence, it is of utmost priority to develop genomic database resources and species-specific markers for studying quantitative genetics in G. robusta. Given this, the present study aimed to develop de novo genome sequencing, robust microsatellites markers, sequence annotation and their validation in different stands of G. robusta in northern India. Library preparation and sequencing were carried out using Illumina paired-end sequencing technology. Approximately, ten gigabases (Gb) sequence data with 70.87 million raw reads assembled into 425,923 contigs (read mapped to 76.48%) comprising 455 Mb genome size (23 × coverage) generated through genome skimming approach. In total, 9421 simple sequence repeat (SSR) primer pairs were successfully designed from 13,335 microsatellite repeats. Afterward, a subset of 161 primer pairs was randomly selected, synthesized and validated. All the tested primers showed successful amplification but only 13 showed polymorphisms. The polymorphic SSRs were further used to estimate the measures of genetic diversity in 12 genotypes each from the states of Punjab, Haryana, Himachal Pradesh and Uttarakhand. Importantly, the average number of alleles (Na), observed heterozygosity (Ho), expected heterozygosity (He), and the polymorphism information content (PIC) were recorded as 2.69, 0.356, 0.557 and 0.388, respectively. The availability of sequence information and newly developed SSR markers could potentially be used in various genetic analyses and improvements through molecular breeding strategies for G. robusta.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01035-w.

An immune-related signature that to improve prognosis prediction of breast cancer

Although the classic molecular subtype of breast cancer (BRCA) has been widely used in clinical diagnosis, as a highly heterogeneous malignant tumor, the classic scheme is not enough to accurately predict the prognosis of breast cancer patients. Immune cells in the tumor microenvironment (TME) are thought to play a paramount role in tumor development and driving poor prognosis. In this study, we aimed to develop a TME-associated, immune-related signature to improve prognosis prediction of BRCA. BRCA_OURS enriched transcriptomic RNA sequencing (RNA-seq) of tumor tissue was acquired from 43 breast cancer patients before any treatment. On the immune gene profiles of 43 patients from BRCA_OURS and 932 BRCA patients from The Cancer Genome Atlas (TCGA), we identified a robust immune-related signature including one positive coefficients gene (IL-10) and other 9 genes (C14orf79, C1orf168, C1orf226, CELSR2, FABP7, FGFBP1, KLRB1, PLEKHO1, and RAC2), of which the negative coefficients suggesting higher expression were correlated with better prognosis. Based on the expression of these genes, patients were grouped into the high- and low-risk group with significant overall survival (OS) (P<0.0001). The high-risk group was likely to have inferior outcomes related to several important cancer-associated pathways, including mobilizing more Golgi vesicle-mediated transport and intensive DNA double-strand breaking, which are closely related to the infiltration of immune cells and holds the key for further growing and metastasizing. Collectively, our results highlight that the immunological value within BRCA is an essential determinant of prognostic factor. Our signature may provide an effective risk stratification tool for clinical prognosis assessment of patients with BRCA.