Genomics, Proteomics & Bioinformatics最新文献

The high host genetic background of tissue biopsies hinders the application of shotgun metagenomic sequencing in characterizing the tissue microbiota. We proposed an optimized method that removed host DNA from colon biopsies and examined the effect on metagenomic analysis. Human or mouse colon biopsies were divided into two groups, with one group undergoing host DNA depletion and the other serving as the control. Host DNA was removed through differential lysis of mammalian and bacterial cells before sequencing. The impact of host DNA depletion on microbiota was compared based on phylogenetic diversity analyses and regression analyses. Removing host DNA enhanced bacterial sequencing depth and improved species discovery, increasing bacterial reads by 2.46 ± 0.20 folds while reducing host reads by 6.80% ± 1.06%. Moreover, 2.40 times more of bacterial species were detected after host DNA depletion. This was confirmed from mouse colon tissues, increasing bacterial reads by 5.46 ± 0.42 folds while decreasing host reads by 10.2% ± 0.83%. Similarly, significantly more bacterial species were detected in the mouse colon tissue upon host DNA depletion (P < 0.001). Furthermore, an increased microbial richness was evident in the host DNA-depleted samples compared with non-depleted controls in human colon biopsies and mouse colon tissues (P < 0.001). Our optimized method of host DNA depletion improves the sensitivity of shotgun metagenomic sequencing in bacteria detection in the biopsy, which may yield a more accurate taxonomic profile of the tissue microbiota and identify bacteria that are important for disease initiation or progression.

In perinatal medicine, intrauterine growth restriction (IUGR) is one of the greatest challenges. The etiology of IUGR is multifactorial, but most cases are thought to arise from placental insufficiency. However, identifying the placental cause of IUGR can be difficult due to numerous confounding factors. Selective IUGR (sIUGR) would be a good model to investigate how impaired placentation affects fetal development, as the growth discordance between monochorionic twins cannot be explained by confounding genetic or maternal factors. Herein, we constructed and analyzed the placental proteomic profiles of IUGR twins and normal cotwins. Specifically, we identified a total of 5481 proteins, of which 233 were differentially expressed (57 up-regulated and 176 down-regulated) in IUGR twins. Bioinformatics analysis indicates that these differentially expressed proteins (DEPs) are mainly associated with cardiovascular system development and function, organismal survival, and organismal development. Notably, 34 DEPs are significantly enriched in angiogenesis, and diminished placental angiogenesis in IUGR twins has been further elaborately confirmed. Moreover, we found decreased expression of metadherin (MTDH) in the placentas of IUGR twins and demonstrated that MTDH contributes to placental angiogenesis and fetal growth in vitro. Collectively, our findings reveal the comprehensive proteomic signatures of placentas for sIUGR twins, and the DEPs identified may provide in-depth insights into the pathogenesis of placental dysfunction and subsequent impaired fetal growth.

The fetal liver (FL) is the key erythropoietic organ during fetal development, but knowledge on human FL erythropoiesis is very limited. In this study, we sorted primary erythroblasts from FL cells and performed RNA sequencing (RNA-seq) analyses. We found that temporal gene expression patterns reflected changes in function during primary human FL terminal erythropoiesis. Notably, the expression of genes enriched in proteolysis and autophagy was up-regulated in orthochromatic erythroblasts (OrthoEs), suggesting the involvement of these pathways in enucleation. We also performed RNA-seq of in vitro cultured erythroblasts derived from FL CD34⁺ cells. Comparison of transcriptomes between the primary and cultured erythroblasts revealed significant differences, indicating impacts of the culture system on gene expression. Notably, the expression of lipid metabolism-related genes was increased in cultured erythroblasts. We further immortalized erythroid cell lines from FL and cord blood (CB) CD34⁺ cells (FL-iEry and CB-iEry, respectively). FL-iEry and CB-iEry were immortalized at the proerythroblast stage and can be induced to differentiate into OrthoEs, but their enucleation ability was very low. Comparison of the transcriptomes between OrthoEs with and without enucleation capability revealed the down-regulation of pathways involved in chromatin organization and mitophagy in OrthoEs without enucleation capacity, indicating that defects in chromatin organization and mitophagy contribute to the inability of OrthoEs to enucleate. Additionally, the expression of HBE1, HBZ, and HBG2 was up-regulated in FL-iEry compared with CB-iEry, and such up-regulation was accompanied by down-regulated expression of BCL11A and up-regulated expression of LIN28B and IGF2BP1. Our study provides new insights into human FL erythropoiesis and rich resources for future studies.

