Par6G suppresses cell proliferation and is targeted by loss-of-function mutations in multiple cancers

Biliary atresia: insights into mechanisms using a toxic model of the disease including Wnt and Hippo signaling pathways and microtubules

BACKGROUND

Mechanisms underlying bile duct injury in biliary atresia (BA) remain unclear and mechanisms of bile duct repair are unknown. This study aimed to explore the roles of microtubule instability and Wnt and Hippo signaling pathways in a biliatresone-induced BA model.

METHODS

Using primary murine neonatal cholangiocytes in both 2D and 3D cultures, and ex-vivo extra hepatic bile ducts (EHBD) which also has peri-cholangiocyte area, we analyzed injury and recovery processes. Injury was induced by the toxin biliatresone and recovery was induced by toxin wash-out.

RESULTS

Microtubule stabilizer paclitaxel prevented biliatresone-induced injury, both to cholangiocytes as well as reduced periductal αSMA stain, this process is mediated by decreased glutathione levels. RhoU and Wnt11 (Wnt signaling) and Pard6g and Amotl1 (Hippo signaling) are involved in both injury and recovery processes, with the latter acting upstream to Wnt signaling.

CONCLUSIONS

Early stages of biliatresone-induced EHBD injury in cholangiocytes and periductal structures are reversible. Wnt and Hippo signaling pathways play crucial roles in injury and recovery, providing insights into BA injury mechanisms and potential recovery avenues.

IMPACT

Microtubule stabilization prevents cholangiocyte injury and lumen obstruction in a toxic model of biliary atresia (biliatresone induced). Early stages of biliatresone-induced injury, affecting both cholangiocytes and periductal structures, are reversible. Both Wnt and Hippo signaling pathways play a crucial role in bile duct injury and recovery, with a noted interplay between the two. Understanding mechanisms of cholangiocyte recovery is imperative to unveil potential therapeutic avenues.

Quantitative relationships of FAM50B and PTCHD3 methylation with reduced intelligence quotients in school aged children exposed to lead: Evidence from epidemiological and in vitro studies

At present, the application of DNA methylation (DNAm) biomarkers in environmental health risk assessment (EHRA) is more challenging due to the unclearly quantitative relationship between them. We aimed to explore the role of FAM50B and PTCHD3 at the level of signaling pathways, and establish the quantitative relationship between them and children's intelligence quotients (IQs). DNAm of target regions was measured in multiple cell models and was compared with the human population data. Then the dose-response relationships of lead exposure with neurotoxicity and DNAm were established by benchmark dose (BMD) model, followed by potential signaling pathway screening. Results showed that there was a quantitative linear relationship between children's IQs and FAM50B/PTCHD3 DNAm (DNAm between 51.40 % - 78.78 % and 31.41 % - 74.19 % for FAM50B and PTCHD3, respectively), and this relationship was more significant when children's IQs > 90. The receiver operating characteristic (ROC) and calibration curves showed that FAM50B/PTCHD3 DNAm had a satisfying accuracy and consistency in predicting children's IQs, which was confirmed by sensitivity analysis of gender and CpG site grouping data. In cell experiments, there was also a quantitative linear relationship between FAM50B DNAm and reactive oxygen species (ROS) production, which was mediated by PI3K-AKT signaling pathway. In addition, the lead BMD of ROS was close to that of FAM50B DNAm, suggesting that FAM50B DNAm was a suitable biomarker for the risk assessments of adverse outcomes induced by lead. Taken collectively, these results suggest that FAM50B/PTCHD3 can be applied to EHRA and the prevention/intervention of adverse effects of lead on children's IQs.

Genome-Wide Integrative Transcriptional Profiling Identifies Age-Associated Signatures in Dogs

