Large-Scale Low-Cost NGS Library Preparation Using a Robust Tn5 Purification and Tagmentation Protocol

Role of PRC2 in the stochastic expression of Aire target genes and development of mimetic cells in the thymus

The transcriptional targets of Aire and the mechanisms controlling their expression in medullary thymic epithelial cells (mTECs) need to be clarified to understand Aire's tolerogenic function. By using a multi-omics single-cell approach coupled with deep scRNA-seq, we examined how Aire controls the transcription of a wide variety of genes in a small fraction of Aire-expressing cells. We found that chromatin repression by PRC2 is an important step for Aire to achieve stochastic gene expression. Aire unleashed the silenced chromatin configuration caused by PRC2, thereby increasing the expression of its functional targets. Besides this preconditioning for Aire's gene induction, we demonstrated that PRC2 also controls the composition of mTECs that mimic the developmental trait of peripheral tissues, i.e., mimetic cells. Of note, this action of PRC2 was independent of Aire and it was more apparent than Aire. Thus, our study uncovered the essential role of polycomb complex for Aire-mediated promiscuous gene expression and the development of mimetic cells.

Medication use is associated with distinct microbial features in anxiety and depression

This study investigated the relationship between gut microbiota and neuropsychiatric disorders (NPDs), specifically anxiety disorder (ANXD) and/or major depressive disorder (MDD), as defined by Diagnostic and Statistical Manual of Mental Disorders (DSM)-IV or V criteria. The study also examined the influence of medication use, particularly antidepressants and/or anxiolytics, classified through the Anatomical Therapeutic Chemical (ATC) Classification System, on the gut microbiota. Both 16S rRNA gene amplicon sequencing (16S) and shallow shotgun sequencing (WGS) were performed on DNA extracted from 666 fecal samples from the Tulsa-1000 and Neurocomputational Mechanisms of Affiliation and Personality Study Center for Biomedical Research Excellence (NeuroMAP CoBRE) cohorts. The results highlight the significant influence of medication use; antidepressant use is associated with significant differences in gut microbiota beta diversity and has a larger effect size than NPD diagnosis. Next, specific microbes were associated with ANXD and MDD, highlighting their potential for non-pharmacological intervention. Finally, the study demonstrated the capability of Random Forest classifiers to predict diagnoses of NPD and medication use from microbial profiles, suggesting a promising direction for the use of gut microbiota as biomarkers for NPD. Though the effect sizes were larger in females than males, similar trends emerged for both sexes. These findings encourage future research on the gut microbiota's role in NPD and its interactions with pharmacological treatments.

Semaphorin 3E-Plexin-D1 Pathway Downstream of the Luteinizing Hormone Surge Regulates Ovulation, Granulosa Cell Luteinization, and Ovarian Angiogenesis in Mice

Ovulation is induced by the luteinizing hormone (LH) surge and accompanied by granulosa cell luteinization and ovarian angiogenesis. Semaphorin 3E (Sema3E)-Plexin-D1 pathway regulates angiogenesis in other tissues, but its role in the ovary is unknown. Evidence indicates that Sema3E-Plexin-D1 pathway plays an important role in the mouse ovary. The expression of Sema3E and its receptor, Plexin-D1, is dynamically regulated in the mouse ovary downstream of the LH surge. This regulation requires the modulation of chromatin accessibility by CCAAT/enhancer-binding proteins α and β. Intraovarian injection of recombinant Sema3E results in reduced ovulation, impaired corpus luteum formation, and aberrant ovarian angiogenesis. These in vivo physiological abnormalities are consistent with altered expression of genes regulating these processes, and with data from in vitro cultured granulosa cells and ovarian stromal tissues treated with Sema3E or neutralizing antibody of Plexin-D1. The findings pinpoint Sema3E-Plexin-D1 pathway as a potential therapeutic target for fertility and infertility management.

CCAAT/Enhancer-Binding Proteins α and β Regulate Ovulation and Gene Expression via Dose- and Stage-Dependent Mechanisms

The preovulatory luteinizing hormone (LH) surge orchestrates complex cellular and molecular events leading to ovulation. CCAAT/enhancer-binding proteins α and β (C/EBPα/β) are transcription factors acutely induced by the LH surge and crucial for ovulation and granulosa cell luteinization. However, biological processes (BPs) and their regulatory mechanisms downstream of C/EBPα/β in the preovulatory ovary are not completely understood. To address this knowledge gap, we generated Cebpa/bfl/fl;Pgr-Cre mutants and compared them with Cebpa/bfl/fl;Cyp19a1-Cre mutant female mice: Cebpa/bfl/fl;Cyp19a1-Cre mutants have undetectable levels of C/EBPα/β throughout the preovulatory stages and do not ovulate, aligning with previous reports; and Cebpa/bfl/fl;Pgr-Cre mutants present gradual depletion of C/EBPα/β during the late preovulatory stage and a reduced ovulation rate. Comparison of these two models indicates that sustained expression of C/EBPα/β throughout the preovulatory stages is necessary for successful ovulation and provides a unique opportunity to interrogate gene regulatory mechanisms by C/EBPα/β during different preovulatory time windows and the effect of dysregulating C/EBPα/β on ovulation-associated BPs. Our study revealed that C/EBPα/β regulate gene expression and distinct biological functions such as vascular remodeling via dose- and preovulatory stage-dependent mechanisms. These findings shed new light on the intricate mechanisms of gene regulation by C/EBPα/β downstream of the LH surge.

Design principles of cell-state-specific enhancers in hematopoiesis

During cellular differentiation, enhancers transform overlapping gradients of transcription factors (TFs) to highly specific gene expression patterns. However, the vast complexity of regulatory DNA impedes the identification of the underlying cis-regulatory rules. Here, we characterized 64,400 fully synthetic DNA sequences to bottom-up dissect design principles of cell-state-specific enhancers in the context of the differentiation of blood stem cells to seven myeloid lineages. Focusing on binding sites for 38 TFs and their pairwise interactions, we found that identical sites displayed both repressive and activating function as a consequence of cell state, site combinatorics, or simply predicted occupancy of a TF on an enhancer. Surprisingly, combinations of activating sites frequently neutralized one another or gained repressive function. These negative synergies convert quantitative imbalances in TF expression into binary activity patterns. We exploit this principle to automatically create enhancers with specificity to user-defined combinations of hematopoietic progenitor cell states from scratch.

Glycolytic activity instructs germ layer proportions through regulation of Nodal and Wnt signaling

