Heat stress affects epigenetic gene silencing in Arabidopsis. To test for a mechanistic involvement of epigenetic regulation in heat-stress responses, we analyzed the heat tolerance of mutants defective in DNA methylation, histone modifications, chromatin-remodeling, or siRNA-based silencing pathways. Plants deficient in NRPD2, the common second-largest subunit of RNA polymerases IV and V, and in the Rpd3-type histone deacetylase HDA6 were hypersensitive to heat exposure. Microarray analysis demonstrated that NRPD2 and HDA6 have independent roles in transcriptional reprogramming in response to temperature stress. The misexpression of protein-coding genes in nrpd2 mutants recovering from heat correlated with defective epigenetic regulation of adjacent transposon remnants which involved the loss of control of heat-stress-induced read-through transcription. We provide evidence that the transcriptional response to temperature stress, at least partially, relies on the integrity of the RNA-dependent DNA methylation pathway.

DNA methylation loss occurs frequently in cancer genomes, primarily within lamina-associated, late-replicating regions termed partially methylated domains (PMDs). We profiled 39 diverse primary tumors and 8 matched adjacent tissues using whole-genome bisulfite sequencing (WGBS) and analyzed them alongside 343 additional human and 206 mouse WGBS datasets. We identified a local CpG sequence context associated with preferential hypomethylation in PMDs. Analysis of CpGs in this context ('solo-WCGWs') identified previously undetected PMD hypomethylation in almost all healthy tissue types. PMD hypomethylation increased with age, beginning during fetal development, and appeared to track the accumulation of cell divisions. In cancer, PMD hypomethylation depth correlated with somatic mutation density and cell cycle gene expression, consistent with its reflection of mitotic history and suggesting its application as a mitotic clock. We propose that late replication leads to lifelong progressive methylation loss, which acts as a biomarker for cellular aging and which may contribute to oncogenesis.

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection can result in severe immune dysfunction, hospitalization, and death. Many patients also develop long-COVID-19, experiencing symptoms months after infection. Although significant progress has been made in understanding the immune response to acute SARS-CoV-2 infection, gaps remain in our knowledge of how innate immunity influences disease kinetics and severity. We hypothesized that cytometry by time-of-flight analysis of PBMCs from healthy and infected subjects would identify novel cell surface markers and innate immune cell subsets associated with COVID-19 severity. In this pursuit, we identified monocyte and dendritic cell subsets that changed in frequency during acute SARS-CoV-2 infection and correlated with clinical parameters of disease severity. Subsets of nonclassical monocytes decreased in frequency in hospitalized subjects, yet increased in the most severe patients and positively correlated with clinical values associated with worse disease severity. CD9, CD163, PDL1, and PDL2 expression significantly increased in hospitalized subjects, and CD9 and 6-Sulfo LacNac emerged as the markers that best distinguished monocyte subsets amongst all subjects. CD9+ monocytes remained elevated, whereas nonclassical monocytes remained decreased, in the blood of hospitalized subjects at 3-4 months postinfection. Finally, we found that CD9+ monocytes functionally released more IL-8 and MCP-1 after LPS stimulation. This study identifies new monocyte subsets present in the blood of COVID-19 patients that correlate with disease severity, and links CD9+ monocytes to COVID-19 progression.

CD8+ T cells drive anti-cancer immunity in response to antigen-presenting cells such as dendritic cells and subpopulations of monocytes and macrophages. While CD14+ classical monocytes modulate CD8+ T cell responses, the contributions of CD16+ nonclassical monocytes to this process remain unclear. Herein we explored the role of nonclassical monocytes in CD8+ T cell activation by utilizing E2-deficient (E2-/-) mice that lack nonclassical monocytes. During early metastatic seeding, modeled by B16F10-OVA cancer cells injected into E2-/- mice, we noted lower CD8+ effector memory and effector T cell frequencies within the lungs as well as in lung-draining mediastinal lymph nodes in the E2-/- mice. Analysis of the myeloid compartment revealed that these changes were associated with depletion of MHC-IIloLy6Clo nonclassical monocytes within these tissues, with little change in other monocyte or macrophage populations. Additionally, nonclassical monocytes preferentially trafficked to primary tumor sites in the lungs, rather than to the lung-draining lymph nodes, and did not cross-present antigen to CD8+ T cells. Examination of the lung microenvironment in E2-/- mice revealed reduced CCL21 expression in endothelial cells, which is chemokine involved in T cell trafficking. Our results highlight the previously unappreciated importance of nonclassical monocytes in shaping the tumor microenvironment via CCL21 production and CD8+ T cell recruitment.