A series of clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR associated protein 9 (Cas9) systems have been engineered for genome editing. The most widely used Cas9 is SpCas9 from Streptococcus pyogenes and SaCas9 from Staphylococcus aureus. However, a comparison of their detailed gene editing outcomes is still lacking. By characterizing the editing outcomes of 11 sites in human induced pluripotent stem cells (iPSCs) and K562 cells, we found that SaCas9 could edit the genome with greater efficiencies than SpCas9. We also compared the effects of spacer lengths of single-guide RNAs (sgRNAs; 18-21 nt for SpCas9 and 19-23 nt for SaCas9) and found that the optimal spacer lengths were 20 nt and 21 nt for SpCas9 and SaCas9, respectively. However, the optimal spacer length for a particular sgRNA was 18-21 nt for SpCas9 and 21-22 nt for SaCas9. Furthermore, SpCas9 exhibited a more substantial bias than SaCas9 for nonhomologous end-joining (NHEJ) +1 insertion at the fourth nucleotide upstream of the protospacer adjacent motif (PAM), indicating a characteristic of a staggered cut. Accordingly, editing with SaCas9 led to higher efficiencies of NHEJ-mediated double-stranded oligodeoxynucleotide (dsODN) insertion or homology-directed repair (HDR)-mediated adeno-associated virus serotype 6 (AAV6) donor knock-in. Finally, GUIDE-seq analysis revealed that SaCas9 exhibited significantly reduced off-target effects compared with SpCas9. Our work indicates the superior performance of SaCas9 to SpCas9 in transgene integration-based therapeutic gene editing and the necessity to identify the optimal spacer length to achieve desired editing results.

The recently developed technologies that allow the analysis of each single omics have provided an unbiased insight into ongoing disease processes. However, it remains challenging to specify the study design for the subsequent integration strategies that can associate sepsis pathophysiology and clinical outcomes. Here, we conducted a time-dependent multi-omics integration (TDMI) in a sepsis-associated liver dysfunction (SALD) model. We successfully deduced the relation of the Toll-like receptor 4 (TLR4) pathway with SALD. Although TLR4 is a critical factor in sepsis progression, it is not specified in single-omics analyses but only in the TDMI analysis. This finding indicates that the TDMI-based approach is more advantageous than single-omics analyses in terms of exploring the underlying pathophysiological mechanism of SALD. Furthermore, TDMI-based approach can be an ideal paradigm for insightful biological interpretations of multi-omics datasets that will potentially reveal novel insights into basic biology, health, and diseases, thus allowing the identification of promising candidates for therapeutic strategies.

The biology of Toxoplasma gondii, the causative pathogen of one of the most widespread parasitic diseases (toxoplasmosis), remains poorly understood. Lactate, which is derived from glucose metabolism, is not only an energy source in a variety of organisms, including T. gondii, but also a regulatory molecule that participates in gene activation and protein function. Lysine lactylation (Kla) is a type of post-translational modifications (PTMs) that has been recently associated with chromatin remodeling; however, Kla of histone and non-histone proteins has not yet been studied in T. gondii. To examine the prevalence and function of lactylation in T. gondii parasites, we mapped the lactylome of proliferating tachyzoite cells and identified 1964 Kla sites on 955 proteins in the T. gondii RH strain. Lactylated proteins were distributed in multiple subcellular compartments and were closely related to a wide variety of biological processes, including mRNA splicing, glycolysis, aminoacyl-tRNA biosynthesis, RNA transport, and many signaling pathways. We also performed a chromatin immunoprecipitation sequencing (ChIP-seq) analysis using a lactylation-specific antibody and found that the histones H4K12la and H3K14la were enriched in the promoter and exon regions of T. gondii associated with microtubule-based movement and cell invasion. We further confirmed the delactylase activity of histone deacetylases TgHDAC2-4, and found that treatment with anti-histone acetyltransferase (TgMYST-A) antibodies profoundly reduced protein lactylation in T. gondii. This study offers the first dataset of the global lactylation proteome and provides a basis for further dissecting the functional biology of T. gondii.