Metabolic pathways can influence cell fate decisions, yet their regulative role during embryonic development remains poorly understood. Here, we demonstrate an instructive role of glycolytic activity in regulating signaling pathways involved in mesoderm and endoderm specification. Using a mouse embryonic stem cell (mESC)-based in vitro model for gastrulation, we found that glycolysis inhibition increases ectodermal cell fates at the expense of mesodermal and endodermal lineages. We demonstrate that this relationship is dose dependent, enabling metabolic control of germ layer proportions through exogenous glucose levels. We further show that glycolysis acts as an upstream regulator of Nodal and Wnt signaling and that its influence on cell fate specification can be decoupled from its effects on growth. Finally, we confirm the generality of our findings using a human gastrulation model. Our work underscores the dependence of signaling pathways on metabolic conditions and provides mechanistic insight into the nutritional regulation of cell fate decision-making.

DIS3L, cytoplasmic exosome catalytic subunit, is essential for development but not cell viability in mice

Among numerous enzymes involved in RNA decay, processive exoribonucleases are the most prominent group responsible for the degradation of entire RNA molecules. The role of mammalian cytoplasmic 3'-5' exonuclease DIS3L at the organismal level remained unknown. Herein, we established knock-in and knockout (KO) mouse models to study DIS3L functions in mice. DIS3L in mice is indeed a subunit of the cytoplasmic exosome complex, the disruption of which leads to severe embryo degeneration and death in mice soon after implantation. These changes could not be prevented by supplementing extraembryonic tissue with functional DIS3L through the construction of chimeric embryos. Preimplantation Dis3l -/- embryos were unaffected in their morphology and ability to produce functional embryonic stem (ES) cells, showing that DIS3L is not essential for cell viability. There were also no major changes at the transcriptome level for both ES cells and blastocysts, as revealed by RNA-seq experiments. Notably, however, Dis3l KO led to inhibition of global protein synthesis. These results point to the essential role of DIS3L in mRNA metabolism, which is crucial for proper protein synthesis during embryo development.

Perturbing nuclear glycosylation in the mouse preimplantation embryo slows down embryonic development

The main form of intracellular protein glycosylation (O-GlcNAc) is reversible and has been mapped on thousands of cytoplasmic and nuclear proteins, including RNA polymerase II, transcription factors, and chromatin modifiers. The O-GlcNAc modification is catalyzed by a single enzyme known as O-GlcNAc Transferase, that is required for mammalian early development. Yet, neither the regulatory function of protein O-GlcNAcylation in the embryo nor the embryonic O-GlcNAc proteome have been documented. Here, we devised a strategy to enzymatically remove O-GlcNAc from preimplantation embryonic nuclei, where this modification accumulates coincidently with embryonic genome activation (EGA). Unexpectedly, the depletion of nuclear O-GlcNAc to undetectable levels has no impact on EGA, but dampens the transcriptional upregulation of the translational machinery, and triggers a spindle checkpoint response. These molecular alterations were phenotypically associated with a developmental delay starting from early cleavage stages and persisting after embryo implantation, establishing a link between nuclear glycosylation and the pace of embryonic development.

Interpreting regulatory mechanisms of Hippo signaling through a deep learning sequence model

Signaling pathway components are well studied, but how they mediate cell-type-specific transcription responses is an unresolved problem. Using the Hippo pathway in mouse trophoblast stem cells as a model, we show that the DNA binding of signaling effectors is driven by cell-type-specific sequence rules that can be learned genome wide by deep learning models. Through model interpretation and experimental validation, we show that motifs for the cell-type-specific transcription factor TFAP2C enhance TEAD4/YAP1 binding in a nucleosome-range and distance-dependent manner, driving synergistic enhancer activation. We also discovered that Tead double motifs are widespread, highly active canonical response elements. Molecular dynamics simulations suggest that TEAD4 binds them cooperatively through surprisingly labile protein-protein interactions that depend on the DNA template. These results show that the response to signaling pathways is encoded in the cis-regulatory sequences and that interpreting the rules reveals insights into the mechanisms by which signaling effectors influence cell-type-specific enhancer activity.

Linking CRISPR-Cas9 double-strand break profiles to gene editing precision with BreakTag

Cas9 can cleave DNA in both blunt and staggered configurations, resulting in distinct editing outcomes, but what dictates the type of Cas9 incisions is largely unknown. In this study, we developed BreakTag, a versatile method for profiling Cas9-induced DNA double-strand breaks (DSBs) and identifying the determinants of Cas9 incisions. Overall, we assessed cleavage by SpCas9 at more than 150,000 endogenous on-target and off-target sites targeted by approximately 3,500 single guide RNAs. We found that approximately 35% of SpCas9 DSBs are staggered, and the type of incision is influenced by DNA:gRNA complementarity and the use of engineered Cas9 variants. A machine learning model shows that Cas9 incision is dependent on the protospacer sequence and that human genetic variation impacts the configuration of Cas9 cuts and the DSB repair outcome. Matched datasets of Cas9 and engineered variant incisions with repair outcomes show that Cas9-mediated staggered breaks are linked with precise, templated and predictable single-nucleotide insertions, demonstrating that a scission-based gRNA design can be used to correct clinically relevant pathogenic single-nucleotide deletions.

Single-nucleus transcriptional profiling of the placenta reveals the syncytiotrophoblast stress response to COVID-19

BACKGROUND

COVID-19 in pregnancy is associated with placental immune activation, inflammation, and vascular malperfusion, but its impact on syncytiotrophoblast biology and function is unclear.

OBJECTIVE

This study aimed to determine the effects of maternal COVID-19 on placental syncytiotrophoblasts using single-nucleus transcriptional profiling and to compare placental stress responses in COVID-19 and preeclampsia.

STUDY DESIGN

For transcriptional characterization of syncytiotrophoblasts, we used the single-nucleus RNA sequencing platform, single-cell combinatorial indexing RNA sequencing (sci-RNA-seq3), to profile placental villi and fetal membranes from unvaccinated patients with symptomatic COVID-19 at birth (n = 4), gestational age-matched controls (n = 4), and a case of critical COVID-19 in the second trimester with delivery at term (n = 1). Clustering of nuclei and differential gene expression analysis was performed in Seurat. Gene ontology analysis was conducted using Enrichr. High-confidence transcriptional target analysis was used to identify key transcription factor nodes governing the syncytiotrophoblast response to maternal SARS-CoV-2 infection. Bioinformatic approaches were further used to compare the COVID-19 dataset to published preeclampsia gene signatures. Tissue analysis, including immunofluorescence, was conducted to validate the transcriptional data and to compare COVID-19 and preeclampsia placental histology for an expanded cohort of placentas: controls (n = 6), asymptomatic COVID-19 (n = 3), symptomatic COVID-19 (n = 5), and preeclampsia with severe features (n = 7).

RESULTS