Most high-grade serous ovarian cancers (HGSCs) likely initiate from fallopian tube (FT) epithelia. While epithelial subtypes have been characterized using single-cell RNA-sequencing (scRNA-Seq), heterogeneity of other compartments and their involvement in tumor progression are poorly defined. Integrated analysis of human FT scRNA-Seq and HGSC-related tissues, including tumors, revealed greater immune and stromal transcriptional diversity than previously reported. We identified abundant monocytes in FTs across two independent cohorts. The ratio of macrophages to monocytes is similar between benign FTs, ovaries, and adjacent normal tissues but significantly greater in tumors. FT-defined monocyte and macrophage signatures, cell-cell communication, and gene set enrichment analyses identified monocyte- and macrophage-specific interactions and functional pathways in paired tumors and adjacent normal tissues. Further reanalysis of HGSC scRNA-Seq identified monocyte and macrophage subsets associated with neoadjuvant chemotherapy. Taken together, our work provides data that an altered FT myeloid cell composition could inform the discovery of early detection markers for HGSC.

Neutrophils are short-lived cells that play important roles in both health and disease. Neutrophils and monocytes originate from the granulocyte monocyte progenitor (GMP) in bone marrow; however, unipotent neutrophil progenitors are not well defined. Here, we use cytometry by time of flight (CyTOF) and single-cell RNA sequencing (scRNA-seq) methodologies to identify a committed unipotent early-stage neutrophil progenitor (NeP) in adult mouse bone marrow. Importantly, we found a similar unipotent NeP (hNeP) in human bone marrow. Both NeP and hNeP generate only neutrophils. NeP and hNeP both significantly increase tumor growth when transferred into murine cancer models, including a humanized mouse model. hNeP are present in the blood of treatment-naive melanoma patients but not of healthy subjects. hNeP can be readily identified by flow cytometry and could be used as a biomarker for early cancer discovery. Understanding the biology of hNeP should allow the development of new therapeutic targets for neutrophil-related diseases, including cancer.

Neutrophils are the most abundant peripheral immune cells and thus, are continually replenished by bone marrow-derived progenitors. Still, how newly identified neutrophil subsets fit into the bone marrow neutrophil lineage remains unclear. Here, we use mass cytometry to show that two recently defined human neutrophil progenitor populations contain a homogeneous progenitor subset we term "early neutrophil progenitors" (eNePs) (Lin-CD66b+CD117+CD71+). Surface marker- and RNA-expression analyses, together with in vitro colony formation and in vivo adoptive humanized mouse transfers, indicate that eNePs are the earliest human neutrophil progenitors. Furthermore, we identified CD71 as a marker associated with the earliest neutrophil developmental stages. Expression of CD71 marks proliferating neutrophils, which were expanded in the blood of melanoma patients and detectable in blood and tumors from lung cancer patients. In summary, we establish CD117+CD71+ eNeP as the inceptive human neutrophil progenitor and propose a refined model of the neutrophil developmental lineage in bone marrow.

