Search

Primary Sjögren’s syndrome (pSS) is a chronic autoimmune disease that is estimated to affect 35 million people worldwide and is characterized by lymphocytic infiltration, elevated circulating autoantibodies, and proinflammatory cytokines. The key immune cell subset changes and the TCR/BCR repertoire alterations in pSS patients remain unclear. In this study, we sought to comprehensively characterize the transcriptional changes in PBMCs of pSS patients by single-cell RNA sequencing and single-cell V(D)J sequencing. Naive CD8+ T cells and mucosal-associated invariant T cells were markedly decreased but regulatory T cells were increased in pSS patients. There were a large number of differentially expressed genes shared by multiple subpopulations of T cells and B cells. Abnormal signaling pathways, including Ag processing and presentation, the BCR signaling pathway, the TCR signaling pathway, and Epstein–Barr virus infection, were highly enriched in pSS patients. Moreover, there were obvious differences in the CD30, FLT3, IFN-II, IL-1, IL-2, IL-6, IL-10, RESISTIN, TGF-β, TNF, and VEGF signaling networks between pSS patients and healthy controls. Single-cell TCR and BCR repertoire analysis showed that there was a lower diversity of T cells in pSS patients than in healthy controls; however, there was no significant difference in the degree of clonal expansion, CDR3 length distribution, or degree of sequence sharing. Notably, our results further emphasize the functional importance of αβ pairing in determining Ag specificity. In conclusion, our analysis provides a comprehensive single-cell map of gene expression and TCR/BCR profiles in pSS patients for a better understanding of the pathogenesis, diagnosis, and treatment of pSS.

BACKGROUND Sternal cleft is a greatly rare congenital thoracic deformity, arising from a failure of the sternal bars fusion process that should be completed in the fetal period, the incidence of which is less than 0.15%. CASE REPORT Herein, we present a case report of a newborn girl having a superior congenital sternal cleft. After the baby was born, scar-like tissue was found in the middle of the chest and extended to the root of the umbilical cord. Based on the imaging data, this newborn was diagnosed with sternal cleft belonging to the superior sternal cleft and not associated with other congenital deformities. CONCLUSIONS As a rare congenital thoracic deformity, postpartum diagnosis of the sternal cleft mainly is currently based on medical imaging, including thoracic computed tomography (CT), three-dimensional (3D) reconstruction CT, and magnetic resonance imaging (MRI). Sternum cleft not only affects the aesthetic appearance but also leads to the destruction of the bone structure of the thorax, resulting in opposing thoracic movements. Therefore, early diagnosis and early treatment play significant roles in the treatment of this congenital sternal deformity. Regardless of whether there are clinical symptoms of sternal cleft, primary repair surgery must be done as soon as possible and during the neonatal period is best, in which simple surgical techniques achieve remarkable effects.

Background Gastric cancer is a common and fatal disease with the highest incidence and mortality in men. Protocadherin beta 6 (PCDHB6) has not been reported much in tumors and the role of PCDHB6 in gastric cancer is not yet clear. Therefore, we aimed to investigate the expression of PCDHB6 in gastric cancer infiltration, invasion and metastasis and to assess the prognostic and functional significance of PCDHB6 in gastric cancer. Methods Based on TCGA, GEO database, bioinformatics analysis of PCDHB6 gene was performed from gene expression, survival analysis, gene mutation, immune infiltration, DNA methylation and enrichment analysis. In addition, Real-Time PCR was used to further confirm the role of PCDHB6 in gastric cancer. Results PCDHB6 gene was highly expressed in gastric cancer and positively correlated with poor patient prognosis. Real-Time PCR results also showed high expression of PCDHB6 in gastric cancer. PCDHB6 expression was positively correlated with tumor-associated fibroblasts, Basophils, CD4 + memory T-cells, CD8 + T- cells, Eosinophils, Macrophages, Regulatory T-cells and Type 2 T-helper cells were positively correlated with immune infiltration. KEGG and GO enrichment analysis showed that PCDHB6 gene was mainly involved in the regulation of Herpes simplex virus 1 infection, Hippo signaling pathway and cell adhesion pathway. Conclusions Our study comprehensive summary and analysis revealed for the first time that PCDHB6 is highly expressed in gastric cancer and is an oncogene. PCDHB6 can be used as a potential prognostic biomarker for gastric cancer and provides a powerful therapeutic target for the treatment of gastric cancer.