The analyzed dataset comprised 15 cell clusters and 47,889 nuclei. We identified 3 clusters of syncytiotrophoblasts representing fusing and mature nuclei with overlapping but distinct transcriptional responses to COVID-19. Bioinformatic analyses indicated that COVID-19 is associated with the following alterations in syncytiotrophoblasts: (1) endoplasmic reticulum stress and activation of stress signaling pathways, including the unfolded protein response and integrated stress response; (2) regulation of gene expression by CCAAT/enhancer-binding protein beta (CEBPB), a master transcription factor of the syncytiotrophoblast lineage; and (3) upregulation of preeclampsia-associated genes. Using complementary methods, we confirmed increased levels of stress response proteins (eg, BiP, G3BP1) in syncytiotrophoblasts, unfolded protein response signaling (spliced XBP1 mRNA), and CEBPB activation (phosphorylation) in COVID-19. Increased cytotrophoblast proliferation (Ki-67) was also detected in COVID-19, consistent with a trophoblast response to injury. Markers of stress detected in preeclampsia demonstrated similarities in the placental stress phenotype of COVID-19 and preeclampsia.

CONCLUSION

Maternal COVID-19 is associated with syncytiotrophoblast endoplasmic reticulum stress and activation of the syncytiotrophoblast lineage transcription factor, CEBPB. Similarities between syncytiotrophoblast stress in COVID-19 and preeclampsia provide insights into their clinical association.

Exportin-1 functions as an adaptor for transcription factor-mediated docking of chromatin at the nuclear pore complex

Nuclear pore proteins (nucleoporins [Nups]) physically interact with hundreds of chromosomal sites, impacting transcription. In yeast, transcription factors mediate interactions between Nups and enhancers and promoters. To define the molecular basis of this mechanism, we exploited a separation-of-function mutation in the Gcn4 transcription factor that blocks its interaction with the nuclear pore complex (NPC). This mutation reduces the interaction of Gcn4 with the highly conserved nuclear export factor Crm1/Xpo1. Crm1 and Nups co-occupy enhancers, and Crm1 inhibition blocks interaction of the nuclear pore protein Nup2 with the genome. In vivo, Crm1 interacts stably with the NPC and in vitro, Crm1 binds directly to both Gcn4 and Nup2. Importantly, the interaction between Crm1 and Gcn4 requires neither Ran-guanosine triphosphate (GTP) nor the nuclear export sequence binding site. Finally, Crm1 and Ran-GTP stimulate DNA binding by Gcn4, suggesting that allosteric coupling between Crm1-Ran-GTP binding and DNA binding facilitates the docking of transcription-factor-bound enhancers at the NPC.

Identifying cross-lineage dependencies of cell-type-specific regulators in mouse gastruloids

Correct gene expression levels are crucial for normal development. Advances in genomics enable the inference of gene regulatory programs active during development but cannot capture the complex multicellular interactions occurring during mammalian embryogenesis in utero. In vitro models of mammalian development, like gastruloids, can overcome this limitation. Using time-resolved single-cell chromatin accessibility analysis, we delineated the regulatory profile during mouse gastruloid development, identifying critical drivers of developmental transitions. Gastruloids develop from bipotent progenitor cells driven by the transcription factors (TFs) OCT4, SOX2, and TBXT, differentiating into the mesoderm (characterized by the mesogenin 1 [MSGN1]) and spinal cord (characterized by CDX2). ΔCDX gastruloids fail to form spinal cord, while Msgn1 ablation inhibits paraxial mesoderm and spinal cord development. Chimeric gastruloids with ΔMSGN1 and wild-type cells formed both tissues, indicating that inter-tissue communication is necessary for spinal cord formation. Our work has important implications for studying inter-tissue communication and gene regulatory programs in development.

Systematic screening of enhancer-blocking insulators in Drosophila identifies their DNA sequence determinants

Long-range transcriptional activation of gene promoters by abundant enhancers in animal genomes calls for mechanisms to limit inappropriate regulation. DNA elements called insulators serve this purpose by shielding promoters from an enhancer when interposed. Unlike promoters and enhancers, insulators have not been systematically characterized due to lacking high-throughput screening assays, and questions regarding how insulators are distributed and encoded in the genome remain. Here, we establish "insulator-seq" as a plasmid-based massively parallel reporter assay in Drosophila cultured cells to perform a systematic insulator screen of selected genomic loci. Screening developmental gene loci showed that not all insulator protein binding sites effectively block enhancer-promoter communication. Deep insulator mutagenesis identified sequences flexibly positioned around the CTCF insulator protein binding motif that are critical for functionality. The ability to screen millions of DNA sequences without positional effect has enabled functional mapping of insulators and provided further insights into the determinants of insulators.

Telomemore enables single-cell analysis of cell cycle and chromatin condensation

Single-cell RNA-seq methods can be used to delineate cell types and states at unprecedented resolution but do little to explain why certain genes are expressed. Single-cell ATAC-seq and multiome (ATAC + RNA) have emerged to give a complementary view of the cell state. It is however unclear what additional information can be extracted from ATAC-seq data besides transcription factor binding sites. Here, we show that ATAC-seq telomere-like reads counter-inituively cannot be used to infer telomere length, as they mostly originate from the subtelomere, but can be used as a biomarker for chromatin condensation. Using long-read sequencing, we further show that modern hyperactive Tn5 does not duplicate 9 bp of its target sequence, contrary to common belief. We provide a new tool, Telomemore, which can quantify nonaligning subtelomeric reads. By analyzing several public datasets and generating new multiome fibroblast and B-cell atlases, we show how this new readout can aid single-cell data interpretation. We show how drivers of condensation processes can be inferred, and how it complements common RNA-seq-based cell cycle inference, which fails for monocytes. Telomemore-based analysis of the condensation state is thus a valuable complement to the single-cell analysis toolbox.

A Senescent Cluster in Aged Human Hematopoietic Stem Cell Compartment as Target for Senotherapy

To identify the differences between aged and young human hematopoiesis, we performed a direct comparison of aged and young human hematopoietic stem and progenitor cells (HSPCs). Alterations in transcriptome profiles upon aging between humans and mice were then compared. Human specimens consist of CD34+ cells from bone marrow, and mouse specimens of hematopoietic stem cells (HSCs; Lin- Kit+ Sca1+ CD150+). Single-cell transcriptomic studies, functional clustering, and developmental trajectory analyses were performed. A significant increase in multipotent progenitor 2A (MPP2A) cluster is found in the early HSC trajectory in old human subjects. This cluster is enriched in senescence signatures (increased telomere attrition, DNA damage, activation of P53 pathway). In mouse models, the accumulation of an analogous subset was confirmed in the aged LT-HSC population. Elimination of this subset has been shown to rejuvenate hematopoiesis in mice. A significant activation of the P53-P21WAF1/CIP1 pathway was found in the MPP2A population in humans. In contrast, the senescent HSCs in mice are characterized by activation of the p16Ink4a pathway. Aging in the human HSC compartment is mainly caused by the clonal evolution and accumulation of a senescent cell cluster. A population with a similar senescence signature in the aged LT-HSCs was confirmed in the murine aging model. Clearance of this senescent population with senotherapy in humans is feasible and potentially beneficial.

