Ground source heat pump (GSHP) system has been installed as the air-conditioning system worldwide due to it has the characteristics of high efficiency, easy access and environmental protection. Since ground heat exchanger (GHE) plays a key role in the performance of GSHPs, many models of GHE have been proposed to simulate temperature distribution around the borehole. However, most of these models depict only the heat conduction process between buried pipes and surrounding soil based on the line source model or cylindrical source model. And these models do not consider water transfer under the action of heat source, which can cause some prediction errors. The objective of this study is to provide a numerical model to simulate the spatiotemporal distribution of temperature and moisture caused by a GHE with constant temperature in unsaturated soils. The numerical model is developed by establishing two tridiagonal matrices and adopting Thomas algorithm to achieve the programming. The experiment is operated at the Taiyuan University of Technology and the comparisons between modeled and experimental data prove the high accuracy of this model. The model shows significant engineering values in designs and operations of GSHP.

Myocardial infarction (MI) induced by acute coronary arterial occlusion is usually secondary to atherosclerotic plaque rupture. Dysregulated response of vascular smooth muscle cells (VSMCs) in atherosclerotic plaques may promote plaque rupture. Cadherins (CDHs) form adherens junctions and are known stabilizers of atherosclerotic plaques. To date, the expression patterns of cadherin have not been well investigated in MI aortic VSMCs. We aimed to investigate the expression of cadherin genes in the aortic wall of patients with and without MI. Laser capture microdissected VSMCs were obtained from aortic tissue samples of patients undergoing coronary artery bypass graft surgery. Integrative bioinformatic analysis of the microarray profiles of the VSMCs revealed that MI is discriminated at the whole transcriptome level by hundreds of differentially expressed genes, including genes involved in cell adhesion, of which the cadherin superfamily genes were among the top structural category. Eleven significantly deregulated candidates of the cadherin superfamily were chosen and formed a new classifier that collectively discriminated MI versus non-MI with approximately 95% accuracy. Significance validation was performed with an independent cohort using quantitative RT-qPCR, confirming overexpression of CDH2, CDH12, PCDH17 and PCDH18 in MI VSMCs. The dysregulation of these cadherin superfamily genes might be related to an MI-induced remote effect on aortic wall VSMCs and to imbalances in signaling pathways and myocardial repair mechanisms. Although pathophysiological significance of our findings requires functional studies, mRNA upregulation of the identified cadherin superfamily members in VSMCs might be associated with the progression of atherosclerosis and angiogenesis activation in MI.

In cohort studies the outcome is often time to a particular event, and subjects are followed at regular intervals. Periodic visits may also monitor a secondary irreversible event influencing the event of primary interest, and a significant proportion of subjects develop the secondary event over the period of follow-up. The status of the secondary event serves as a time-varying covariate, but is recorded only at the times of the scheduled visits, generating incomplete time-varying covariates. While information on a typical time-varying covariate is missing for entire follow-up period except the visiting times, the status of the secondary event are unavailable only between visits where the status has changed, thus interval-censored. One may view interval-censored covariate of the secondary event status as missing time-varying covariates, yet missingness is partial since partial information is provided throughout the follow-up period. Current practice of using the latest observed status produces biased estimators, and the existing missing covariate techniques cannot accommodate the special feature of missingness due to interval censoring. To handle interval-censored covariates in the Cox proportional hazards model, we propose an available-data estimator, a doubly robust-type estimator as well as the maximum likelihood estimator via EM algorithm and present their asymptotic properties. We also present practical approaches that are valid. We demonstrate the proposed methods using our motivating example from the Northern Manhattan Study.

