Intratumoral radiation dose heterogeneity augments antitumor immunity in mice and primes responses to checkpoint blockade.

Radiation therapy (RT) activates multiple immunologic effects in the tumor microenvironment (TME), with diverse dose-response relationships observed. We hypothesized that, in contrast with homogeneous RT, a heterogeneous RT dose would simultaneously optimize activation of multiple immunogenic effects in a single TME, resulting in a more effective antitumor immune response. Using high-dose-rate brachytherapy, we treated mice bearing syngeneic tumors with a single fraction of heterogeneous RT at a dose ranging from 2 to 30 gray. When combined with dual immune checkpoint inhibition in murine models, heterogeneous RT generated more potent antitumor responses in distant, nonirradiated tumors compared with any homogeneous dose. The antitumor effect after heterogeneous RT required CD4 and CD8 T cells and low-dose RT to a portion of the tumor. At the 3-day post-RT time point, dose heterogeneity imprinted the targeted TME with spatial differences in immune-related gene expression, antigen presentation, and susceptibility of tumor cells to immune-mediated destruction. At a later 10-day post-RT time point, high-, moderate-, or low-RT-dose regions demonstrated distinct infiltrating immune cell populations. This was associated with an increase in the expression of effector-associated cytokines in circulating CD8 T cells. Consistent with enhanced adaptive immune priming, heterogeneous RT promoted clonal expansion of effector CD8 T cells. These findings illuminate the breadth of dose-dependent effects of RT on the TME and the capacity of heterogeneous RT to promote antitumor immunity when combined with immune checkpoint inhibitors.

Elevated GRHL2 Imparts Plasticity in ER-Positive Breast Cancer Cells.

Estrogen receptor (ER)-positive breast cancer is characterized by late recurrences following initial treatment. The epithelial cell fate transcription factor Grainyhead-like protein 2 (GRHL2) is overexpressed in ER-positive breast cancers and is linked to poorer prognosis as compared to ER-negative breast cancers. To understand how GRHL2 contributes to progression, GRHL2 was overexpressed in ER-positive cells. We demonstrated that elevated GRHL2 imparts plasticity with stem cell- and dormancy-associated traits. RNA sequencing and immunocytochemistry revealed that high GRHL2 not only strengthens the epithelial identity but supports a hybrid epithelial to mesenchymal transition (EMT). Proliferation and tumor studies exhibited a decrease in growth and an upregulation of dormancy markers, such as NR2F1 and CDKN1B. Mammosphere assays and flow cytometry revealed enrichment of stem cell markers CD44 and ALDH1, and increased self-renewal capacity. Cistrome analyses revealed a change in transcription factor motifs near GRHL2 sites from developmental factors to those associated with disease progression. Together, these data support the idea that the plasticity and properties induced by elevated GRHL2 may provide a selective advantage to explain the association between GRHL2 and breast cancer progression.

Assessing immune factors in maternal milk and paired infant plasma antibody binding to human rhinoviruses.

Introduction: Before they can produce their own antibodies, newborns are protected from infections by transplacental transfer of maternal IgG antibodies and after birth through breast milk IgA antibodies. Rhinovirus (RV) infections are extremely common in early childhood, and while RV infections often result in only mild upper respiratory illnesses, they can also cause severe lower respiratory illnesses such as bronchiolitis and pneumonia. Methods: We used high-density peptide arrays to profile infant and maternal antibody reactivity to capsid and full proteome sequences of three human RVs - A16, B52, and C11. Results: Numerous plasma IgG and breast milk IgA RV epitopes were identified that localized to regions of the RV capsid surface and interior, and also to several non-structural proteins. While most epitopes were bound by both IgG and IgA, there were several instances where isotype-specific and RV-specific binding were observed. We also profiled 62 unique RV-C protein loop sequences characteristic of this species' capsid VP1 protein. Discussion: Many of the RV-C loop sequences were highly bound by IgG from one-year-old infants, indicating recent or ongoing active infections, or alternatively, a level of cross-reactivity among homologous RV-C sites.

NK cells propagate T cell immunity following in situ tumor vaccination.