Mammals experience similar stages of embryonic development, birth, infancy, youth, adolescence, maturity, and senescence. While embryonic developmental processes have been extensively researched, many molecular mechanisms regulating the different life stages after birth, such as aging, remain unresolved. We investigated the conserved and global molecular transitions in transcriptional remodeling with age in dogs of 15 breeds, which revealed that genes underlying hormone level regulation and developmental programs were differentially regulated during aging. Subsequently, we show that the candidate genes associated with tumorigenesis also exhibit age-dependent DNA methylation patterns, which might have contributed to the tumor state through inhibiting the plasticity of cell differentiation processes during aging, and ultimately suggesting the molecular events that link the processes of aging and cancer. These results highlight that the rate of age-related transcriptional remodeling is influenced not only by the lifespan, but also by the timing of critical physiological milestones.

Orchestration of tissue-scale mechanics and fate decisions by polarity signalling

Eukaryotic development relies on dynamic cell shape changes and segregation of fate determinants to achieve coordinated compartmentalization at larger scale. Studies in invertebrates have identified polarity programmes essential for morphogenesis; however, less is known about their contribution to adult tissue maintenance. While polarity-dependent fate decisions in mammals utilize molecular machineries similar to invertebrates, the hierarchies and effectors can differ widely. Recent studies in epithelial systems disclosed an intriguing interplay of polarity proteins, adhesion molecules and mechanochemical pathways in tissue organization. Based on major advances in biophysics, genome editing, high-resolution imaging and mathematical modelling, the cell polarity field has evolved to a remarkably multidisciplinary ground. Here, we review emerging concepts how polarity and cell fate are coupled, with emphasis on tissue-scale mechanisms, mechanobiology and mammalian models. Recent findings on the role of polarity signalling for tissue mechanics, micro-environmental functions and fate choices in health and disease will be summarized.

Recurrent HBV Integration Targets as Potential Drivers in Hepatocellular Carcinoma

Chronic hepatitis B virus (HBV) infection is the major etiology of hepatocellular carcinoma (HCC), frequently with HBV integrating into the host genome. HBV integration, found in 85% of HBV-associated HCC (HBV–HCC) tissue samples, has been suggested to be oncogenic. Here, we investigated the potential of HBV–HCC driver identification via the characterization of recurrently targeted genes (RTGs). A total of 18,596 HBV integration sites from our in-house study and others were analyzed. RTGs were identified by applying three criteria: at least two HCC subjects, reported by at least two studies, and the number of reporting studies. A total of 396 RTGs were identified. Among the 28 most frequent RTGs, defined as affected in at least 10 HCC patients, 23 (82%) were associated with carcinogenesis and 5 (18%) had no known function. Available breakpoint positions from the three most frequent RTGs, TERT, MLL4/KMT2B, and PLEKHG4B, were analyzed. Mutual exclusivity of TERT promoter mutation and HBV integration into TERT was observed. We present an RTG consensus through comprehensive analysis to enable the potential identification and discovery of HCC drivers for drug development and disease management.

DNA Methylation Biomarkers of IQ Reduction are Associated with Long-term Lead Exposure in School Aged Children in Southern China

Although lead associated with intelligence decline in children has long been reported, studies combining intelligence determination, molecular mechanisms exploration, and biomarker screen are quite rare. In this study, based on 333 children aged 9-11, we determined the role of DNA methylation (DNAm) in the relationship of lead exposure with children's intelligence. DNAm was measured from children's blood DNA specimens, and mediation analysis was performed to identify DNAm biomarkers mediating the lead-intelligence relationship. We identified forty-three differentially methylated regions (DMRs), and two fragments (FAM50B1 and PTCHD3) significantly mediated the lead-intelligence relationship, with contribution rates of 30.36% (p = 0.010) and 60.36% (p < 0.001), respectively. In addition, blood lead levels (BLLs) lower than 100 μg/L still adversely affected children's IQs and DNAm of the two fragments. Our data suggests that DNAm mediates lead-associated cognitive delay in children and blood lead reference value for school-aged children (100 μg/L) should be revised, and the candidate biomarkers can be used in related neurological diseases in future.

Time-Varying Gene Network Analysis of Human Prefrontal Cortex Development

