ggCyto: next generation open-source visualization software for cytometry

Optimization of a Cas12a-Driven Synthetic Gene Regulatory Network System

Gene regulatory networks, which control gene expression patterns in development and in response to stimuli, use regulatory logic modules to coordinate inputs and outputs. One example of a regulatory logic module is the gene regulatory cascade (GRC), where a series of transcription factor genes turn on in order. Synthetic biologists have derived artificial systems that encode regulatory rules, including GRCs. Furthermore, the development of single-cell approaches has enabled the discovery of gene regulatory modules in a variety of experimental settings. However, the tools available for validating these observations remain limited. Based on a synthetic GRC using DNA cutting-defective Cas9 (dCas9), we designed and implemented an alternative synthetic GRC utilizing DNA cutting-defective Cas12a (dCas12a). Comparing the ability of these two systems to express a fluorescent reporter, the dCas9 system was initially more active, while the dCas12a system was more streamlined. Investigating the influence of individual components of the systems identified nuclear localization as a major driver of differences in activity. Improving nuclear localization for the dCas12a system resulted in 1.5-fold more reporter-positive cells and a 15-fold increase in reporter intensity relative to the dCas9 system. We call this optimized system the "Synthetic Gene Regulatory Network" (SGRN, pronounced "sojourn").

MftG is crucial for ethanol metabolism of mycobacteria by linking mycofactocin oxidation to respiration

Mycofactocin is a redox cofactor essential for the alcohol metabolism of mycobacteria. While the biosynthesis of mycofactocin is well established, the gene mftG, which encodes an oxidoreductase of the glucose-methanol-choline superfamily, remained functionally uncharacterized. Here, we show that MftG enzymes are almost exclusively found in genomes containing mycofactocin biosynthetic genes and are present in 75% of organisms harboring these genes. Gene deletion experiments in Mycolicibacterium smegmatis demonstrated a growth defect of the ∆mftG mutant on ethanol as a carbon source, accompanied by an arrest of cell division reminiscent of mild starvation. Investigation of carbon and cofactor metabolism implied a defect in mycofactocin reoxidation. Cell-free enzyme assays and respirometry using isolated cell membranes indicated that MftG acts as a mycofactocin dehydrogenase shuttling electrons toward the respiratory chain. Transcriptomics studies also indicated remodeling of redox metabolism to compensate for a shortage of redox equivalents. In conclusion, this work closes an important knowledge gap concerning the mycofactocin system and adds a new pathway to the intricate web of redox reactions governing the metabolism of mycobacteria.

Doramectin Induces Apoptosis in B16 Melanoma Cells

INTRODUCTION/OBJECTIVE

Metastatic melanoma resists current pharmacological regimens that act through apoptosis. This indicates that therapies acting via non-apoptotic cell-death pathways could be pursued. Doramectin has shown promising results in another cancer of neural crest origin, neuroblastoma, through the inhibition of growth via autophagy. Our research hypothesis is that doramectin induces autophagy in B16F10 melanoma cells.

METHODS

Cells were treated with doramectin (15 uM) or a combination of both doramectin and a cell-death inhibitor, compared to untreated control cells (media), and then analyzed with MTT analysis. Likewise, MDC analysis was completed to detect autophagy involvement with doramectin treatment. Flow cytometry and TUNEL Assay were conducted to observe cell death-related effects.

RESULTS

MTT analysis of doramectin-treated cells displayed a decrease in cell growth compared to control. Apoptotic morphology was prominent in melanoma cells treated with doramectin. Increased autophagy was not detected by fluorometric microscopic analysis. Flow cytometry analysis of doramectin-treated cells showed apoptosis as a major mode of cell death with some necrosis.

CONCLUSION

Doramectin induces a novel cell-death mechanism in melanoma compared to other forms of cancer and should be studied as an effective anti-cancer agent for melanoma treatment.

Overcoming fixation and permeabilization challenges in flow cytometry by optical barcoding and multi-pass acquisition

The fixation and permeabilization of cells are essential for labeling intracellular biomarkers in flow cytometry. However, these chemical treatments often alter fragile targets, such as cell surface and fluorescent proteins (FPs), and can destroy chemically-sensitive fluorescent labels. This reduces measurement accuracy and introduces compromises into sample workflows, leading to losses in data quality. Here, we demonstrate a novel multi-pass flow cytometry approach to address this long-standing problem. Our technique utilizes individual cell barcoding with laser particles, enabling sequential analysis of the same cells with single-cell resolution maintained. Chemically-fragile protein markers and their fluorochrome conjugates are measured prior to destructive sample processing and adjoined to subsequent measurements of intracellular markers after fixation and permeabilization. We demonstrate the effectiveness of our technique in accurately measuring intracellular FPs and methanol-sensitive antigens and fluorophores, along with various surface and intracellular markers. This approach significantly enhances assay flexibility, enabling accurate and comprehensive cellular analysis without the constraints of conventional one-time measurement flow cytometry. This innovation paves new avenues in flow cytometry for a wide range of applications in immuno-oncology, stem cell research, and cell biology.

Epstein-Barr virus induces host shutoff extensively via BGLF5-independent mechanisms

Epstein-Barr virus (EBV) is a ubiquitous oncogenic virus associated with multiple cancers and autoimmune diseases. Unlike most herpesviruses, EBV reactivation from latency occurs asymptomatically, allowing it to spread efficiently to other hosts. However, available models are limited by the inefficient and asynchronous reactivation from latency into lytic replication. To address this problem, we develop a dual-fluorescent lytic reporter (DFLR) EBV that specifically labels cells in the early and late stages of replication. Using lymphoblastoid cell lines transformed by DFLR EBV as a model for EBV reactivation in B cells, we observe extensive reprogramming of the host cell transcriptome during lytic cycle progression. This includes widespread shutoff of host gene expression and disruption of mRNA processing. Unexpectedly, host shutoff remains extensive even in cells infected with DFLR EBV deleted for the BGLF5 nuclease. These findings implicate BGLF5-independent mechanisms as the primary drivers of host transcriptome remodeling during EBV lytic replication.

Long non-coding RNA GRASLND links melanoma differentiation and interferon-gamma response

Melanoma is a highly malignant tumor, that stands as the most lethal form of skin cancer and is characterized by notable phenotypic plasticity and intratumoral heterogeneity. Melanoma plasticity is involved in tumor growth, metastasis and therapy resistance. Long non-coding RNAs (lncRNAs) could influence plasticity due to their regulatory function. However, their role and mode of action are poorly studied. Here, we show a relevance of lncRNA GRASLND in melanoma differentiation and IFNγ signaling. GRASLND knockdown revealed switching of differentiated, melanocytic melanoma cells towards a dedifferentiated, slow-proliferating and highly-invasive cell state. Interestingly, GRASLND is overexpressed in differentiated melanomas and associated with poor prognosis. Accordingly, we found GRASLND expressed in immunological "cold" tumors and it negatively correlates with gene signatures of immune response activation. In line, silencing of GRASLND under IFNγ enhanced the expression of IFNγ-stimulated genes, including HLA-I antigen presentation, demonstrating suppressive activity of GRASLND on IFNγ signaling. Our findings demonstrate that in differentiated melanomas elevated expression of GRASLND interferes with anti-tumor effects of IFNγ, suggesting a role of GRASLND in tumor immune evasion.

Spontaneously regenerative corticospinal neurons in mice

The spinal cord receives inputs from the cortex via corticospinal neurons (CSNs). While predominantly a contralateral projection, a less-investigated minority of its axons terminate in the ipsilateral spinal cord. We analyzed the spatial and molecular properties of these ipsilateral axons and their post-synaptic targets in mice and found they project primarily to the ventral horn, including directly to motor neurons. Barcode-based reconstruction of the ipsilateral axons revealed a class of primarily bilaterally-projecting CSNs with a distinct cortical distribution. The molecular properties of these ipsilaterally-projecting CSNs (IP-CSNs) are strikingly similar to the previously described molecular signature of embryonic-like regenerating CSNs. Finally, we show that IP-CSNs are spontaneously regenerative after spinal cord injury. The discovery of a class of spontaneously regenerative CSNs may prove valuable to the study of spinal cord injury. Additionally, this work suggests that the retention of juvenile-like characteristics may be a widespread phenomenon in adult nervous systems.

