The horizontal gene transfer of perchlorate reduction genomic island in three bacteria from an ecological niche

Three new strains of dissimilatory perchlorate-reducing bacteria (DPRB), QD19-16, QD1-5, and P3-1, were isolated from an active sludge. Phylogenetic trees based on 16S rRNA genes indicated that QD19-16, QD1-5, and P3-1 belonged to Brucella, Acidovorax, and Citrobacter, respectively, expanding the distribution of DPRB in the Proteobacteria. The three strains were gram-negative and facultative anaerobes with rod-shaped cells without flagella, which were 1.0-1.6 μm long and 0.5-0.6 μm wide. The three DPRB strains utilized similar broad spectrum of electron donors and acceptors and demonstrated a similar capability to reduce perchlorate within 6 days. The enzyme activity of perchlorate reductase in QD19-16 toward chlorate was higher than that toward perchlorate. The high sequence similarity of the perchlorate reductase operon and chlorite dismutase genes in the perchlorate reduction genomic islands (PRI) of the three strains implied that they were monophyletic origin from a common ancestral PRI. Two transposase genes (tnp1 and tnp2) were found in the PRIs of strain QD19-16 and QD1-5, but were absent in the strain P3-1 PRI. The presence of fragments of IR sequences in the P3-1 PRI suggested that P3-1 PRI had previously contained these two tnp genes. Therefore, it is plausible to suggest that a common ancestral PRI transferred across the strains Brucella sp. QD19-16, Acidovorax sp. QD1-5, and Citrobacter sp. P3-1 through horizontal gene transfer, facilitated by transposases. These results provided a direct evidence of horizontal gene transfer of PRI that could jump across phylogenetically unrelated bacteria through transposase. KEY POINTS: • Three new DPRB strains can effectively remove high concentration of perchlorate. • The PRIs of three DPRB strains are acquired from a single ancestral PRI. • PRIs are incorporated into different bacteria genome through HGT by transposase.

An optimized polymeric delivery system for piggyBac transposition

The piggyBac transposon/transposase system has been explored for long-term, stable gene expression to execute genomic integration of therapeutic genes, thus emerging as a strong alternative to viral transduction. Most studies with piggyBac transposition have employed physical methods for successful delivery of the necessary components of the piggyBac system into the cells. Very few studies have explored polymeric gene delivery systems. In this short communication, we report an effective delivery system based on low molecular polyethylenimine polymer with lipid substitution (PEI-L) capable of delivering three components, (i) a piggyBac transposon plasmid DNA carrying a gene encoding green fluorescence protein (PB-GFP), (ii) a piggyBac transposase plasmid DNA or mRNA, and (iii) a 2 kDa polyacrylic acid as additive for transfection enhancement, all in a single complex. We demonstrate an optimized formulation for stable GFP expression in two model cell lines, MDA-MB-231 and SUM149 recorded till day 108 (3.5 months) and day 43 (1.4 months), respectively, following a single treatment with very low cell number as starting material. Moreover, the stability of the transgene (GFP) expression mediated by piggyBac/PEI-L transposition was retained following three consecutive cryopreservation cycles. The success of this study highlights the feasibility and potential of employing a polymeric delivery system to obtain piggyBac-based stable expression of therapeutic genes.

Knockdown of PGBD5 inhibits the malignant progression of glioma through upregulation of the PPAR pathway

Glioma is the most common type of primary intracranial malignant tumor, and because of its high invasiveness and recurrence, its prognosis remains poor. The present study investigated the biological function of piggyBac transportable element derived 5 (PGBD5) in glioma. Glioma and para-cancerous tissues were obtained from five patients. Reverse transcription-quantitative PCR and western blotting were used to detect the expression levels of PGBD5. Transwell assay and flow cytometry were used to evaluate cell migration, invasion, apoptosis and cell cycle distribution. In addition, a nude mouse tumor transplantation model was established to study the downstream pathways of PGBD5 and the molecular mechanism was analyzed using transcriptome sequencing. The mRNA and protein expression levels of PGBD5 were increased in glioma tissues and cells. Notably, knockdown of PGBD5 in vitro could inhibit the migration and invasion of glioma cells. In addition, the knockdown of PGBD5 expression promoted apoptosis and caused cell cycle arrest in the G2/M phase, thus inhibiting cell proliferation. Furthermore, in vivo experiments revealed that knockdown of PGBD5 expression could inhibit Ki67 expression and slow tumor growth. Changes in PGBD5 expression were also shown to be closely related to the peroxisome proliferator-activated receptor (PPAR) signaling pathway. In conclusion, interference with PGBD5 could inhibit the malignant progression of glioma through the PPAR pathway, suggesting that PGBD5 may be a potential molecular target of glioma.

Flexible TAM requirement of TnpB enables efficient single-nucleotide editing with expanded targeting scope

TnpBs encoded by the IS200/IS605 family transposon are among the most abundant prokaryotic proteins from which type V CRISPR-Cas nucleases may have evolved. Since bacterial TnpBs can be programmed for RNA-guided dsDNA cleavage in the presence of a transposon-adjacent motif (TAM), these nucleases hold immense promise for genome editing. However, the activity and targeting specificity of TnpB in homology-directed gene editing remain unknown. Here we report that a thermophilic archaeal TnpB enables efficient gene editing in the natural host. Interestingly, the TnpB has different TAM requirements for eliciting cell death and for facilitating gene editing. By systematically characterizing TAM variants, we reveal that the TnpB recognizes a broad range of TAM sequences for gene editing including those that do not elicit apparent cell death. Importantly, TnpB shows a very high targeting specificity on targets flanked by a weak TAM. Taking advantage of this feature, we successfully leverage TnpB for efficient single-nucleotide editing with templated repair. The use of different weak TAM sequences not only facilitates more flexible gene editing with increased cell survival, but also greatly expands targeting scopes, and this strategy is probably applicable to diverse CRISPR-Cas systems.

Deciphering cell types by integrating scATAC-seq data with genome sequences

The single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) technology provides insight into gene regulation and epigenetic heterogeneity at single-cell resolution, but cell annotation from scATAC-seq remains challenging due to high dimensionality and extreme sparsity within the data. Existing cell annotation methods mostly focus on the cell peak matrix without fully utilizing the underlying genomic sequence. Here we propose a method, SANGO, for accurate single-cell annotation by integrating genome sequences around the accessibility peaks within scATAC data. The genome sequences of peaks are encoded into low-dimensional embeddings, and then iteratively used to reconstruct the peak statistics of cells through a fully connected network. The learned weights are considered as regulatory modes to represent cells, and utilized to align the query cells and the annotated cells in the reference data through a graph transformer network for cell annotations. SANGO was demonstrated to consistently outperform competing methods on 55 paired scATAC-seq datasets across samples, platforms and tissues. SANGO was also shown to be able to detect unknown tumor cells through attention edge weights learned by the graph transformer. Moreover, from the annotated cells, we found cell-type-specific peaks that provide functional insights/biological signals through expression enrichment analysis, cis-regulatory chromatin interaction analysis and motif enrichment analysis.

Iterative assay for transposase-accessible chromatin by sequencing to isolate functionally relevant neuronal subtypes

The Drosophila brain contains tens of thousands of distinct cell types. Thousands of different transgenic lines reproducibly target specific neuron subsets, yet most still express in several cell types. Furthermore, most lines were developed without a priori knowledge of where the transgenes would be expressed. To aid in the development of cell type-specific tools for neuronal identification and manipulation, we developed an iterative assay for transposase-accessible chromatin (ATAC) approach. Open chromatin regions (OCRs) enriched in neurons, compared to whole bodies, drove transgene expression preferentially in subsets of neurons. A second round of ATAC-seq from these specific neuron subsets revealed additional enriched OCR2s that further restricted transgene expression within the chosen neuron subset. This approach allows for continued refinement of transgene expression, and we used it to identify neurons relevant for sleep behavior. Furthermore, this approach is widely applicable to other cell types and to other organisms.