Validation of Golden Gate Assemblies Using Highly Multiplexed Nanopore Amplicon Sequencing

Golden Gate cloning has revolutionized synthetic biology. Its concept of modular, highly characterized libraries of parts that can be combined into higher order assemblies allows engineering principles to be applied to biological systems. The basic parts, typically stored in Level 0 plasmids, are sequence validated by the method of choice and can be combined into higher order assemblies on demand. Higher order assemblies are typically transcriptional units, and multiple transcriptional units can be assembled into multi-gene constructs. Higher order Golden Gate assembly based on defined and validated parts usually does not introduce sequence changes. Therefore, simple validation of the assemblies, e.g., by colony polymerase chain reaction (PCR) or restriction digest pattern analysis is sufficient. However, in many experimental setups, researchers do not use defined parts, but rather part libraries, resulting in assemblies of high combinatorial complexity where sequencing again becomes mandatory. Here, we present a detailed protocol for the use of a highly multiplexed dual barcode amplicon sequencing using the Nanopore sequencing platform for in-house sequence validation. The workflow, called DuBA.flow, is a start-to-finish procedure that provides all necessary steps from a single colony to the final easy-to-interpret sequencing report.

Targeted barcoding of variable antibody domains and individual transcriptomes of the human B-cell repertoire using Link-Seq

Here, we present Link-Seq, a highly efficient droplet microfluidic method for combined sequencing of antibody-encoding genes and the transcriptome of individual B cells at large scale. The method is based on 3' barcoding of the transcriptome and subsequent single-molecule PCR in droplets, which freely shift the barcode along specific gene regions, such as the antibody heavy- and light-chain genes. Using the immune repertoire of COVID-19 patients and healthy donors as a model system, we obtain up to 91.7% correctly paired immunoglobulin heavy and light chains. Furthermore, we map the V(D)J usage and obtain sensitivities comparable with the current gold-standard 10× Genomics commercial systems while offering full flexibility in experimental setup and significant cost savings. A further unique feature of Link-Seq is the possibility of barcoding multiple target genes in a site-specific manner. Based on the open character of the platform and its conceptual advantages, we expect Link-Seq to become a versatile tool for single-cell analysis, especially for applications requiring additional processing steps that cannot be implemented on commercially available platforms.

A protocol for time-resolved transcriptomics through metabolic labeling and combinatorial indexing

The snapshot nature of single-cell transcriptomics presents a challenge for studying the dynamics of gene expression. Metabolic labeling, where nascent RNA is labeled with 4-thiouridine (4sU), captures temporal information at the single-cell level, providing greater insight into expression dynamics. Here, we present an optimized, automation-friendly protocol for the metabolic labeling of RNA alongside single-cell RNA sequencing through combinatorial indexing. We describe steps for 4sU labeling, cell fixation and chemical treatment, and automated two-level combinatorial indexing. For complete details on the use and execution of this protocol, please refer to Maizels et al.1.

Protocol for detecting genomic insulators in Drosophila using insulator-seq, a massively parallel reporter assay

Genomic insulators are DNA elements that prevent transcriptional activation of a promoter by an enhancer when interposed. We present a protocol for insulator-seq that enables high-throughput screening of genomic insulators using a plasmid-based massively parallel reporter assay in Drosophila cultured cells. We describe steps for insulator reporter plasmid library generation, transient transfection into cultured cells, and sequencing library preparation and provide a pipeline for data analysis. For complete details on the use and execution of this protocol, please refer to Tonelli et al.1.

Dextromethorphan inhibits collagen and collagen-like cargo secretion to ameliorate lung fibrosis

Excessive deposition of fibrillar collagen in the interstitial extracellular matrix (ECM) of human lung tissue causes fibrosis, which can ultimately lead to organ failure. Despite our understanding of the molecular mechanisms underlying the disease, no cure for pulmonary fibrosis has yet been found. We screened a drug library and found that dextromethorphan (DXM), a cough expectorant, reduced the amount of excess fibrillar collagen deposited in the ECM in cultured primary human lung fibroblasts, a bleomycin mouse model, and a cultured human precision-cut lung slice model of lung fibrosis. The reduced extracellular fibrillar collagen upon DXM treatment was due to reversible trafficking inhibition of collagen type I (COL1) in the endoplasmic reticulum (ER) in TANGO1- and HSP47-positive structures. Mass spectrometric analysis showed that DXM promoted hyperhydroxylation of proline and lysine residues on various collagens (COL1, COL3, COL4, COL5, COL7, and COL12) and latent transforming growth factor-β-binding protein (LTBP1 and LTBP2) peptides, coinciding with their secretion block. Additionally, proteome profiling of DXM-treated cells showed increased thermal stability of prolyl-hydroxylases P3H2, P3H3, P3H4, P4HA1, and P4HA2, suggesting a change in their activity. Transcriptome analysis of profibrotic stimulated primary human lung fibroblasts and human ex vivo lung slices after DXM treatment showed activation of an antifibrotic program through regulation of multiple pathways, including the MMP-ADAMTS axis, WNT signaling, and fibroblast-to-myofibroblast differentiation. Together, these data obtained from in vitro, in vivo, and ex vivo models of lung fibrogenesis show that DXM has the potential to limit fibrosis through inhibition of COL1 membrane trafficking in the ER.

HIPSD&R-seq enables scalable genomic copy number and transcriptome profiling

Single-cell DNA sequencing (scDNA-seq) enables decoding somatic cancer variation. Existing methods are hampered by low throughput or cannot be combined with transcriptome sequencing in the same cell. We propose HIPSD&R-seq (HIgh-throughPut Single-cell Dna and Rna-seq), a scalable yet simple and accessible assay to profile low-coverage DNA and RNA in thousands of cells in parallel. Our approach builds on a modification of the 10X Genomics platform for scATAC and multiome profiling. In applications to human cell models and primary tissue, we demonstrate the feasibility to detect rare clones and we combine the assay with combinatorial indexing to profile over 17,000 cells.

Single-cell multi-omic analysis of fibrolamellar carcinoma reveals rewired cell-to-cell communication patterns and unique vulnerabilities