Transplantation is currently the best treatment for patients with end-stage renal disease. However, acute rejection (AR) is the major source of failure in renal transplantation. The current best practice for the diagnosis of AR involves renal biopsy, but it is invasive, time-consuming, costly and inconvenient. Sensitive and less invasive detection of AR episodes in renal transplant patients is essential to preserve allograft function. The present study applied isobaric tags for relative and absolute quantitation (iTRAQ) mass spectrometry to analyze serum protein expression in patients with AR and healthy controls. Overall, 1,399 proteins were identified. Using a cut-off of Q1.2 for the variation in expression, 109 proteins were identified to be differentially expressed between the AR and control groups, 72 of which were upregulated and 37 were downregulated. Several proteins, including properdin, keratin 1, lipoprotein(a) and vitamin D-binding protein, may have roles in the pathogenesis of AR. The present study focused on iTRAQ-based proteomic profiling of serum samples in AR. Insight from the present study may help advance the understanding of the molecular mechanisms of AR and identify potential novel biomarkers of AR for further characterization.

Posttranslational modifications (PTMs) to histones such as lysine crotonylation are classified as epigenetic changes. Lysine crotonylation participates in various cellular processes and occurs in active promoters, directly accelerating transcription. The present study performed a proteomics analysis of crotonylation between healthy controls and patients with immunoglobulin A (IgA) nephropathy using tandem mass spectrometry and high‑resolution liquid chromatography. The present results identified 353 crotonylated proteins and 770 modification sites, including 155 upregulated and 198 downregulated crotonylated proteins. In total, seven conserved motifs were identified in the present study. The present bioinformatics analysis results suggested a number of the crotonylated proteins exhibited various subcellular localization patterns, such as in the cytoplasm. Protein domains, including thioredoxin, moesin tail and myosin like IQ motif domains were markedly enriched in crotonylated proteins. Kyoto Encyclopedia of Genes and Genomes and functional enrichment analyses suggested significant enrichment of crotonylated proteins in complement and coagulation cascades, and antigen processing and presentation pathways displaying important relationships with IgA nephropathy. The present results suggested that crotonylation occurred in numerous proteins and may play key regulatory roles in IgA nephropathy.

The dynamic immunity of kidney transplant patients has not been fully elucidated. In this study, we explored the repertoire features of B/T cell receptor (BCR/TCR) of kidney transplant patients.Using combined multiplex PCR amplification and high-throughput sequencing technique, we analyzed the uremic patients' BCR H chain and TCR beta chain repertoire which obtained 1 day before kidney transplantation (PRE-1), 1 day and 7 day after kidney transplantation (POST-1 and POST-7).Our analysis results showed the diversity of TCRβ CDR3 in POST-7 group was highest. In addition, there were specific skewed usage of TRBV gene subfamilies, and V-J combinations in different time points during kidney transplantation. Moreover, the overlap degrees of BCR-H (TCR-β) CDR3 repertoire among each group were identified. Notably, the abundance of some TCR-β CDR3 sequences changed regularly in the time point of kidney transplantation.The BCR-H (TCR-β) CDR3 repertoire of kidney transplant patients changed dynamically.

The current pandemic of coronavirus disease 2019 (COVID-19), transmitted person-to-person by the severe acute respiratory syndrome of coronavirus 2 (SARS-CoV-2), poses a threat to global public health.In this study, we performed the comprehensive analysis of the T cell receptor (TCR) repertoire may contribute to a more in-depth understanding of the pathogenesis of COVID-19.A comprehensive immunological analysis was performed to explore the features of the TCR repertoire and identified TCR sequences correlated with SARS-CoV-2 viral antigens.we analyzed the COVID-19 patients' TCR repertoires in peripheral blood mononuclear cells (PBMC) which obtained before (baseline), during (acute), and after rehabilitation (convalescent) by ImmunoSEQ-technology, and found that repertoire features of TCRβ-chain (TCRβ) complementary-determining region 3 (CDR3) in COVID-19 patients were remarkable difference, including decreased TCR diversity, abnormal CDR3 length, difference of TRBV/J gene usage and higher TCR sequence overlap. Besides, we identified some COVID-19 disease-associated TCRβ clones, and the abundance of them changed with the progression of the disease. Importantly, these disease-associated TCRβ clones could be used to distinguish COVID-19 patients from healthy controls with high accuracy.We provide a clear understanding of the TCR repertoire of COVID-19 patients, which lays the foundation for better diagnosis and treatment of COVID-19 patients.