OBJECTIVES: Diabetic nephropathy (DN) is one of the most common complications in patients with diabetes, and the discovery of novel targeted therapeutic approaches for DN treatment still faces severe challenges. In the current study, we aimed to discover a novel natural product for potential DN treatment and determine its molecular mechanisms. MATERIALS AND METHODS: Methylthiazoltetrazolium (MTT) assay was employed to evaluate cell viability. Transmission electron microscopy, GFP-LC3 fluorescence fusion plasmid, and Annexin V/PI apoptosis assay were carried out to determine cellular autophagy and apoptosis. Moreover, quantitative proteomics and bioinformatics analysis, Western blotting, and RNA interference were performed to investigate potential molecular mechanisms. RESULTS: Hispidulin displayed protective capacity on the high-glucose-induced podocyte injury models by activating autophagy and inhibiting apoptosis. The mechanism for hispidulin-induced autophagy was associated to Pim1 inhibition and the regulation of Pim1-p21-mTOR signaling axis. Moreover, quantitative proteomics and bioinformatics analysis revealed that the hispidulin-regulated Pim1 inhibition was associated to RAB18, NRas, PARK7, and FIS1. CONCLUSIONS: These results indicate that hispidulin induces autophagy and inhibits apoptosis induced by high glucose in murine podocytes. This study will illuminate future developments in DN-targeted therapy.

We report two Egyptian male genomes (EGP1 and EGP2) sequenced at ~ 30× sequencing depths. EGP1 had 4.7 million variants, where 198,877 were novel variants while EGP2 had 209,109 novel variants out of 4.8 million variants. The mitochondrial haplogroup of the two individuals were identified to be H7b1 and L2a1c, respectively. We also identified the Y haplogroup of EGP1 (R1b) and EGP2 (J1a2a1a2 > P58 > FGC11). EGP1 had a mutation in the NADH gene of the mitochondrial genome ND4 (m.11778 G > A) that causes Leber's hereditary optic neuropathy. Some SNPs shared by the two genomes were associated with an increased level of cholesterol and triglycerides, probably related with Egyptians obesity. Comparison of these genomes with African and Western-Asian genomes can provide insights on Egyptian ancestry and genetic history. This resource can be used to further understand genomic diversity and functional classification of variants as well as human migration and evolution across Africa and Western-Asia.

Studies in E. sinensis have shown that ubiquitination mediated by Cullin-RING E3 ligases (CRLs) plays important roles in spermatogenesis. In other species, CRLs are also essential in cell cycle progression, DNA replication, signal transduction, gene transcription, and development. The catalytic RING component, the RING box protein, is an important part of CRLs. However, there have been few studies on CRLs in crustaceans. In this study, we cloned two RING box protein genes from the Chinese mitten crab, Eriocheir sinensis, termed Es-RBX1 and Es-RBX2 The full length Es-RBX1 cDNA comprises 741 nucleotides, and encodes a protein of 124 amino acid residues, whereas the Es-RBX2 cDNA comprises 1325 nucleotides, and encodes a protein of 110 amino acid residues. Bioinformatics analysis showed that the domains and structure of the RBX proteins have been highly conserved during evolution. Quantitative real-time polymerase chain reaction and western blotting showed that Es-RBX1 is highly expressed in the testis, particularly during the spermatocyte stage, whereas Es-RBX2 did not show specific expression in the male reproductive system. Furthermore, Es-RBX1 is mainly distributed in the nucleus, and changed its location with the development of the nucleus. Co-immunoprecipitation showed that Es-RBX1 could bind Cullin4. These results suggested that Es-RBX1 plays a key role in spermatogenesis of E. sinensis though forming a complex with Cullin4.

MicroRNAs (miRNAs) are non-coding RNAs that regulate gene expression at the posttranscriptional level. In this study, the function of microRNA-7 (miR-7) in host-virus interaction was investigated. Replication of White spot syndrome virus (WSSV) was enhanced with the overexpression of miR-7 and inhibited with the downregulation of miR-7 by using anti-miRNA oligonucleotide AMO-miR-7. The target gene of miR-7 was predicted using bioinformatics methods. Results showed that crab myeloid differentiation factor 88 (Myd88) could be targeted by miR-7. When the expression of Myd88 was knocked down by sequence-specific siRNA, WSSV copies in crabs were significantly increased. Further findings revealed that knockdown of Myd88, Tube, or Pelle inhibited the expressions of interleukin enhancer-binding factor 2 homolog (ILF2) and interleukin-16-like gene (IL-16L). While ILF2 was silenced, IL-16L expression was inhibited. The overexpression of miR-7 inhibited the expressions of ILF2 and IL-16L. Moreover, when ILF2 or IL-16L was silenced, WSSV copies in crabs were increased. Thus, the upregulated expression of miR-7 during WSSV challenge suppressed the host Myd88-ILF2-(IL-16L) signaling pathway in crabs and enhanced WSSV replication. Our study indicated that WSSV utilized crab miR-7 to enhance virus replication during infection.