The prefrontal cortex (PFC) constitutes a large part of the human central nervous system and is essential for the normal social affection and executive function of humans and other primates. Despite ongoing research in this region, the development of interactions between PFC genes over the lifespan is still unknown. To investigate the conversion of PFC gene interaction networks and further identify hub genes, we obtained time-series gene expression data of human PFC tissues from the Gene Expression Omnibus (GEO) database. A statistical model, loggle, was used to construct time-varying networks and several common network attributes were used to explore the development of PFC gene networks with age. Network similarity analysis showed that the development of human PFC is divided into three stages, namely, fast development period, deceleration to stationary period, and recession period. We identified some genes related to PFC development at these different stages, including genes involved in neuronal differentiation or synapse formation, genes involved in nerve impulse transmission, and genes involved in the development of myelin around neurons. Some of these genes are consistent with findings in previous reports. At the same time, we explored the development of several known KEGG pathways in PFC and corresponding hub genes. This study clarified the development trajectory of the interaction between PFC genes, and proposed a set of candidate genes related to PFC development, which helps further study of human brain development at the genomic level supplemental to regular anatomical analyses. The analytical process used in this study, involving the loggle model, similarity analysis, and central analysis, provides a comprehensive strategy to gain novel insights into the evolution and development of brain networks in other organisms.

Secreted breast tumor interstitial fluid microRNAs and their target genes are associated with triple-negative breast cancer, tumor grade, and immune infiltration

Background

Studies on tumor-secreted microRNAs point to a functional role of these in cellular communication and reprogramming of the tumor microenvironment. Uptake of tumor-secreted microRNAs by neighboring cells may result in the silencing of mRNA targets and, in turn, modulation of the transcriptome. Studying miRNAs externalized from tumors could improve cancer patient diagnosis and disease monitoring and help to pinpoint which miRNA-gene interactions are central for tumor properties such as invasiveness and metastasis.

Methods

Using a bioinformatics approach, we analyzed the profiles of secreted tumor and normal interstitial fluid (IF) microRNAs, from women with breast cancer (BC). We carried out differential abundance analysis (DAA), to obtain miRNAs, which were enriched or depleted in IFs, from patients with different clinical traits. Subsequently, miRNA family enrichment analysis was performed to assess whether any families were over-represented in the specific sets. We identified dysregulated genes in tumor tissues from the same cohort of patients and constructed weighted gene co-expression networks, to extract sets of co-expressed genes and co-abundant miRNAs. Lastly, we integrated miRNAs and mRNAs to obtain interaction networks and supported our findings using prediction tools and cancer gene databases.

Results

Network analysis showed co-expressed genes and miRNA regulators, associated with tumor lymphocyte infiltration. All of the genes were involved in immune system processes, and many had previously been associated with cancer immunity. A subset of these, BTLA, CXCL13, IL7R, LAMP3, and LTB, was linked to the presence of tertiary lymphoid structures and high endothelial venules within tumors. Co-abundant tumor interstitial fluid miRNAs within this network, including miR-146a and miR-494, were annotated as negative regulators of immune-stimulatory responses. One co-expression network encompassed differences between BC subtypes. Genes differentially co-expressed between luminal B and triple-negative breast cancer (TNBC) were connected with sphingolipid metabolism and predicted to be co-regulated by miR-23a. Co-expressed genes and TIF miRNAs associated with tumor grade were BTRC, CHST1, miR-10a/b, miR-107, miR-301a, and miR-454.

Conclusion

Integration of IF miRNAs and mRNAs unveiled networks associated with patient clinicopathological traits, and underlined molecular mechanisms, specific to BC sub-groups. Our results highlight the benefits of an integrative approach to biomarker discovery, placing secreted miRNAs within a biological context.

Identification and Classification of Differentially Expressed Genes and Network Meta-Analysis Reveals Potential Molecular Signatures Associated With Tuberculosis