To achieve the enormous potential of gene-editing technology in clinical therapies, one needs to evaluate both the on-target efficiency and unintended editing consequences comprehensively. However, there is a lack of a pipelined, large-scale, and economical workflow for detecting genome editing outcomes, in particular insertion or deletion of a large fragment. Here, we describe an approach for efficient and accurate detection of multiple genetic changes after CRISPR/Cas9 editing by pooled nanopore sequencing of barcoded long-range PCR products. Recognizing the high error rates of Oxford nanopore sequencing, we developed a novel pipeline to capture the barcoded sequences by grepping reads of nanopore amplicon sequencing (GREPore-seq). GREPore-seq can assess nonhomologous end-joining (NHEJ)-mediated double-stranded oligodeoxynucleotide (dsODN) insertions with comparable accuracy to Illumina next-generation sequencing (NGS). GREPore-seq also reveals a full spectrum of homology-directed repair (HDR)-mediated large gene knock-in, correlating well with the fluorescence-activated cell sorting (FACS) analysis results. Of note, we discovered low-level fragmented and full-length plasmid backbone insertion at the CRISPR cutting site. Therefore, we have established a practical workflow to evaluate various genetic changes, including quantifying insertions of short dsODNs, knock-ins of long pieces, plasmid insertions, and large fragment deletions after CRISPR/Cas9-mediated editing. GREPore-seq is freely available at GitHub (https://github.com/lisiang/GREPore-seq) and the National Genomics Data Center (NGDC) BioCode (https://ngdc.cncb.ac.cn/biocode/tools/BT007293).

Acid-base homeostasis is a fundamental property of living cells, and its persistent disruption in human cells can lead to a wide range of diseases. In this study, we conducted a computational modeling analysis of transcriptomic data of 4750 human tissue samples of 9 cancer types in The Cancer Genome Atlas (TCGA) database. Built on our previous study, we quantitatively estimated the average production rate of OH^- by cytosolic Fenton reactions, which continuously disrupt the intracellular pH (pH_i) homeostasis. Our predictions indicate that all or at least a subset of 43 reprogrammed metabolisms (RMs) are induced to produce net protons (H⁺) at comparable rates of Fenton reactions to keep the pH_i stable. We then discovered that a number of well-known phenotypes of cancers, including increased growth rate, metastasis rate, and local immune cell composition, can be naturally explained in terms of the Fenton reaction level and the induced RMs. This study strongly suggests the possibility to have a unified framework for studies of cancer-inducing stressors, adaptive metabolic reprogramming, and cancerous behaviors. In addition, strong evidence is provided to demonstrate that a popular view that Na⁺/H⁺ exchangers along with lactic acid exporters and carbonic anhydrases are responsible for the intracellular alkalization and extracellular acidification in cancer may not be justified.

Streptomycetes possess numerous gene clusters and the potential to produce a large amount of natural products. Histone deacetylase (HDAC) inhibitors play an important role in the regulation of histone modifications in fungi, but their roles in prokaryotes remain poorly understood. Here, we investigated the global effects of the HDAC inhibitor, sodium butyrate (SB), on marine-derived Streptomycesolivaceus FXJ 8.021, particularly focusing on the activation of secondary metabolite biosynthesis. The antiSMASH analysis revealed 33 secondary metabolite biosynthetic gene clusters (BGCs) in strain FXJ 8.021, among which the silent lobophorin BGC was activated by SB. Transcriptomic data showed that the expression of genes involved in lobophorin biosynthesis (ge00097-ge00139) and CoA-ester formation (e.g., ge02824), as well as the glycolysis/gluconeogenesis pathway (e.g., ge01661), was significantly up-regulated in the presence of SB. Intracellular CoA-ester analysis confirmed that SB triggered the biosynthesis of CoA-ester, thereby increasing the precursor supply for lobophorin biosynthesis. Further acetylomic analysis revealed that the acetylation levels on 218 sites of 190 proteins were up-regulated and those on 411 sites of 310 proteins were down-regulated. These acetylated proteins were particularly enriched in transcriptional and translational machinery components (e.g., elongation factor GE04399), and their correlations with the proteins involved in lobophorin biosynthesis were established by protein-protein interaction network analysis, suggesting that SB might function via a complex hierarchical regulation to activate the expression of lobophorin BGC. These findings provide solid evidence that acetylated proteins triggered by SB could affect the expression of genes involved in the biosynthesis of primary and secondary metabolites in prokaryotes.