Tertiary lymphoid structure-related immune infiltrates in NSCLC tumor lesions correlate with low tumor-reactivity of TIL products

Adoptive transfer of tumor infiltrating lymphocytes (TIL therapy) has proven highly effective for treating solid cancers, including non-small cell lung cancer (NSCLC). However, not all patients benefit from this therapy for yet unknown reasons. Defining markers that correlate with high tumor-reactivity of the autologous TIL products is thus key for achieving better tailored immunotherapies. We questioned whether the composition of immune cell infiltrates correlated with the tumor-reactivity of expanded TIL products. Unbiased flow cytometry analysis of immune cell infiltrates of 26 early-stage and 20 late-stage NSCLC tumor lesions was used for correlations with the T cell differentiation and activation status, and with the expansion rate and anti-tumor response of generated TIL products. The composition of tumor immune infiltrates was highly variable between patients. Spearman's Rank Correlation revealed that high B cell infiltration negatively correlated with the tumor-reactivity of the patient's expanded TIL products, as defined by cytokine production upon exposure to autologous tumor digest. In-depth analysis revealed that tumor lesions with high B cell infiltrates contained tertiary lymphoid structure (TLS)-related immune infiltrates, including BCL6+ antibody-secreting B cells, IgD+BCL6+ B cells and CXCR5+BLC6+ CD4+ T cells, and higher percentages of naïve CD8+ T cells. In conclusion, the composition of immune cell infiltrates in NSCLC tumors associates with the functionality of the expanded TIL product. Our findings may thus help improve patient selection for TIL therapy.

Overcoming fixation and permeabilization challenges in flow cytometry by optical barcoding and multi-pass acquisition

The fixation and permeabilization of cells are essential for labeling intracellular biomarkers in flow cytometry. However, these chemical treatments often alter fragile targets, such as cell surface and fluorescent proteins, and can destroy chemically-sensitive fluorescent labels. This reduces measurement accuracy and introduces compromises into sample workflows, leading to losses in data quality. Here, we demonstrate a novel multi-pass flow cytometry approach to address this long-standing problem. Our technique utilizes individual cell barcoding with laser particles, enabling sequential analysis of the same cells with single-cell resolution maintained. Chemically-fragile protein markers and their fluorochrome conjugates are measured prior to destructive sample processing and adjoined to subsequent measurements of intracellular markers after fixation and permeabilization. We demonstrate the effectiveness of our technique in accurately measuring intracellular fluorescent proteins and methanol-sensitive antigens and fluorophores, along with various surface and intracellular markers. This approach significantly enhances assay flexibility, enabling accurate and comprehensive cell analysis without the constraints of conventional one-time measurement flow cytometry. This innovation paves new avenues in flow cytometry for a wide range of applications in immuno-oncology, stem cell research, and cell biology.

IGHMBP2 deletion suppresses translation and activates the integrated stress response

IGHMBP2 is a nonessential, superfamily 1 DNA/RNA helicase that is mutated in patients with rare neuromuscular diseases SMARD1 and CMT2S. IGHMBP2 is implicated in translational and transcriptional regulation via biochemical association with ribosomal proteins, pre-rRNA processing factors, and tRNA-related species. To uncover the cellular consequences of perturbing IGHMBP2, we generated full and partial IGHMBP2 deletion K562 cell lines. Using polysome profiling and a nascent protein synthesis assay, we found that IGHMBP2 deletion modestly reduces global translation. We performed Ribo-seq and RNA-seq and identified diverse gene expression changes due to IGHMBP2 deletion, including ATF4 up-regulation. With recent studies showing the integrated stress response (ISR) can contribute to tRNA metabolism-linked neuropathies, we asked whether perturbing IGHMBP2 promotes ISR activation. We generated ATF4 reporter cell lines and found IGHMBP2 knockout cells demonstrate basal, chronic ISR activation. Our work expands upon the impact of IGHMBP2 in translation and elucidates molecular mechanisms that may link mutant IGHMBP2 to severe clinical phenotypes.

Chromosome evolution screens recapitulate tissue-specific tumor aneuploidy patterns

Whole chromosome and arm-level copy number alterations occur at high frequencies in tumors, but their selective advantages, if any, are poorly understood. Here, utilizing unbiased whole chromosome genetic screens combined with in vitro evolution to generate arm- and subarm-level events, we iteratively selected the fittest karyotypes from aneuploidized human renal and mammary epithelial cells. Proliferation-based karyotype selection in these epithelial lines modeled tissue-specific tumor aneuploidy patterns in patient cohorts in the absence of driver mutations. Hi-C-based translocation mapping revealed that arm-level events usually emerged in multiples of two via centromeric translocations and occurred more frequently in tetraploids than diploids, contributing to the increased diversity in evolving tetraploid populations. Isogenic clonal lineages enabled elucidation of pro-tumorigenic mechanisms associated with common copy number alterations, revealing Notch signaling potentiation as a driver of 1q gain in breast cancer. We propose that intrinsic, tissue-specific proliferative effects underlie tumor copy number patterns in cancer.

Differential Response of Human Dendritic Cells upon Stimulation with Encapsulated or Non-Encapsulated Isogenic Strains of Porphyromonas gingivalis

During periodontitis, the extracellular capsule of Porphyromonas gingivalis favors alveolar bone loss by inducing Th1 and Th17 patterns of lymphocyte response in the infected periodontium. Dendritic cells recognize bacterial antigens and present them to T lymphocytes, defining their activation and polarization. Thus, dendritic cells could be involved in the Th1 and Th17 response induced against the P. gingivalis capsule. Herein, monocyte-derived dendritic cells were obtained from healthy individuals and then stimulated with different encapsulated strains of P. gingivalis or two non-encapsulated isogenic mutants. Dendritic cell differentiation and maturation were analyzed by flow cytometry. The mRNA expression levels for distinct Th1-, Th17-, or T-regulatory-related cytokines and transcription factors, as well as TLR2 and TLR4, were assessed by qPCR. In addition, the production of IL-1β, IL-6, IL-23, and TNF-α was analyzed by ELISA. The encapsulated strains and non-encapsulated mutants of P. gingivalis induced dendritic cell maturation to a similar extent; however, the pattern of dendritic cell response was different. In particular, the encapsulated strains of P. gingivalis induced higher expression of IRF4 and NOTCH2 and production of IL-1β, IL-6, IL-23, and TNF-α compared with the non-encapsulated mutants, and thus, they showed an increased capacity to trigger Th1 and Th17-type responses in human dendritic cells.

Dynamic hydrogen peroxide levels reveal a rate-dependent sensitivity in B-cell lymphoma signaling

Although in vivo extracellular microenvironments are dynamic, most in vitro studies are conducted under static conditions. Here, we exposed diffuse large B-cell lymphoma (DLBCL) cells to gradient increases in the concentration of hydrogen peroxide (H2O2), thereby capturing some of the dynamics of the tumour microenvironment. Subsequently, we measured the phosphorylation response of B-cell receptor (BCR) signalling proteins CD79a, SYK and PLCγ2 at a high temporal resolution via single-cell phospho-specific flow cytometry. We demonstrated that the cells respond bimodally to static extracellular H2O2, where the percentage of cells that respond is mainly determined by the concentration. Computational analysis revealed that the bimodality results from a combination of a steep dose-response relationship and cell-to-cell variability in the response threshold. Dynamic gradient inputs of varying durations indicated that the H2O2 concentration is not the only determinant of the signalling response, as cells exposed to more shallow gradients respond at lower H2O2 levels. A minimal model of the proximal BCR network qualitatively reproduced the experimental findings and uncovered a rate-dependent sensitivity to H2O2, where a lower rate of increase correlates to a higher sensitivity. These findings will bring us closer to understanding how cells process information from their complex and dynamic in vivo environments.