Childhood cancer mutagenesis caused by transposase-derived PGBD5

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

Mage transposon: a novel gene delivery system for mammalian cells

Transposons, as non-viral integration vectors, provide a secure and efficient method for stable gene delivery. In this study, we have discovered Mage (MG), a novel member of the piggyBac(PB) family, which exhibits strong transposability in a variety of mammalian cells and primary T cells. The wild-type MG showed a weaker insertion preference for near genes, transcription start sites (TSS), CpG islands, and DNaseI hypersensitive sites in comparison to PB, approaching the random insertion pattern. Utilizing in silico virtual screening and feasible combinatorial mutagenesis in vitro, we effectively produced the hyperactive MG transposase (hyMagease). This variant boasts a transposition rate 60% greater than its native counterpart without significantly altering its insertion pattern. Furthermore, we applied the hyMagease to efficiently deliver chimeric antigen receptor (CAR) into T cells, leading to stable high-level expression and inducing significant anti-tumor effects both in vitro and in xenograft mice models. These findings provide a compelling tool for gene transfer research, emphasizing its potential and prospects in the domains of genetic engineering and gene therapy.

Increasing transposase abundance with ocean depth correlates with a particle-associated lifestyle

Transposases are mobile genetic elements that move within and between genomes, promoting genomic plasticity in microorganisms. In marine microbial communities, the abundance of transposases increases with depth, but the reasons behind this trend remain unclear. Our analysis of metagenomes from the Tara Oceans and Malaspina Expeditions suggests that a particle-associated lifestyle is the main covariate for the high occurrence of transposases in the deep ocean, and this trend holds true for individual genomes as well as in a community-wide sense. We observed a strong and depth-independent correlation between transposase abundance and the presence of biofilm-associated genes, as well as the prevalence of secretory enzymes. This suggests that mobile genetic elements readily propagate among microbial communities within crowded biofilms. Furthermore, we show that particle association positively correlates with larger genome size, which is in turn associated with higher transposase abundance. Cassette sequences associated with transposons are enriched with genes related to defense mechanisms, which are more highly expressed in the deep sea. Thus, while transposons spread at the expense of their microbial hosts, they also introduce novel genes and potentially benefit the hosts in helping to compete for limited resources. Overall, our results suggest a new understanding of deep ocean particles as highways for gene sharing among defensively oriented microbial genomes.IMPORTANCEGenes can move within and between microbial genomes via mobile genetic elements, which include transposases and transposons. In the oceans, there is a puzzling increase in transposase abundance in microbial genomes as depth increases. To gain insight into this trend, we conducted an extensive analysis of marine microbial metagenomes and metatranscriptomes. We found a significant correlation between transposase abundance and a particle-associated lifestyle among marine microbes at both the metagenome and genome-resolved levels. We also observed a link between transposase abundance and genes related to defense mechanisms. These results suggest that as microbes become densely packed into crowded particles, mobile genes are more likely to spread and carry genetic material that provides a competitive advantage in crowded habitats. This may enable deep sea microbes to effectively compete in such environments.

Transposition mechanism of ISApl1-the determinant of colistin resistance dissemination

Multidrug-resistant Enterobacteriaceae, a prominent family of gram-negative pathogenic bacteria, causes a wide range of severe diseases. Strains carrying the mobile colistin resistance (mcr-1) gene show resistance to polymyxin, the last line of defense against multidrug-resistant gram-negative bacteria. However, the transmission of mcr-1 is not well understood. In this study, genomes of mcr-1-positive strains were obtained from the NCBI database, revealing their widespread distribution in China. We also showed that ISApl1, a crucial factor in mcr-1 transmission, is capable of self-transposition. Moreover, the self-cyclization of ISApl1 is mediated by its own encoded transposase. The electrophoretic mobility shift assay experiment validated that the transposase can bind to the inverted repeats (IRs) on both ends, facilitating the cyclization of ISApl1. Through knockout or shortening of IRs at both ends of ISApl1, we demonstrated that the cyclization of ISApl1 is dependent on the sequences of the IRs at both ends. Simultaneously, altering the ATCG content of the bases at both ends of ISApl1 can impact the excision rate by modifying the binding ability between IRs and ISAPL1. Finally, we showed that heat-unstable nucleoid protein (HU) can inhibit ISApl1 transposition by binding to the IRs and preventing ISAPL1 binding and expression. In conclusion, the regulation of ISApl1-self-circling is predominantly controlled by the inverted repeat (IR) sequence and the HU protein. This molecular mechanism deepens our comprehension of mcr-1 dissemination.

Bacterial genome engineering using CRISPR-associated transposases

Clustered regularly interspaced short palindromic repeats (CRISPR)-associated transposases have the potential to transform the technology landscape for kilobase-scale genome engineering, by virtue of their ability to integrate large genetic payloads with high accuracy, easy programmability and no requirement for homologous recombination machinery. These transposons encode efficient, CRISPR RNA-guided transposases that execute genomic insertions in Escherichia coli at efficiencies approaching ~100%. Moreover, they generate multiplexed edits when programmed with multiple guides, and function robustly in diverse Gram-negative bacterial species. Here we present a detailed protocol for engineering bacterial genomes using CRISPR-associated transposase (CAST) systems, including guidelines on the available vectors, customization of guide RNAs and DNA payloads, selection of common delivery methods, and genotypic analysis of integration events. We further describe a computational CRISPR RNA design algorithm to avoid potential off-targets, and a CRISPR array cloning pipeline for performing multiplexed DNA insertions. The method presented here allows the isolation of clonal strains containing a novel genomic integration event of interest within 1-2 weeks using available plasmid constructs and standard molecular biology techniques.

Chromosome recombination and modification by LoxP-mediated evolution in Vibrio natriegens using CRISPR-associated transposases

Chromosome rearrangement by LoxP-mediated evolution has emerged as a powerful approach to studying how chromosome architecture impacts phenotypes. However, it relies on the in vitro synthesis of artificial chromosomes. The recently reported CRISPR-associated transposases (CASTs) held great promise for the efficient insertion of abundant LoxP sites directly onto the genome of wild-type strains. In this study, with the fastest-growing bacterium Vibrio natrigens (V. natriegens) as an object, a multiplex genome integration tool derived from CASTs was employed to achieve the insertion of cargo genes at eight specific genomic loci within 2 days. Next, we introduced 30 LoxP sites onto chromosome 2 (Chr2) of V. natriegens. Rigorously induced Cre recombinase was used to demonstrate Chromosome Rearrangement and Modification by LoxP-mediated Evolution (CRaMbLE). Growth characterization and genome sequencing showed that the ~358 kb fragment on Chr2 was accountable for the rapid growth of V. natriegens. The enabling tools we developed can help identify genomic regions that influence the rapid growth of V. natriegens without a prior understanding of genome mechanisms. This groundbreaking demonstration may also be extended to other organisms such as Escherichia coli, Pseudomonas putida, Bacillus subtilis, and so on.

Application of RNA interference and protein localization to investigate housekeeping and developmentally regulated genes in the emerging model protozoan Paramecium caudatum

Unicellular eukaryotes represent tremendous evolutionary diversity. However, the molecular mechanisms underlying this diversity remain largely unexplored, partly due to a limitation of genetic tools to only a few model species. Paramecium caudatum is a well-known unicellular eukaryote with an unexpectedly large germline genome, of which only two percent is retained in the somatic genome following sexual processes, revealing extensive DNA elimination. However, further progress in understanding the molecular mechanisms governing this process is hampered by a lack of suitable genetic tools. Here, we report the successful application of gene knockdown and protein localization methods to interrogate the function of both housekeeping and developmentally regulated genes in P. caudatum. Using these methods, we achieved the expected phenotypes upon RNAi by feeding, and determined the localization of these proteins by microinjection of fusion constructs containing fluorescent protein or antibody tags. Lastly, we used these methods to reveal that P. caudatum PiggyMac, a domesticated piggyBac transposase, is essential for sexual development, and is likely to be an active transposase directly involved in DNA cleavage. The application of these methods lays the groundwork for future studies of gene function in P. caudatum and can be used to answer important biological questions in the future.