Fibrolamellar carcinoma (FLC) is a rare malignancy disproportionately affecting adolescents and young adults with no standard of care. FLC is characterized by thick stroma, which has long suggested an important role of the tumor microenvironment. Over the past decade, several studies have revealed aberrant markers and pathways in FLC. However, a significant drawback of these efforts is that they were conducted on bulk tumor samples. Consequently, identities and roles of distinct cell types within the tumor milieu, and the patterns of intercellular communication, have yet to be explored. In this study we unveil cell-type specific gene signatures, transcription factor networks, and super-enhancers in FLC using a multi-omics strategy that leverages both single-nucleus ATAC-seq and single-nucleus RNA-seq. We also infer completely rewired cell-to-cell communication patterns in FLC including signaling mediated by SPP1-CD44, MIF-ACKR3, GDF15-TGFBR2, and FGF7-FGFR. Finally, we validate findings with loss-of-function studies in several models including patient tissue slices, identifying vulnerabilities that merit further investigation as candidate therapeutic targets in FLC.

Uridylation regulates mRNA decay directionality in fission yeast

Cytoplasmic mRNA decay is effected by exonucleolytic degradation in either the 5' to 3' or 3' to 5' direction. Pervasive terminal uridylation is implicated in mRNA degradation, however, its functional relevance for bulk mRNA turnover remains poorly understood. In this study, we employ genome-wide 3'-RACE (gw3'-RACE) in the model system fission yeast to elucidate the role of uridylation in mRNA turnover. We observe widespread uridylation of shortened poly(A) tails, promoting efficient 5' to 3' mRNA decay and ensuring timely and controlled mRNA degradation. Inhibition of this uridylation process leads to excessive deadenylation and enhanced 3' to 5' mRNA decay accompanied by oligouridylation. Strikingly we found that uridylation of poly(A) tails and oligouridylation of non-polyadenylated substrates are catalysed by different terminal uridyltransferases Cid1 and Cid16 respectively. Our study sheds new light on the intricate regulatory mechanisms underlying bulk mRNA turnover, demonstrating the role of uridylation in modulating mRNA decay pathways.

Protocol for multimodal profiling of human kidneys with simultaneous high-throughput ATAC and RNA expression with sequencing

Simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq) profiles transcriptomics and chromatin accessibility in the same cells at high throughput. Here, we present a protocol for multimodal profiling of human kidneys with SHARE-seq. We describe steps for processing fixed nuclei for SHARE-seq split-pool barcoding and library preparation. We also detail how to determine the optimal working concentration of Tn5 transposase for transposition and tagmentation. This protocol allows researchers to generate large-scale single-cell multiomics data at low reagent cost. For complete details on the use and execution of this protocol, please refer to Li et al.1.

Modular control of vertebrate axis segmentation in time and space

How the timing of development is linked to organismal size is a longstanding question. Although numerous studies have reported a correlation of temporal and spatial traits, the developmental or selective constraints underlying this link remain largely unexplored. We address this question by studying the periodic process of embryonic axis segmentation in-vivo in Oryzias fish. Interspecies comparisons reveal that the timing of segmentation correlates to segment, tissue and organismal size. Segment size in turn scales according to tissue and organism size. To probe for underlying causes, we genetically hybridised two closely related species. Quantitative analysis in ~600 phenotypically diverse F2 embryos reveals a decoupling of timing from size control, while spatial scaling is preserved. Using developmental quantitative trait loci (devQTL) mapping we identify distinct genetic loci linked to either the control of segmentation timing or tissue size. This study demonstrates that a developmental constraint mechanism underlies spatial scaling of axis segmentation, while its spatial and temporal control are dissociable modules.

Genotype imputation in F2 crosses of inbred lines

Motivation

Crosses among inbred lines are a fundamental tool for the discovery of genetic loci associated with phenotypes of interest. In organisms for which large reference panels or SNP chips are not available, imputation from low-pass whole-genome sequencing is an effective method for obtaining genotype data from a large number of individuals. To date, a structured analysis of the conditions required for optimal genotype imputation has not been performed.

Results

We report a systematic exploration of the effect of several design variables on imputation performance in F2 crosses of inbred medaka lines using the imputation software STITCH. We determined that, depending on the number of samples, imputation performance reaches a plateau when increasing the per-sample sequencing coverage. We also systematically explored the trade-offs between cost, imputation accuracy, and sample numbers. We developed a computational pipeline to streamline the process, enabling other researchers to perform a similar cost-benefit analysis on their population of interest.

Availability and implementation

The source code for the pipeline is available at https://github.com/birneylab/stitchimpute. While our pipeline has been developed and tested for an F2 population, the software can also be used to analyse populations with a different structure.

OpenTn5: Open-Source Resource for Robust and Scalable Tn5 Transposase Purification and Characterization

Tagmentation combines DNA fragmentation and sequencing adapter addition by leveraging the transposition activity of the bacterial cut-and-paste Tn5 transposase, to enable efficient sequencing library preparation. Here we present an open-source protocol for the generation of multi-purpose hyperactive Tn5 transposase, including its benchmarking in CUT&Tag, bulk and single-cell ATAC-seq. The OpenTn5 protocol yields multi-milligram quantities of pG-Tn5E54K, L372P protein per liter of E. coli culture, sufficient for thousands of tagmentation reactions and the enzyme retains activity in storage for more than a year.

Transposase-assisted tagmentation: an economical and scalable strategy for single-worm whole-genome sequencing

AlphaMissense identifies 23 million human missense variants as likely pathogenic, but only 0.1% have been clinically classified. To experimentally validate these predictions, chemical mutagenesis presents a rapid, cost-effective method to produce billions of mutations in model organisms. However, the prohibitive costs and limitations in the throughput of whole-genome sequencing (WGS) technologies, crucial for variant identification, constrain its widespread application. Here, we introduce a Tn5 transposase-assisted tagmentation technique for conducting WGS in Caenorhabditis elegans, Escherichia coli, Saccharomyces cerevisiae, and Chlamydomonas reinhardtii. This method, demands merely 20 min of hands-on time for a single-worm or single-cell clones and incurs a cost below 10 US dollars. It effectively pinpoints causal mutations in mutants defective in cilia or neurotransmitter secretion and in mutants synthetically sterile with a variant analogous to the B-Raf Proto-oncogene, Serine/Threonine Kinase (BRAF) V600E mutation. Integrated with chemical mutagenesis, our approach can generate and identify missense variants economically and efficiently, facilitating experimental investigations of missense variants in diverse species.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

TENT5-mediated polyadenylation of mRNAs encoding secreted proteins is essential for gametogenesis in mice