Fallopian tube secretory epithelial cells (FTSECs) are likely the main precursor cell type of high-grade serous ovarian cancers (HGSOCs), but these tumors may also arise from ovarian surface epithelial cells (OSECs). We profiled global landscapes of gene expression and active chromatin to characterize molecular similarities between OSECs (n = 114), FTSECs (n = 74), and HGSOCs (n = 394). A one-class machine learning algorithm predicts that most HGSOCs derive from FTSECs, with particularly high FTSEC scores in mesenchymal-type HGSOCs (padj < 8 × 10-4). However, a subset of HGSOCs likely derive from OSECs, particularly HGSOCs of the proliferative type (padj < 2 × 10-4), suggesting a dualistic model for HGSOC origins. Super-enhancer (SE) landscapes were also more similar between FTSECs and HGSOCs than between OSECs and HGSOCs (p < 2.2 × 10-16). The SOX18 transcription factor (TF) coincided with a HGSOC-specific SE, and ectopic overexpression of SOX18 in FTSECs caused epithelial-to-mesenchymal transition, indicating that SOX18 plays a role in establishing the mesenchymal signature of fallopian-derived HGSOCs.

We investigated whether in human head and neck squamous cell carcinoma (HNSCC) high levels of expression of stress keratin 17 (K17) are associated with poor survival and resistance to immunotherapy.

Complement factor H (CFH) negatively regulates consumption of complement component 3 (C3), thereby restricting complement activation. Genetic variants in CFH predispose to chronic inflammatory disease. Here, we examined the impact of CFH on atherosclerosis development. In a mouse model of atherosclerosis, CFH deficiency limited plaque necrosis in a C3-dependent manner. Deletion of CFH in monocyte-derived inflammatory macrophages propagated uncontrolled cell-autonomous C3 consumption without downstream C5 activation and heightened efferocytotic capacity. Among leukocytes, Cfh expression was restricted to monocytes and macrophages, increased during inflammation, and coincided with the accumulation of intracellular C3. Macrophage-derived CFH was sufficient to dampen resolution of inflammation, and hematopoietic deletion of CFH in atherosclerosis-prone mice promoted lesional efferocytosis and reduced plaque size. Furthermore, we identified monocyte-derived inflammatory macrophages expressing C3 and CFH in human atherosclerotic plaques. Our findings reveal a regulatory axis wherein CFH controls intracellular C3 levels of macrophages in a cell-autonomous manner, evidencing the importance of on-site complement regulation in the pathogenesis of inflammatory diseases.

Esophageal squamous cell carcinoma (ESCC) is among the most common malignancies, but little is known about its spatial intratumoral heterogeneity (ITH) and temporal clonal evolutionary processes. To address this, we performed multiregion whole-exome sequencing on 51 tumor regions from 13 ESCC cases and multiregion global methylation profiling for 3 of these 13 cases. We found an average of 35.8% heterogeneous somatic mutations with strong evidence of ITH. Half of the driver mutations located on the branches of tumor phylogenetic trees targeted oncogenes, including PIK3CA, NFE2L2 and MTOR, among others. By contrast, the majority of truncal and clonal driver mutations occurred in tumor-suppressor genes, including TP53, KMT2D and ZNF750, among others. Interestingly, phyloepigenetic trees robustly recapitulated the topological structures of the phylogenetic trees, indicating a possible relationship between genetic and epigenetic alterations. Our integrated investigations of spatial ITH and clonal evolution provide an important molecular foundation for enhanced understanding of tumorigenesis and progression in ESCC.

The regulation of gene expression is mediated at the transcriptional level by enhancer regions that are bound by sequence-specific transcription factors (TFs). Recent studies have shown that the in vivo binding sites of single TFs differ between developmental or cellular contexts. How this context-specific binding is encoded in the cis-regulatory DNA sequence has, however, remained unclear. We computationally dissect context-specific TF binding sites in Drosophila, Caenorhabditis elegans, mouse, and human and find distinct combinations of sequence motifs for partner factors, which are predictive and reveal specific motif requirements of individual binding sites. We predict that TF binding in the early Drosophila embryo depends on motifs for the early zygotic TFs Vielfaltig (also known as Zelda) and Tramtrack. We validate experimentally that the activity of Twist-bound enhancers and Twist binding itself depend on Vielfaltig motifs, suggesting that Vielfaltig is more generally important for early transcription. Our finding that the motif content can predict context-specific binding and that the predictions work across different Drosophila species suggests that characteristic motif combinations are shared between sites, revealing context-specific motif codes (cis-regulatory signatures), which appear to be conserved during evolution. Taken together, this study establishes a novel approach to derive predictive cis-regulatory motif requirements for individual TF binding sites and enhancers. Importantly, the method is generally applicable across different cell types and organisms to elucidate cis-regulatory sequence determinants and the corresponding trans-acting factors from the increasing number of tissue- and cell-type-specific TF binding studies.