Zinc finger C3H1 domain-containing protein (ZFC3H1) is differentially expressed between primary tumor and the normal in most cancers. Additionally, a recent study has suggested that ZFC3H1 could serve as a novel marker for the prognosis of prostate adenocarcinoma (PRAD). However, the relationship between ZFC3H1 expression and the prognostic values in most tumors remains unclear. Our study is mainly for exploring the prognosis of ZFC3H1 in pan-cancer and for further discovering a potential therapeutics target.Based on the clinical big data, we performed a pan-cancer analysis of ZFC3H1, including gene expression, survival prognosis, genetic alteration, protein phosphorylation, immune infiltration and enrichment analysis. In addition, Real-Time PCR and Western Blot were used to further confirm the role of ZFC3H1 in the colorectal cancer.We found that ZFC3H1 expression was connected with the prognosis of multiple malignant tumors. Furthermore, we also observed that ZFC3H1 was highly expressed in colorectal cancer through Real-Time PCR and Western Blot. The primary tumors presented higher phosphorylation level of the S655 site in lung adenocarcinoma, colon adenocarcinoma and uterine corpus endometrial carcinoma. ZFC3H1 expression was positively correlated with the immune infiltration of Cancer-associated fibroblasts (CAFs) in some tumors, such as liver hepatocellular carcinoma. And RNA surveillance pathways may be closely associated with the occurrence of tumors.Our study first reveals that ZFC3H1 could serve as a novel prognostic biomarker of pan-cancer, especially colorectal cancer.

Background Autosomal dominant osteopetrosis type II (ADO2) is a genetically and phenotypically metabolic bone disease, caused by osteoclast abnormalities. The pathways dysregulated in ADO2 could lead to the defects in osteoclast formation and function. However, the mechanism remains elusive. Materials and methods To systematically explore the molecular characterization of ADO2, we performed a multi-omics profiling from the autosomal dominant osteopetrosis type II iPSCs (ADO2-iPSCs) and healthy normal control iPSCs (NC-iPSCs) using whole genome re-sequencing, DNA methylation and N6-methyladenosine (m6A) analysis in this study. Results Totally, we detected 7,095,817 single nucleotide polymorphisms (SNPs) and 1,179,573 insertion and deletions (InDels), 1,001,943 differentially methylated regions (DMRs) and 2984 differential m6A peaks, and the comprehensive multi-omics profile was generated from the two cells. Interestingly, the ISG15 m6A level in ADO2-iPSCs is higher than NC-iPSCs by IGV software, and the differentially expressed m6A-modified genes (DEMGs) were highly enriched in the osteoclast differentiation and p53 signaling pathway, which associated with the development of osteopetrosis. In addition, combining our previously published transcriptome and proteome datasets, we found that the change in DNA methylation levels correlates inversely with some gene expression levels. Conclusion Our results indicate that the global multi-omics landscape not only provides a high-quality data resource but also reveals a dynamic pattern of gene expression, and found that the pathogenesis of ADO2 may begin early in life.

Background Autosomal dominant osteopetrosis type II (ADO2) is a rare human genetic disease that has been broadly studied as an important osteopetrosis model; however, there are no disease-specific induced pluripotent stem cells (ADO2-iPSCs) that may be valuable for understanding the pathogenesis and may be a potential source of cells for autologous cell-based therapies. Methods To generate the first human ADO2-iPSCs from a Chinese family with ADO2 and to identify their characteristics, blood samples were collected from the proband and his parents and were used for genotyping by whole-exome sequencing (WES); the urine-derived cells of the proband were reprogrammed with episomal plasmids that contained transcription factors, such as KLF4, OCT4, c-MYC, and SOX2. The proteome-wide protein quantification and lysine 2-hydroxyisobutyrylation detection of the ADO2-iPSCs and normal control iPSCs (NC-iPSCs) were performed by high-resolution LC-MS/MS and bioinformatics analysis. Results WES with filtering strategies identified a mutation in CLCN7 (R286W) in the proband and his father, which was absent in the proband’s mother and the healthy controls; this was confirmed by Sanger sequencing. The ADO2-iPSCs were successfully generated, which carried a normal male karyotype (46, XY) and the mutation of CLCN7 (R286W); the ADO2-iPSCs positively expressed alkaline phosphatase and other surface markers; and no vector and transgene were detected. The ADO2-iPSCs could differentiate into all three germ cell layers, both in vitro and in vivo. The proteomic profiling revealed similar expression of pluripotency markers in the two cell lines and identified 7405 proteins and 3664 2-hydroxyisobutyrylated peptides in 1036 proteins in the ADO2-iPSCs. Conclusions Our data indicated that the mutation CLCN7 (R286W) may be a cause of the osteopetrosis family. The generated vector-free and transgene-free ADO2-iPSCs with known proteomic characteristics may be valuable for personalized and cell-based regenerative medicine in the future. Electronic supplementary material The online version of this article (10.1186/s13287-019-1369-8) contains supplementary material, which is available to authorized users.