Cytoplasmic polyadenylation plays a vital role in gametogenesis; however, the participating enzymes and substrates in mammals remain unclear. Using knockout and knock-in mouse models, we describe the essential role of four TENT5 poly(A) polymerases in mouse fertility and gametogenesis. TENT5B and TENT5C play crucial yet redundant roles in oogenesis, with the double knockout of both genes leading to oocyte degeneration. Additionally, TENT5B-GFP knock-in females display a gain-of-function infertility effect, with multiple chromosomal aberrations in ovulated oocytes. TENT5C and TENT5D both regulate different stages of spermatogenesis, as shown by the sterility in males following the knockout of either gene. Finally, Tent5a knockout substantially lowers fertility, although the underlying mechanism is not directly related to gametogenesis. Through direct RNA sequencing, we discovered that TENT5s polyadenylate mRNAs encoding endoplasmic reticulum-targeted proteins essential for gametogenesis. Sequence motif analysis and reporter mRNA assays reveal that the presence of an endoplasmic reticulum-leader sequence represents the primary determinant of TENT5-mediated regulation.

Specialized Tfh cell subsets driving type-1 and type-2 humoral responses in lymphoid tissue

Effective antibody responses are essential to generate protective humoral immunity. Different inflammatory signals polarize T cells towards appropriate effector phenotypes during an infection or immunization. Th1 and Th2 cells have been associated with the polarization of humoral responses. However, T follicular helper cells (Tfh) have a unique ability to access the B cell follicle and support the germinal center (GC) responses by providing B cell help. We investigated the specialization of Tfh cells induced under type-1 and type-2 conditions. We first studied homogenous Tfh cell populations generated by adoptively transferred TCR-transgenic T cells in mice immunized with type-1 and type-2 adjuvants. Using a machine learning approach, we established a gene expression signature that discriminates Tfh cells polarized towards type-1 and type-2 response, defined as Tfh1 and Tfh2 cells. The distinct signatures of Tfh1 and Tfh2 cells were validated against datasets of Tfh cells induced following lymphocytic choriomeningitis virus (LCMV) or helminth infection. We generated single-cell and spatial transcriptomics datasets to dissect the heterogeneity of Tfh cells and their localization under the two immunizing conditions. Besides a distinct specialization of GC Tfh cells under the two immunizations and in different regions of the lymph nodes, we found a population of Gzmk+ Tfh cells specific for type-1 conditions. In human individuals, we could equally identify CMV-specific Tfh cells that expressed Gzmk. Our results show that Tfh cells acquire a specialized function under distinct types of immune responses and with particular properties within the B cell follicle and the GC.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Reconstructing developmental trajectories using latent dynamical systems and time-resolved transcriptomics

The snapshot nature of single-cell transcriptomics presents a challenge for studying the dynamics of cell fate decisions. Metabolic labeling and splicing can provide temporal information at single-cell level, but current methods have limitations. Here, we present a framework that overcomes these limitations: experimentally, we developed sci-FATE2, an optimized method for metabolic labeling with increased data quality, which we used to profile 45,000 embryonic stem (ES) cells differentiating into neural tube identities. Computationally, we developed a two-stage framework for dynamical modeling: VelvetVAE, a variational autoencoder (VAE) for velocity inference that outperforms all other tools tested, and VelvetSDE, a neural stochastic differential equation (nSDE) framework for simulating trajectory distributions. These recapitulate underlying dataset distributions and capture features such as decision boundaries between alternative fates and fate-specific gene expression. These methods recast single-cell analyses from descriptions of observed data to models of the dynamics that generated them, providing a framework for investigating developmental fate decisions.

Transcriptomic, epigenomic, and spatial metabolomic cell profiling redefines regional human kidney anatomy

A large-scale multimodal atlas that includes major kidney regions is lacking. Here, we employed simultaneous high-throughput single-cell ATAC/RNA sequencing (SHARE-seq) and spatially resolved metabolomics to profile 54 human samples from distinct kidney anatomical regions. We generated transcriptomes of 446,267 cells and chromatin accessibility profiles of 401,875 cells and developed a package to analyze 408,218 spatially resolved metabolomes. We find that the same cell type, including thin limb, thick ascending limb loop of Henle and principal cells, display distinct transcriptomic, chromatin accessibility, and metabolomic signatures, depending on anatomic location. Surveying metabolism-associated gene profiles revealed non-overlapping metabolic signatures between nephron segments and dysregulated lipid metabolism in diseased proximal tubule (PT) cells. Integrating multimodal omics with clinical data identified PLEKHA1 as a disease marker, and its in vitro knockdown increased gene expression in PT differentiation, suggesting possible pathogenic roles. This study highlights previously underrepresented cellular heterogeneity underlying the human kidney anatomy.

The rise of single-cell transcriptomics in yeast

The field of single-cell omics has transformed our understanding of biological processes and is constantly advancing both experimentally and computationally. One of the most significant developments is the ability to measure the transcriptome of individual cells by single-cell RNA-seq (scRNA-seq), which was pioneered in higher eukaryotes. While yeast has served as a powerful model organism in which to test and develop transcriptomic technologies, the implementation of scRNA-seq has been significantly delayed in this organism, mainly because of technical constraints associated with its intrinsic characteristics, namely the presence of a cell wall, a small cell size and little amounts of RNA. In this review, we examine the current technologies for scRNA-seq in yeast and highlight their strengths and weaknesses. Additionally, we explore opportunities for developing novel technologies and the potential outcomes of implementing single-cell transcriptomics and extension to other modalities. Undoubtedly, scRNA-seq will be invaluable for both basic and applied yeast research, providing unique insights into fundamental biological processes.

High-throughput sequencing of insect specimens with sub-optimal DNA preservation using a practical, plate-based Illumina-compatible Tn5 transposase library preparation method

Entomological sampling and storage conditions often prioritise efficiency, practicality and conservation of morphological characteristics, and may therefore be suboptimal for DNA preservation. This practice can impact downstream molecular applications, such as the generation of high-throughput genomic libraries, which often requires substantial DNA input amounts. Here, we use a practical Tn5 transposase tagmentation-based library preparation method optimised for 96-well plates and low yield DNA extracts from insect legs that were stored under sub-optimal conditions for DNA preservation. The samples were kept in field vehicles for extended periods of time, before long-term storage in ethanol in the freezer, or dry at room temperature. By reducing DNA input to 6ng, more samples with sub-optimal DNA yields could be processed. We matched this low DNA input with a 6-fold dilution of a commercially available tagmentation enzyme, significantly reducing library preparation costs. Costs and workload were further suppressed by direct post-amplification pooling of individual libraries. We generated medium coverage (>3-fold) genomes for 88 out of 90 specimens, with an average of approximately 10-fold coverage. While samples stored in ethanol yielded significantly less DNA compared to those which were stored dry, these samples had superior sequencing statistics, with longer sequencing reads and higher rates of endogenous DNA. Furthermore, we find that the efficiency of tagmentation-based library preparation can be improved by a thorough post-amplification bead clean-up which selects against both short and large DNA fragments. By opening opportunities for the use of sub-optimally preserved, low yield DNA extracts, we broaden the scope of whole genome studies of insect specimens. We therefore expect these results and this protocol to be valuable for a range of applications in the field of entomology.