We report an in situ vaccination, adaptable to nearly any type of cancer, that combines radiotherapy targeting one tumor and intratumoral injection of this site with tumor-specific antibody and interleukin-2 (IL-2; 3xTx). In a phase I clinical trial, administration of 3xTx (with an immunocytokine fusion of tumor-specific antibody and IL-2, hu14.18-IL2) to subjects with metastatic melanoma increases peripheral CD8+ T cell effector polyfunctionality. This suggests the potential for 3xTx to promote antitumor immunity against metastatic tumors. In poorly immunogenic syngeneic murine melanoma or head and neck carcinoma models, 3xTx stimulates CD8+ T cell-mediated antitumor responses at targeted and non-targeted tumors. During 3xTx treatment, natural killer (NK) cells promote CTLA4+ regulatory T cell (Treg) apoptosis in non-targeted tumors. This is dependent on NK cell expression of CD86, which is upregulated downstream of KLRK1. NK cell depletion increases Treg infiltration, diminishing CD8+ T cell-dependent antitumor response. These findings demonstrate that NK cells sustain and propagate CD8+ T cell immunity following 3xTx.

Assessing Immune Factors in Maternal Milk and Paired Infant Plasma Antibody Binding to Human Rhinoviruses.

Before they can produce their own antibodies, newborns are protected from infections by transplacental transfer of maternal IgG antibodies and after birth through breast milk IgA antibodies. Rhinovirus (RV) infections are extremely common in early childhood, and while RV infections often result in only mild upper respiratory illnesses, they can also cause severe lower respiratory illnesses such as bronchiolitis and pneumonia. We used high-density peptide arrays to profile infant and maternal antibody reactivity to capsid and full proteome sequences of three human RVs - A16, B52, and C11. Numerous plasma IgG and breast milk IgA RV epitopes were identified that localized to regions of the RV capsid surface and interior, and also to several non-structural proteins. While most epitopes were bound by both IgG and IgA, there were several instances where isotype-specific and RV-specific binding were observed. We also profiled 62 unique RV-C dominant protein loop sequences characteristic of this species' capsid VP1 protein. Many of these RV-C sites were highly bound by IgG from one-year-old infants, indicating recent or ongoing active infections, or alternatively, a level of cross-reactivity among homologous RV-C sites.

A majority of Rhodobacter sphaeroides promoters lack a crucial RNA polymerase recognition feature, enabling coordinated transcription activation.

Using an in vitro transcription system with purified RNA polymerase (RNAP) to investigate rRNA synthesis in the photoheterotrophic α-proteobacterium Rhodobacter sphaeroides, we identified a surprising feature of promoters recognized by the major holoenzyme. Transcription from R. sphaeroides rRNA promoters was unexpectedly weak, correlating with absence of -7T, the very highly conserved thymine found at the last position in -10 elements of promoters in most bacterial species. Thymine substitutions for adenine at position -7 in the three rRNA promoters strongly increased intrinsic promoter activity, indicating that R. sphaeroides RNAP can utilize -7T when present. rRNA promoters were activated by purified R. sphaeroides CarD, a transcription factor found in many bacterial species but not in ß- and γ-proteobacteria. Overall, CarD increased the activity of 15 of 16 native R. sphaeroides promoters tested in vitro that lacked -7T, whereas it had no effect on three of the four native promoters that contained -7T. Genome-wide bioinformatic analysis of promoters from R. sphaeroides and two other α-proteobacterial species indicated that 30 to 43% contained -7T, whereas 90 to 99% of promoters from non-α-proteobacteria contained -7T. Thus, promoters lacking -7T appear to be widespread in α-proteobacteria and may have evolved away from consensus to enable their coordinated regulation by transcription factors like CarD. We observed a strong reduction in R. sphaeroides CarD levels when cells enter stationary phase, suggesting that reduced activation by CarD may contribute to inhibition of rRNA transcription when cells enter stationary phase, the stage of growth when bacterial ribosome synthesis declines.

Genome-Wide Identification of Transcription Start Sites in Two Alphaproteobacteria, Rhodobacter sphaeroides 2.4.1 and Novosphingobium aromaticivorans DSM 12444.

Here, we report the genome-wide identification of transcription start sites (TSSs) from two Alphaproteobacteria grown under conditions that result in significant changes in gene expression. TSSs that were identified as present in one condition or both will be an important resource for future studies of these, and possibly other, Alphaproteobacteria.

Heterologous expression of a glycosyl hydrolase and cellular reprogramming enable Zymomonas mobilis growth on cellobiose.