Comparative Proteomic Analysis Reveals Candidate Pathways Related to the Effect of Different Light Qualities on the Development of Mycelium and Fruiting Body of Pleurotus ostreatus

Light affects the morphology and physiology of Pleurotus ostreatus. However, the underlying molecular mechanism of this effect remains unclear. In this study, a label-free comparative proteomic analysis was conducted to investigate the global protein expression profile of the mycelia and fruiting bodies of P. ostreatus PH11 growing under four different light quality treatments. Among all the 2234 P. ostreatus proteins, 1349 were quantifiable under all tested conditions. A total of 1100 differentially expressed proteins were identified by comparing the light group data with those of the darkness group. GO and KEGG enrichment analyses indicated that the oxidative phosphorylation, proteasome, and mRNA surveillance pathways were the most related pathways under the light condition. qRT-PCR verified that the expression of the white collar 1 protein was significantly enhanced under white light. Additionally, glutamine synthetase and aldehyde dehydrogenase played important roles during light exposure. This study provides valuable insight into the P. ostreatus light response mechanism, which will lay the foundation for improved cultivation.

A flow-cytometry-based pipeline for the rapid quantification of C2C12 cell differentiation

The C2C12 cell line represents a simple in vitro model for cell differentiation. Here, we present a flow-cytometry-based pipeline to quantitate C2C12 cell differentiation based on myosin heavy-chain marker expression. We describe steps for cell seeding, transfection, drug treatment, differentiation, and labeling. We then detail procedures for flow cytometry acquisition and introduce the R script FlowFate for automated analysis, including the study of dose-dependent effects of GFP-tagged genes on differentiation. For complete details on the use and execution of this protocol, please refer to Chippalkatti et al. (2023).1.

Variation in transcription regulator expression underlies differences in white-opaque switching between the SC5314 reference strain and the majority of Candida albicans clinical isolates

Candida albicans, a normal member of the human microbiome and an opportunistic fungal pathogen, undergoes several morphological transitions. One of these transitions is white-opaque switching, where C. albicans alternates between 2 stable cell types with distinct cellular and colony morphologies, metabolic preferences, mating abilities, and interactions with the innate immune system. White-to-opaque switching is regulated by mating type; it is repressed by the a1/α2 heterodimer in a/α cells, but this repression is lifted in a/a and α/α mating type cells (each of which are missing half of the repressor). The widely used C. albicans reference strain, SC5314, is unusual in that white-opaque switching is completely blocked when the cells are a/α; in contrast, most other C. albicans a/α strains can undergo white-opaque switching at an observable level. In this paper, we uncover the reason for this difference. We show that, in addition to repression by the a1/α2 heterodimer, SC5314 contains a second block to white-opaque switching: 4 transcription regulators of filamentous growth are upregulated in this strain and collectively suppress white-opaque switching. This second block is missing in the majority of clinical strains, and, although they still contain the a1/α2 heterodimer repressor, they exhibit a/α white-opaque switching at an observable level. When both blocks are absent, white-opaque switching occurs at very high levels. This work shows that white-opaque switching remains intact across a broad group of clinical strains, but the precise way it is regulated and therefore the frequency at which it occurs varies from strain to strain.

Streamlining N-terminally anchored yeast surface display via structural insights into S. cerevisiae Pir proteins

Surface display co-opts yeast's innate ability to embellish its cell wall with mannoproteins, thus converting the yeast's outer surface into a growing and self-sustaining catalyst. However, the efficient toolbox for converting the enzyme of interest into its surface-displayed isoform is currently lacking, especially if the isoform needs to be anchored to the cell wall near the isoform's N-terminus, e.g., through a short GPI-independent protein anchor. Aiming to advance such N-terminally anchored surface display, we employed in silico and machine-learning strategies to study the 3D structure, function, genomic organisation, and evolution of the Pir protein family, whose members evolved to covalently attach themselves near their N-terminus to the β-1,3-glucan of the cell wall. Through the newly-gained insights, we rationally engineered 14 S. cerevisiae Hsp150 (Pir2)-based fusion proteins. We quantified their performance, uncovering guidelines for efficient yeast surface display while developing a construct that promoted a 2.5-fold more efficient display of a reporter protein than the full-length Hsp150. Moreover, we developed a Pir-tag, i.e., a peptide spanning only 4.5 kDa but promoting as efficient surface display of a reporter protein as the full-length Hsp150. These constructs fortify the existing surface display toolbox, allowing for a prompt and routine refitting of intracellular proteins into their N-terminally anchored isoforms.

SGRN: A Cas12a-driven Synthetic Gene Regulatory Network System

Gene regulatory networks, which control gene expression patterns in development and in response to stimuli, use regulatory logic modules to coordinate inputs and outputs. One example of a regulatory logic module is the gene regulatory cascade (GRC), where a series of transcription factor genes turn on in order. Synthetic biologists have derived artificial systems that encode regulatory rules, including GRCs. Furthermore, the development of single-cell approaches has enabled the discovery of gene regulatory modules in a variety of experimental settings. However, the tools available for validating these observations remain limited. Based on a synthetic GRC using DNA cutting-defective Cas9 (dCas9), we designed and implemented an alternative synthetic GRC utilizing DNA cutting-defective Cas12a (dCas12a). Comparing the ability of these two systems to express a fluorescent reporter, the dCas9 system was initially more active, while the dCas12a system was more streamlined. Investigating the influence of individual components of the systems identified nuclear localization as a major driver of differences in activity. Improving nuclear localization for the dCas12a system resulted in 1.5-fold more reporter-positive cells and a 15-fold increase in reporter intensity relative to the dCas9 system. We call this optimized system the "Synthetic Gene Regulatory Network" (SGRN, pronounced "sojourn").

Transcriptome Dynamics of Pseudomonas aeruginosa during Transition from Overlapping To Non-Overlapping Cell Cycles

Bacteria either duplicate their chromosome once per cell division or a new round of replication is initiated before the cells divide, thus cell cycles overlap. Here, we show that the opportunistic pathogen Pseudomonas aeruginosa switches from fast growth with overlapping cell cycles to sustained slow growth with only one replication round per cell division when cultivated under standard laboratory conditions. The transition was characterized by fast-paced, sequential changes in transcriptional activity along the ori-ter axis of the chromosome reflecting adaptation to the metabolic needs during both growth phases. Quorum sensing (QS) activity was highest at the onset of the slow growth phase with non-overlapping cell cycles. RNA sequencing of subpopulations of these cultures sorted based on their DNA content, revealed a strong gene dosage effect as well as specific expression patterns for replicating and nonreplicating cells. Expression of flagella and mexE, involved in multidrug efflux was restricted to cells that did not replicate, while those that did showed a high activity of the cell division locus and recombination genes. A possible role of QS in the formation of these subpopulations upon switching to non-overlapping cell cycles could be a subject of further research. IMPORTANCE The coordination of gene expression with the cell cycle has so far been studied only in a few bacteria, the bottleneck being the need for synchronized cultures. Here, we determined replication-associated effects on transcription by comparing Pseudomonas aeruginosa cultures that differ in their growth mode and number of replicating chromosomes. We further show that cell cycle-specific gene regulation can be principally identified by RNA sequencing of subpopulations from cultures that replicate only once per cell division and that are sorted according to their DNA content. Our approach opens the possibility to study asynchronously growing bacteria from a wide phylogenetic range and thereby enhance our understanding of the evolution of cell cycle control on the transcriptional level.

Comparative Proteomic Analyses within Three Developmental Stages of the Mushroom White Hypsizygus marmoreus