txci-ATAC-seq: a massive-scale single-cell technique to profile chromatin accessibility

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16-/-) model. Our study of CC16-/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model.

Medication Use is Associated with Distinct Microbial Features in Anxiety and Depression

This study investigated the relationship between gut microbiota and neuropsychiatric disorders (NPDs), specifically anxiety disorder (ANXD) and/or major depressive disorder (MDD), as defined by DSM-IV or V criteria. The study also examined the influence of medication use, particularly antidepressants and/or anxiolytics, classified through the Anatomical Therapeutic Chemical (ATC) Classification System, on the gut microbiota. Both 16S rRNA gene amplicon sequencing and shallow shotgun sequencing were performed on DNA extracted from 666 fecal samples from the Tulsa-1000 and NeuroMAP CoBRE cohorts. The results highlight the significant influence of medication use; antidepressant use is associated with significant differences in gut microbiota beta diversity and has a larger effect size than NPD diagnosis. Next, specific microbes were associated with ANXD and MDD, highlighting their potential for non-pharmacological intervention. Finally, the study demonstrated the capability of Random Forest classifiers to predict diagnoses of NPD and medication use from microbial profiles, suggesting a promising direction for the use of gut microbiota as biomarkers for NPD. The findings suggest that future research on the gut microbiota's role in NPD and its interactions with pharmacological treatments are needed.

ChEC-seq2: an improved chromatin endogenous cleavage sequencing method and bioinformatic analysis pipeline for mapping in vivo protein-DNA interactions

Defining the in vivo DNA binding specificity of transcription factors (TFs) has relied nearly exclusively on chromatin immunoprecipitation (ChIP). While ChIP reveals TF binding patterns, its resolution is low. Higher resolution methods employing nucleases such as ChIP-exo, chromatin endogenous cleavage (ChEC-seq) and CUT&RUN resolve both TF occupancy and binding site protection. ChEC-seq, in which an endogenous TF is fused to micrococcal nuclease, requires neither fixation nor antibodies. However, the specificity of DNA cleavage during ChEC has been suggested to be lower than the specificity of the peaks identified by ChIP or ChIP-exo, perhaps reflecting non-specific binding of transcription factors to DNA. We have simplified the ChEC-seq protocol to minimize nuclease digestion while increasing the yield of cleaved DNA. ChEC-seq2 cleavage patterns were highly reproducible between replicates and with published ChEC-seq data. Combined with DoubleChEC, a new bioinformatic pipeline that removes non-specific cleavage sites, ChEC-seq2 identified high-confidence cleavage sites for three different yeast TFs that are strongly enriched for their known binding sites and adjacent to known target genes.

DuBA.flow─A Low-Cost, Long-Read Amplicon Sequencing Workflow for the Validation of Synthetic DNA Constructs

Modern biological science, especially synthetic biology, relies heavily on the construction of DNA elements, often in the form of plasmids. Plasmids are used for a variety of applications, including the expression of proteins for subsequent purification, the expression of heterologous pathways for the production of valuable compounds, and the study of biological functions and mechanisms. For all applications, a critical step after the construction of a plasmid is its sequence validation. The traditional method for sequence determination is Sanger sequencing, which is limited to approximately 1000 bp per reaction. Here, we present a highly scalable in-house method for rapid validation of amplified DNA sequences using long-read Nanopore sequencing. We developed two-step amplicon and transposase strategies to provide maximum flexibility for dual barcode sequencing. We also provide an automated analysis pipeline to quickly and reliably analyze sequencing results and provide easy-to-interpret results for each sample. The user-friendly DuBA.flow start-to-finish pipeline is widely applicable. Furthermore, we show that construct validation using DuBA.flow can be performed by barcoded colony PCR amplicon sequencing, thus accelerating research.

The UBP5 histone H2A deubiquitinase counteracts PRCs-mediated repression to regulate Arabidopsis development

Polycomb Repressive Complexes (PRCs) control gene expression through the incorporation of H2Aub and H3K27me3. In recent years, there is increasing evidence of the complexity of PRCs' interaction networks and the interplay of these interactors with PRCs in epigenome reshaping, which is fundamental to understand gene regulatory mechanisms. Here, we identified UBIQUITIN SPECIFIC PROTEASE 5 (UBP5) as a chromatin player able to counteract the deposition of the two PRCs' epigenetic hallmarks in Arabidopsis thaliana. We demonstrated that UBP5 is a plant developmental regulator based on functional analyses of ubp5-CRISPR Cas9 mutant plants. UBP5 promotes H2A monoubiquitination erasure, leading to transcriptional de-repression. Furthermore, preferential association of UBP5 at PRC2 recruiting motifs and local H3K27me3 gaining in ubp5 mutant plants suggest the existence of functional interplays between UBP5 and PRC2 in regulating epigenome dynamics. In summary, acting as an antagonist of the pivotal epigenetic repressive marks H2Aub and H3K27me3, UBP5 provides novel insights to disentangle the complex regulation of PRCs' activities.

A modular dCas9-based recruitment platform for combinatorial epigenome editing

Targeted epigenome editing tools allow precise manipulation and investigation of genome modifications, however they often display high context dependency and variable efficacy between target genes and cell types. While systems that simultaneously recruit multiple distinct 'effector' chromatin regulators can improve efficacy, they generally lack control over effector composition and spatial organisation. To overcome this we have created a modular combinatorial epigenome editing platform, called SSSavi. This system is an interchangeable and reconfigurable docking platform fused to dCas9 that enables simultaneous recruitment of up to four different effectors, allowing precise control of effector composition and spatial ordering. We demonstrate the activity and specificity of the SSSavi system and, by testing it against existing multi-effector targeting systems, demonstrate its comparable efficacy. Furthermore, we demonstrate the importance of the spatial ordering of the recruited effectors for effective transcriptional regulation. Together, the SSSavi system enables exploration of combinatorial effector co-recruitment to enhance manipulation of chromatin contexts previously resistant to targeted editing.

Selection of Flax Genotypes for Pan-Genomic Studies by Sequencing Tagmentation-Based Transcriptome Libraries