Liver cancer is the second leading cause of cancer mortality worldwide. Safer and more effective diagnostic methods for liver cancer are desirable, and biomarkers represent a potentially alternative method for diagnosis. The present study was designed to identify liver cancer biomarkers. We quantified the changes in serum protein levels between liver transplantation and healthy (control) females using isobaric tags for relative and absolute quantitation (iTRAQ) as well as proteomic analysis. A total of 1399 proteins were identified; of these, three proteins showed significantly different concentrations between the before transplantation group and the control group. These proteins may thus be relevant to liver cancer and constitute potential liver cancer biomarkers.

BACKGROUND There are many situations of abnormal metabolism influencing liver graft function. This study aims to provide data for the development of liver function recovery after liver transplantation by dynamically analyzing metabolites of bile acids pathway in serum. MATERIAL AND METHODS A comprehensive metabolomics profiling of serum of 9 liver transplantation patients before transplantation, on the 1st, 3rd, and 7th days after liver transplantation, and healthy individuals were performed by ultra-performance liquid chromatography-mass spectrometry (UPLC-MS). Multivariate data and dynamic analysis were used to search for biomarkers between the metabolomics profiles present in perioperative liver transplantation and normal controls. RESULTS Thirty-three differential endogenous metabolites were screened by the threshold of variable importance in the projection (VIP) from an orthogonal partial least square discriminant analysis (OPLS-DA) greater than 1.0, q-value

Background Dysregulation of lipid metabolism has been implicated in the progression of hepatocellular carcinoma (HCC). We therefore investigated the molecular characteristics of lipid-metabolism-related genes for the prognostic prediction of HCC. Methods Multi-dimensional bioinformatics analysis was conducted to comprehensively analyzed the lipid metabolism-related genes (IMRG) and construct the prognostic prediction signature. Results A total of 770 HCC patients and their corresponding 776 IMRGs were downloaded from three databases. The HCC patients were classified into 2 molecular clusters, which were associated with overall survival, clinical characteristics, and immune cells. The biological function of the differentially expressed IMRGs in the 2 clusters showed that the genes were associated with tumor-related metabolism pathways. A 6 IMRGs signature (6-IS), including FMO3, SLC11A1, RNF10, KCNH2, ME1, and ZIC2 were established for HCC prognostic prediction, which was found to be an independent prognostic factor. Performance of the 6-IS prognostic signature was verified in a validation set and compared with an external data set. Results revealed that the 6-IS signature could effectively predict the prognosis of patients with HCC. Conclusion This study provides new insights into the role of IMRG in the pathogenesis of HCC and presents a novel signature 6-IS to predict the prognosis of HCC.

Embryonic development is a complex process that is regulated by a series of precise cellular behaviours. The limited number of cells in the early stages of embryonic development represents a challenge for studying early gene regulation and maintaining cell sternness. Single-cell sequencing is a new technology for high-throughput sequencing analysis at the single-cell level that not only reflects the heterogeneity between cells but also reveals gene expression characteristics in different cells from limited samples. Currently, the widespread application of single-cell RNA sequencing technology is gradually changing our understanding of disease pathogenesis. This article reviews the application of single-cell RNA sequencing in embryonic development in recent years and provides innovative ideas for research on embryonic development and the treatment of diseases related to embryonic development.

This research aims to investigate the immune-associated gene signature from databases to improve the prognostic value in hepatocellular carcinoma (HCC) by multidimensional methods using various bioinformatic methods. Fifty-one immune-associated genes were mined out, which were associated with clinical characters through univariate and multivariate Cox regression analyses, and 51 immune-associated genes could be well-divided HCC samples into high-risk and low-risk clusters. Next, we performed least absolute shrinkage and selection operator (LASSO) Cox regression method to reveal 18 immune-associated genes' signature and calculate risk score of each gene for receiver operating characteristic (ROC) analysis. Comparing with low-risk cluster, high-risk cluster had higher risk score with unfavorable prognosis. Then, multivariate Cox regression analysis showed that risk score of 18 immune-associated genes' signature was associated with tumor invasion and tumor-node-metastasis (TNM) stage. ROC analysis indicated combined TNM stage, and risk score performed more sensitive and specific than single TNM stage or risk score in survival prediction. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the pathways enriched in tumorigenesis were related to risk score, and those pathways could separate HCC samples into high and low clusters. In addition, the survival prediction of 18 immune-associated genes' signature was well validated in independent test data set, external data set, and Real-time Quantitative PCR (RT-qPCR) experiment. The 18 immune-associated genes' signature was constructed, which could be used in effective prediction of HCC prognosis.