AAV-mediated delivery of a Sleeping Beauty transposon and an mRNA-encoded transposase for the engineering of therapeutic immune cells

Engineering cells for adoptive therapy requires overcoming limitations in cell viability and, in the efficiency of transgene delivery, the duration of transgene expression and the stability of genomic integration. Here we report a gene-delivery system consisting of a Sleeping Beauty (SB) transposase encoded into a messenger RNA delivered by an adeno-associated virus (AAV) encoding an SB transposon that includes the desired transgene, for mediating the permanent integration of the transgene. Compared with lentiviral vectors and with the electroporation of plasmids of transposon DNA or minicircle DNA, the gene-delivery system, which we named MAJESTIC (for 'mRNA AAV-SB joint engineering of stable therapeutic immune cells'), offers prolonged transgene expression, as well as higher transgene expression, therapeutic-cell yield and cell viability. MAJESTIC can deliver chimeric antigen receptors (CARs) into T cells (which we show lead to strong anti-tumour activity in vivo) and also transduce natural killer cells, myeloid cells and induced pluripotent stem cells with bi-specific CARs, kill-switch CARs and synthetic T-cell receptors.

Insertion sequence excision is enhanced by a protein that catalyzes branch migration and promotes microhomology-mediated end joining

Insertion sequence (IS)-excision enhancer (IEE) promotes the excision of ISs in the genome of enterohemorrhagic Escherichia coli O157. Because IEE-dependent IS excision occurs in the presence of transposase, the process of IS transposition may be involved in IS excision; however, little is understood about the molecular mechanisms of IS excision. Our in vitro analysis revealed that IEE exhibits DNA-dependent ATPase activity, which is activated by branched DNA. IEE also catalyzes the branch migration of fork-structured DNA. These results suggest that IEE remodels branched structures of the IS transposition intermediate. Sequence analysis of recombination sites in IS-excision products suggested that microhomologous sequences near the ends of the IS are involved in IS excision. IEE promoted microhomology-mediated end joining (MMEJ), in which base pairing between 6-nucleotides complementary ends of two 3'-protruding DNAs and subsequent elongation of the paired DNA strand occurred. IS-excision frequencies were significantly decreased in cells producing IEE mutants that had lost either branch migration or MMEJ activity, which suggests that these activities of IEE are required for IS excision. Based on our results, we propose a model for IS excision triggered by IEE and transposase.

Application of Single-Cell Assay for Transposase-Accessible Chromatin with High Throughput Sequencing in Plant Science: Advances, Technical Challenges, and Prospects

The Single-cell Assay for Transposase-Accessible Chromatin with high throughput sequencing (scATAC-seq) has gained increasing popularity in recent years, allowing for chromatin accessibility to be deciphered and gene regulatory networks (GRNs) to be inferred at single-cell resolution. This cutting-edge technology now enables the genome-wide profiling of chromatin accessibility at the cellular level and the capturing of cell-type-specific cis-regulatory elements (CREs) that are masked by cellular heterogeneity in bulk assays. Additionally, it can also facilitate the identification of rare and new cell types based on differences in chromatin accessibility and the charting of cellular developmental trajectories within lineage-related cell clusters. Due to technical challenges and limitations, the data generated from scATAC-seq exhibit unique features, often characterized by high sparsity and noise, even within the same cell type. To address these challenges, various bioinformatic tools have been developed. Furthermore, the application of scATAC-seq in plant science is still in its infancy, with most research focusing on root tissues and model plant species. In this review, we provide an overview of recent progress in scATAC-seq and its application across various fields. We first conduct scATAC-seq in plant science. Next, we highlight the current challenges of scATAC-seq in plant science and major strategies for cell type annotation. Finally, we outline several future directions to exploit scATAC-seq technologies to address critical challenges in plant science, ranging from plant ENCODE(The Encyclopedia of DNA Elements) project construction to GRN inference, to deepen our understanding of the roles of CREs in plant biology.

Targeted DNA integration in human cells without double-strand breaks using CRISPR-associated transposases

Conventional genome engineering with CRISPR-Cas9 creates double-strand breaks (DSBs) that lead to undesirable byproducts and reduce product purity. Here we report an approach for programmable integration of large DNA sequences in human cells that avoids the generation of DSBs by using Type I-F CRISPR-associated transposases (CASTs). We optimized DNA targeting by the QCascade complex through protein design and developed potent transcriptional activators by exploiting the multi-valent recruitment of the AAA+ ATPase TnsC to genomic sites targeted by QCascade. After initial detection of plasmid-based integration, we screened 15 additional CAST systems from a wide range of bacterial hosts, identified a homolog from Pseudoalteromonas that exhibits improved activity and further increased integration efficiencies. Finally, we discovered that bacterial ClpX enhances genomic integration by multiple orders of magnitude, likely by promoting active disassembly of the post-integration CAST complex, akin to its known role in Mu transposition. Our work highlights the ability to reconstitute complex, multi-component machineries in human cells and establishes a strong foundation to exploit CRISPR-associated transposases for eukaryotic genome engineering.

Dynamics of chromatin accessibility and genome wide control of desiccation tolerance in the resurrection plant Haberlea rhodopensis

BACKGROUND

Drought is one of the main consequences of global climate change and this problem is expected to intensify in the future. Resurrection plants evolved the ability to withstand the negative impact of long periods of almost complete desiccation and to recover at rewatering. In this respect, many physiological, transcriptomic, proteomic and genomic investigations have been performed in recent years, however, few epigenetic control studies have been performed on these valuable desiccation-tolerant plants so far.

RESULTS

In the present study, for the first time for resurrection plants we provide evidences about the differential chromatin accessibility of Haberlea rhodopensis during desiccation stress by ATAC-seq (Assay for Transposase Accessible Chromatin with high-throughput sequencing). Based on gene similarity between species, we used the available genome of the closely related resurrection plant Dorcoceras hygrometricum to identify approximately nine hundred transposase hypersensitive sites (THSs) in H. rhodopensis. The majority of them corresponds to proximal and distal regulatory elements of different genes involved in photosynthesis, carbon metabolism, synthesis of secondary metabolites, cell signalling and transcriptional regulation, cell growth, cell wall, stomata conditioning, chaperons, oxidative stress, autophagy and others. Various types of binding motifs recognized by several families of transcription factors have been enriched from the THSs found in different stages of drought. Further, we used the previously published RNA-seq data from H. rhodopensis to evaluate the expression of transcription factors putatively interacting with the enriched motifs, and the potential correlation between the identified THS and the expression of their corresponding genes.

CONCLUSIONS

These results provide a blueprint for investigating the epigenetic regulation of desiccation tolerance in resurrection plant H. rhodopensis and comparative genomics between resurrection and non-resurrection species with available genome information.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

BACKGROUND

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility.

RESULTS

Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships.

CONCLUSIONS

Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Mechanism of target site selection by type V-K CRISPR-associated transposases

CRISPR-associated transposases (CASTs) repurpose nuclease-deficient CRISPR effectors to catalyze RNA-guided transposition of large genetic payloads. Type V-K CASTs offer potential technology advantages but lack accuracy, and the molecular basis for this drawback has remained elusive. Here, we reveal that type V-K CASTs maintain an RNA-independent, "untargeted" transposition pathway alongside RNA-dependent integration, driven by the local availability of TnsC filaments. Using cryo-electron microscopy, single-molecule experiments, and high-throughput sequencing, we found that a minimal, CRISPR-less transpososome preferentially directs untargeted integration at AT-rich sites, with additional local specificity imparted by TnsB. By exploiting this knowledge, we suppressed untargeted transposition and increased type V-K CAST specificity up to 98.1% in cells without compromising on-target integration efficiency. These findings will inform further engineering of CAST systems for accurate, kilobase-scale genome engineering applications.