Plant-derived fuels and chemicals from renewable biomass have significant potential to replace reliance on petroleum and improve global carbon balance. However, plant biomass contains significant fractions of oligosaccharides that are not usable natively by many industrial microorganisms, including Escherichia coli, Saccharomyces cerevisiae, and Zymomonas mobilis. Even after chemical or enzymatic hydrolysis, some carbohydrate remains as non-metabolizable oligosaccharides (e.g., cellobiose or longer cellulose-derived oligomers), thus reducing the efficiency of conversion to useful products. To begin to address this problem for Z. mobilis, we engineered a strain (Z. mobilis GH3) that expresses a glycosyl hydrolase (GH) with ß-glucosidase activity from a related α-proteobacterial species, Caulobacter crescentus, and subjected it to an adaptation in cellobiose medium. Growth on cellobiose was achieved after a prolonged lag phase in cellobiose medium that induced changes in gene expression and cell composition, including increased expression and extracellular release of GH. These changes were reversible upon growth in glucose-containing medium, meaning they did not result from genetic mutation but could be retained upon transfer of cells to fresh cellobiose medium. After adaptation to cellobiose, our GH-expressing strain was able to convert about 50% of cellobiose to glucose within 24 h and use it for growth and ethanol production. Alternatively, pre-growth of Z. mobilis GH3 in sucrose medium enabled immediate growth on cellobiose. Proteomic analysis of cellobiose- and sucrose-adapted strains revealed upregulation of secretion-, transport-, and outer membrane-related proteins, which may aid release or surface display of GHs, entry of cellobiose into the periplasm, or both. Our two key findings are that Z. mobilis can be reprogrammed to grow on cellobiose as a sole carbon source and that this reprogramming is related to a natural response of Z. mobilis to sucrose that promotes sucrase production.

Genome-Scale Transcription-Translation Mapping Reveals Features of Zymomonas mobilis Transcription Units and Promoters.

Zymomonas mobilis is an ethanologenic alphaproteobacterium with promise for the industrial conversion of renewable plant biomass into fuels and chemical bioproducts. Limited functional annotation of the Z. mobilis genome is a current barrier to both fundamental studies of Z. mobilis and its development as a synthetic biology chassis. To gain insight, we collected sample-matched multiomics data, including RNA sequencing (RNA-seq), transcription start site (TSS) sequencing (TSS-seq), termination sequencing (term-seq), ribosome profiling, and label-free shotgun proteomic mass spectrometry, across different growth conditions and used these data to improve annotation and assign functional sites in the Z. mobilis genome. Proteomics and ribosome profiling informed revisions of protein-coding genes, which included 44 start codon changes and 42 added proteins. We developed statistical methods for annotating transcript 5' and 3' ends, enabling the identification of 3,940 TSSs and their corresponding promoters and 2,091 transcription termination sites, which were distinguished from RNA processing sites by the lack of an adjacent RNA 5' end. Our results revealed that Z. mobilis σA -35 and -10 promoter elements closely resemble canonical Escherichia coli -35 and -10 elements, with one notable exception: the Z. mobilis -10 element lacks the highly conserved -7 thymine observed in E. coli and other previously characterized σA promoters. The σA promoters of another alphaproteobacterium, Caulobacter crescentus, similarly lack the conservation of -7 thymine in their -10 elements. Our results anchor the development of Z. mobilis as a platform for synthetic biology and establish strategies for empirical genome annotation that can complement purely computational methods.IMPORTANCE Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Multiomic Fermentation Using Chemically Defined Synthetic Hydrolyzates Revealed Multiple Effects of Lignocellulose-Derived Inhibitors on Cell Physiology and Xylose Utilization in Zymomonas mobilis.

Utilization of both C5 and C6 sugars to produce biofuels and bioproducts is a key goal for the development of integrated lignocellulosic biorefineries. Previously we found that although engineered Zymomonas mobilis 2032 was able to ferment glucose to ethanol when fermenting highly concentrated hydrolyzates such as 9% glucan-loading AFEX-pretreated corn stover hydrolyzate (9% ACSH), xylose conversion after glucose depletion was greatly impaired. We hypothesized that impaired xylose conversion was caused by lignocellulose-derived inhibitors (LDIs) in hydrolyzates. To investigate the effects of LDIs on the cellular physiology of Z. mobilis during fermentation of hydrolyzates, including impacts on xylose utilization, we generated synthetic hydrolyzates (SynHs) that contained nutrients and LDIs at concentrations found in 9% ACSH. Comparative fermentations of Z. mobilis 2032 using SynH with or without LDIs were performed, and samples were collected for end product, transcriptomic, metabolomic, and proteomic analyses. Several LDI-specific effects were observed at various timepoints during fermentation including upregulation of sulfur assimilation and cysteine biosynthesis, upregulation of RND family efflux pump systems (ZMO0282-0285) and ZMO1429-1432, downregulation of a Type I secretion system (ZMO0252-0255), depletion of reduced glutathione, and intracellular accumulation of mannose-1P and mannose-6P. Furthermore, when grown in SynH containing LDIs, Z. mobilis 2032 only metabolized ∼50% of xylose, compared to ∼80% in SynH without LDIs, recapitulating the poor xylose utilization observed in 9% ACSH. Our metabolomic data suggest that the overall flux of xylose metabolism is reduced in the presence of LDIs. However, the expression of most genes involved in glucose and xylose assimilation was not affected by LDIs, nor did we observe blocks in glucose and xylose metabolic pathways. Accumulations of intracellular xylitol and xylonic acid was observed in both SynH with and without LDIs, which decreased overall xylose-to-ethanol conversion efficiency. Our results suggest that xylose metabolism in Z. mobilis 2032 may not be able to support the cellular demands of LDI mitigation and detoxification during fermentation of highly concentrated lignocellulosic hydrolyzates with elevated levels of LDIs. Together, our findings identify several cellular responses to LDIs and possible causes of impaired xylose conversion that will enable future strain engineering of Z. mobilis.