(1) Background: The Hypsizygus marmoreus is a popular edible mushroom in East Asian markets. In a previous study, we reported the proteomic analyses of different developmental stages of H. marmoreus, from primordium to mature fruiting body. However, the growth and protein expression changes from scratching to primordium are unclear. (2) Methods: A label-free LC-MS/MS quantitative proteomic analysis technique was adopted to obtain the protein expression profiles of three groups of samples collected in different growth stages from scratching to the tenth day after scratching. The Pearson's correlation coefficient analysis and principal component analysis were performed to reveal the correlation among samples. The differentially expressed proteins (DEPs) were organized. Gene Ontology (GO) analysis was performed to divide the DEPs into different metabolic processes and pathways. (3) Results: From the 3rd day to the 10th day after scratching, mycelium recovered gradually and formed primordia. Compared with the Rec stage, 218 highly expressed proteins were identified in the Knot stage. Compared with the Pri stage, 217 highly expressed proteins were identified in the Rec stage. Compared with the Pri stage, 53 highly expressed proteins were identified in the Knot stage. A variety of the same highly expressed proteins were identified in these three developmental stages, including: glutathione S-transferase, acetyltransferase, importin, dehydrogenase, heat-shock proteins, ribosomal proteins, methyltransferase, etc. The key pathways in the development of H. marmoreus are metabolic process, catabolic process, oxidoreductase activity and hydrolase activity. DEPs in the Knot or Pri stages compared with the Rec stage were significantly decreased in the metabolic-, catabolic- and carbohydrate-related process; and the oxidoreductase, peptidase, and hydrolase activity, which can serve as targets for selectable molecular breeding in H. marmoreus. A total of 2000 proteins were classified into eight different modules by WGCNA, wherein 490 proteins were classified into the turquoise module. (4) Conclusions: Generally, from the 3rd day to the 10th day after scratching, mycelium recovered gradually and formed primordia. Importin, dehydrogenase, heat-shock proteins, ribosomal proteins, transferases were all highly expressed in these three developmental stages. DEPs in the Rec stage compared with the Knot or Pri stages were significantly enriched in the metabolic-, catabolic- and carbohydrate-related process; and in oxidoreductase, peptidase and hydrolase activities. This research contributes to the understanding of the mechanisms of the development changes before primordium of H. marmoreus.

Vaccination induces HIV broadly neutralizing antibody precursors in humans

Broadly neutralizing antibodies (bnAbs) can protect against HIV infection but have not been induced by human vaccination. A key barrier to bnAb induction is vaccine priming of rare bnAb-precursor B cells. In a randomized, double-blind, placebo-controlled phase 1 clinical trial, the HIV vaccine-priming candidate eOD-GT8 60mer adjuvanted with AS01B had a favorable safety profile and induced VRC01-class bnAb precursors in 97% of vaccine recipients with median frequencies reaching 0.1% among immunoglobulin G B cells in blood. bnAb precursors shared properties with bnAbs and gained somatic hypermutation and affinity with the boost. The results establish clinical proof of concept for germline-targeting vaccine priming, support development of boosting regimens to induce bnAbs, and encourage application of the germline-targeting strategy to other targets in HIV and other pathogens.

The memory B cell response to influenza vaccination is impaired in older persons

Influenza infection imparts an age-related increase in mortality and morbidity. The most effective countermeasure is vaccination; however, vaccines offer modest protection in older adults. To investigate how aging impacts the memory B cell response, we track hemagglutinin-specific B cells by indexed flow sorting and single-cell RNA sequencing (scRNA-seq) in 20 healthy adults that were administered the trivalent influenza vaccine. We demonstrate age-related skewing in the memory B cell compartment 6 weeks after vaccination, with younger adults developing hemagglutinin-specific memory B cells with an FcRL5+ "atypical" phenotype, showing evidence of somatic hypermutation and positive selection, which happened to a lesser extent in older persons. We use publicly available scRNA-seq from paired human lymph node and blood samples to corroborate that FcRL5+ atypical memory B cells can derive from germinal center (GC) precursors. Together, this study shows that the aged human GC reaction and memory B cell response following vaccination is defective.

Construction and Evaluation of a Risk Score Model for Lymph Node Metastasis-Associated Circadian Clock Genes in Esophageal Squamous Carcinoma

Background: Studies suggested that circadian clock genes (CCGs) in human esophageal squamous carcinoma (ESCC) samples are dysregulated. However, the relevance of CCGs to lymph node metastasis (LNM) and prognosis of ESCC remains unclear. Methods: The differentially expressed genes (DEGs) between normal and ESCC samples in The Cancer Genome Atlas database (TCGA) database were intersected with the genes associated with LNM (LNMGs) in ESCC samples and 300 CCGs to obtain the differentially expressed LNM-associated CCGs (DE-LNM-CCGs). The risk model was constructed by Cox regression analysis in the TCGA-ESCC training set, and the accuracy of the risk model was verified by risk profile and overall survival profile. Furthermore, differences of 23 immune cells, 13 immune functions, and immune checkpoint molecules between the high- and low-risk groups were assessed using the single-sample gene set enrichment analysis (ssGSEA) algorithm. Gene set enrichment analysis (GSEA) was conducted to investigate the functional differences between low- and high-risk groups. Finally, we validated the mRNA expression levels of prognostic model genes by quantitative real-time polymerase chain reaction (qRT-PCR). Results: A total of six DE-LNM-CCGs were identified in TCGA-ESCC. TP53 and NAGLU were selected by Cox regression analysis to construct the risk model. Risk profile plots, overall survival plots, and validation results of the risk model in the validation set indicated that the constructed risk model was reliable. The result of ssGSEA showed that the percentages of activated B cells, activated dendritic cells, effector memory CD8 T cells, immune function in neutrophils, plasmacytoid dendritic cells, T cell co-inhibition, and Type 17 T helper cells were different between the high- and low-risk groups. In addition, the expression of CD274, PDCD1, TNFRSF18, and TNFRSF9 was dysregulated between the high- and low-risk groups. GSEA revealed that the high-risk group was associated with cell differentiation, oxidative phosphorylation, and steroid biosynthesis pathways, while the low-risk group was associated with chromosome, ECM–receptor interaction, and other pathways. Finally, qRT-PCR results showed that the mRNA expression levels of two prognostic genes were consistent with TCGA. Conclusion: In conclusion, the risk model constructed based on TP53 and NAGLU could accurately predict the prognosis.

Using colony size to measure fitness in Saccharomyces cerevisiae

Competitive fitness assays in liquid culture have been a mainstay for characterizing experimental evolution of microbial populations. Growth of microbial strains has also been extensively characterized by colony size and could serve as a useful alternative if translated to per generation measurements of relative fitness. To examine fitness based on colony size, we established a relationship between cell number and colony size for strains of Saccharomyces cerevisiae robotically pinned onto solid agar plates in a high-density format. This was used to measure growth rates and estimate relative fitness differences between evolved strains and their ancestors. After controlling for edge effects through both normalization and agar-trimming, we found that colony size is a sensitive measure of fitness, capable of detecting 1% differences. While fitnesses determined from liquid and solid mediums were not equivalent, our results demonstrate that colony size provides a sensitive means of measuring fitness that is particularly well suited to measurements across many environments.

Comprehensive Analysis of GDF10 Methylation Site-Associated Genes as Prognostic Markers for Endometrial Cancer