Hiker, a new family of DNA transposons encoding transposases with DD35E motifs, displays a distinct phylogenetic relationship with most known DNA transposon families of IS630-Tc1-mariner (ITm)

DNA transposons play a crucial role in determining the size and structure of eukaryotic genomes. In this study, a new family of IS630-Tc1-mariner (ITm) DNA transposons, named Hiker (HK), was identified. HK is characterized by a DD35E catalytic domain and is distinct from all previously known families of the ITm group. Phylogenetic analyses showed that DD35E/Hiker forms a monophyletic clade with DD34E/Gambol, indicating that they may represent a separate superfamily of ITm. A total of 178 Hiker species were identified, with 170 found mainly in Actinopterygii, one in Chondrichthyes, six in Anura and one in Mollusca. Gambol (GM), on the other hand, are found in invertebrates, with 18 in Arthropoda and one in Platyhelminthes. Hiker transposons have a total length ranging from 2.14 to 3.67 kb and contain a single open reading frame that encodes a protein of approximately 370 amino acids (range 311-413 aa). They are flanked by short terminal inverted repeats (TIRs) of 16-30 base pairs and two base pair (TA) target-site duplications. In contrast, most transposons of the Gambol family have a total length of 1.35-5.96 kb, encode a transposase protein of approximately 350 amino acids (range 306-374 aa), and are flanked by TIRs that range from 32 to 1097 bp in length. Both Hiker and Gambol transposases have several conserved motifs, including helix-turn-helix (HTH) motifs and a DDE domain. Our study observed multiple amplification waves and repeated horizontal transfer (HT) events of HK transposons in vertebrate genomes, indicating their role in diversifying and shaping the genomes of Actinopterygii, Chondrichthyes, and Anura. Conversely, GM transposons showed few Horizontal transfer events. According to cell-based transposition assays, most HK transposons are likely inactive due to the truncated DNA binding domains of their transposases. We present an updated classification of the ITm group based on these findings, which will enhance the understanding of both the evolution of ITm transposons and that of their hosts.

Improved piggyBac Transformation with Capped Transposase mRNA in Pest Insects

Creating transgenic insects is a key technology in insect genetics and molecular biology. A widely used instrument in insect transgenesis is the piggyBac transposase, resulting in essentially random genomic integrations. In contrast, site-specific recombinases allow the targeted integration of the transgene construct into a specific genomic target site. Both strategies, however, often face limitations due to low transgenesis efficiencies. We aimed to enhance transgenesis efficiencies by utilizing capped mRNA as a source of transposase or recombinase instead of a helper plasmid. A systematic comparison of transgenesis efficiencies in Aedes mosquitoes, as models for hard-to-transform insects, showed that suppling piggyBac transposase as mRNA increased the average transformation efficiency in Aedes aegypti from less than 5% with the plasmid source to about 50% with mRNA. Similar high activity was observed in Ae. albopictus with pBac mRNA. No efficiency differences between plasmid and mRNA were observed in recombination experiments. Furthermore, a hyperactive version of piggyBac transposase delivered as a plasmid did not improve the transformation efficiency in Ae. aegypti or the agricultural pest Drosophila suzukii. We believe that the use of mRNA has strong potential for enhancing piggyBac transformation efficiencies in other mosquitoes and important agricultural pests, such as tephritids.

Transposon-encoded nucleases use guide RNAs to promote their selfish spread

Insertion sequences are compact and pervasive transposable elements found in bacteria, which encode only the genes necessary for their mobilization and maintenance1. IS200- and IS605-family transposons undergo 'peel-and-paste' transposition catalysed by a TnpA transposase2, but they also encode diverse, TnpB- and IscB-family proteins that are evolutionarily related to the CRISPR-associated effectors Cas12 and Cas9, respectively3,4. Recent studies have demonstrated that TnpB and IscB function as RNA-guided DNA endonucleases5,6, but the broader biological role of this activity has remained enigmatic. Here we show that TnpB and IscB are essential to prevent permanent transposon loss as a consequence of the TnpA transposition mechanism. We selected a family of related insertion sequences from Geobacillus stearothermophilus that encode several TnpB and IscB orthologues, and showed that a single TnpA transposase was broadly active for transposon mobilization. The donor joints formed upon religation of transposon-flanking sequences were efficiently targeted for cleavage by RNA-guided TnpB and IscB nucleases, and co-expression of TnpB and TnpA led to substantially greater transposon retention relative to conditions in which TnpA was expressed alone. Notably, TnpA and TnpB also stimulated recombination frequencies, surpassing rates observed with TnpB alone. Collectively, this study reveals that RNA-guided DNA cleavage arose as a primal biochemical activity to bias the selfish inheritance and spread of transposable elements, which was later co-opted during the evolution of CRISPR-Cas adaptive immunity for antiviral defence.

Genomic features of recombinant CHO clones arising from transposon-based and randomized integration

The use of transposase in cell line development (CLD) programs has experienced increased popularity over the past decade. However, few studies have described the mechanism of action and the genomic and phenotypic characteristics of clones derived from transposase. Additionally, how these traits impact long-term bioproduction is unknown. Here, we use chromosome painting, deep sequencing, and ddPCR to characterize the unique fingerprints associated with transposase-derived clones. Transposase reduces the cellular pool of transient vector as early as three days post transfection following transfection and expedites stable pool establishment by up to two weeks. Furthermore, recombinant DNA expression is significantly improved up to ∼3 fold along with a greater balance of antibody heavy and light chain transcripts, resulting in higher titers in transposase generated pools. Transposase derived pools contained an often innumerable number of integration sites, representing a vast increase in integration site diversity over randomly generated pools, which were bottlenecked at 1-3 integration sites per pool. These transposase mediated integrations typically occurred in clean singlets, free of genomic scars such as deletions, inversions, and other modifications associated with legacy transfection methods which exhibited higher copy numbers per integration site. Relative declines in gene expression occur with copy number increase in the randomly generated, but not the transposase derived clones. Furthermore, transposase-derived clones were more likely to exhibit enhanced a long term stability profile, including product quality attributes such as mannose-5. This improved stability may result from circumventing mechanisms associated with the silencing of tandem repeats. Thus, transposase-mediated approaches can provide multifaceted molecular and phenotypic advantages in cell line development when compared to legacy random-integration methods.

Independent evolution of transposase and TIRs facilitated by recombination between Mutator transposons from divergent clades in maize

Nearly all eukaryotes carry DNA transposons of the Robertson's Mutator (Mu) superfamily, a widespread source of genome instability and genetic variation. Despite their pervasive impact on host genomes, much remains unknown about the evolution of these transposons. Transposase recognition of terminal inverted repeats (TIRs) is thought to drive and constrain coevolution of MuDR transposase genes and TIRs. To address the extent of this relationship and its impact, we compared separate phylogenies of TIRs and MuDR gene sequences from Mu elements in the maize genome. Five major clades were identified. As expected, most Mu elements were bound by highly similar TIRs from the same clade (homomorphic type). However, a subset of elements contained dissimilar TIRs derived from divergent clades. These "heteromorphs" typically occurred in multiple copies indicating active transposition in the genome. In addition, analysis of internal sequences showed that exchanges between elements having divergent TIRs produced new mudra and mudrb gene combinations. In several instances, TIR homomorphs had been regenerated within a heteromorph clade with retention of distinctive internal MuDR sequence combinations. Results reveal that recombination between divergent clades facilitates independent evolution of transposase (mudra), transposase-binding targets (TIRs), and capacity for insertion (mudrb) of active Mu elements. This mechanism would be enhanced by the preference of Mu insertions for recombination-rich regions near the 5' ends of genes. We suggest that cycles of recombination give rise to alternating homo- and heteromorph forms that enhance the diversity on which selection for Mu fitness can operate.