Tuberculosis (TB) is one of deadly transmissible disease that causes death worldwide; however, only 10% of people infected with Mycobacteriumtuberculosis develop disease, indicating that host genetic factors may play key role in determining susceptibility to TB disease. In this way, the analysis of gene expression profiling of TB infected individuals can give us a snapshot of actively expressed genes and transcripts under various conditions. In the present study, we have analyzed microarray data set and compared the gene expression profiles of patients with different datasets of healthy control, latent infection, and active TB. We observed the transition of genes from normal condition to different stages of the TB and identified and annotated those genes/pathways/processes that have important roles in TB disease during its cyclic interventions in the human body. We identified 488 genes that were differentially expressed at various stages of TB and allocated to pathways and gene set enrichment analysis. These pathways as well as GSEA’s importance were evaluated according to the number of DEGs presents in both. In addition, we studied the gene regulatory networks that may help to further understand the molecular mechanism of immune response against the TB infection and provide us a new angle for future biomarker and therapeutic targets. In this study, we identified 26 leading hubs which are deeply rooted from top to bottom in the gene regulatory network and work as the backbone of the network. These leading hubs contains 31 key regulator genes, of which 14 genes were up-regulated and 17 genes were down-regulated. The proposed approach is based on gene-expression profiling, and network analysis approaches predict some unknown TB-associated genes, which can be considered (or can be tested) as reliable candidates for further (in vivo/in vitro) studies.

Molecular subtype classification of papillary renal cell cancer using miRNA expression

Background

Renal papillary cell carcinoma (KIRP) is a relatively rare renal malignancy. Although KIRP subtyping about clinical relevance has been defined, there have been scarce number of studies on the molecular characteristics of KIRP subtypes.

Method

In this study, a independent samples t-test was used to identify differentially expressed (DE) miRNAs between tumor and normal samples of KIRP. Meanwhile, we performed unsupervised clustering using DE miRNA expression data to analyze molecular characteristics of KIRP. The Partitioning Around Medoids clustering approach was used to identify molecular subtypes. The cumulative distribution function (CDF), proportion of ambiguously clustered pairs (PAC), principal component analysis (PCA) and consensus heatmaps were used to assess the optimal subtypes. In the differential molecular subtypes, we performed an integrated analysis of survival, DE genes, biological function and somatic mutations on the cohort of KIRP patients from The Cancer Genome Atlas.

Results

From solutions with 2, 3, 4, 5, 6 and 7 clusters we selected three KIRP molecular subtypes after assessing PCA, PAC, CDF and consensus heatmaps. We found that the three subtypes are associated with different overall survival and molecular characteristics. Compared with subtype1 and subtype3, subtype2 had a better prognosis in KIRP patients. After exploring their signaling pathways and biological characteristics, we identified the significantly enriched KEGG pathways and Gene Ontology terms for the three subtypes. The distribution of PARD6B, SETD2, STAG2, CUL3, TNRC18, LRBA, IGSF9B and DUNC1H1 mutations differed between the subtypes.

Conclusion

We performed unsupervised clustering using differentially expressed miRNA expression data and described the three KIRP molecular subtypes. The three subtypes differed in overall survival, molecular characteristics and gene mutation frequency.

Roles of partitioning-defective protein 6 (Par6) and its complexes in the proliferation, migration and invasion of cancer cells

A pivotal regulator of cell polarity and homeostasis, partitioning-defective protein 6 (Par6) forms multicomponent complexes that not only regulate cell polarity and stabilize cell morphology, but have also been demonstrated to participate in the proliferation, migration and invasion of cancer cells. The transforming growth factor (TGF)-β and extracellular signal-regulated kinase (Erk) 1/2 pathways are the most thoroughly studied pathways involving Par6 in many cancers. Aurothiomalate has been used to disrupt the interaction between Par6 and atypical protein kinase C within the multicomponent complexes, and has been shown to effectively block transformed growth and metastasis in vitro and/or in vivo in a variety of cancers, including pancreatic, prostate and lung cancers, as well as alveolar rhabdomyosarcoma. It is likely that with further revelations regarding the critical roles of Par6 in cancer initiation, progression and metastasis, targeted therapies against Par6 will be discovered and prove effective preclinically, and hopefully clinically, in cancer treatment.

Phenotype-driven identification of epithelial signalling clusters