Complete genome sequence and the expression pattern of plasmids of the model ethanologen Zymomonas mobilis ZM4 and its xylose-utilizing derivatives 8b and 2032.

BACKGROUND: Zymomonas mobilis is a natural ethanologen being developed and deployed as an industrial biofuel producer. To date, eight Z. mobilis strains have been completely sequenced and found to contain 2-8 native plasmids. However, systematic verification of predicted Z. mobilis plasmid genes and their contribution to cell fitness has not been hitherto addressed. Moreover, the precise number and identities of plasmids in Z. mobilis model strain ZM4 have been unclear. The lack of functional information about plasmid genes in ZM4 impedes ongoing studies for this model biofuel-producing strain. RESULTS: In this study, we determined the complete chromosome and plasmid sequences of ZM4 and its engineered xylose-utilizing derivatives 2032 and 8b. Compared to previously published and revised ZM4 chromosome sequences, the ZM4 chromosome sequence reported here contains 65 nucleotide sequence variations as well as a 2400-bp insertion. Four plasmids were identified in all three strains, with 150 plasmid genes predicted in strain ZM4 and 2032, and 153 plasmid genes predicted in strain 8b due to the insertion of heterologous DNA for expanded substrate utilization. Plasmid genes were then annotated using Blast2GO, InterProScan, and systems biology data analyses, and most genes were found to have apparent orthologs in other organisms or identifiable conserved domains. To verify plasmid gene prediction, RNA-Seq was used to map transcripts and also compare relative gene expression under various growth conditions, including anaerobic and aerobic conditions, or growth in different concentrations of biomass hydrolysates. Overall, plasmid genes were more responsive to varying hydrolysate concentrations than to oxygen availability. Additionally, our results indicated that although all plasmids were present in low copy number (about 1-2 per cell), the copy number of some plasmids varied under specific growth conditions or due to heterologous gene insertion. CONCLUSIONS: The complete genome of ZM4 and two xylose-utilizing derivatives is reported in this study, with an emphasis on identifying and characterizing plasmid genes. Plasmid gene annotation, validation, expression levels at growth conditions of interest, and contribution to host fitness are reported for the first time.

Survey of cryptic unstable transcripts in yeast.

BACKGROUND: Cryptic unstable transcripts (CUTs) are a largely unexplored class of nuclear exosome degraded, non-coding RNAs in budding yeast. It is highly debated whether CUT transcription has a functional role in the cell or whether CUTs represent noise in the yeast transcriptome. We sought to ascertain the extent of conserved CUT expression across a variety of Saccharomyces yeast strains to further understand and characterize the nature of CUT expression. RESULTS: We sequenced the WT and rrp6Δ transcriptomes of three S.cerevisiae strains: S288c, Σ1278b, JAY291 and the S.paradoxus strain N17 and utilized a hidden Markov model to annotate CUTs in these four strains. Utilizing a four-way genomic alignment we identified a large population of CUTs with conserved syntenic expression across all four strains. By identifying configurations of gene-CUT pairs, where CUT expression originates from the gene 5' or 3' nucleosome free region, we observed distinct gene expression trends specific to these configurations which were most prevalent in the presence of conserved CUT expression. Divergent pairs correlate with higher expression of genes, and convergent pairs correlate with reduced gene expression. CONCLUSIONS: Our RNA-seq based method has greatly expanded upon previous CUT annotations in S.cerevisiae underscoring the extensive and pervasive nature of unstable transcription. Furthermore we provide the first assessment of conserved CUT expression in yeast and globally demonstrate possible modes of CUT-based regulation of gene expression.