Solid-phase capture and profiling of open chromatin by spatial ATAC

Current methods for epigenomic profiling are limited in their ability to obtain genome-wide information with spatial resolution. We introduce spatial ATAC, a method that integrates transposase-accessible chromatin profiling in tissue sections with barcoded solid-phase capture to perform spatially resolved epigenomics. We show that spatial ATAC enables the discovery of the regulatory programs underlying spatial gene expression during mouse organogenesis, lineage differentiation and in human pathology.

A subclass of the IS1202 family of bacterial insertion sequences targets XerCD recombination sites

We describe here a new family of IS which are related to IS1202, originally isolated from Streptococcus pneumoniae in the mid-1990s and previously tagged as an emerging IS family in the ISfinder database. Members of this family have impacted some important properties of their hosts. We describe here another potentially important property of certain family members: specific targeting of xrs recombination sites. The family could be divided into three subgroups based on their transposase sequences and the length on the target repeats (DR) they generate on insertion: subgroup IS1202 (24-29 bp); ISTde1 (15-18 bp); and ISAba32 (5-6 bp). Members of the ISAba32 subgroup were repeatedly found abutting Xer recombinase recombination sites (xrs), separated by an intervening copy of a DR. These xrs sites, present in multiple copies in a number of Acinetobacter plasmids flanking antibiotic resistance genes, were proposed to form a new type of mobile genetic element using the chromosomally-encoded XerCD recombinase for mobility. Transposase alignments identified subgroup-specific indels which may be responsible for the differences in the transposition properties of the three subgroups (i.e. DR length and target specificity). We propose that this collection of IS be classed as a new insertion sequence family: the IS1202 family composed of three subgroups, only one of which specifically targets plasmid-borne xrs. We discuss the implications of xrs targeting for gene mobility.

AFM-based force spectroscopy unravels stepwise formation of the DNA transposition complex in the widespread Tn3 family mobile genetic elements

Transposon Tn4430 belongs to a widespread family of bacterial transposons, the Tn3 family, which plays a prevalent role in the dissemination of antibiotic resistance among pathogens. Despite recent data on the structural architecture of the transposition complex, the molecular mechanisms underlying the replicative transposition of these elements are still poorly understood. Here, we use force-distance curve-based atomic force microscopy to probe the binding of the TnpA transposase of Tn4430 to DNA molecules containing one or two transposon ends and to extract the thermodynamic and kinetic parameters of transposition complex assembly. Comparing wild-type TnpA with previously isolated deregulated TnpA mutants supports a stepwise pathway for transposition complex formation and activation during which TnpA first binds as a dimer to a single transposon end and then undergoes a structural transition that enables it to bind the second end cooperatively and to become activated for transposition catalysis, the latter step occurring at a much faster rate for the TnpA mutants. Our study thus provides an unprecedented approach to probe the dynamic of a complex DNA processing machinery at the single-particle level.

Optimization of piggyBac transposon-mediated gene transfer method in common marmoset embryos

Generating non-human primate models of human diseases is important for the development of therapeutic strategies especially for neurodegenerative diseases. The common marmoset has attracted attention as a new experimental animal model, and many transgenic marmosets have been produced using lentiviral vector-mediated transgenesis. However, lentiviral vectors have a size limitation of up to 8 kb in length for transgene applications. Therefore, the present study aimed to optimize a piggyBac transposon-mediated gene transfer method in which transgenes longer than 8 kb were injected into the perivitelline space of marmoset embryos, followed by electroporation. We constructed a long piggyBac vector carrying the gene responsible for Alzheimer's disease. The optimal weight ratio of the piggyBac transgene vector to the piggyBac transposase mRNA was examined using mouse embryos. Transgene integration into the genome was confirmed in 70.7% of embryonic stem cells established from embryos injected with 1000 ng of transgene and transposase mRNA. Under these conditions, long transgenes were introduced into marmoset embryos. All embryos survived after transgene introduction treatment, and transgenes were detected in 70% of marmoset embryos. The transposon-mediated gene transfer method developed in this study can be applied to the genetic modification of non-human primates, as well as large animals.

Novel molecular requirements for CRISPR RNA-guided transposition

CRISPR-associated transposases (CASTs) direct DNA integration downstream of target sites using the RNA-guided DNA binding activity of nuclease-deficient CRISPR-Cas systems. Transposition relies on several key protein-protein and protein-DNA interactions, but little is known about the explicit sequence requirements governing efficient transposon DNA integration activity. Here, we exploit pooled library screening and high-throughput sequencing to reveal novel sequence determinants during transposition by the Type I-F Vibrio cholerae CAST system (VchCAST). On the donor DNA, large transposon end libraries revealed binding site nucleotide preferences for the TnsB transposase, as well as an additional conserved region that encoded a consensus binding site for integration host factor (IHF). Remarkably, we found that VchCAST requires IHF for efficient transposition, thus revealing a novel cellular factor involved in CRISPR-associated transpososome assembly. On the target DNA, we uncovered preferred sequence motifs at the integration site that explained previously observed heterogeneity with single-base pair resolution. Finally, we exploited our library data to design modified transposon variants that enable in-frame protein tagging. Collectively, our results provide new clues about the assembly and architecture of the paired-end complex formed between TnsB and the transposon DNA, and inform the design of custom payload sequences for genome engineering applications with CAST systems.

GTAC enables parallel genotyping of multiple genomic loci with chromatin accessibility profiling in single cells

Understanding clonal evolution and cancer development requires experimental approaches for characterizing the consequences of somatic mutations on gene regulation. However, no methods currently exist that efficiently link high-content chromatin accessibility with high-confidence genotyping in single cells. To address this, we developed Genotyping with the Assay for Transposase-Accessible Chromatin (GTAC), enabling accurate mutation detection at multiple amplified loci, coupled with robust chromatin accessibility readout. We applied GTAC to primary acute myeloid leukemia, obtaining high-quality chromatin accessibility profiles and clonal identities for multiple mutations in 88% of cells. We traced chromatin variation throughout clonal evolution, showing the restriction of different clones to distinct differentiation stages. Furthermore, we identified switches in transcription factor motif accessibility associated with a specific combination of driver mutations, which biased transformed progenitors toward a leukemia stem cell-like chromatin state. GTAC is a powerful tool to study clonal heterogeneity across a wide spectrum of pre-malignant and neoplastic conditions.

Genome sequence of Vibrio anguillarum isolates carrying a novel class A β-lactamase VAN-1: do migratory fish transport novel resistance factors?

OBJECTIVES

The aim of the study was to understand the genetic basis of resistance of five β-lactam resistant Vibrio anguillarum isolates obtained from the gut content of Atlantic mackerel (Scomber scomberus), using whole genome sequencing and to characterize a novel β-lactamase (VAN-1) from these isolates.

METHOD

Antibiotic sensitivity pattern was determined using Sensititre™ plates and whole genome sequencing was carried out using Illumina MiSeq-based sequencing. The bla_VAN-1 gene was synthesized and expressed in Escherichia coli Top10 cells.

RESULTS

Five isolates obtained (out of 73) from the gut content of Atlantic mackerel were identified as Vibrio anguillarum. Whole genome assemblies ranged from 3.894 to 3.906 million bases in length with an average of 50 contigs. A novel β-lactamase bla_VAN-1, sharing 77.7% nucleotide identity with a known mobile β-lactamase from Vibrio species was detected. The bla_VAN-1 gene in these isolates is flanked by a truncated IS5 family transposase on one end and a hypothetical protein and outer membrane protein followed by another IS5 family transposase on the other end, suggesting its potential for mobility. The bla_VAN-1 gene was absent in V. anguillarum type strain (ATCC 14181) and V. anguillarum isolates from bivalves and sea water in Norway. VAN-1 conferred ampicillin resistance when expressed in E. coli, thus confirming the functionality of this gene.