Growth differentiation factor-10 (GDF10) with its methylation trait has recently been found to play a crucial regulatory and communication role in cancers. This investigation aims to identify GDF10 methylation site-associated genes that are closely associated with endometrial cancer (EC) patients' survival based on normal and UCEC samples from the UCSC Xena database. Our study revealed for the first time that EC exhibited significantly higher levels of GDF10 promoter methylation in comparison with normal tissues. Multiple differentiated methylation sites, which have prognostic value due to their apparent survival differences, were found in the GDF10 promoter region. We performed weighted gene coexpression network analysis (WGCNA) on EC tissues and paraneoplastic tissues while using these differentially methylated sites as phenotypes for selecting the most correlated key modules and their internal genes. To obtain a gene set, the key module genes and differentially expressed genes (DEGs) of EC were intersected. The least absolute shrinkage and selection operator (LASSO) regression along with multivariate Cox regression were performed from the gene set and we screened out the key genes B4GALNT3, DNAJC22, and GREB1. Finally, a prognostic model was validated for effectiveness based on these genes. Additionally, Kaplan-Meier analysis and time-dependent receiver operating characteristics (ROC) were applied to assess and verify the model, and they showed good prognosis prediction. Moreover, the differences in risk scores were statistically significant with age, tumor stage, and grade. They may be related to the immune infiltration of tumors as well. In conclusion, based on the methylation-related genes associated with GDF10, we developed a prognosis model for EC patients. It might provide a fresh view for further research and treatment of EC.

Identification and analysis of mitochondria-related key genes of heart failure

Mitochondria-induced cell death is a vital mechanism of heart failure (HF). Thus, identification of mitochondria-related genes (Mito-RGs) based on transcriptome sequencing data of HF might provide novel diagnostic markers and therapeutic targets for HF. First, bioinformatics analysis was conducted on the GSE57338, GSE76701, GSE136547, and GSE77399 datasets in the Gene Expression Omnibus. Next, we analyzed HF-Mito differentially expressed genes (DEGs) using the protein-protein interaction (PPI) network for obtaining critical genes and exploring their functions. Subsequently, immune cell scores of the HF and normal groups were compared. The potential alteration mechanisms of the key genes were investigated by constructing a competing endogenous RNA network. Finally, we predicted potential therapeutic agents and validated the expression levels of the key genes. Twenty-three HF-Mito DEGs were acquired in the GSE57338 dataset, and the PPI network obtained four key genes, including IFIT3, XAF1, RSAD2, and MX1. According to gene set enrichment analysis, the key genes showed high enrichment in myogenesis and hypoxia. Immune cell analysis demonstrated that aDCs, B cells, and 20 other immune cell types varied between the HF and normal groups. Moreover, we observed that H19 might affect the expression of IFIT3, AXF1, and RSAD2. PCGEM1 might regulate RSAD2 expression. A total of 515 potential therapeutic drugs targeting the key genes, such as tretinoin, silicon dioxide, and bisphenol A, were acquired. Finally, IFIT3, RSAD2, and MX1 expression increased in HF samples compared with normal samples in the GSE76701 dataset, conforming to the GSE57338 dataset analysis. This work screened four key genes, namely, IFIT3, XAF1, RSAD2, and MX1, which can be further explored in subsequent studies for their specific molecular mechanisms in HF.

Integrative Analysis of Angiogenesis-Related Long Non-Coding RNA and Identification of a Six-DEARlncRNA Signature Associated with Prognosis and Therapeutic Response in Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is a lethal gastrointestinal malignancy worldwide. We aimed to identify an angiogenesis-related lncRNAs (ARlncRNAs) signature that could predict the prognosis in ESCC. The GSE53624 and GSE53622 datasets were derived from the GEO database. The differently expressed ARlncRNAs (DEARlncRNAs) were retrieved by the weighted gene co-expression network analysis (WGCNA), differential expression analysis, and correlation analysis. Optimal lncRNA biomarkers were screened from the training set and the six-DEARlncRNA signature comprising AP000696.2, LINC01711, RP11-70C1.3, AP000487.5, AC011997.1, and RP11-225N10.1 could separate patients into high- and low-risk groups with markedly different survival. The validation of the reliability of the risk model was performed by the Kaplan-Meier test, ROC curves, and risk curves in the test set and validation set. Predictive independence analysis indicated that risk score is an independent prognostic biomarker for predicting the prognosis of ESCC patients. Subsequently, a ceRNA regulatory network and functional enrichment analysis were performed. The IC50 test revealed that patients in the high-risk group were resistant to Gefitinib and Lapatinib. Finally, the six DEARlncRNAs were detected by qRT-PCR. In conclusion, we demonstrated a novel ARlncRNA signature as an independent prognostic factor to distinguish the risk of ESCC patients and benefit the personalized clinical applications.

An R-based reproducible and user-friendly preprocessing pipeline for CyTOF data

Mass cytometry (CyTOF) has become a method of choice for in-depth characterization of tissue heterogeneity in health and disease, and is currently implemented in multiple clinical trials, where higher quality standards must be met. Currently, preprocessing of raw files is commonly performed in independent standalone tools, which makes it difficult to reproduce. Here, we present an R pipeline based on an updated version of CATALYST that covers all preprocessing steps required for downstream mass cytometry analysis in a fully reproducible way. This new version of CATALYST is based on Bioconductor’s SingleCellExperiment class and fully unit tested. The R-based pipeline includes file concatenation, bead-based normalization, single-cell deconvolution, spillover compensation and live cell gating after debris and doublet removal. Importantly, this pipeline also includes different quality checks to assess machine sensitivity and staining performance while allowing also for batch correction. This pipeline is based on open source R packages and can be easily be adapted to different study designs. It therefore has the potential to significantly facilitate the work of CyTOF users while increasing the quality and reproducibility of data generated with this technology.

BCyto: A shiny app for flow cytometry data analysis

Flow cytometry (FCM) is a widely used technique to simultaneously measure various characteristics of a cell or particle. Most of reliable software for FCM data analysis are commercially available at high costs. A few R packages are available; however, programming skills are required. In addition, currently, no package offers an interactive user interface (UI) with efficient tools for conventional FCM data analysis. BCyto is an open-source R package/shiny app that allows cell biologists with no coding skills to reliably analyse FCM data in R. BCyto provides intuitive axis transformation for better data visualization, easy visualization of compensation plots and modification of compensation matrices, fast generation of backgating and overlay plots, as well as built-in proliferation and dimensionality reduction tools. BCyto will not only improve accessibility to FCM data analysis but might also serve as an environment for the development of further R-based computational FCM algorithms.

Responses of the Mushroom Pleurotus ostreatus under Different CO2 Concentration by Comparative Proteomic Analyses

BACKGROUND

Pleurotus ostreatus is a popular edible mushroom in East Asian markets. Research on the responses of P. ostreatus under different carbon dioxide concentrations is limited.

METHODS

Label-free LC-MS/MS quantitative proteomics analysis technique was adopted to obtain the protein expression profiles of P. ostreatus fruiting body pileus collected under different carbon dioxide concentrations. The Pearson correlation coefficient analysis and principal component analysis were performed to reveal the correlation among samples. The differentially expressed proteins (DEPs) were organized. Gene ontology analysis was performed to divide the DEPs into different metabolic processes and pathways.

RESULTS

The expansion of stipes was inhibited in the high CO2 group compared with that in the low CO2 group. There were 415 DEPs (131 up- and 284 down-regulated) in P. ostreatus PH11 treated with 1% CO2 concentration compared with P. ostreatus under atmospheric conditions. Proteins related to hydrolase activity, including several amidohydrolases and cell wall synthesis proteins, were highly expressed under high CO2 concentration. Most of the kinases and elongation factors were significantly down-regulated under high CO2 concentration. The results suggest that the metabolic regulation and development processes were inhibited under high CO2 concentrations. In addition, the sexual differentiation process protein Isp4 was inhibited under high CO2 concentrations, indicating that the sexual reproductive process was also inhibited under high CO2 concentrations, which is inconsistent with the small fruiting body pileus under high CO2 concentrations.

CONCLUSIONS

This research reports the proteome analysis of commercially relevant edible fungi P. ostreatus under different carbon dioxide concentrations. This study deepens our understanding of the mechanism for CO2-induced morphological change in the P. ostreatus fruiting body, which will facilitate the artificial cultivation of edible mushrooms.

