A multi-omic single-cell landscape of human gynecologic malignancies.

Deconvolution of regulatory mechanisms that drive transcriptional programs in cancer cells is key to understanding tumor biology. Herein, we present matched transcriptome (scRNA-seq) and chromatin accessibility (scATAC-seq) profiles at single-cell resolution from human ovarian and endometrial tumors processed immediately following surgical resection. This dataset reveals the complex cellular heterogeneity of these tumors and enabled us to quantitatively link variation in chromatin accessibility to gene expression. We show that malignant cells acquire previously unannotated regulatory elements to drive hallmark cancer pathways. Moreover, malignant cells from within the same patients show substantial variation in chromatin accessibility linked to transcriptional output, highlighting the importance of intratumoral heterogeneity. Finally, we infer the malignant cell type-specific activity of transcription factors. By defining the regulatory logic of cancer cells, this work reveals an important reliance on oncogenic regulatory elements and highlights the ability of matched scRNA-seq/scATAC-seq to uncover clinically relevant mechanisms of tumorigenesis in gynecologic cancers.

Single-Cell Transcriptional and Epigenetic Profiles of Male Breast Cancer Nominate Salient Cancer-Specific Enhancers.

Male breast cancer represents about 1% of all breast cancer diagnoses and, although there are some similarities between male and female breast cancer, the paucity of data available on male breast cancer makes it difficult to establish targeted therapies. To date, most male breast cancers (MBCs) are treated according to protocols established for female breast cancer (FBC). Thus, defining the transcriptional and epigenetic landscape of MBC with improved resolution is critical for developing better avenues for therapeutic intervention. In this study, we present matched transcriptional (scRNA-seq) and epigenetic (scATAC-seq) profiles at single-cell resolution of two treatment naïve MBC tumors processed immediately after surgical resection. These data enable the detection of differentially expressed genes between male and female breast tumors across immune, stromal, and malignant cell types, to highlight several genes that may have therapeutic implications. Notably, MYC target genes and mTORC1 signaling genes were significantly upregulated in the malignant cells of MBC compared to the female counterparts. To understand how the regulatory landscape of MBC gives rise to these male-specific gene expression patterns, we leveraged the scATAC-seq data to systematically link changes in chromatin accessibility to changes in gene expression within each cell type. We observed cancer-specific rewiring of several salient enhancers and posit that these enhancers have a higher regulatory load than lineage-specific enhancers. We highlight two examples of previously unannotated cancer-cell-specific enhancers of ANXA2 and PRDX4 gene expression and show evidence for super-enhancer regulation of LAMB3 and CD47 in male breast cancer cells. Overall, this dataset annotates clinically relevant regulatory networks in male breast tumors, providing a useful resource that expands our current understanding of the gene expression programs that underlie the biology of MBC.

CRISPR-Cas9 editing of CAFFEOYL SHIKIMATE ESTERASE 1 and 2 shows their importance and partial redundancy in lignification in Populus tremula × P. alba.

Lignins are cell wall-located aromatic polymers that provide strength and hydrophobicity to woody tissues. Lignin monomers are synthesized via the phenylpropanoid pathway, wherein CAFFEOYL SHIKIMATE ESTERASE (CSE) converts caffeoyl shikimate into caffeic acid. Here, we explored the role of the two CSE homologs in poplar (Populus tremula × P. alba). Reporter lines showed that the expression conferred by both CSE1 and CSE2 promoters is similar. CRISPR-Cas9-generated cse1 and cse2 single mutants had a wild-type lignin level. Nevertheless, CSE1 and CSE2 are not completely redundant, as both single mutants accumulated caffeoyl shikimate. In contrast, the cse1 cse2 double mutants had a 35% reduction in lignin and associated growth penalty. The reduced-lignin content translated into a fourfold increase in cellulose-to-glucose conversion upon limited saccharification. Phenolic profiling of the double mutants revealed large metabolic shifts, including an accumulation of p-coumaroyl, 5-hydroxyferuloyl, feruloyl and sinapoyl shikimate, in addition to caffeoyl shikimate. This indicates that the CSEs have a broad substrate specificity, which was confirmed by in vitro enzyme kinetics. Taken together, our results suggest an alternative path within the phenylpropanoid pathway at the level of the hydroxycinnamoyl-shikimates, and show that CSE is a promising target to improve plants for the biorefinery.

2H and 13C metabolic flux analysis elucidates in vivo thermodynamics of the ED pathway in Zymomonas mobilis.