CONCLUSIONS

Our study highlights the importance of the marine environment as a reservoir of new antibiotic resistance genes. Our results suggest that migratory fish may transport novel antibiotic resistance determinants over long distances.

Structural basis for topological regulation of Tn3 resolvase

Site-specific DNA recombinases play a variety of biological roles, often related to the dissemination of antibiotic resistance, and are also useful synthetic biology tools. The simplest site-specific recombination systems will recombine any two cognate sites regardless of context. Other systems have evolved elaborate mechanisms, often sensing DNA topology, to ensure that only one of multiple possible recombination products is produced. The closely related resolvases from the Tn3 and γδ transposons have historically served as paradigms for the regulation of recombinase activity by DNA topology. However, despite many proposals, models of the multi-subunit protein-DNA complex (termed the synaptosome) that enforces this regulation have been unsatisfying due to a lack of experimental constraints and incomplete concordance with experimental data. Here, we present new structural and biochemical data that lead to a new, detailed model of the Tn3 synaptosome, and discuss how it harnesses DNA topology to regulate the enzymatic activity of the recombinase.

Phylogenetic Relationships among TnpB-Containing Mobile Elements in Six Bacterial Species

Some families of mobile elements in bacterial genomes encode not only a transposase but also an accessory TnpB gene. This gene has been shown to encode an RNA-guided DNA endonuclease, co-evolving with Y1 transposase and serine recombinase in mobile elements IS605 and IS607. In this paper, we reveal the evolutionary relationships among TnpB-containing mobile elements (TCMEs) in well-assembled genomes of six bacterial species: Bacillus cereus, Clostridioides difficile, Deinococcus radiodurans, Escherichia coli, Helicobacter pylori and Salmonella enterica. In total, 9996 TCMEs were identified in 4594 genomes. They belonged to 39 different insertion sequences (ISs). Based on their genetic structures and sequence identities, the 39 TCMEs were classified into three main groups and six subgroups. According to our phylogenetic analysis, TnpBs include two main branches (TnpB-A and TnpB-B) and two minor branches (TnpB-C and TnpB-D). The key TnpB motifs and the associated Y1 and serine recombinases were highly conserved across species, even though their overall sequence identities were low. Substantial variation was observed for the rate of invasion across bacterial species and strains. Over 80% of the genomes of B. cereus, C. difficile, D. radiodurans and E. coli contained TCMEs; however, only 64% of the genomes of H. pylori and 44% of S. enterica genomes contained TCMEs. IS605 showed the largest rate of invasion in these species, while IS607 and IS1341 had a relatively narrow distribution. Co-invasions of IS605, IS607 and IS1341 elements were observed in various genomes. The largest average copy number was observed for IS605b elements in C. difficile. The average copy numbers of most other TCMEs were smaller than four. Our findings have important implications for understanding the co-evolution of TnpB-containing mobile elements and their biological roles in host genome evolution.

Imaging Chromatin Accessibility by Assay of Transposase-Accessible Chromatin with Visualization

Chromatin accessibility is one of the fundamental structures regulating genome functions including transcription and DNA repair. Recent technological advantages to analyze chromatin accessibility begun to explore the dynamics of local chromatin structures. Here I describe protocols for Assay of Transposase-Accessible Chromatin with Visualization (ATAC-see), which allows us to analyze subnuclear localization of accessible chromatin and quantify accessible chromatin at single-cell level.

ATAC-Seq Analysis of Accessible Chromatin: From Experimental Steps to Data Analysis

Accessible chromatin often represents gene regulatory elements, including promoters and enhancers, essential for gene expression. Assay for Transposase Accessible Chromatin sequencing (ATAC-seq) is one of the most popular techniques to investigate chromatin accessibility across the genome. Here we describe, step by step, a series of optimized experimental methods and bioinformatics pipelines for ATAC-seq analysis. As an example, we present an analysis of murine spermatogenic cells: a method to isolate germ cells, a reaction step using Tn5 transposase to insert sequencing adapters into accessible DNA, a library preparation method for high-throughput sequencing, and bioinformatics analysis of sequencing data. Overall, we introduce a framework of ATAC-seq analysis that can be applied to any cell population to identify cell-type-specific gene regulatory elements and their cis-regulatory networks.

Structural basis for the assembly of the type V CRISPR-associated transposon complex

CRISPR-Cas systems have been co-opted by Tn7-like transposable elements to direct RNA-guided transposition. Type V-K CRISPR-associated transposons rely on the concerted activities of the pseudonuclease Cas12k, the AAA+ ATPase TnsC, the Zn-finger protein TniQ, and the transposase TnsB. Here we present a cryo-electron microscopic structure of a target DNA-bound Cas12k-transposon recruitment complex comprised of RNA-guided Cas12k, TniQ, a polymeric TnsC filament and, unexpectedly, the ribosomal protein S15. Complex assembly, mediated by a network of interactions involving the guide RNA, TniQ, and S15, results in R-loop completion. TniQ contacts two TnsC protomers at the Cas12k-proximal filament end, likely nucleating its polymerization. Transposition activity assays corroborate our structural findings, implying that S15 is a bona fide component of the type V crRNA-guided transposon machinery. Altogether, our work uncovers key mechanistic aspects underpinning RNA-mediated assembly of CRISPR-associated transposons to guide their development as programmable tools for site-specific insertion of large DNA payloads.

Para-chloro-meta-xylenol reshaped the fates of antibiotic resistance genes during sludge fermentation: Insights of cell membrane permeability, bacterial structure and biological pathways

The occurrence of para-chloro-meta-xylenol (PCMX, as largely consumed antimicrobial chemicals) in waste activated sludge (WAS) would pose environmental risks for WAS utilization. This study revealed that PCMX principally prompted the abundances and diversity of antibiotic resistance genes (ARGs), particularly for the multidrug- genes (i.e., acrB and mexW), and reshaped the resistance mechanism categories during WAS fermentation process. The genotype and phenotype results indicated that PCMX upregulated abundances of transposase and increased cell permeability via disrupting WAS structure, which further facilitated the horizontal transfer of ARGs. The network and correlation analysis among ARGs, mobile genetic elements (MGEs) and genera (i.e., Sphingopyxis and Pseudoxanthomonas) verified that PCMX enriched the potential ARGs hosts associated with multidrug resistance mechanism. Also, PCMX upregulated the genes involved in ARGs-associated metabolic pathways, such as two-component (i.e., phoP and vcaM) and quorum sensing systems (i.e., lasR and cciR), which determined the ARGs proliferation via multidrug efflux pump and outer membrane proteins, and facilitated the recognition between ARGs hosts. Variance partitioning analysis (VPA) implied that the shift of microbial community contributed predominantly to the dissemination of ARGs. These findings unveiled the environmental behaviors and risks of exogenous pollutants in WAS with insightful understanding, which could guide the WAS utilization for resource recovery.

Microbiome of the freshwater sponge Ephydatia muelleri shares compositional and functional similarities with those of marine sponges

Sponges are known for hosting diverse communities of microbial symbionts, but despite persistent interest in the sponge microbiome, most research has targeted marine sponges; freshwater sponges have been the focus of less than a dozen studies. Here, we used 16 S rRNA gene amplicon sequencing and shotgun metagenomics to characterize the microbiome of the freshwater sponge Ephydatia muelleri and identify potential indicators of sponge-microbe mutualism. Using samples collected from the Sooke, Nanaimo, and Cowichan Rivers on Vancouver Island, British Columbia, we show that the E. muelleri microbiome is distinct from the ambient water and adjacent biofilms and is dominated by Sediminibacterium, Comamonas, and unclassified Rhodospirillales. We also observed phylotype-level differences in sponge microbiome taxonomic composition among different rivers. These differences were not reflected in the ambient water, suggesting that other environmental or host-specific factors may drive the observed geographic variation. Shotgun metagenomes and metagenome-assembled genomes further revealed that freshwater sponge-associated bacteria share many genomic similarities with marine sponge microbiota, including an abundance of defense-related proteins (CRISPR, restriction-modification systems, and transposases) and genes for vitamin B12 production. Overall, our results provide foundational information on the composition and function of freshwater sponge-associated microbes, which represent an important yet underappreciated component of the global sponge microbiome.