Optochemical Control of TET Dioxygenases Enables Kinetic Insights into the Domain-Dependent Interplay of TET1 and MBD1 while Oxidizing and Reading 5-Methylcytosine

Methyl-CpG binding domain (MBD) proteins and ten-eleven-translocation (TET) dioxygenases are the readers and erasers of 5-methylcytosine (5mC), the central epigenetic mark of mammalian DNA. We employ light-activatable human TET1 controlled by a genetically encoded photocaged serine to enable in vivo kinetic studies of their interplay at the common substrate methylated cytosine–guanine (mCpG). We identify the multidomain reader MBD1 to negatively regulate TET1-catalyzed 5mC oxidation kinetics via its mCpG-binding MBD domain. However, we also identify the third Cys-x-x-Cys (CXXC3) domain of MBD1 to promote oxidation kinetics by TET1, dependent on its ability to bind nonmethylated CpG, the final product of TET-mediated mCpG oxidation and active demethylation. In contrast, we do not observe differences in TET1 regulation for MBD1 variants with or without the transcriptional repressor domain. Our approach reveals a complex, domain-dependent interplay of these readers and erasers of 5mC with different domain-specific contributions of MBD1 to the overall kinetics of TET1-catalyzed global 5mC oxidation kinetics that contribute to a better understanding of dynamic methylome shaping.

Cas9-induced large deletions and small indels are controlled in a convergent fashion

Repair of Cas9-induced double-stranded breaks results primarily in formation of small insertions and deletions (indels), but can also cause potentially harmful large deletions. While mechanisms leading to the creation of small indels are relatively well understood, very little is known about the origins of large deletions. Using a library of clonal NGS-validated mouse embryonic stem cells deficient for 32 DNA repair genes, we have shown that large deletion frequency increases in cells impaired for non-homologous end joining and decreases in cells deficient for the central resection gene Nbn and the microhomology-mediated end joining gene Polq. Across deficient clones, increase in large deletion frequency was closely correlated with the increase in the extent of microhomology and the size of small indels, implying a continuity of repair processes across different genomic scales. Furthermore, by targeting diverse genomic sites, we identified examples of repair processes that were highly locus-specific, discovering a role for exonuclease Trex1. Finally, we present evidence that indel sizes increase with the overall efficiency of Cas9 mutagenesis. These findings may have impact on both basic research and clinical use of CRISPR-Cas9, in particular in conjunction with repair pathway modulation.

Boron clusters as broadband membrane carriers

The membrane translocation of hydrophilic substances constitutes a challenge for their application as therapeutic compounds and labelling probes^1–4. To remedy this, charged amphiphilic molecules have been classically used as carriers^3,5. However, such amphiphilic carriers may cause aggregation and non-specific membrane lysis^6,7. Here we show that globular dodecaborate clusters, and prominently B₁₂Br₁₂²⁻, can function as anionic inorganic membrane carriers for a broad range of hydrophilic cargo molecules (with molecular mass of 146–4,500 Da). We show that cationic and neutral peptides, amino acids, neurotransmitters, vitamins, antibiotics and drugs can be carried across liposomal membranes. Mechanistic transport studies reveal that the carrier activity is related to the superchaotropic nature of these cluster anions^8–12. We demonstrate that B₁₂Br₁₂²⁻ affects cytosolic uptake of different small bioactive molecules, including the antineoplastic monomethyl auristatin F, the proteolysis targeting chimera dBET1 and the phalloidin toxin, which has been successfully delivered in living cells for cytoskeleton labelling. We anticipate the broad and distinct delivery spectrum of our superchaotropic carriers to be the starting point of conceptually distinct cell-biological, neurobiological, physiological and pharmaceutical studies.

The superchaotropic nature of globular boron cluster anions enables direct passage of a wide range of molecules across lipid membranes.

Heterogeneity of Starved Yeast Cells in IF1 Levels Suggests the Role of This Protein in vivo

In mitochondria, a small protein IF1 suppresses the hydrolytic activity of ATP synthase and presumably prevents excessive ATP hydrolysis under conditions of energy deprivation. In yeast Saccharomyces cerevisiae, IF1 homologs are encoded by two paralogous genes: INH1 and STF1. INH1 expression is known to aggravate the deleterious effects of mitochondrial DNA (mtDNA) depletion. Surprisingly, no beneficial effects of INH1 and STF1 were documented for yeast so far, and the functions of INH1 and STF1 in wild type cells are unclear. Here, we put forward a hypothesis that INH1 and STF1 bring advantage during the fast start of proliferation after reentry into exponential growth from post-diauxic or stationary phases. We found that yeast cells increase the concentration of both proteins in the post-diauxic phase. Post-diauxic phase yeast cells formed two subpopulations distinct in Inh1p and Stf1p concentrations. Upon exit from the post-diauxic phase cells with high level of Inh1-GFP started growing earlier than cells devoid of Inh1-GFP. However, double deletion of INH1 and STF1 did not increase the lag period necessary for stationary phase yeast cells to start growing after reinoculation into the fresh medium. These results point to a redundancy of the mechanisms preventing uncontrolled ATP hydrolysis during energy deprivation.

Advances in automated real-time flow cytometry for monitoring of bioreactor processes

Flow cytometry and its technological possibilities have greatly advanced in the past decade as analysis tool for single cell properties and population distributions of different cell types in bioreactors. Along the way, some solutions for automated real-time flow cytometry (ART-FCM) were developed for monitoring of bioreactor processes without operator interference over extended periods with variable sampling frequency. However, there is still great potential for ART-FCM to evolve and possibly become a standard application in bioprocess monitoring and process control. This review first addresses different components of an ART-FCM, including the sampling device, the sample-processing unit, the unit for sample delivery to the flow cytometer and the settings for measurement of pre-processed samples. Also, available algorithms are presented for automated data analysis of multi-parameter fluorescence datasets derived from ART-FCM experiments. Furthermore, challenges are discussed for integration of fluorescence-activated cell sorting into an ART-FCM setup for isolation and separation of interesting subpopulations that can be further characterized by for instance omics-methods. As the application of ART-FCM is especially of interest for bioreactor process monitoring, including investigation of population heterogeneity and automated process control, a summary of already existing setups for these purposes is given. Additionally, the general future potential of ART-FCM is addressed.

CombiFlow: Flow cytometry-based identification and characterization of genetically and functionally distinct AML subclones

Many cancers, including leukemias, are dynamic oligoclonal diseases. Tools to identify and prospectively isolate genetically distinct clones for functional studies are needed. We describe our CombiFlow protocol, which is a combinatorial flow cytometry-based approach to identify and isolate such distinct clones. CombiFlow enables the visualization of clonal evolution during disease progression and the identification of potential relapse-inducing cells at minimal residual disease (MRD) time points. The protocol can be adapted to various research questions and allows functional studies on live sorted cell populations.

For complete details on the use and execution of this protocol, please refer to de Boer et al. (2018).

•

Subclones can be identified and sorted based on aberrant marker expression

•

Functional in vitro and in vivo studies can be performed on identified subclones

•

Combiflow can visualize the phenotype and track clonal evolution in AML

Comparative Proteomic Analysis within the Developmental Stages of the Mushroom White Hypsizygus marmoreus