Zymomonas mobilis is an industrially relevant bacterium notable for its ability to rapidly ferment simple sugars to ethanol using the Entner-Doudoroff (ED) glycolytic pathway, an alternative to the well-known Embden-Meyerhof-Parnas (EMP) pathway used by most organisms. Recent computational studies have predicted that the ED pathway is substantially more thermodynamically favorable than the EMP pathway, a potential factor explaining the high glycolytic rate in Z. mobilis. Here, to investigate the in vivo thermodynamics of the ED pathway and central carbon metabolism in Z. mobilis, we implemented a network-level approach that integrates quantitative metabolomics with 2H and 13C metabolic flux analysis to estimate reversibility and Gibbs free energy (ΔG) of metabolic reactions. This analysis revealed a strongly thermodynamically favorable ED pathway in Z. mobilis that is nearly twice as favorable as the EMP pathway in E. coli or S. cerevisiae. The in vivo step-by-step thermodynamic profile of the ED pathway was highly similar to previous in silico predictions, indicating that maximizing ΔG for each pathway step likely constitutes a cellular objective in Z. mobilis. Our analysis also revealed novel features of Z. mobilis metabolism, including phosphofructokinase-like enzyme activity, tricarboxylic acid cycle anaplerosis via PEP carboxylase, and a metabolic imbalance in the pentose phosphate pathway resulting in excretion of shikimate pathway intermediates. The integrated approach we present here for in vivo ΔG quantitation may be applied to the thermodynamic profiling of pathways and metabolic networks in other microorganisms and will contribute to the development of quantitative models of metabolism.

Roles of Water Molecules and Counterion on HS- Sensing Reaction Utilizing a Pyrylium Derivative: A Computational Study.

In this paper, we present a comprehensive computational study on the hydrogen sulfide ion (HS-) sensing mechanism in aqueous solution using pyrylium-thiopyrylium transformation. Explicit water molecules up to three water molecules are considered using supramolecular models. The effect of water bulk solvent is also taken into account according to the polarizable continuum model. Our results demonstrate that water molecules are directly involved in the sensing reactions by altering reaction mechanisms and dramatically lower the activation energies. The most favorable HS- sensing mechanism involves a 10-membered ring transition structure formed by three water molecules and one hydronium. The catalytic effects of water molecule(s) due to the alleviation of ring strain and the stabilization from deprotonated hydronium significantly lower the activation energy. The activation energies in aqueous solution decrease from 40.2 kcal/mol for the hydronium-only-catalyzed reaction to 15.7, 14.8, and 7.4 kcal/mol for one-water-, two-water-, and three-water-catalyzed mechanisms, respectively. In addition, the effect of the counterion tetrafluoroborate (BF4-) on the reaction mechanisms was also investigated. Our results demonstrate that the counterion BF4- most likely behaves as a spectator and has minor influence on the reaction mechanism.

Tricin, a flavonoid monomer in monocot lignification.