Liver cirrhosis of hepatitis B is an immune-related disease in which liver cells die during the body’s immune system activation to clear the virus, and the progress is closely related to T lymphocytes. T lymphocyte cells recognise antigens, specifically by major histocompatibility complex (MHC), through a membrane protein T cell receptor (TCR). Here, we used high throughput immune repertoire sequencing technique to study the characteristics and diversity of the TCR repertoire between patients who underwent liver transplantation and healthy controls (NC). We sequenced the TCR β-chain complementary-determining region 3 (CDR3) repertoire in peripheral blood mononuclear cells (PBMCs) from 6 liver transplantation patients before transplantation (Pre) and on the first (Post1) and seventh days (Post7) after transplantation along with 6 NC. We observed that the distributions of CDR3, VD indel, and DJ indel lengths were similar among the Pre, Post1, Post7 and NC groups. We found that the TCR repertoire diversity of transplantation groups was relatively lower compared to NC group. The Pre-group had more highly expanded T cell clones compared to Post1, Post7 and NC groups, and the diversity of the T cell repertoire of the Post7 group was significantly decreased compared to the Pre, Post1 and NC groups. In addition, we found our results also show that various TRBV expression increased and some public sequences at different time points after liver transplantation, and the expression levels of 3 TRBV segments and 2 TRBJ segments were also significantly different in Pre, Post1, Post7 and NC groups. Moreover, 1 aa sequence shared by all liver transplantation patients and 2 aa sequences shared by at least two groups, which may serve as biomarkers to monitor the immune status of liver transplant patients.

Recent studies have shown that circular RNAs (circRNAs) exhibit differential expression in certain diseases. However, to the best of our knowledge, maternal fetal-derived circRNAs and mRNAs associated with Down's syndrome (DS) have not yet been investigated. A total of 12 umbilical cord blood samples were collected from pregnant women, including six women carrying fetuses with DS (diagnosed by G-banding karyotype analysis), and six women carrying fetuses without DS. In addition, 12 peripheral blood samples were obtained from children, including six children with DS and six children without DS. Gene chip technology was used to screen for differentially expressed circRNAs and mRNAs in the cord blood samples, and were subsequently verified by reverse transcription-quantitative polymerase chain reaction in peripheral blood from the children to identify potential biomarkers. Furthermore, circRNA/microRNA (miRNA) interactions were predicted using Arraystar miRNA target prediction software. There was a significant difference in the expression of hsa_circRNA_103127, hsa_circRNA_103112 and hsa_circRNA_104907 between cord blood obtained from the women carrying fetuses with and without DS, and between peripheral blood obtained from children with and without DS (P

Thrombophilia is a multifactorial disorder involving environmental and genetic factors. Well-known factors that result in predisposition to congenital disorders associated with thrombophilia include antithrombin deficiency, protein C and S deficiency, Factor V Leiden mutation, abnormal prothrombin and antiphospholipid syndrome. The present study revealed an association between a mutation of the F2 gene, which codes for coagulation factor II, thrombin, and the risk of thrombophilia in a Han Chinese family, of which four members (I-2, II-2, II-3 and III-1) had a history of deep venous thromboembolism. The disease was measured in this family using laboratory measurements and computed tomography angiography. To identify the abnormality underlying the increased thrombophilia risk, whole-exome sequencing technology was used to analyze two affected individuals (II-2 and III-1). An exonic missense F2 mutation, T165M ({"type":"entrez-nucleotide","attrs":{"text":"NM_000506","term_id":"1732746292","term_text":"NM_000506"}}NM_000506:c.C494T:p.T165M;rs5896), was identified from a total of 2,222 and 2,203 genetic variations observed in the two affected individuals, respectively, which were subsequently filtered and confirmed using Sanger sequencing. I-2, II-3 and III-1 shared this mutation with the proband (II-2), and II-6 had a heterozygous form of the mutation. This deleterious mutation was not identified in normal controls. The present study may improve understanding of the function of the F2 gene.

The aim of the present study was to screen gastric cancer (GC) tissue and adjacent tissue for differences in mRNA and circular (circRNA) expression, to analyze the differences in circRNA and mRNA expression, and to investigate the circRNA expression in gastric carcinoma and its mechanism. circRNA and mRNA differential expression profiles generated using Agilent microarray technology were analyzed in the GC tissues and adjacent tissues. qRT-PCR was used to verify the differential expression of circRNAs and mRNAs according to the interactions between circRNAs and miRNAs as well as the possible existence of miRNA and mRNA interactions. We found that: i) the circRNA expression profile revealed 1,285 significant differences in circRNA expression, with circRNA expression downregulated in 594 samples and upregulated in 691 samples via interactions with miRNAs. The qRT-PCR validation experiments showed that hsa_circRNA_400071, hsa_circRNA_000543 and hsa_circRNA_001959 expression was consistent with the microarray analysis results. ii) 29,112 genes were found in the GC tissues and adjacent tissues, including 5,460 differentially expressed genes. Among them, 2,390 differentially expressed genes were upregulated and 3,070 genes were downregulated. Gene Ontology (GO) analysis of the differentially expressed genes revealed these genes involved in biological process classification, cellular component classification and molecular function classification. Pathway analysis of the differentially expressed genes identified 83 significantly enriched genes, including 28 upregulated genes and 55 downregulated genes. iii) 69 differentially expressed circRNAs were found that might adsorb specific miRNAs to regulate the expression of their target gene mRNAs. The conclusions are: i) differentially expressed circRNAs had corresponding miRNA binding sites. These circRNAs regulated the expression of target genes through interactions with miRNAs and might become new molecular biomarkers for GC in the future. ii) Differentially expressed genes may be involved in the occurrence of GC via a variety of mechanisms. iii) CD44, CXXC5, MYH9, MALAT1 and other genes may have important implications for the occurrence and development of GC through the regulation, interaction, and mutual influence of circRNA-miRNA-mRNA via different mechanisms.