In metazoans, epithelial architecture provides a context that dynamically modulates most if not all epithelial cell responses to intrinsic and extrinsic signals, including growth or survival signalling and transforming oncogene action. Three-dimensional (3D) epithelial culture systems provide tractable models to interrogate the function of human genetic determinants in establishment of context-dependency. We performed an arrayed genetic shRNA screen in mammary epithelial 3D cultures to identify new determinants of epithelial architecture, finding that the key phenotype impacting shRNAs altered not only the data population average but even more noticeably the population distribution. The broad distributions were attributable to sporadic gene silencing actions by shRNA in unselected populations. We employed Maximum Mean Discrepancy concept to capture similar population distribution patterns and demonstrate here the feasibility of the test in identifying an impact of shRNA in populations of 3D structures. Integration of the clustered morphometric data with protein-protein interactions data enabled hypothesis generation of novel biological pathways underlying similar 3D phenotype alterations. The results present a new strategy for 3D phenotype-driven pathway analysis, which is expected to accelerate discovery of context-dependent gene functions in epithelial biology and tumorigenesis.

Downregulation of PKCζ/Pard3/Pard6b is responsible for lung adenocarcinoma cell EMT and invasion

Atypical protein kinase C ζ (PKCζ) forms an apico-basal polarity complex with Partitioning Defective (Pard) 3 and Pard6 to regulate normal epithelial cell apico-basolateral polarization. The dissociation of the PKCζ/Pard3/Pard6 complex is essential for the disassembly of the tight/adherens junction and epithelial-mesenchymal transition (EMT) that is critical for tumor spreading. Loss of cell polarity and epithelial organization is strongly correlated with malignancy and tumor progression in some other cancer types. However, it is unclear whether the PKCζ/Pard3/Pard6 complex plays a role in the progression of non-small-cell lung cancer (NSCLC). We found that hypoxia downregulated the PKCζ/Pard3/Pard6 complex, correlating with induction of lung cancer cell migration and invasion. Silencing of the PKCζ/Pard3/Pard6 polarity complex components induced lung cancer cell EMT, invasion, and colonization in vivo. Suppression of Pard3 was associated with altered expression of genes regulating wound healing, cell apoptosis/death and cell motility, and particularly upregulation of MAP3K1 and fibronectin which are known to contribute to lung cancer progression. Human lung adenocarcinoma tissues expressed less Pard6b and PKCζ than the adjacent normal tissues and in experimental mouse lung adenocarcinoma, the levels of Pard3 and PKCζ were also decreased. In addition, we showed that a methylation locus in the gene body of Pard3 is positively associated with the expression of Pard3 and that methylation of the Pard3 gene increased cellular sensitivity to carboplatin, a common chemotherapy drug. Suppression of Pard3 increased chemoresistance in lung cancer cells. Together, these results suggest that reduced expression of PKCζ/Pard3/Pard6 contributes to NSCLC EMT, invasion, and chemoresistance.

"Stealth dissemination" of macrophage-tumor cell fusions cultured from blood of patients with pancreatic ductal adenocarcinoma

Here we describe isolation and characterization of macrophage-tumor cell fusions (MTFs) from the blood of pancreatic ductal adenocarcinoma (PDAC) patients. The MTFs were generally aneuploidy, and immunophenotypic characterizations showed that the MTFs express markers characteristic of PDAC and stem cells, as well as M2-polarized macrophages. Single cell RNASeq analyses showed that the MTFs express many transcripts implicated in cancer progression, LINE1 retrotransposons, and very high levels of several long non-coding transcripts involved in metastasis (such as MALAT1). When cultured MTFs were transplanted orthotopically into mouse pancreas, they grew as obvious well-differentiated islands of cells, but they also disseminated widely throughout multiple tissues in "stealth" fashion. They were found distributed throughout multiple organs at 4, 8, or 12 weeks after transplantation (including liver, spleen, lung), occurring as single cells or small groups of cells, without formation of obvious tumors or any apparent progression over the 4 to 12 week period. We suggest that MTFs form continually during PDAC development, and that they disseminate early in cancer progression, forming "niches" at distant sites for subsequent colonization by metastasis-initiating cells.