Tricin was recently discovered in lignin preparations from wheat (Triticum aestivum) straw and subsequently in all monocot samples examined. To provide proof that tricin is involved in lignification and establish the mechanism by which it incorporates into the lignin polymer, the 4'-O-ß-coupling products of tricin with the monolignols (p-coumaryl, coniferyl, and sinapyl alcohols) were synthesized along with the trimer that would result from its 4'-O-ß-coupling with sinapyl alcohol and then coniferyl alcohol. Tricin was also found to cross couple with monolignols to form tricin-(4'-O-ß)-linked dimers in biomimetic oxidations using peroxidase/hydrogen peroxide or silver (I) oxide. Nuclear magnetic resonance characterization of gel permeation chromatography-fractionated acetylated maize (Zea mays) lignin revealed that the tricin moieties are found in even the highest molecular weight fractions, ether linked to lignin units, demonstrating that tricin is indeed incorporated into the lignin polymer. These findings suggest that tricin is fully compatible with lignification reactions, is an authentic lignin monomer, and, because it can only start a lignin chain, functions as a nucleation site for lignification in monocots. This initiation role helps resolve a long-standing dilemma that monocot lignin chains do not appear to be initiated by monolignol homodehydrodimerization as they are in dicots that have similar syringyl-guaiacyl compositions. The term flavonolignin is recommended for the racemic oligomers and polymers of monolignols that start from tricin (or incorporate other flavonoids) in the cell wall, in analogy with the existing term flavonolignan that is used for the low-molecular mass compounds composed of flavonoid and lignan moieties.

Defining the Regulatory Logic of Breast Cancer Using Single-Cell Epigenetic and Transcriptome Profiling.

Annotation of the cis-regulatory elements that drive transcriptional dysregulation in cancer cells is critical to improving our understanding of tumor biology. Herein, we present a compendium of matched chromatin accessibility (scATAC-seq) and transcriptome (scRNA-seq) profiles at single-cell resolution from human breast tumors and healthy mammary tissues processed immediately following surgical resection. We identify the most likely cell-of-origin for luminal breast tumors and basal breast tumors and then introduce a novel methodology that implements linear mixed-effects models to systematically quantify associations between regions of chromatin accessibility (i.e. regulatory elements) and gene expression in malignant cells versus normal mammary epithelial cells. These data unveil regulatory elements with that switch from silencers of gene expression in normal cells to enhancers of gene expression in cancer cells, leading to the upregulation of clinically relevant oncogenes. To translate the utility of this dataset into tractable models, we generated matched scATAC-seq and scRNA-seq profiles for breast cancer cell lines, revealing, for each subtype, a conserved oncogenic gene expression program between in vitro and in vivo cells. Together, this work highlights the importance of non-coding regulatory mechanisms that underlie oncogenic processes and the ability of single-cell multi-omics to define the regulatory logic of BC cells at single-cell resolution.

Enhancer RNA Transcription Is Essential for a Novel CSF1 Enhancer in Triple-Negative Breast Cancer.

Enhancers are critical regulatory elements in the genome that help orchestrate spatiotemporal patterns of gene expression during development and normal physiology. In cancer, enhancers are often rewired by various genetic and epigenetic mechanisms for the activation of oncogenes that lead to initiation and progression. A key feature of active enhancers is the production of non-coding RNA molecules called enhancer RNAs, whose functions remain unknown but can be used to specify active enhancers de novo. Using a combination of eRNA transcription and chromatin modifications, we have identified a novel enhancer located 30 kb upstream of Colony Stimulating Factor 1 (CSF1). Notably, CSF1 is implicated in the progression of breast cancer, is overexpressed in triple-negative breast cancer (TNBC) cell lines, and its enhancer is primarily active in TNBC patient tumors. Genomic deletion of the enhancer (via CRISPR/Cas9) enabled us to validate this regulatory element as a bona fide enhancer of CSF1 and subsequent cell-based assays revealed profound effects on cancer cell proliferation, colony formation, and migration. Epigenetic silencing of the enhancer via CRISPR-interference assays (dCas9-KRAB) coupled to RNA-sequencing, enabled unbiased identification of additional target genes, such as RSAD2, that are predictive of clinical outcome. Additionally, we repurposed the RNA-guided RNA-targeting CRISPR-Cas13 machinery to specifically degrade the eRNAs transcripts produced at this enhancer to determine the consequences on CSF1 mRNA expression, suggesting a post-transcriptional role for these non-coding transcripts. Finally, we test our eRNA-dependent model of CSF1 enhancer function and demonstrate that our results are extensible to other forms of cancer. Collectively, this work describes a novel enhancer that is active in the TNBC subtype, which is associated with cellular growth, and requires eRNA transcripts for proper enhancer function. These results demonstrate the significant impact of enhancers in cancer biology and highlight their potential as tractable targets for therapeutic intervention.

A multi-omic dissection of super-enhancer driven oncogenic gene expression programs in ovarian cancer.

The human genome contains regulatory elements, such as enhancers, that are often rewired by cancer cells for the activation of genes that promote tumorigenesis and resistance to therapy. This is especially true for cancers that have little or no known driver mutations within protein coding genes, such as ovarian cancer. Herein, we utilize an integrated set of genomic and epigenomic datasets to identify clinically relevant super-enhancers that are preferentially amplified in ovarian cancer patients. We systematically probe the top 86 super-enhancers, using CRISPR-interference and CRISPR-deletion assays coupled to RNA-sequencing, to nominate two salient super-enhancers that drive proliferation and migration of cancer cells. Utilizing Hi-C, we construct chromatin interaction maps that enable the annotation of direct target genes for these super-enhancers and confirm their activity specifically within the cancer cell compartment of human tumors using single-cell genomics data. Together, our multi-omic approach examines a number of fundamental questions about how regulatory information encoded into super-enhancers drives gene expression networks that underlie the biology of ovarian cancer.

Systems-Level Analysis of Oxygen Exposure in Zymomonas mobilis: Implications for Isoprenoid Production.

Zymomonas mobilis is an aerotolerant anaerobe and prolific ethanologen with attractive characteristics for industrial bioproduct generation. However, there is currently insufficient knowledge of the impact that environmental factors have on flux through industrially relevant biosynthetic pathways. Here, we examined the effect of oxygen exposure on metabolism and gene expression in Z. mobilis by combining targeted metabolomics, mRNA sequencing, and shotgun proteomics. We found that exposure to oxygen profoundly influenced metabolism, inducing both transient metabolic bottlenecks and long-term metabolic remodeling. In particular, oxygen induced a severe but temporary metabolic bottleneck in the methyl erythritol 4-phosphate pathway for isoprenoid biosynthesis caused by oxidative damage to the iron-sulfur cofactors of the final two enzymes in the pathway. This bottleneck was resolved with minimal changes in expression of isoprenoid biosynthetic enzymes. Instead, it was associated with pronounced upregulation of enzymes related to iron-sulfur cluster maintenance and biogenesis (i.e., flavodoxin reductase and the suf operon). We also detected major changes in glucose utilization in the presence of oxygen. Specifically, we observed increased gluconate production following exposure to oxygen, accounting for 18% of glucose uptake. Our results suggest that under aerobic conditions, electrons derived from the oxidation of glucose to gluconate are diverted to the electron transport chain, where they can minimize oxidative damage by reducing reactive oxygen species such as H2O2. This model is supported by the simultaneous upregulation of three membrane-bound dehydrogenases, cytochrome c peroxidase, and a cytochrome bd oxidase following exposure to oxygen. IMPORTANCE Microbially generated biofuels and bioproducts have the potential to provide a more environmentally sustainable alternative to fossil-fuel-derived products. In particular, isoprenoids, a diverse class of natural products, are chemically suitable for use as high-grade transport fuels and other commodity molecules. However, metabolic engineering for increased production of isoprenoids and other bioproducts is limited by an incomplete understanding of factors that control flux through biosynthetic pathways. Here, we examined the native regulation of the isoprenoid biosynthetic pathway in the biofuel producer Zymomonas mobilis. We leveraged oxygen exposure as a means to perturb carbon flux, allowing us to observe the formation and resolution of a metabolic bottleneck in the pathway. Our multi-omics analysis of this perturbation enabled us to identify key auxiliary enzymes whose expression correlates with increased production of isoprenoid precursors, which we propose as potential targets for future metabolic engineering.

Enhanced delignification of steam-pretreated poplar by a bacterial laccase.

The recalcitrance of woody biomass, particularly its lignin component, hinders its sustainable transformation to fuels and biomaterials. Although the recent discovery of several bacterial ligninases promises the development of novel biocatalysts, these enzymes have largely been characterized using model substrates: direct evidence for their action on biomass is lacking. Herein, we report the delignification of woody biomass by a small laccase (sLac) from Amycolatopsis sp. 75iv3. Incubation of steam-pretreated poplar (SPP) with sLac enhanced the release of acid-precipitable polymeric lignin (APPL) by ~6-fold, and reduced the amount of acid-soluble lignin by ~15%. NMR spectrometry revealed that the APPL was significantly syringyl-enriched relative to the original material (~16:1 vs. ~3:1), and that sLac preferentially oxidized syringyl units and altered interunit linkage distributions. sLac's substrate preference among monoaryls was also consistent with this observation. In addition, sLac treatment reduced the molar mass of the APPL by over 50%, as determined by gel-permeation chromatography coupled with multi-angle light scattering. Finally, sLac acted synergistically with a commercial cellulase cocktail to increase glucose production from SPP ~8%. Overall, this study establishes the lignolytic activity of sLac on woody biomass and highlights the biocatalytic potential of bacterial enzymes.

Ionic liquids with thioether motifs as synthetic cationic lipids for gene delivery.

This study introduces a novel class of imidazolium- and ammonium-based ionic liquids possessing two C12 and C14 tails and thioether linkers designed for lipoplex-mediated DNA delivery. Imidazolium-based ionic liquids displayed efficient gene delivery properties with low toxicity. Thiol-yne click chemistry was employed for the facile and robust synthesis of these thioether-based cationic lipioids with enhanced lipophilicity and low fluidity.