The human fallopian tube harbors the cell of origin for the majority of high-grade serous "ovarian" cancers (HGSCs), but its cellular composition, particularly the epithelial component, is poorly characterized. We perform single-cell transcriptomic profiling of around 53,000 individual cells from 12 primary fallopian specimens to map their major cell types. We identify 10 epithelial subpopulations with diverse transcriptional programs. Based on transcriptional signatures, we reconstruct a trajectory whereby secretory cells differentiate into ciliated cells via a RUNX3high intermediate. Computational deconvolution of advanced HGSCs identifies the "early secretory" population as a likely precursor state for the majority of HGSCs. Its signature comprises both epithelial and mesenchymal features and is enriched in mesenchymal-type HGSCs (p = 6.7 × 10-27), a group known to have particularly poor prognoses. This cellular and molecular compendium of the human fallopian tube in cancer-free women is expected to advance our understanding of the earliest stages of fallopian epithelial neoplasia.

Homeostatic trafficking to lymph nodes allows T cells to efficiently survey the host for cognate antigen. Nonmammalian jawed vertebrates lack lymph nodes but maintain diverse T cell pools. Here, we exploit in vivo imaging of transparent zebrafish to investigate how T cells organize and survey for antigen in an animal devoid of lymph nodes. We find that naïve-like T cells in zebrafish organize into a previously undescribed whole-body lymphoid network that supports streaming migration and coordinated trafficking through the host. This network has the cellular hallmarks of a mammalian lymph node, including naïve T cells and CCR7-ligand expressing nonhematopoietic cells, and facilitates rapid collective migration. During infection, T cells transition to a random walk that supports antigen-presenting cell interactions and subsequent activation. Our results reveal that T cells can toggle between collective migration and individual random walks to prioritize either large-scale trafficking or antigen search in situ. This lymphoid network thus facilitates whole-body T cell trafficking and antigen surveillance in the absence of a lymph node system.

Deep sequencing after bisulfite conversion (BS-Seq) is the method of choice to generate whole genome maps of cytosine methylation at single base-pair resolution. Its application to genomic DNA of Arabidopsis flower bud tissue resulted in the first complete methylome, determining a methylation rate of 6.7% in this tissue. BS-Seq reads were mapped onto an in silico converted reference genome, applying the so-called 3-letter genome method. Here, we present BiSS (Bisufite Sequencing Scorer), a new method applying Smith-Waterman alignment to map bisulfite-converted reads to a reference genome. In addition, we introduce a comprehensive adaptive error estimate that accounts for sequencing errors, erroneous bisulfite conversion and also wrongly mapped reads. The re-analysis of the Arabidopsis methylome data with BiSS mapped substantially more reads to the genome. As a result, it determines the methylation status of an extra 10% of cytosines and estimates the methylation rate to be 7.7%. We validated the results by individual traditional bisulfite sequencing for selected genomic regions. In addition to predicting the methylation status of each cytosine, BiSS also provides an estimate of the methylation degree at each genomic site. Thus, BiSS explores BS-Seq data more extensively and provides more information for downstream analysis.