PURPOSE: Male lower urinary tract symptoms (LUTS) secondary to benign prostatic hyperplasia (BPH) is common in men and can have negative effects on quality of life (QOL). It is the hope that this Guideline becomes a reference on the effective evidence-based surgical management of LUTS/BPH. MATERIALS AND METHODS: The evidence team searched Ovid MEDLINE, the Cochrane Library, and the Agency for Healthcare Research and Quality (AHRQ) database to identify studies indexed between January 2007 and September 2017. When sufficient evidence existed, the body of evidence was assigned a strength rating of A (high), B (moderate), or C (low) for support of Strong, Moderate, or Conditional Recommendations. In the absence of sufficient evidence, additional information is provided as Clinical Principles and Expert Opinions (table 1 in supplementary unabridged guideline, http://jurology.com/). RESULTS: This Guideline provides updated, evidence-based recommendations regarding management of LUTS/BPH utilizing surgery and minimally invasive surgical therapies; additional statements are made regarding diagnostic and pre-operative tests. Clinical statements are made in comparison to what is generally accepted as the gold standard (i.e. transurethral resection of the prostate [TURP]- monopolar and/or bipolar). This guideline is designed to be used in conjunction with the associated treatment algorithm. CONCLUSIONS: The prevalence and the severity of LUTS increases as men age and is an important diagnosis in the healthcare of patients and the welfare of society. This document will undergo additional literature reviews and updating as the knowledge regarding current treatments and future surgical options continues to expand.

Many signaling proteins consist of globular domains connected by flexible linkers that allow for substantial domain motion. Because these domains often serve as complementary functional modules, the possibility of functionally important domain motions arises. To explore this possibility, we require knowledge of the ensemble of protein conformations sampled by interdomain motion. Measurements of NMR residual dipolar couplings (RDCs) of backbone HN bonds offer a per-residue characterization of interdomain dynamics, as the couplings are sensitive to domain orientation. A challenge in reaching this potential is the need to interpret the RDCs as averages over dynamic ensembles of domain conformations. Here, we address this challenge by introducing an efficient protocol for generating conformational ensembles appropriate for flexible, multi-domain proteins. The protocol uses map-restrained self-guided Langevin dynamics simulations to promote collective, interdomain motion while restraining the internal domain motion to near rigidity. Critically, the simulations retain an all-atom description for facile inclusion of site-specific NMR RDC restraints. The result is the rapid generation of conformational ensembles consistent with the RDC data. We illustrate this protocol on human Pin1, a two-domain peptidyl-prolyl isomerase relevant for cancer and Alzheimer's disease. The results include the ensemble of domain orientations sampled by Pin1, as well as those of a dysfunctional variant, I28A-Pin1. The differences between the ensembles corroborate our previous spin relaxation results that showed weakened interdomain contact in the I28A variant relative to wild type. Our protocol extends our abilities to explore the functional significance of protein domain motions.

Loading of small RNAs into Argonaute, the core protein in RNA silencing, requires the Hsp70/Hsp90 chaperone machinery. This machinery also activates many other clients, including steroid hormone receptors and kinases, but how their structures change during chaperone-dependent activation remains unclear. Here, we utilized single-molecule Förster resonance energy transfer (smFRET) to probe the conformational changes of Drosophila Ago2 mediated by the chaperone machinery. We found that empty Ago2 exists in various closed conformations. The Hsp70 system (Hsp40 and Hsp70) and the Hsp90 system (Hop, Hsp90, and p23) together render Ago2 into an open, active form. The Hsp70 system, but not the Hsp90 system alone, is sufficient for Ago2 to partially populate the open form. Instead, the Hsp90 system is required to extend the dwell time of Ago2 in the open state, which must be transiently primed by the Hsp70 system. Our data uncover distinct and coordinated actions of the chaperone machinery, where the Hsp70 system expands the structural ensembles of Ago2 and the Hsp90 system captures and stabilizes the active form.

Animal cells have a remarkable capacity to adopt durable and heritable gene expression programs or epigenetic states that define the physical properties and diversity of somatic cell types. The maintenance of epigenetic programs depends on poorly understood pathways that prevent gain or loss of inherited signals. In the germline, epigenetic factors are enriched in liquid-like perinuclear condensates called nuage. Here, we identify the deeply conserved helicase-domain protein, ZNFX-1, as an epigenetic regulator and component of nuage that interacts with Argonaute systems to balance epigenetic inheritance. Our findings suggest that ZNFX-1 promotes the 3' recruitment of machinery that propagates the small RNA epigenetic signal and thus counteracts a tendency for Argonaute targeting to shift 5' along the mRNA. These functional insights support the idea that recently identified subdomains of nuage, including ZNFX-1 granules or "Z-granules," may define spatial and temporal zones of molecular activity during epigenetic regulation.

Ligand binding sites in proteins are often localized to deeply buried cavities, inaccessible to bulk solvent. Yet, in many cases binding of cognate ligands occurs rapidly. An intriguing system is presented by the L99A cavity mutant of T4 Lysozyme (T4L L99A) that rapidly binds benzene (~106 M-1s-1). Although the protein has long served as a model system for protein thermodynamics and crystal structures of both free and benzene-bound T4L L99A are available, the kinetic pathways by which benzene reaches its solvent-inaccessible binding cavity remain elusive. The current work, using extensive molecular dynamics simulation, achieves this by capturing the complete process of spontaneous recognition of benzene by T4L L99A at atomistic resolution. A series of multi-microsecond unbiased molecular dynamics simulation trajectories unequivocally reveal how benzene, starting in bulk solvent, diffuses to the protein and spontaneously reaches the solvent inaccessible cavity of T4L L99A. The simulated and high-resolution X-ray derived bound structures are in excellent agreement. A robust four-state Markov model, developed using cumulative 60 μs trajectories, identifies and quantifies multiple ligand binding pathways with low activation barriers. Interestingly, none of these identified binding pathways required large conformational changes for ligand access to the buried cavity. Rather, these involve transient but crucial opening of a channel to the cavity via subtle displacements in the positions of key helices (helix4/helix6, helix7/helix9) leading to rapid binding. Free energy simulations further elucidate that these channel-opening events would have been unfavorable in wild type T4L. Taken together and via integrating with results from experiments, these simulations provide unprecedented mechanistic insights into the complete ligand recognition process in a buried cavity. By illustrating the power of subtle helix movements in opening up multiple pathways for ligand access, this work offers an alternate view of ligand recognition in a solvent-inaccessible cavity, contrary to the common perception of a single dominant pathway for ligand binding.

Extracellular matrix (ECM) proteins have been shown to play important roles regulating multiple biological processes in an array of organ systems, including the cardiovascular system. By using a novel bioinformatics text-mining tool, we studied six categories of cardiovascular disease (CVD), namely ischemic heart disease (IHD), cardiomyopathies (CM), cerebrovascular accident (CVA), congenital heart disease (CHD), arrhythmias (ARR), and valve disease (VD), anticipating novel ECM protein-disease and protein-protein relationships hidden within vast quantities of textual data. We conducted a phrase-mining analysis, delineating the relationships of 709 ECM proteins with the six groups of CVDs reported in 1,099,254 abstracts. The technology pipeline known as Context-aware Semantic Online Analytical Processing (CaseOLAP) was applied to semantically rank the association of proteins to each and all six CVDs, performing analyses to quantify each protein-disease relationship. We performed principal component analysis and hierarchical clustering of the data, where each protein is visualized as a six dimensional vector. We found that ECM proteins display variable degrees of association with the six CVDs; certain CVDs share groups of associated proteins whereas others have divergent protein associations. We identified 82 ECM proteins sharing associations with all six CVDs. Our bioinformatics analysis ascribed distinct ECM pathways (via Reactome) from this subset of proteins, namely insulin-like growth factor regulation and interleukin-4 and interleukin-13 signaling, suggesting their contribution to the pathogenesis of all six CVDs. Finally, we performed hierarchical clustering analysis and identified protein clusters associated with a targeted CVD; analyses revealed unexpected insights underlying ECM-pathogenesis of CVDs.

Neural precursor cell expressed, developmentally downregulated 9 (NEDD9) supports oncogenic signaling in a number of solid and hematologic tumors. Little is known about the role of NEDD9 in ovarian carcinoma (OC), but available data suggest elevated mRNA and protein expression in advanced stage high-grade cancers. We used a transgenic MISIIR-TAg mouse OC model combined with genetic ablation of Nedd9 to investigate its action in the development and progression of OC. A Nedd9-/- genotype delayed tumor growth rate, reduced incidence of ascites, and reduced expression and activation of signaling proteins including SRC, STAT3, E-cadherin, and AURKA. Cell lines established from MISIIR-TAg;Nedd9-/- and MISIIR-TAg;Nedd9+/+ mice exhibited altered migration and invasion. Growth of these cells in a syngeneic allograft model indicated that systemic Nedd9 loss in the microenvironment had little impact on tumor allograft growth, but in a Nedd9 wild-type background Nedd9-/- allografts exhibited significantly reduced growth, dissemination, and oncogenic signaling compared to Nedd9+/+ allografts. Gene expression analysis revealed that Nedd9+/+ tumors exhibited more mesenchymal "stem-like" transcriptional program, including increased expression of Aldh1a1 and Aldh1a2. Conversely, loss of Nedd9 resulted in increased expression of differentiation genes, including fallopian tube markers Foxj1, Ovgp1, and Pax8. Collectively, these data suggest that tumor cell-intrinsic Nedd9 expression promotes OC development and progression by broad induction of oncogenic protein signaling and stem/mesenchymal gene expression.

AIM: Non-small-cell lung cancer (NSCLC) is a leading global health threat that impairs patient health outcomes. Health state utilities are fundamental values in economic evaluation and significantly vary across countries. Given the scarce data on the Chinese population, the current study measured utility values in the Chinese patients with NSCLC. METHODS: This study was conducted as a cross-sectional survey of patients with advanced NSCLC at the Shanghai Chest Hospital. Utility values were assessed using the EuroQol five-dimension (EQ-5D) instrument and scored based on the Chinese-specific value algorithm. Predictors of utility values were examined using a subgroup analysis and a multiple regression model. RESULTS: The mean EQ-5D utility value of recruited patients was 0.814. The regression analysis revealed that tumor stage, treatment regimen and line of therapy were the potential predictors of utility values. CONCLUSION: This study provides the Chinese-specific health utility data for advanced NSCLC using the EQ-5D.

Mass-spectrometry-based proteomics has become the standard approach for identifying and quantifying proteins. A vital step consists of analyzing experimentally generated mass spectra to identify the underlying peptide sequences for later mapping to the originating proteins. We here present the latest developments in SearchGUI, a common open-source interface for the most frequently used freely available proteomics search and de novo engines that has evolved into a central component in numerous bioinformatics workflows.

PURPOSE: MicroRNAs (miRNAs) have been shown to be involved in the initiation and progression of glioma. However, the underlying molecular mechanisms are still unclear. METHODS: We performed microarray analysis to evaluate miRNA expression levels in 158 glioma tissue samples, and examined miR-1231 levels in glioma samples and healthy brain tissues using qRT-PCR. In vitro analyses were performed using miR-1231 mimics, inhibitors, and siRNA targeting EGFR. We used flow cytometry, CCK-8 assays, and colony formation assays to examine glioma proliferation and cell cycle analysis. A dual luciferase reporter assay was performed to examine miR-1231 regulation of EGFR, and the effect of upregulated miR-1231 was investigated in a subcutaneous GBM model. RESULTS: We found that miR-1231 expression was decreased in human glioma tissues and negatively correlated with EGFR levels. Moreover, the downregulation of miR-1231 negatively correlated with the clinical stage of human glioma patients. miR-1231 overexpression dramatically downregulated glioma cell proliferation, and suppressed tumor growth in a nude mouse model. Bioinformatics prediction and a luciferase assay confirmed EGFR as a direct target of miR-1231. EGFR overexpression abrogated the suppressive effect of miR-1231 on the PI3K/AKT pathway and G1 arrest. CONCLUSIONS: Taken together, these results demonstrated that EGFR is a direct target of miR-1231. Our findings suggest that the miR-1231/EGFR axis may be a helpful future diagnostic target for malignant glioma.

AIMS: To investigate whether PDCD1 gene polymorphisms are functional, and their associations with T1D risk and related clinical characteristics. METHODS: A total of 3060 Chinese Han individuals (1019 T1D patients and 2041 healthy controls) were genotyped for 4 tag single nucleotide polymorphisms (SNPs) within the PDCD1 region (rs2227982, rs7421861, rs10204525, and rs6710479) and another most studied synonymous SNP, rs2227981. In addition, 251 healthy individuals underwent an oral glucose tolerance test (OGTT); measures of insulin release and sensitivity were estimated from insulinogenic, BIGTT, Matsuda. Further, we performed in silico bioinformatics analysis to explore potential functional annotation of the investigated SNPs in PDCD1 gene. RESULTS: Both rs2227982 and rs2227981 polymorphisms were associated with T1D risk in Chinese Han population under additive model (OR = 0.84, 95% CI 0.75-0.93 and OR = 1.23, 95% CI 1.08-1.40, respectively), but not the other three SNPs in PDCD1 gene. Our meta-analysis revealed that rs2227982 and rs2227981 polymorphisms also have significant associations with T1D risk in East Asians (OR = 0.82, 95% CI 0.74-0.90 and OR = 1.23, 95% CI 1.12-1.36, respectively), but not Europeans. And the T allele of rs2227982 polymorphism is associated with increased 30 min post OGTT glucose level (P = 0.023) and 120 min post OGTT insulin level (P = 0.033). Furthermore, the genetic and regulatory architecture suggested all the 5 investigated SNPs in PDCD1 are putatively functional. CONCLUSIONS: Both rs2227982 and rs2227981 polymorphisms were associated with T1D risk in East Asians, and rs2227982 also had a significant association with glycemic traits, which suggested PDCD1 gene polymorphisms might participate in facilitating T1D risk.

Phytoplankton growth is limited in vast oceanic regions by the low bioavailability of iron. Iron fertilization often results in diatom blooms, yet the physiological underpinnings for how diatoms survive in chronically iron-limited waters and outcompete other phytoplankton when iron becomes available are unresolved. We show that some diatoms can use siderophore-bound iron, and exhibit a species-specific recognition for siderophore types. In Phaeodactylum tricornutum, hydroxamate siderophores are taken up without previous reduction by a high-affinity mechanism that involves binding to the cell surface followed by endocytosis-mediated uptake and delivery to the chloroplast. The affinity recorded is the highest ever described for an iron transport system in any eukaryotic cell. Collectively, our observations suggest that there are likely a variety of iron uptake mechanisms in diatoms besides the well-established reductive mechanism. We show that iron starvation-induced protein 1 (ISIP1) plays an important role in the uptake of siderophores, and through bioinformatics analyses we deduce that this protein is largely diatom-specific. We quantify expression of ISIP1 in the global ocean by querying the Tara Oceans atlas of eukaryotic genes and show a link between the abundance and distribution of diatom-associated ISIP1 with ocean provinces defined by chronic iron starvation.

PCR-based RNA splicing assays are commonly used in diagnostic and research settings to assess the potential effects of variants of uncertain clinical significance in BRCA1 and BRCA2. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium completed a multicentre investigation to evaluate differences in assay design and the integrity of published data, raising a number of methodological questions associated with cell culture conditions and PCR-based protocols. We utilized targeted RNA-seq to re-assess BRCA1 and BRCA2 mRNA isoform expression patterns in lymphoblastoid cell lines (LCLs) previously used in the multicentre ENIGMA study. Capture of the targeted cDNA sequences was carried out using 34 BRCA1 and 28 BRCA2 oligonucleotides from the Illumina Truseq Targeted RNA Expression platform. Our results show that targeted RNA-seq analysis of LCLs overcomes many of the methodology limitations associated with PCR-based assays leading us to make the following observations and recommendations: (1) technical replicates (n > 2) of variant carriers to capture methodology induced variability associated with RNA-seq assays, (2) LCLs can undergo multiple freeze/thaw cycles and can be cultured up to 2 weeks without noticeably influencing isoform expression levels, (3) nonsense-mediated decay inhibitors are essential prior to splicing assays for comprehensive mRNA isoform detection, (4) quantitative assessment of exon:exon junction levels across BRCA1 and BRCA2 can help distinguish between normal and aberrant isoform expression patterns. Experimentally derived recommendations from this study will facilitate the application of targeted RNA-seq platforms for the quantitation of BRCA1 and BRCA2 mRNA aberrations associated with sequence variants of uncertain clinical significance.

The Polymerase Associated Factor 1 complex (PAF1c) is an epigenetic co-modifying complex that directly contacts RNA polymerase II (RNAPII) and several epigenetic regulating proteins. Mutations, overexpression and loss of expression of subunits of the PAF1c are observed in various forms of cancer suggesting proper regulation is needed for cellular development. However, the biochemical interactions with the PAF1c that allow dynamic gene regulation are unclear. We and others have shown that the PAF1c makes a direct interaction with MLL fusion proteins, which are potent oncogenic drivers of acute myeloid leukemia (AML). This interaction is critical for the maintenance of MLL translocation driven AML by targeting MLL fusion proteins to the target genes Meis1 and Hoxa9. Here, we use a proteomics approach to identify protein-protein interactions with the PAF1c subunit CDC73 that regulate the function of the PAF1c. We identified a novel interaction with a histone H3 lysine 9 (H3K9) methyltransferase protein, SETDB1. This interaction is stabilized with a mutant CDC73 that is incapable of supporting AML cell growth. Importantly, transcription of Meis1 and Hoxa9 is reduced and promoter H3K9 trimethylation (H3K9me3) increased by overexpression of SETDB1 or stabilization of the PAF1c-SETDB1 interaction in AML cells. These findings were corroborated in human AML patients where increased SETDB1 expression was associated with reduced HOXA9 and MEIS1. To our knowledge, this is the first proteomics approach to search for CDC73 protein-protein interactions in AML, and demonstrates that the PAF1c may play a role in H3K9me3-mediated transcriptional repression in AML.

Exosomes are a potential source of cancer biomarkers. Probing tumor-derived exosomes can offer a potential non-invasive way to diagnose cancer, assess cancer progression, and monitor treatment responses. Novel molecular methods would facilitate exosome analysis and accelerate basic and clinical exosome research. Methods: A standard gold-coated glass microscopy slide was used to develop a miniaturized affinity-based device to capture exosomes in a target-specific manner with the assistance of low-cost 3-D printing technology. Gold nanorods coated with QSY21 Raman reporters were used as the label agent to quantitatively detect the target proteins based on surface enhanced Raman scattering spectroscopy. The expressions of several surface protein markers on exosomes from conditioned culture media of breast cancer cells and from HER2-positive breast cancer patients were quantitatively measured. The data was statistically analyzed and compared with healthy controls. Results: A miniaturized 17 × 5 Au array device with 2-mm well size was fabricated to capture exosomes in a target-specific manner and detect the target proteins on exosomes with surface enhanced Raman scattering gold nanorods. This assay can specifically detect exosomes with a limit of detection of 2×106 exosomes/mL and analyze over 80 purified samples on a single device within 2 h. Using the assay, we have showed that exosomes derived from MDA-MB-231, MDA-MB-468, and SKBR3 breast cancer cells give distinct protein profiles compared to exosomes derived from MCF12A normal breast cells. We have also showed that exosomes in the plasma from HER2-positive breast cancer patients exhibit significantly (P ≤ 0.01) higher level of HER2 and EpCAM than those from healthy donors. Conclusion: We have developed a simple, inexpensive, highly efficient, and portable Raman exosome assay for detection and protein profiling of exosomes. Using the assay and model exosomes from breast cancer cells, we have showed that exosomes exhibit diagnostic surface protein markers, reflecting the protein profile of their donor cells. Through proof-of-concept studies, we have identified HER2 and EpCAM biomarkers on exosomes in plasma from HER2-positive breast cancer patients, suggesting the diagnostic potential of these markers for breast cancer diagnostics. This assay would accelerate exosome research and pave a way to the development of novel cancer liquid biopsy for cancer detection and monitoring.

Background and Aims: Improving liver regeneration (LR) capacity and thereby liver function reserve is a critical bridging strategy for managing liver failure patients. Since estrogen signaling may participate in LR, our aim was to characterize the roles of ERα and ERβ in LR. Methods: LR capacity and estradiol levels following 2/3rd partial hepatectomy (PHx) were compared in ERα-KO or ERβ-KO vs. wildtype mice. The ERα- or ERβ-related transcriptome and interactome were analyzed from regenerating livers, and then bioinformatics was used for pathway discovery and analysis of interactome-transcriptome relationships. Human hepatic progenitors (HepRG cells) and mouse Hepa1-6 hepatocytes were used to elucidate molecular interactions and functions. Results: This paper demonstrated that estrogen signals orchestrate hepatic repopulation and differentiation via distinct transcriptome patterns governed by ERα or ERβ. Cell repopulation pathway was associated with the ERα-transcriptome, but cell differentiation and metabolic function were associated with the ERβ transcriptome. Mechanistic studies linking ERs interactomes and transcriptomes discovered that ERα-Chd1 interaction promoted cell growth by upregulating Ssxb6, Crygc, and Cst1; and, ERβ-Ube3a interaction facilitated hepatic progenitor cell differentiation to hepatocytes and cholangiocytes, specifically by upregulating Ifna5. Conclusions: ERα and ERβ orchestrate liver cell proliferation and differentiation respectively, thereby promoting LR.

Comparing metagenomic samples is crucial for understanding microbial communities. For different groups of microbial communities, such as human gut metagenomic samples from patients with a certain disease and healthy controls, identifying group-specific sequences offers essential information for potential biomarker discovery. A sequence that is present, or rich, in one group, but absent, or scarce, in another group is considered "group-specific" in our study. Our main purpose is to discover group-specific sequence regions between control and case groups as disease-associated markers. We developed a long k-mer (k ≥ 30 bps)-based computational pipeline to detect group-specific sequences at strain resolution free from reference sequences, sequence alignments, and metagenome-wide de novo assembly. We called our method MetaGO: Group-specific oligonucleotide analysis for metagenomic samples. An open-source pipeline on Apache Spark was developed with parallel computing. We applied MetaGO to one simulated and three real metagenomic datasets to evaluate the discriminative capability of identified group-specific markers. In the simulated dataset, 99.11% of group-specific logical 40-mers covered 98.89% disease-specific regions from the disease-associated strain. In addition, 97.90% of group-specific numerical 40-mers covered 99.61 and 96.39% of differentially abundant genome and regions between two groups, respectively. For a large-scale metagenomic liver cirrhosis (LC)-associated dataset, we identified 37,647 group-specific 40-mer features. Any one of the features can predict disease status of the training samples with the average of sensitivity and specificity higher than 0.8. The random forests classification using the top 10 group-specific features yielded a higher AUC (from ∼0.8 to ∼0.9) than that of previous studies. All group-specific 40-mers were present in LC patients, but not healthy controls. All the assembled 11 LC-specific sequences can be mapped to two strains of Veillonella parvula: UTDB1-3 and DSM2008. The experiments on the other two real datasets related to Inflammatory Bowel Disease and Type 2 Diabetes in Women consistently demonstrated that MetaGO achieved better prediction accuracy with fewer features compared to previous studies. The experiments showed that MetaGO is a powerful tool for identifying group-specific k-mers, which would be clinically applicable for disease prediction. MetaGO is available at https://github.com/VVsmileyx/MetaGO.

Clinical interpretation of germline missense variants represents a major challenge, including those in the TP53 Li-Fraumeni syndrome gene. Bioinformatic prediction is a key part of variant classification strategies. We aimed to optimize the performance of the Align-GVGD tool used for p53 missense variant prediction, and compare its performance to other bioinformatic tools (SIFT, PolyPhen-2) and ensemble methods (REVEL, BayesDel). Reference sets of assumed pathogenic and assumed benign variants were defined using functional and/or clinical data. Area under the curve and Matthews correlation coefficient (MCC) values were used as objective functions to select an optimized protein multisequence alignment with best performance for Align-GVGD. MCC comparison of tools using binary categories showed optimized Align-GVGD (C15 cut-off) combined with BayesDel (0.16 cut-off), or with REVEL (0.5 cut-off), to have the best overall performance. Further, a semi-quantitative approach using multiple tiers of bioinformatic prediction, validated using an independent set of nonfunctional and functional variants, supported use of Align-GVGD and BayesDel prediction for different strength of evidence levels in ACMG/AMP rules. We provide rationale for bioinformatic tool selection for TP53 variant classification, and have also computed relevant bioinformatic predictions for every possible p53 missense variant to facilitate their use by the scientific and medical community.