(1) Background: The white Hypsizygus marmoreus is a popular edible mushroom in East Asia markets. Research on the systematic investigation of the protein expression changes in the cultivation process of this mushroom are few. (2) Methods: Label-free LC-MS/MS quantitative proteomics analysis technique was adopted to obtain the protein expression profiles of six groups of samples collected in different growth stages. A total of 3468 proteins were identified. The UpSetR plot analysis, Pearson correlation coefficient (PCC) analysis, and principal component (PC) analysis were performed to reveal the correlation among the six groups of samples. The differentially expressed proteins (DEPs) were organised by One-way ANOVA test and divided into four clusters. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to divide the DEPs into different metabolic processes and pathways in each cluster. (3) Results: The DEPs in cluster 1 are of the highest abundance in the mycelium and are mainly involved in protein biosynthesis, biosynthesis of cofactors, lipid metabolism, spliceosome, cell cycle regulation, and MAPK signaling pathway. The DEPs in cluster 2 are enriched in the stem and are mainly associated with protein biosynthesis, biosynthesis of cofactors, carbon, and energy metabolism. The DEPs in cluster 3 are highly expressed in the primordia and unmatured fruiting bodies and are related to amino acids metabolism, carbon and carbohydrate metabolism, protein biosynthesis and processing, biosynthesis of cofactors, cell cycle regulation, MAPK signaling pathway, ubiquitin-mediated proteolysis, and proteasome. The DEPs in cluster 4 are of the highest abundance in the cap and are mainly associated with spliceosome, endocytosis, nucleocytoplasmic transport, protein processing, oxidative phosphorylation, biosynthesis of cofactors, amino acids metabolism, and lipid metabolism. (4) Conclusions: This research reports the proteome analysis of different developmental stages during the cultivation of the commercially relevant edible fungi the white H. marmoreus. In the mycelium stage, most of the DEPs are associated with cell proliferation, signal response, and mycelium growth. In the primordia and unmatured fruiting bodies stage, the DEPs are mainly involved in biomass increase, cell proliferation, signal response, and differentiation. In the mature fruiting body stage, the DEPs in the stem are largely associated with cell elongation and increase in biomass, and most of the DEPs in the cap are mainly related to pileus expansion. Several carbohydrate-active enzymes, transcription factors, heat shock proteins, and some DEPs involved in MAPK and cAMP signaling pathways were determined. These proteins might play vital roles in metabolic processes and activities. This research can add value to the understanding of mechanisms concerning mushroom development during commercial production.

Comprehensive Analysis of the Immune and Prognostic Implication of MMP14 in Lung Cancer

More and more studies have indicated an association between immune infiltration in lung cancer and clinical outcomes. Matrix metalloproteinase 14 (MMP14) has been reported to be dysregulated in many types of tumors and involved in the development and progression of tumors. However, its contribution to cancer immunity was rarely reported. In the study, we found that MMP14 expression was distinctly upregulated in lung cancer specimens compared with nontumor lung specimens. High MMP14 expression predicted a poor prognosis of lung squamous cell carcinoma (LUSC) patients. Increased MMP14 expressions were observed to be positively related to high immune infiltration levels in most of the immune cells. A pathway enrichment analysis of 32 MMP14-associated immunomodulators indicated the involvement of T cell receptor signaling pathway and Toll-like receptor signaling pathway. Based on MMP14-associated immunomodulators, we applied multivariate assays to construct multiple-gene risk prediction signatures. We observed that risk scores were independently associated with overall survival. These data highlighted that MMP14 was involved in tumor immunity, indicating that MMP14 could serve as a novel prognostic biomarker and therapeutic target for lung cancer. Our data suggest that the four genes identified in this study may serve as valuable biomarkers of lung cancer patient outcomes.

Phage-delivered CRISPR-Cas9 for strain-specific depletion and genomic deletions in the gut microbiome

Mechanistic insights into the role of the human microbiome in the predisposition to and treatment of disease are limited by the lack of methods to precisely add or remove microbial strains or genes from complex communities. Here, we demonstrate that engineered bacteriophage M13 can be used to deliver DNA to Escherichia coli within the mouse gastrointestinal (GI) tract. Delivery of a programmable exogenous CRISPR-Cas9 system enables the strain-specific depletion of fluorescently marked isogenic strains during competitive colonization and genomic deletions that encompass the target gene in mice colonized with a single strain. Multiple mechanisms allow E. coli to escape targeting, including loss of the CRISPR array or even the entire CRISPR-Cas9 system. These results provide a robust and experimentally tractable platform for microbiome editing, a foundation for the refinement of this approach to increase targeting efficiency, and a proof of concept for the extension to other phage-bacterial pairs of interest.

Longitudinal immune dynamics of mild COVID-19 define signatures of recovery and persistence

SARS-CoV-2 has infected over 200 million and caused more than 4 million deaths to date. Most individuals (>80%) have mild symptoms and recover in the outpatient setting, but detailed studies of immune responses have focused primarily on moderate to severe COVID-19. We deeply profiled the longitudinal immune response in individuals with mild COVID-19 beginning with early time points post-infection (1-15 days) and proceeding through convalescence to >100 days after symptom onset. We correlated data from single cell analyses of peripheral blood cells, serum proteomics, virus-specific cellular and humoral immune responses, and clinical metadata. Acute infection was characterized by vigorous coordinated innate and adaptive immune activation that differed in character by age (young vs. old). We then characterized signals associated with recovery and convalescence to define and validate a new signature of inflammatory cytokines, gene expression, and chromatin accessibility that persists in individuals with post-acute sequelae of SARS-CoV-2 infection (PASC).

FYN and CD247: key Genes for Septic Shock Based on Bioinformatics and Meta-Analysis

BACKGROUND

Septic shock is sepsis accompanied by hemodynamic instability and high clinical mortality.

MATERIAL AND METHODS

GSE95233, GSE57065, GSE131761 gene-expression profiles of healthy control subjects and septic shock patients were downloaded from the Gene-Expression Omnibus (GEO) database, and differences of expression profiles and their intersection were analysed using GEO2R. Function and pathway enrichment analysis was performed on common differentially expressed genes (DEG), and key genes for septic shock were screened using a protein-protein interaction network created with STRING. Also, data from the GEO database were used for survival analysis for key genes, and a meta-analysis was used to explore expression trends of core genes. Finally, high-throughput sequencing using the blood of a murine sepsis model was performed to analyse the expression of CD247 and FYN in mice.

RESULTS

A total of 539 DEGs were obtained (p < 0.05). Gene ontology analysis showed that key genes were enriched in functions, such as immune response and T cell activity, and DEGs were enriched in signal pathways, such as T cell receptors. FYN and CD247 are in the centre of the protein-protein interaction network, and survival analysis found that they are positively correlated with survival from sepsis. Further, meta-analysis results showed that FYN could be useful for the prognosis of patients, and CD247 might distinguish between sepsis and systemic inflammatory response syndrome patients. Finally, RNA sequencing using a mouse septic shock model showed low expression of CD247 and FYN in this model.

CONCLUSION

FYN and CD247 are expected to become new biomarkers of septic shock.

Sickle-trait hemoglobin reduces adhesion to both CD36 and EPCR by Plasmodium falciparum-infected erythrocytes

Sickle-trait hemoglobin protects against severe Plasmodium falciparum malaria. Severe malaria is governed in part by the expression of the Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP1) that are encoded by var genes, specifically those variants that bind Endothelial Protein C Receptor (EPCR). In this study, we investigate the effect of sickle-trait on parasite var gene expression and function in vitro and in field-collected parasites. We mapped var gene reads generated from RNA sequencing in parasite cultures in normal and sickle-cell trait blood throughout the asexual lifecycle. We investigated sickle-trait effect on PfEMP1 interactions with host receptors CD36 and EPCR using static adhesion assays and flow cytometry. Var expression in vivo was compared by assembling var domains sequenced from total RNA in parasites infecting Malian children with HbAA and HbAS. Sickle-trait did not alter the abundance or type of var gene transcripts in vitro, nor the abundance of overall transcripts or of var functional domains in vivo. In adhesion assays using recombinant host receptors, sickle-trait reduced adhesion by 73-86% to CD36 and 83% to EPCR. Similarly, sickle-trait reduced the surface expression of EPCR-binding PfEMP1. In conclusion, Sickle-cell trait does not directly affect var gene transcription but does reduce the surface expression and function of PfEMP1. This provides a direct mechanism for protection against severe malaria conferred by sickle-trait hemoglobin. Trial Registration: ClinicalTrials.gov Identifier: NCT02645604.