Humans are infected with two types of EBV (Type 1 (T1) and Type 2 (T2)) that differ substantially in their EBNA2 and EBNA 3A/B/C latency proteins and have different phenotypes in B cells. T1 EBV transforms B cells more efficiently than T2 EBV in vitro, and T2 EBV-infected B cells are more lytic. We previously showed that both increased NFATc1/c2 activity, and an NFAT-binding motif within the BZLF1 immediate-early promoter variant (Zp-V3) contained in all T2 strains, contribute to lytic infection in T2 EBV-infected B cells. Here we compare cellular and viral gene expression in early-passage lymphoblastoid cell lines (LCLs) infected with either T1 or T2 EBV strains. Using bulk RNA-seq, we show that T2 LCLs are readily distinguishable from T1 LCLs, with approximately 600 differentially expressed cellular genes. Gene Set Enrichment Analysis (GSEA) suggests that T2 LCLs have increased B-cell receptor (BCR) signaling, NFAT activation, and enhanced expression of epithelial-mesenchymal-transition-associated genes. T2 LCLs also have decreased RNA and protein expression of a cellular gene required for survival of T1 LCLs, IRF4. In addition to its essential role in plasma cell differentiation, IRF4 decreases BCR signaling. Knock-down of IRF4 in a T1 LCL (infected with the Zp-V3-containing Akata strain) induced lytic reactivation whereas over-expression of IRF4 in Burkitt lymphoma cells inhibited both NFATc1 and NFATc2 expression and lytic EBV reactivation. Single-cell RNA-seq confirmed that T2 LCLs have many more lytic cells compared to T1 LCLs and showed that lytically infected cells have both increased NFATc1, and decreased IRF4, compared to latently infected cells. These studies reveal numerous differences in cellular gene expression in B cells infected with T1 versus T2 EBV and suggest that decreased IRF4 contributes to both the latent and lytic phenotypes in cells with T2 EBV.

Analogous to genetically distinct alleles, epialleles represent heritable states of different gene expression from sequence-identical genes. Alleles and epialleles both contribute to phenotypic heterogeneity. While alleles originate from mutation and recombination, the source of epialleles is less well understood. We analyze active and inactive epialleles that were found at a transgenic insert with a selectable marker gene in Arabidopsis. Both converse expression states are stably transmitted to progeny. The silent epiallele was previously shown to change its state upon loss-of-function of trans-acting regulators and drug treatments. We analyzed the composition of the epialleles, their chromatin features, their nuclear localization, transcripts, and homologous small RNA. After mutagenesis by T-DNA transformation of plants carrying the silent epiallele, we found new active alleles. These switches were associated with different, larger or smaller, and non-overlapping deletions or rearrangements in the 3' regions of the epiallele. These cis-mutations caused different degrees of gene expression stability depending on the nature of the sequence alteration, the consequences for transcription and transcripts, and the resulting chromatin organization upstream. This illustrates a tight dependence of epigenetic regulation on local structures and indicates that sequence alterations can cause epigenetic changes at some distance in regions not directly affected by the mutation. Similar effects may also be involved in gene expression and chromatin changes in the vicinity of transposon insertions or excisions, recombination events, or DNA repair processes and could contribute to the origin of new epialleles.

Alternative splicing (AS) coupled to nonsense-mediated decay (NMD) is a post-transcriptional mechanism for regulating gene expression. We have used a high-resolution AS RT-PCR panel to identify endogenous AS isoforms which increase in abundance when NMD is impaired in the Arabidopsis NMD factor mutants, upf1-5 and upf3-1. Of 270 AS genes (950 transcripts) on the panel, 102 transcripts from 97 genes (32%) were identified as NMD targets. Extrapolating from these data around 13% of intron-containing genes in the Arabidopsis genome are potentially regulated by AS/NMD. This cohort of naturally occurring NMD-sensitive AS transcripts also allowed the analysis of the signals for NMD in plants. We show the importance of AS in introns in 5' or 3'UTRs in modulating NMD-sensitivity of mRNA transcripts. In particular, we identified upstream open reading frames overlapping the main start codon as a new trigger for NMD in plants and determined that NMD is induced if 3'-UTRs were >350 nt. Unexpectedly, although many intron retention transcripts possess NMD features, they are not sensitive to NMD. Finally, we have shown that AS/NMD regulates the abundance of transcripts of many genes important for plant development and adaptation including transcription factors, RNA processing factors and stress response genes.

RNA expression at isoform level is biologically more informative than at gene level and can potentially reveal cellular subsets and corresponding biomarkers that are not visible at gene level. However, due to the strong 3' bias sequencing protocol, mRNA quantification for high-throughput single-cell RNA sequencing such as Chromium Single Cell 3' 10× Genomics is currently performed at the gene level.