Klinefelter syndrome (KS) is one of the most common congenital disorders of male infertility. Given its high heterogeneity in clinical and genetic presentation, the relationship between transcriptome, clinical phenotype, and associated co-morbidities seen in KS has not been fully clarified. Here, we report a 47,XXY Chinese male with infertility and analyzed the differences in gene expression patterns of peripheral blood mononuclear cells (PBMCs) with regard to a Chinese male and a female control with normal karyotype by single-cell sequencing. A total of 24,439 cells were analyzed and divided into 5 immune cell types (including B cells, T cells, macrophage cells, dendritic cells, and natural killer cells) according to marker genes. Using unsupervised dimensionality reduction and clustering algorithms, we identified molecularly distinct subpopulations of cells between the KS patient and both controls. Gene ontology enrichment analyses yielded terms associated with well-known comorbidities seen in KS as well as an affected immune system and type I diabetes mellitus. Based on our data, we identified several candidate genes which may be implicated in regulating the phenotype of KS. Overall, this analysis provides a comprehensive map of the cell types of PBMCs in a KS patient at the single-cell level, which will contribute to the prevention of comorbidity and improvement of the life quality of KS patients.

Hepatocellular carcinoma (HCC) is one of the most lethal malignancies that is closely associated with the Hepatitis B virus (HBV). HBV integration into host genomes can induce instability and the aberrant expression of human genomic DNA. To directly assess HBV integration breakpoints at whole genome level, four small sequencing libraries were constructed and the HBV integration profiles of four patients with HCC were characterized. In total, the current study identified 11,800,974, 11,216,998, 11,026,546 and 11,607,842 clean reads for patients 1-3 and 4, respectively, of which 92.82, 95.95, 97.21 and 97.29% were properly aligned to the hybrid reference genome. In addition, 220 HBV integration events were detected from the tumor tissues of four patients with HCC and an average of 55 breakpoints per sample was calculated. The results indicated that HBV integration events may be implicated in HCC physiologies and diseases. The results acquired may also provide insight into the pathogenesis of HCC, which may be valuable for future HCC therapy.

Bladder cancer (BCa) is a tumor associated with high morbidity and mortality and its incidence is increasing worldwide. However, the pathogenesis of bladder cancer is not well understood. To further illustrate the molecular mechanisms involved in the pathogenesis of BCa and identify potential therapeutic targets, we combined the transcriptomic analysis with RNA sequencing and tandem mass tags (TMT)-based proteomic methods to quantitatively screen the differentially expressed genes and proteins between bladder cancer tissues (BC) and adjacent normal tissues (AN). Transcriptome and proteome studies indicated 7094 differentially expressed genes (DEGs) and 596 differentially expressed proteins (DEPs) between BC and AN, respectively. GO enrichment analyses revealed that cell adhesion, calcium ion transport, and regulation of ATPase activity were highly enriched in BCa. Moreover, several key signaling pathway were identified as of relevance to BCa, in particular the ECM-receptor interaction, cell adhesion molecules (CAMs), and PPAR signaling pathway. Interestingly, 367 genes were shared by DEGs and DEPs, and a significant positive correlation between mRNA and translation profiles was found. In summary, this joint analysis of transcript and protein profiles provides a comprehensive reference map of gene activity regarding the disease status of BCa.