Integration, exploration, and analysis of high-dimensional single-cell cytometry data using Spectre

As the size and complexity of high-dimensional (HD) cytometry data continue to expand, comprehensive, scalable, and methodical computational analysis approaches are essential. Yet, contemporary clustering and dimensionality reduction tools alone are insufficient to analyze or reproduce analyses across large numbers of samples, batches, or experiments. Moreover, approaches that allow for the integration of data across batches or experiments are not well incorporated into computational toolkits to allow for streamlined workflows. Here we present Spectre, an R package that enables comprehensive end-to-end integration and analysis of HD cytometry data from different batches or experiments. Spectre streamlines the analytical stages of raw data pre-processing, batch alignment, data integration, clustering, dimensionality reduction, visualization, and population labelling, as well as quantitative and statistical analysis. Critically, the fundamental data structures used within Spectre, along with the implementation of machine learning classifiers, allow for the scalable analysis of very large HD datasets, generated by flow cytometry, mass cytometry, or spectral cytometry. Using open and flexible data structures, Spectre can also be used to analyze data generated by single-cell RNA sequencing or HD imaging technologies, such as Imaging Mass Cytometry. The simple, clear, and modular design of analysis workflows allow these tools to be used by bioinformaticians and laboratory scientists alike. Spectre is available as an R package or Docker container. R code is available on Github (https://github.com/immunedynamics/spectre).

Long insert clone experimental evidence for assembly improvement and chimeric chromosomes detection in an allopentaploid beer yeast

Lager beer is made with the hybrid Saccharomyces pastorianus. Many publicly available S. pastorianus genome assemblies are highly fragmented due to the difficulties of assembling hybrid genomes, such as the presence of homeologous chromosomes from both parental types, and translocations between them. To improve the assembly of a previously sequenced lager yeast hybrid Saccharomyces sp. 790 and elucidate its genome structure, we proposed the use of alternative experimental evidence. We determined the phylogenetic position of Saccharomyces sp. 790 and established it as S. pastorianus 790. Then, we obtained from this yeast a bacterial artificial chromosome (BAC) genomic library with its BAC-end sequences (BESs). To analyze these data, we developed a pipeline (applicable to other assemblies) that classifies BES pairs alignments according to their orientation. For the case of S. pastorianus 790, paired-end BESs alignments validated parts of the assembly and unpaired-end ones suggested contig joins or misassemblies. Importantly, the BACs library was preserved and used for verification experiments. Unpaired-end alignments were used to upgrade the previous assembly and provided an improved detection of translocations. With this, we proposed a genome structure of S. pastorianus 790, which was similar to that of other lager yeasts; however, when we estimated chromosome copy number and experimentally measured its genome size, we discovered that one key difference is the outstanding S. pastorianus 790 ploidy level (allopentaploid). Altogether, our results show the value of combining bioinformatic analyses with experimental data such as long-insert clone information to improve a short-read assembly of a hybrid genome.

Discovery of Bivalent GalNAc-Conjugated Betulin as a Potent ASGPR-Directed Agent against Hepatocellular Carcinoma

Herein, we describe the design, synthesis, and biological evaluation of novel betulin and N-acetyl-d-galactosamine (GalNAc) glycoconjugates and suggest them as targeted agents against hepatocellular carcinoma. We prepared six conjugates derived via the C-3 and C-28 positions of betulin with one or two saccharide ligands. These molecules demonstrate high affinity to the asialoglycoprotein receptor (ASGPR) of hepatocytes assessed by in silico modeling and surface plasmon resonance tests. Cytotoxicity studies in vitro revealed a bivalent conjugate with moderate activity, selectivity of action, and cytostatic properties against hepatocellular carcinoma cells HepG2. An additional investigation confirmed the specific engagement with HepG2 cells by the enhanced generation of reactive oxygen species. Stability tests demonstrated its lability to acidic media and to intracellular enzymes. Therefore, the selected bivalent conjugate represents a new potential agent targeted against hepatocellular carcinoma. Further extensive studies of the cellular uptake in vitro and the real-time microdistribution in the murine liver in vivo for fluorescent dye-labeled analogue showed its selective internalization into hepatocytes due to the presence of GalNAc ligand in comparison with reference compounds. The betulin and GalNAc glycoconjugates can therefore be considered as a new strategy for developing therapeutic agents based on natural triterpenoids.

A Semi-automated Organoid Screening Method Demonstrates Epigenetic Control of Intestinal Epithelial Differentiation

Intestinal organoids are an excellent model to study epithelial biology. Yet, the selection of analytical tools to accurately quantify heterogeneous organoid cultures remains limited. Here, we developed a semi-automated organoid screening method, which we applied to a library of highly specific chemical probes to identify epigenetic regulators of intestinal epithelial biology. The role of epigenetic modifiers in adult stem cell systems, such as the intestinal epithelium, is still undefined. Based on this resource dataset, we identified several targets that affected epithelial cell differentiation, including HDACs, EP300/CREBBP, LSD1, and type I PRMTs, which were verified by complementary methods. For example, we show that inhibiting type I PRMTs, which leads enhanced epithelial differentiation, blocks the growth of adenoma but not normal organoid cultures. Thus, epigenetic probes are powerful tools to study intestinal epithelial biology and may have therapeutic potential.

Computational Analysis of Microbial Flow Cytometry Data

Flow cytometry is an important technology for the study of microbial communities. It grants the ability to rapidly generate phenotypic single-cell data that are both quantitative, multivariate and of high temporal resolution. The complexity and amount of data necessitate an objective and streamlined data processing workflow that extends beyond commercial instrument software. No full overview of the necessary steps regarding the computational analysis of microbial flow cytometry data currently exists. In this review, we provide an overview of the full data analysis pipeline, ranging from measurement to data interpretation, tailored toward studies in microbial ecology. At every step, we highlight computational methods that are potentially useful, for which we provide a short nontechnical description. We place this overview in the context of a number of open challenges to the field and offer further motivation for the use of standardized flow cytometry in microbial ecology research.

An R-based reproducible and user-friendly preprocessing pipeline for CyTOF data

Mass cytometry (CyTOF) has become a method of choice for in-depth characterization of tissue heterogeneity in health and disease, and is currently implemented in multiple clinical trials, where higher quality standards must be met. Currently, preprocessing of raw files is commonly performed in independent standalone tools, which makes it difficult to reproduce. Here, we present an R pipeline based on an updated version of CATALYST that covers all preprocessing steps required for downstream mass cytometry analysis in a fully reproducible way. This new version of CATALYST is based on Bioconductor's SingleCellExperiment class and fully unit tested. The R-based pipeline includes file concatenation, bead-based normalization, single-cell deconvolution, spillover compensation and live cell gating after debris and doublet removal. Importantly, this pipeline also includes different quality checks to assess machine sensitivity and staining performance while allowing also for batch correction. This pipeline is based on open source R packages and can be easily be adapted to different study designs. It therefore has the potential to significantly facilitate the work of CyTOF users while increasing the quality and reproducibility of data generated with this technology.

An Opaque Cell-Specific Expression Program of Secreted Proteases and Transporters Allows Cell-Type Cooperation in Candida albicans

An unusual feature of the opportunistic pathogen Candida albicans is its ability to switch stochastically between two distinct, heritable cell types called white and opaque. Here, we show that only opaque cells, in response to environmental signals, massively upregulate a specific group of secreted proteases and peptide transporters, allowing exceptionally efficient use of proteins as sources of nitrogen. We identify the specific proteases [members of the secreted aspartyl protease (SAP) family] needed for opaque cells to proliferate under these conditions, and we identify four transcriptional regulators of this specialized proteolysis and uptake program. We also show that, in mixed cultures, opaque cells enable white cells to also proliferate efficiently when proteins are the sole nitrogen source. Based on these observations, we suggest that one role of white-opaque switching is to create mixed populations where the different phenotypes derived from a single genome are shared between two distinct cell types.