Global changes in chromatin accessibility and transcription in growth hormone-secreting pituitary adenoma

PURPOSE

Growth hormone-secreting pituitary adenoma (GHPA) is an insidious disease with persistent hypersecretion of growth hormone and insulin-like growth factor 1, causing increased morbidity and mortality. Previous studies have investigated the transcription of GHPA. However, the gene regulatory landscape has not been fully characterized. The objective of our study was to unravel the changes in chromatin accessibility and transcription in GHPA.

METHODS

Six patients diagnosed with GHPA in the Department of Neurosurgery at Huashan Hospital were enrolled in our study. Primary pituitary adenoma tissues and adjacent normal pituitary specimens with no morphologic abnormalities from these six patients were obtained at surgery. RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) were applied to investigate the underlying relationship between gene expression and chromatin accessibility changes in GHPA.

RESULTS

Totally, 1528 differential expression genes (DEGs) were identified by transcriptomics analyses, including 725 up-regulated and 803 down-regulated. Further, we obtained 64 significantly DEGs including 10 DEGs were elevated and 54 DEGs were negligibly expressed in tumors tissues. The up-regulated DEGs were mainly involved in terms related to synapse formation, nervous system development and secretory pathway. In parallel, 3916 increased and 2895 decreased chromatin-accessible regions were mapped by ATAC-seq. Additionally, the chromatin accessible changes were frequently located adjacent to transcription factor CTCF and Rfx2 binding site.

CONCLUSIONS

Our results are the first to demonstrate the landscape of chromatin accessibility in GHPA, which may contribute to illustrate the underlying transcriptional regulation mechanism of this disease.

The complete and closed genome of the facultative generalist Candidatus Endoriftia persephone from deep-sea hydrothermal vents

The mutualistic interactions between Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone (short Endoriftia) have been extensively researched. However, the closed Endoriftia genome is still lacking. Here, by employing single-molecule real-time sequencing we present the closed chromosomal sequence of Endoriftia. In contrast to theoretical predictions of enlarged and mobile genetic element-rich genomes related to facultative endosymbionts, the closed Endoriftia genome is streamlined with fewer than expected coding sequence regions, insertion-, prophage-sequences and transposase-coding sequences. Automated and manually curated functional analyses indicated that Endoriftia is more versatile regarding sulphur metabolism than previously reported. We identified the presence of two identical rRNA operons and two long CRISPR regions in the closed genome. Additionally, pangenome analyses revealed the presence of three types of secretion systems (II, IV and VI) in the different Endoriftia populations indicating lineage-specific adaptations. The in depth mobilome characterization identified the presence of shared genomic islands in the different Endoriftia drafts and in the closed genome, suggesting that the acquisition of foreign DNA predates the geographical dispersal of the different endosymbiont populations. Finally, we found no evidence of epigenetic regulation in Endoriftia, as revealed by gene screenings and absence of methylated modified base motifs in the genome. As a matter of fact, the restriction-modification system seems to be dysfunctional in Endoriftia, pointing to a higher importance of molecular memory-based immunity against phages via spacer incorporation into CRISPR system. The Endoriftia genome is the first closed tubeworm endosymbiont to date and will be valuable for future gene oriented and evolutionary comparative studies.

Congenital corneal opacities as a new feature in an unusual case of White-Sutton syndrome

A 2-week-old girl presented with bilateral congenital corneal opacities. Additional systemic manifestations included microcephaly, patent foramen ovale, and poor feeding. Patient and parents underwent whole exome sequencing trio analysis that revealed a de novo pathogenic variant in POGZ (p.Val1150GlyfsX8), which is causative of the White-Sutton syndrome. This rare genetic condition is usually associated with intellectual and developmental delay, facial dysmorphism, strabismus, refractive error, and retinal changes. To our knowledge, this is the first reported case of White-Sutton syndrome presenting with congenital corneal opacities.

Single-Cell Chromatin Accessibility Data Combined with GWAS Improves Detection of Relevant Cell Types in 59 Complex Phenotypes

Several disease risk variants reside on non-coding regions of DNA, particularly in open chromatin regions of specific cell types. Identifying the cell types relevant to complex traits through the integration of chromatin accessibility data and genome-wide association studies (GWAS) data can help to elucidate the mechanisms of these traits. In this study, we created a collection of associations between the combinations of chromatin accessibility data (bulk and single-cell) with an array of 201 complex phenotypes. We integrated the GWAS data of these 201 phenotypes with bulk chromatin accessibility data from 137 cell types measured by DNase-I hypersensitive sequencing and found significant results (FDR adjusted p-value ≤ 0.05) for at least one cell type in 21 complex phenotypes, such as atopic dermatitis, Graves’ disease, and body mass index. With the integration of single-cell chromatin accessibility data measured by an assay for transposase-accessible chromatin with high-throughput sequencing (scATAC-seq), taken from 111 adult and 111 fetal cell types, the resolution of association was magnified, enabling the identification of further cell types. This resulted in the identification of significant correlations (FDR adjusted p-value ≤ 0.05) between 15 categories of single-cell subtypes and 59 phenotypes ranging from autoimmune diseases like Graves’ disease to cardiovascular traits like diastolic/systolic blood pressure.

Integrating ATAC-seq and RNA-seq Reveals the Dynamics of Chromatin Accessibility and Gene Expression in Apple Response to Drought

Drought resistance in plants is influenced by multiple signaling pathways that involve various transcription factors, many target genes, and multiple types of epigenetic modifications. Studies on epigenetic modifications of drought focus on DNA methylation and histone modifications, with fewer on chromatin remodeling. Changes in chromatin accessibility can play an important role in abiotic stress in plants by affecting RNA polymerase binding and various regulatory factors. However, the changes in chromatin accessibility during drought in apples are not well understood. In this study, the landscape of chromatin accessibility associated with the gene expression of apple (GL3) under drought conditions was analyzed by Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq. Differential analysis between drought treatment and control identified 23,466 peaks of upregulated chromatin accessibility and 2447 peaks of downregulated accessibility. The drought-induced chromatin accessibility changed genes were mainly enriched in metabolism, stimulus, and binding pathways. By combining results from differential analysis of RNA-seq and ATAC-seq, we identified 240 genes with higher chromatin accessibility and increased gene expression under drought conditions that may play important functions in the drought response process. Among them, a total of nine transcription factor genes were identified, including ATHB7, HAT5, and WRKY26. These transcription factor genes are differentially expressed with different chromatin accessibility motif binding loci that may participate in apple response to drought by regulating downstream genes. Our study provides a reference for chromatin accessibility under drought stress in apples and the results will facilitate subsequent studies on chromatin remodelers and transcription factors.

Hypercompact adenine base editors based on transposase B guided by engineered RNA

Transposon-associated transposase B (TnpB) is deemed an ancestral protein for type V, Cas12 family members, and the closest ancestor to UnCas12f1. Previously, we reported a set of engineered guide RNAs supporting high indel efficiency for Cas12f1 in human cells. Here we suggest a new technology whereby the engineered guide RNAs also manifest high-efficiency programmable endonuclease activity for TnpB. We have termed this technology TaRGET (TnpB-augment RNA-based Genome Editing Technology). Having this feature in mind, we established TnpB-based adenine base editors (ABEs). A Tad-Tad mutant (V106W, D108Q) dimer fused to the C terminus of dTnpB (D354A) showed the highest levels of A-to-G conversion. The limited targetable sites for TaRGET-ABE were expanded with engineered variants of TnpB or optimized deaminases. Delivery of TaRGET-ABE also ensured potent A-to-G conversion rates in mammalian genomes. Collectively, the TaRGET-ABE will contribute to improving precise genome-editing tools that can be delivered by adeno-associated viruses, thereby harnessing the development of clustered regularly interspaced short palindromic repeats (CRISPR)-based gene therapy.