The present study aimed to investigate the complexity and diversity of the T lymphocyte immune repertoire in patients with bladder cancer. To do so, the immune state of patients was assessed. The study also aimed to elucidate the aetiology and pathogenesis of bladder cancer to provide a novel theoretical basis for disease prevention, diagnosis, treatment and prognosis monitoring. Cancerous and paracancerous (control) tissue samples were collected from five patients diagnosed with muscle-invasive bladder cancer. Multiplex polymerase chain reaction and Illumina high-throughput sequencing were used to determine the characteristics and clonal diversity of the T-cell receptor (TCR) β-chain complementarity-determining region 3 (CDR3) gene in the cancerous and paracancerous tissues of patients with bladder cancer. The degree of clonal expansion in malignant samples was significantly higher than in adjacent samples. Furthermore, ΤCRβ variable (TRBV), ΤCRβ diversity (TRBD) and ΤCRβ joining (TRBJ) repertoires were significantly different in cancerous samples compared with adjacent samples. In addition, 13 identified V-J pairs were highly expressed in cancerous samples whereas they had low expression in control samples. In conclusion, the degree of T-cell clonal expansion in bladder cancerous tissue was higher than in paracancerous tissue, whereas the immune diversity of the tissues of patients with bladder cancer was significantly lower. The DNA sequence and amino acid sequences, and V-J combination level may be used to comprehensively understand the diversity and characteristics of TCR CDR3 in bladder cancer and paracancerous tissues, and to evaluate the immune status of bladder cancer to develop therapeutic targets and biomarkers for prognosis monitoring.

The forkhead box A1 (FOXA1), one of the forkhead class of DNA-binding proteins, functions as a transcription factor and plays a vital role in cellular control of embryonic development and cancer progression. Downregulation of FOXA1 has reported in several types of cancer, which contributes to cancer cell survival and chemoresistance. However, the mechanism for FOXA1 downregulation in cancer remains unclear. Here, we report that the ubiquitination enzyme zinc finger protein 91 (ZFP91) ubiquitinates and destabilizes FOXA1, which promotes cancer cell growth. High level of ZFP91 expression correlates with low level of FOXA1 protein in human gastric cancer (GC) cell lines and patient samples. Furthermore, ZFP91 knockdown reduces FOXA1 polyubiquitination, which decreases FOXA1 turnover and enhances cellular sensitivity to chemotherapy. Taken together, our findings reveal ZFP91-FOXA1 axis plays an important role in promoting GC progression and provides us a potential therapeutic intervention in the treatment of GC.

Full proteomic analysis methods. (DOCX 21 kb)

With the advent of next-generation sequencing technology, the cost of sequencing has significantly decreased. However, sequencing costs remain high for large-scale studies. In the present study, DNA pooling was applied as a cost-effective strategy for sequencing. The sequencing results for 100 healthy individuals obtained via whole-genome resequencing and using DNA pooling are presented in the present study. In order to minimise the likelihood of systematic bias in sampling, paired-end libraries with an insert size of 500 bp were prepared for all samples and then subjected to whole-genome sequencing using four lanes for each library and resulting in at least a 30-fold haploid coverage for each sample. The NCBI human genome build37 (hg19) was used as a reference genome for the present study and the short reads were aligned to the reference genome achieving 99.84% coverage. In addition, the average sequencing depth was 32.76. In total, ~3 million single-nucleotide polymorphisms were identified, of which 99.88% were in the NCBI dbSNP database. Furthermore, ~600,000 small insertion/deletions, 500,000 structure variants, 5,000 copy number variations and 13,000 single nucleotide variants were identified. According to the present study, the whole genome has been sequenced for a small sample subjects from southern China for the first time. Furthermore, new variation sites were identified by comparing with the reference sequence, and new knowledge of the human genome variation was added to the human genomic databases. Furthermore, the particular distribution regions of variation were illustrated by analyzing various sites of variation, such as single-nucleotide polymorphisms.

Down syndrome (DS) is caused by trisomy of human chromosome 21 and is associated with a number of deleterious phenotypes. To investigate the role of microRNA (miRNA) in the regulation of DS, high-throughput Illumina sequencing technology and isobaric tagging for relative and absolute protein quantification analysis were utilized for simultaneous expression profiling of miRNA and protein in fetuses with DS and normal fetuses. A total of 344 miRNAs were associated with DS. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were used to investigate the proteins found to be differentially expressed. Functionally important miRNAs were determined by identifying enriched or depleted targets in the transcript and the protein expression levels were consistent with miRNA regulation. The results indicated that GRB2, TMSB10, RUVBL2, the hsa-miR-329 and hsa-miR-27b, hsa-miR-27a targets, and MAPK1, PTPN11, ACTA2 and PTK2 or other differentially expressed proteins were connected with each other directly or indirectly. Integrative analysis of miRNAs and proteins provided an expansive view of the molecular signaling pathways in DS.