Expression and prognostic value of the cell polarity PAR complex members in thyroid cancer

Establishment and maintenance of the apical-basal cell polarity, required for proper replication, migration, specialized functions and tissue morphogenesis, relies on three evolutionary conserved complexes: PAR, CRUMBS and SCRIBBLE. Loss of cell polarity/cohesiveness (LOP/C) is implicated in cancer progression, and members of the polarity complex have been described as either oncogenes or oncosuppressors. However, no information on their role in thyroid cancer (TC) progression is available. In the present study, we evaluated the gene expression of the PAR complex members aPKCι, PARD3α/β and PARD6α/β/γ in 95 papillary TC (PTC), compared to their normal matched tissues and in 12 anaplastic TC (ATC). The mRNA and protein levels of investigated genes were altered in the majority of PTC and ATC tissues. In PTC, univariate analysis showed that reduced expression of aPKCι, PARD3β and PARD6γ mRNAs is associated with increased tumor size, and the reduced expression of PARD3β mRNA is associated also with recurrences. Multivariate analysis demonstrated that the presence of lymph node metastasis at diagnosis and the reduced expression of PARD3β are independent risk factors for recurrences, with hazard ratio, respectively, of 8.21 (p=0.006) and 3.04 (p=0.029). The latter result was confirmed by the Kaplan-Meier analysis, which evidenced the association between decreased PARD3β mRNA levels and shorter disease-free interval. In conclusion, we demonstrated that the expression of PAR complex components is deregulated in the majority of PTC and there is a general trend towards their reduction in ATC tissues. Moreover, a prognostic value for the PARD3β gene in PTCs is suggested.

YAP Subcellular Localization and Hippo Pathway Transcriptome Analysis in Pediatric Hepatocellular Carcinoma

Pediatric hepatocellular carcinoma (HCC) is a rare tumor which is associated with an extremely high mortality rate due to lack of effective chemotherapy. Recently, the Hippo pathway and its transcriptional co-activator Yes-associated protein (YAP) have been shown to play a role in hepatocyte proliferation and development of HCC in animal models. Therefore, we sought to examine the activity of YAP and the expression of Hippo pathway components in tumor and non-neoplastic liver tissue from 7 pediatric patients with moderately differentiated HCC. None of the patients had underlying cirrhosis or viral hepatitis, which is commonly seen in adults with HCC. This highlights a major difference in the pathogenesis of HCC between children and adults. We found a statistically significant increase in YAP nuclear localization in 100% of tumors. YAP target gene (CCNE1, CTGF, Cyr61) mRNA expression was also increased in the tumors that had the most significant increase in YAP nuclear localization. Based on Ki67 co-localization studies YAP nuclear localization was not simply a marker of proliferation. Our results demonstrate a clear increase in YAP activity in moderately differentiated pediatric HCC, providing evidence that it may play an important role in tumor survival and propagation.

Cell adhesion and polarity in squamous cell carcinoma of the lung

Lung cancer is a deadly disease that can roughly be classified into three histopathological groups: lung adenocarcinomas, lung squamous cell carcinomas (LSCCs), and small cell carcinomas. These types of lung cancer are molecularly, phenotypically, and regionally diverse neoplasms, reflecting differences in their cells of origin. LSCCs commonly arise in the airway epithelium of a main or lobar bronchus, which is an important line of defence against the external environment. Furthermore, most LSCCs are characterized histopathologically by the presence of keratinization and/or intercellular bridges, consistent with the molecular features of these tumours, characterized by high levels of transcripts encoding keratins and proteins relevant to intercellular junctions and cell polarity. In this review, the relationships between the molecular features of LSCCs and the types of cell and epithelia of origin are discussed. Recurrent alterations in genes involved in intercellular adhesion and cell polarity in LSCCs are also reviewed, emphasizing the importance of the disruption of PAR3 and the PAR complex. Finally, the possible functional effects of these alterations on epithelial homeostasis, and how they contribute to the development of LSCC, are discussed.