Embryologic Origin Influences Smooth Muscle Cell Phenotypic Modulation Signatures in Murine Marfan Syndrome Aortic Aneurysm

BACKGROUND

Aortic root smooth muscle cells (SMC) develop from both the second heart field (SHF) and neural crest. Disparate responses to disease-causing Fbn1 variants by these lineages are proposed to promote focal aortic root aneurysm formation in Marfan syndrome (MFS), but lineage-stratified SMC analysis in vivo is lacking.

METHODS

We generated SHF lineage-traced MFS mice and performed integrated multiomic (single-cell RNA and assay for transposase-accessible chromatin sequencing) analysis stratified by embryological origin. SMC subtypes were spatially identified via RNA in situ hybridization. Response to TWIST1 overexpression was determined via lentiviral transduction in human aortic SMCs.

RESULTS

Lineage stratification enabled nuanced characterization of aortic root cells. We identified heightened SHF-derived SMC heterogeneity including a subset of Tnnt2 (cardiac troponin T)-expressing cells distinguished by altered proteoglycan expression. MFS aneurysm-associated SMC phenotypic modulation was identified in both SHF-traced and nontraced (neural crest-derived) SMCs; however, transcriptomic responses were distinct between lineages. SHF-derived modulated SMCs overexpressed collagen synthetic genes and small leucine-rich proteoglycans while nontraced SMCs activated chondrogenic genes. These modulated SMCs clustered focally in the aneurysmal aortic root at the region of SHF/neural crest lineage overlap. Integrated RNA-assay for transposase-accessible chromatin analysis identified enriched Twist1 and Smad2/3/4 complex binding motifs in SHF-derived modulated SMCs. TWIST1 overexpression promoted collagen and SLRP gene expression in vitro, suggesting TWIST1 may drive SHF-enriched collagen synthesis in MFS aneurysm.

CONCLUSIONS

SMCs derived from both SHF and neural crest lineages undergo phenotypic modulation in MFS aneurysm but are defined by subtly distinct transcriptional responses. Enhanced TWIST1 transcription factor activity may contribute to enriched collagen synthetic pathways SHF-derived SMCs in MFS.

Mechanistic details of CRISPR-associated transposon recruitment and integration revealed by cryo-EM

CRISPR-associated transposons (CASTs) are Tn7-like elements that are capable of RNA-guided DNA integration. Although structural data are known for nearly all core transposition components, the transposase component, TnsB, remains uncharacterized. Using cryo-electron microscopy (cryo-EM) structure determination, we reveal the conformation of TnsB during transposon integration for the type V-K CAST system from Scytonema hofmanni (ShCAST). Our structure of TnsB is a tetramer, revealing strong mechanistic relationships with the overall architecture of RNaseH transposases/integrases in general, and in particular the MuA transposase from bacteriophage Mu. However, key structural differences in the C-terminal domains indicate that TnsB's tetrameric architecture is stabilized by a different set of protein-protein interactions compared with MuA. We describe the base-specific interactions along the TnsB binding site, which explain how different CAST elements can function on cognate mobile elements independent of one another. We observe that melting of the 5' nontransferred strand of the transposon end is a structural feature stabilized by TnsB and furthermore is crucial for donor-DNA integration. Although not observed in the TnsB strand-transfer complex, the C-terminal end of TnsB serves a crucial role in transposase recruitment to the target site. The C-terminal end of TnsB adopts a short, structured 15-residue "hook" that decorates TnsC filaments. Unlike full-length TnsB, C-terminal fragments do not appear to stimulate filament disassembly using two different assays, suggesting that additional interactions between TnsB and TnsC are required for redistributing TnsC to appropriate targets. The structural information presented here will help guide future work in modifying these important systems as programmable gene integration tools.

A combinatorial indexing strategy for low-cost epigenomic profiling of plant single cells

Understanding how cis-regulatory elements facilitate gene expression is a key question in biology. Recent advances in single-cell genomics have led to the discovery of cell-specific chromatin landscapes that underlie transcription programs in animal models. However, the high equipment and reagent costs of commercial systems limit their applications for many laboratories. In this study, we developed a combinatorial index and dual PCR barcode strategy to profile the Arabidopsis thaliana root single-cell epigenome without any specialized equipment. We generated chromatin accessibility profiles for 13 576 root nuclei with an average of 12 784 unique Tn5 integrations per cell. Integration of the single-cell assay for transposase-accessible chromatin sequencing and RNA sequencing data sets enabled the identification of 24 cell clusters with unique transcription, chromatin, and cis-regulatory signatures. Comparison with single-cell data generated using the commercial microfluidic platform from 10X Genomics revealed that this low-cost combinatorial index method is capable of unbiased identification of cell-type-specific chromatin accessibility. We anticipate that, by removing cost, instrumentation, and other technical obstacles, this method will be a valuable tool for routine investigation of single-cell epigenomes and provide new insights into plant growth and development and plant interactions with the environment.

The Chromatin Accessibility Landscape of Nonalcoholic Fatty Liver Disease Progression

The advent of the assay for transposase-accessible chromatin using sequencing (ATAC-seq) has shown great potential as a leading method for analyzing the genome-wide profiling of chromatin accessibility. A comprehensive reference to the ATAC-seq dataset for disease progression is important for understanding the regulatory specificity caused by genetic or epigenetic changes. In this study, we present a genome-wide chromatin accessibility profile of 44 liver samples spanning the full histological spectrum of nonalcoholic fatty liver disease (NAFLD). We analyzed the ATAC-seq signal enrichment, fragment size distribution, and correlation coefficients according to the histological severity of NAFLD (healthy control vs steatosis vs fibrotic nonalcoholic steatohepatitis), demonstrating the high quality of the dataset. Consequently, 112,303 merged regions (genomic regions containing one or multiple overlapping peak regions) were identified. Additionally, we found differentially accessible regions (DARs) and performed transcription factor binding motif enrichment analysis and de novo motif analysis to determine new biomarker candidates. These data revealed the generegulatory interactions and noncoding factors that can affect NAFLD progression. In summary, our study provides a valuable resource for the human epigenome by applying an advanced approach to facilitate diagnosis and treatment by understanding the non-coding genome of NAFLD.

Harnessing changes in open chromatin determined by ATAC-seq to generate insulin-responsive reporter constructs

BACKGROUND

Gene regulation is critical for proper cellular function. Next-generation sequencing technology has revealed the presence of regulatory networks that regulate gene expression and essential cellular functions. Studies investigating the epigenome have begun to uncover the complex mechanisms regulating transcription. Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is quickly becoming the assay of choice for many epigenomic investigations. However, whether intervention-mediated changes in accessible chromatin determined by ATAC-seq can be harnessed to generate intervention-inducible reporter constructs has not been systematically assayed.

RESULTS

We used the insulin signaling pathway as a model to investigate chromatin regions and gene expression changes using ATAC- and RNA-seq in insulin-treated Drosophila S2 cells. We found correlations between ATAC- and RNA-seq data, especially when stratifying differentially-accessible chromatin regions by annotated feature type. In particular, our data demonstrated a weak but significant correlation between chromatin regions annotated to enhancers (1-2 kb from the transcription start site) and downstream gene expression. We cloned candidate enhancer regions upstream of luciferase and demonstrate insulin-inducibility of several of these reporters.

CONCLUSIONS

Insulin-induced chromatin accessibility determined by ATAC-seq reveals enhancer regions that drive insulin-inducible reporter gene expression.