There is increasing evidence that several genes are associated with an increased risk of type 2 diabetes (T2D); genome-wide association investigations and whole-genome re-sequencing investigations offer a useful approach for the identification of genes involved in common human diseases. To further investigate which polymorphisms confer susceptibility to T2D, the present study screened for high-contribution susceptibility gene variants Chinese patients with T2D using whole-genome re-sequencing with DNA pooling. In total, 100 Chinese individuals with T2D and 100 healthy Chinese individuals were analyzed using whole-genome re-sequencing using DNA pooling. To minimize the likelihood of systematic bias in sampling, paired-end libraries with an insert size of 500 bp were prepared for in T2D in all samples, which were then subjected to whole-genome sequencing. Each library contained four lanes. The average sequencing depth was 35.70. In the present study, 1.36 GB of clean sequence data were generated, and the resulting calculated T2D genome consensus sequence covered 99.88% of the hg19 sequence. A total of 3,974,307 single nucleotide polymorphisms were identified, of which 99.88% were in the dbSNP database. The present study also found 642,189 insertions and deletions, 5,590 structure variants (SVs), 4,713 copy number variants (CNVs) and 13,049 single nucleotide variants. A total of 1,884 somatic CNVs and 74 somatic SVs were significantly different between the cases and controls. Therefore, the present study provided validation of whole-genome re-sequencing using the DNA pooling approach. It also generated a whole-genome re-sequencing genotype database for future investigations of T2D.

Objective: The objective of this study was to measure the 3′-untranslated region (3′-UTR) polymorphism lengths in peripheral blood mononuclear cells (PBMCs) from uremia patients.Method: We sequenced the alternative polyadenylation sites in the 3′-UTR of PBMCs from 10 uremic patients and 10 healthy people to detect different gene expression levels between uremia patients and healthy controls. Quantitative reverse transcription polymerase chain reaction was used as validation. Result: Compared with the healthy control group, 691 genes in uremic patients had significantly different 3′-UTR lengths. Of these genes, 475 genes showed shortened 3′-UTRs, and the 3′-UTRs of 216 genes were lengthened. The verification results matched the original sequencing results. Conclusion: There were significant differences in 3′-UTR lengths between uremic patients and healthy controls, and analysis of the differential genes may contribute to the understanding of uremia pathogenesis.

Proteus syndrome is a rare complex syndrome involving clinical presentation with atypical skeletal growth. Only a limited number of cases with this syndrome have been reported in the literature to date. We herein report the case of a Chinese male patient with Proteus syndrome and review the clinical and molecular characteristics of this disease. The patient was a 34-year-old man with clinical manifestations suggestive of the Proteus syndrome, including mosaic distribution of the lesions, sporadic occurrence, progressive course, disproportionate overgrowth of the legs, epidermal nevi, lipomas, venous malformations and characteristic facial phenotype. Genetic mosaicism, such as mutations involving the phosphoinositide 3 kinase-AKT pathway in the affected tissues, may be important causes of Proteus syndrome. In the present case, samples from the affected tissues were collected from the patient and were further analyzed using whole-exome sequencing. However, no mutation of the genes reportedly associated with Proteus syndrome was identified in the affected tissues. Proteus syndrome is a complex mosaic disorder with a number of variable characteristics. Although activating AKT1 mutations have been found to be associated with this disorder, the molecular etiology remains to be fully elucidated and diagnostic criteria must be established in the clinical setting.

Liver cirrhosis of hepatitis B is an immune-related disease in which liver cells die during the body's immune system activation to clear the virus, and the progress is closely related to T lymphocytes. T lymphocyte cells recognise antigens, specifically by major histocompatibility complex (MHC), through a membrane protein T cell receptor (TCR). Here, we used high throughput immune repertoire sequencing technique to study the characteristics and diversity of the TCR repertoire between patients who underwent liver transplantation and healthy controls (NC). We sequenced the TCR β-chain complementary-determining region 3 (CDR3) repertoire in peripheral blood mononuclear cells (PBMCs) from 6 liver transplantation patients before transplantation (Pre) and on the first (Post1) and seventh days (Post7) after transplantation along with 6 NC. We observed that the distributions of CDR3, VD indel, and DJ indel lengths were similar among the Pre, Post1, Post7 and NC groups. We found that the TCR repertoire diversity of transplantation groups was relatively lower compared to NC group. The Pre-group had more highly expanded T cell clones compared to Post1, Post7 and NC groups, and the diversity of the T cell repertoire of the Post7 group was significantly decreased compared to the Pre, Post1 and NC groups. In addition, we found our results also show that various TRBV expression increased and some public sequences at different time points after liver transplantation, and the expression levels of 3 TRBV segments and 2 TRBJ segments were also significantly different in Pre, Post1, Post7 and NC groups. Moreover, 1 aa sequence shared by all liver transplantation patients and 2 aa sequences shared by at least two groups, which may serve as biomarkers to monitor the immune status of liver transplant patients.