Flax (Linum usitatissimum L.) products are used in the food, pharmaceutical, textile, polymer, medical, and other industries. The creation of a pan-genome will be an important advance in flax research and breeding. The selection of flax genotypes that sufficiently cover the species diversity is a crucial step for the pan-genomic study. For this purpose, we have adapted a method based on Illumina sequencing of transcriptome libraries prepared using the Tn5 transposase (tagmentase). This approach reduces the cost of sample preparation compared to commercial kits and allows the generation of a large number of cDNA libraries in a short time. RNA-seq data were obtained for 192 flax plants (3-6 individual plants from 44 flax accessions of different morphology and geographical origin). Evaluation of the genetic relationship between flax plants based on the sequencing data revealed incorrect species identification for five accessions. Therefore, these accessions were excluded from the sample set for the pan-genomic study. For the remaining samples, typical genotypes were selected to provide the most comprehensive genetic diversity of flax for pan-genome construction. Thus, high-throughput sequencing of tagmentation-based transcriptome libraries showed high efficiency in assessing the genetic relationship of flax samples and allowed us to select genotypes for the flax pan-genomic analysis.

ChEC-seq2: an improved Chromatin Endogenous Cleavage sequencing method and bioinformatic analysis pipeline for mapping in vivo protein-DNA interactions

Defining the in vivo DNA binding specificity of transcription factors (TFs) has relied nearly exclusively on chromatin immunoprecipitation (ChIP). While ChIP reveals TF binding patterns, its resolution is low. Higher resolution methods employing nucleases such as ChIP-exo, chromatin endogenous cleavage (ChEC-seq) and CUT&RUN resolve both TF occupancy and binding site protection. ChEC-seq, in which an endogenous TF is fused to micrococcal nuclease, requires neither fixation nor antibodies. However, the specificity of DNA cleavage during ChEC has been suggested to be lower than the specificity of the peaks identified by ChIP or ChIP-exo, perhaps reflecting non-specific binding of transcription factors to DNA. We have simplified the ChEC-seq protocol to minimize nuclease digestion while increasing the yield of cleaved DNA. ChEC-seq2 cleavage patterns were highly reproducible between replicates and with published ChEC-seq data. Combined with DoubleChEC, a new bioinformatic pipeline that removes non-specific cleavage sites, ChEC-seq2 identified high-confidence cleavage sites for three different yeast TFs that are strongly enriched for their known binding sites and adjacent to known target genes.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

AP profiling resolves co-translational folding pathway and chaperone interactions in vivo

Natural proteins have evolved to fold robustly along specific pathways. Folding begins during synthesis, guided by interactions of the nascent protein with the ribosome and molecular chaperones. However, the timing and progression of co-translational folding remain largely elusive, in part because the process is difficult to measure in the natural environment of the cytosol. We developed a high-throughput method to quantify co-translational folding in live cells that we term Arrest Peptide profiling (AP profiling). We employed AP profiling to delineate co-translational folding for a set of GTPase domains with very similar structures, defining how topology shapes folding pathways. Genetic ablation of major nascent chain-binding chaperones resulted in localized folding changes that suggest how functional redundancies among chaperones are achieved by distinct interactions with the nascent protein. Collectively, our studies provide a window into cellular folding pathways of complex proteins and pave the way for systematic studies on nascent protein folding at unprecedented resolution and throughput.

Novel mechanisms of cadmium tolerance and Cd-induced fungal stress in wheat: Transcriptomic and metagenomic insights

Although several studies on the effects of cadmium (Cd) on wheat have been reported, the gene expression profiles of different wheat tissues in response to gradient concentrations of Cd, and whether soil microorganisms are involved in the damage to wheat remain to be discovered. To gain further insight into the molecular mechanisms of Cd-resistance in wheat, we sowed bread wheat (Triticum aestivum) in artificially Cd-contaminated soil and investigated the transcriptomic response of the wheat roots, stems, and leaves to gradient concentrations of Cd, as well as the alteration of the soil microbiome. Results indicated that the root bioaccumulation factors increased with Cd when concentrations were < 10 mg/kg, but at even higher concentrations, the bioaccumulation factors decreased, which is consistent with the overexpression of metal transporters and other genes related to Cd tolerance. In the Cd-contaminated soil, the abundance of fungal pathogens increased, and the antimicrobial response in wheat root was observed. Most of the differentially expressed genes (DEGs) of wheat changed significantly when the Cd concentration increased above 10 mg/kg, and the transcriptional response is much greater in roots than in stems and leaves. The DEGs are mainly involved in Cd transport and chelation, antioxidative stress, antimicrobial responses, and growth regulation. COPT3 and ZnT1 were identified for the first time as the major transporters responding to Cd in wheat. Overexpression of the nicotianamine synthase and pectinesterase genes suggested that nicotianamine and pectin are the key chelators in Cd detoxification. endochitinase, chitinase, and snakin2 were involved in the anti-fungal stress caused by Cd-induced cell damage. Several phytohormone-related DEGs are involved in the root's growth and repair. Overall, this study presents the novel Cd tolerance mechanisms in wheat and the changes in soil fungal pathogens that increase plant damage.

Clonally resolved single-cell multi-omics identifies routes of cellular differentiation in acute myeloid leukemia

Inter-patient variability and the similarity of healthy and leukemic stem cells (LSCs) have impeded the characterization of LSCs in acute myeloid leukemia (AML) and their differentiation landscape. Here, we introduce CloneTracer, a novel method that adds clonal resolution to single-cell RNA-seq datasets. Applied to samples from 19 AML patients, CloneTracer revealed routes of leukemic differentiation. Although residual healthy and preleukemic cells dominated the dormant stem cell compartment, active LSCs resembled their healthy counterpart and retained erythroid capacity. By contrast, downstream myeloid progenitors constituted a highly aberrant, disease-defining compartment: their gene expression and differentiation state affected both the chemotherapy response and leukemia's ability to differentiate into transcriptomically normal monocytes. Finally, we demonstrated the potential of CloneTracer to identify surface markers misregulated specifically in leukemic cells. Taken together, CloneTracer reveals a differentiation landscape that mimics its healthy counterpart and may determine biology and therapy response in AML.

High-resolution single-molecule long-fragment rRNA gene amplicon sequencing of bacterial and eukaryotic microbial communities

Sequencing of hypervariable regions as well as internal transcribed spacer regions of ribosomal RNA genes (rDNA) is broadly used to identify bacteria and fungi, but taxonomic and phylogenetic resolution is hampered by insufficient sequencing length using high throughput, cost-efficient second-generation sequencing. We developed a method to obtain nearly full-length rDNA by assembling single DNA molecules combining DNA co-barcoding with single-tube long fragment read technology and second-generation sequencing. Benchmarking was performed using mock bacterial and fungal communities as well as two forest soil samples. All mock species rDNA were successfully recovered with identities above 99.5% compared to the reference sequences. From the soil samples we obtained good coverage with identification of more than 20,000 unknown species, as well as high abundance correlation between replicates. This approach provides a cost-effective method for obtaining extensive and accurate information on complex environmental microbial communities.