Technical considerations for cost-effective transposon directed insertion-site sequencing (TraDIS).

Transposon directed insertion-site sequencing (TraDIS), a variant of transposon insertion sequencing commonly known as Tn-Seq, is a high-throughput assay that defines essential bacterial genes across diverse growth conditions. However, the variability between laboratory environments often requires laborious, time-consuming modifications to its protocol. In this technical study, we aimed to refine the protocol by identifying key parameters that can impact the complexity of mutant libraries. Firstly, we discovered that adjusting electroporation parameters including transposome concentration, transposome assembly conditions, and cell densities can significantly improve the recovery of viable mutants for different Escherichia coli strains. Secondly, we found that post-electroporation conditions, such as recovery time and the use of different mediums for selecting mutants may also impact the complexity of viable mutants in the library. Finally, we developed a simplified sequencing library preparation workflow based on a Nextera-TruSeq hybrid design where ~ 80% of sequenced reads correspond to transposon-DNA junctions. The technical improvements presented in our study aim to streamline TraDIS protocols, making this powerful technique more accessible for a wider scientific audience.

Antipsychotic quetiapine alters the mouse fecal resistome by impacting antibiotic efflux, cell membrane, and cell wall synthesis genes.

IMPORTANCE: This study significantly contributes to our understanding of how certain medications can unintentionally contribute to a major global health issue, i.e., antibiotic resistance. Quetiapine, a widely used antipsychotic medication, was found to increase key resistance mechanisms of gut bacteria to antibiotics in mice. Specifically, these data suggest that quetiapine may target elements of the bacterial cell membrane. If similar effects are found in humans, this medicine could unexpectedly make it harder to treat certain infections. This research emphasizes the importance of being mindful about not just antibiotics themselves, but also about other medications that could inadvertently contribute to this problem. Ultimately, these findings underline the necessity for more in-depth research on the broader impact of pharmaceuticals.

Using a modular massively parallel reporter assay to discover context-specific regulatory grammars in type 2 diabetes.

Recent genome-wide association studies have established that most complex disease-associated loci are found in noncoding regions where defining their function is nontrivial. In this study, we leverage a modular massively parallel reporter assay (MPRA) to uncover sequence features linked to context-specific regulatory activity. We screened enhancer activity across a panel of 198-bp fragments spanning over 10k type 2 diabetes- and metabolic trait-associated variants in the 832/13 rat insulinoma cell line, a relevant model of pancreatic beta cells. We explored these fragments' context sensitivity by comparing their activities when placed up-or downstream of a reporter gene, and in combination with either a synthetic housekeeping promoter (SCP1) or a more biologically relevant promoter corresponding to the human insulin gene ( INS ). We identified clear effects of MPRA construct design on measured fragment enhancer activity. Specifically, a subset of fragments (n = 702/11,656) displayed positional bias, evenly distributed across up- and downstream preference. A separate set of fragments exhibited promoter bias (n = 698/11,656), mostly towards the cell-specific INS promoter (73.4%). To identify sequence features associated with promoter preference, we used Lasso regression with 562 genomic annotations and discovered that fragments with INS promoter-biased activity are enriched for HNF1 motifs. HNF1 family transcription factors are key regulators of glucose metabolism disrupted in maturity onset diabetes of the young (MODY), suggesting genetic convergence between rare coding variants that cause MODY and common T2D-associated regulatory variants. We designed a follow-up MPRA containing HNF1 motif-enriched fragments and observed several instances where deletion or mutation of HNF1 motifs disrupted the INS promoter-biased enhancer activity, specifically in the beta cell model but not in a skeletal muscle cell line, another diabetes-relevant cell type. Together, our study suggests that cell-specific regulatory activity is partially influenced by enhancer-promoter compatibility and indicates that careful attention should be paid when designing MPRA libraries to capture context-specific regulatory processes at disease-associated genetic signals.

Metamorphic gene regulation programs in Xenopus tropicalis tadpole brain.

Amphibian metamorphosis is controlled by thyroid hormone (TH), which binds TH receptors (TRs) to regulate gene expression programs that underlie morphogenesis. Gene expression screens using tissues from premetamorphic tadpoles treated with TH identified some TH target genes, but few studies have analyzed genome-wide changes in gene regulation during spontaneous metamorphosis. We analyzed RNA sequencing data at four developmental stages from the beginning to the end of spontaneous metamorphosis, conducted on the neuroendocrine centers of Xenopus tropicalis tadpole brain. We also conducted chromatin immunoprecipitation sequencing (ChIP-seq) for TRs, and we compared gene expression changes during metamorphosis with those induced by exogenous TH. The mRNA levels of 26% of protein coding genes changed during metamorphosis; about half were upregulated and half downregulated. Twenty four percent of genes whose mRNA levels changed during metamorphosis had TR ChIP-seq peaks. Genes involved with neural cell differentiation, cell physiology, synaptogenesis and cell-cell signaling were upregulated, while genes involved with cell cycle, protein synthesis, and neural stem/progenitor cell homeostasis were downregulated. There is a shift from building neural structures early in the metamorphic process, to the differentiation and maturation of neural cells and neural signaling pathways characteristic of the adult frog brain. Only half of the genes modulated by treatment of premetamorphic tadpoles with TH for 16 h changed expression during metamorphosis; these represented 33% of the genes whose mRNA levels changed during metamorphosis. Taken together, our results provide a foundation for understanding the molecular basis for metamorphosis of tadpole brain, and they highlight potential caveats for interpreting gene regulation changes in premetamorphic tadpoles induced by exogenous TH.

Psychotropic Drugs in the Discussion of Antimicrobial-Resistant Microorganisms.

Psychotropic drugs have long been known to possess antimicrobial activity against several groups of microorganisms. Although this property has been extensively studied both alone and when combined with antibiotics against antimicrobial-resistant bacterial and fungal species, relatively little attention has been given to their ability to contribute to the emergence of antimicrobial resistance (AMR). We have recently reported the acquisition of multidrug resistance in Escherichia coli after exposure to gut-relevant concentrations of the antipsychotic quetiapine. Considering these observations, this review attempts to establish if a relationship between psychotropics and AMR in microorganisms has been defined in the scientific literature.

Ceftazidime-avibactam based combinations against carbapenemase producing Klebsiella pneumoniae harboring hypervirulence plasmids.

The combination of carbapenem resistance and hypervirulence in Klebsiella pneumoniae is an emerging and urgent threat due to its potential to resist common antibiotics and cause life-threatening infections in healthy hosts. This study aimed to evaluate the activity of clinically relevant antibiotic regimens against carbapenem-resistant K. pneumoniae with hypervirulence plasmids and to identify pathways associated with antibiotic tolerance using transcriptomics. We studied two carbapenem-resistant K. pneumoniae isolates, CDI694 and CDI231, both harboring hypervirulence plasmids. Time-kill and dynamic one-compartment pharmacokinetic/pharmacodynamic assays were used to assess ceftazidime/avibactam-based therapies. RNAseq was performed following 48 h of antibiotic exposure. Closed genomes of CDI694 and CDI231 were obtained; each isolate harbored carbapenem-resistance and hypervirulence (containing rmpA/rmpA2 and iut genes) plasmids. Ceftazidime/avibactam-based regimens were bactericidal, though both isolates continued to grow in the presence of antibiotics despite no shifts in MIC. Transcriptomic analyses suggested that perturbations to cell respiration, carbohydrate transport, and stress-response pathways contributed to the antibiotic tolerance in CDI231. Genes associated with hypervirulence and antibiotic resistance were not strongly impacted by drug exposure except for ompW, which was significantly downregulated. Treatment of carbapenem-resistant K. pneumoniae harboring hypervirulence plasmids with ceftazidime/avibactam-based regimens may yield a tolerant population due to altered transcription of multiple key pathways.

The Atypical Antipsychotic Quetiapine Promotes Multiple Antibiotic Resistance in Escherichia coli.

Atypical antipsychotic (AAP) medication is a critical tool for treating symptoms of psychiatric disorders. While AAPs primarily target dopamine (D2) and serotonin (5HT2A and 5HT1A) receptors, they also exhibit intrinsic antimicrobial activity as an off-target effect. Because AAPs are often prescribed to patients for many years, a potential risk associated with long-term AAP use is the unintended emergence of bacteria with antimicrobial resistance (AMR). Here, we show that exposure to the AAP quetiapine at estimated gut concentrations promotes AMR in Escherichia coli after 6 weeks. Quetiapine-exposed isolates exhibited an increase in MICs for ampicillin, tetracycline, ceftriaxone, and levofloxacin. By whole-genome sequencing analysis, we identified mutations in genes that confer AMR, including the repressor for the multiple antibiotic resistance mar operon (marR), and real-time reverse transcription-quantitative PCR (RT-qPCR) analysis showed increased levels of marA, acrA, and tolC mRNAs and reduced levels of ompF mRNA in the isolates carrying marR mutations. To determine the contribution of each marR mutation to AMR, we constructed isogenic strains carrying individual mutant marR alleles in the parent background and reevaluated their resistance phenotypes using MIC and RT-qPCR assays. While marR mutations induced robust activity of the mar operon, they resulted in only modest increases in MICs. Interestingly, although these marR mutations did not fully recapitulate the AMR phenotype of the quetiapine-exposed isolates, we show that marR mutations promote growth fitness in the presence of quetiapine. Our findings revealed an important link between the use of AAPs and AMR development in E. coli. IMPORTANCE AAP medication is a cornerstone in the treatment of serious psychiatric disease. The AAPs are known to exhibit antimicrobial activity; therefore, a potential unintended risk of long-term AAP use may be the emergence of AMR, although such risk has received little attention. In this study, we describe the development of multidrug antibiotic resistance in Escherichia coli after 6 weeks of exposure to the AAP quetiapine. Investigation of mutations in the marR gene, which encodes a repressor for the multiple antibiotic resistance (mar) operon, reveals a potential mechanism that increases the fitness of E. coli in the presence of quetiapine. Our findings establish a link between the use of AAPs and AMR development in bacteria.

A Transcription Start Site Map in Human Pancreatic Islets Reveals Functional Regulatory Signatures.

Identifying the tissue-specific molecular signatures of active regulatory elements is critical to understand gene regulatory mechanisms. Here, we identify transcription start sites (TSS) using cap analysis of gene expression (CAGE) across 57 human pancreatic islet samples. We identify 9,954 reproducible CAGE tag clusters (TCs), ∼20% of which are islet specific and occur mostly distal to known gene TSS. We integrated islet CAGE data with histone modification and chromatin accessibility profiles to identify epigenomic signatures of transcription initiation. Using a massively parallel reporter assay, we validated the transcriptional enhancer activity for 2,279 of 3,378 (∼68%) tested islet CAGE elements (5% false discovery rate). TCs within accessible enhancers show higher enrichment to overlap type 2 diabetes genome-wide association study (GWAS) signals than existing islet annotations, which emphasizes the utility of mapping CAGE profiles in disease-relevant tissue. This work provides a high-resolution map of transcriptional initiation in human pancreatic islets with utility for dissecting active enhancers at GWAS loci.

Chromatin information content landscapes inform transcription factor and DNA interactions.

Interactions between transcription factors and chromatin are fundamental to genome organization and regulation and, ultimately, cell state. Here, we use information theory to measure signatures of organized chromatin resulting from transcription factor-chromatin interactions encoded in the patterns of the accessible genome, which we term chromatin information enrichment (CIE). We calculate CIE for hundreds of transcription factor motifs across human samples and identify two classes: low and high CIE. The 10-20% of common and tissue-specific high CIE transcription factor motifs, associate with higher protein-DNA residence time, including different binding site subclasses of the same transcription factor, increased nucleosome phasing, specific protein domains, and the genetic control of both chromatin accessibility and gene expression. These results show that variations in the information encoded in chromatin architecture reflect functional biological variation, with implications for cell state dynamics and memory.

STARRPeaker: uniform processing and accurate identification of STARR-seq active regions.

STARR-seq technology has employed progressively more complex genomic libraries and increased sequencing depths. An issue with the increased complexity and depth is that the coverage in STARR-seq experiments is non-uniform, overdispersed, and often confounded by sequencing biases, such as GC content. Furthermore, STARR-seq readout is confounded by RNA secondary structure and thermodynamic stability. To address these potential confounders, we developed a negative binomial regression framework for uniformly processing STARR-seq data, called STARRPeaker. Moreover, to aid our effort, we generated whole-genome STARR-seq data from the HepG2 and K562 human cell lines and applied STARRPeaker to comprehensively and unbiasedly call enhancers in them.

Thyroid Hormone Induces DNA Demethylation in Xenopus Tadpole Brain.

Thyroid hormone (T3) plays pivotal roles in vertebrate development, acting via nuclear T3 receptors (TRs) that regulate gene transcription by promoting post-translational modifications to histones. Methylation of cytosine residues in deoxyribonucleic acid (DNA) also modulates gene transcription, and our recent finding of predominant DNA demethylation in the brain of Xenopus tadpoles at metamorphosis, a T3-dependent developmental process, caused us to hypothesize that T3 induces these changes in vivo. Treatment of premetamorphic tadpoles with T3 for 24 or 48 hours increased immunoreactivity in several brain regions for the DNA demethylation intermediates 5-hydroxymethylcytosine (5-hmC) and 5-carboxylcytosine, and the methylcytosine dioxygenase ten-eleven translocation 3 (TET3). Thyroid hormone treatment induced locus-specific DNA demethylation in proximity to known T3 response elements within the DNA methyltransferase 3a and Krüppel-like factor 9 genes, analyzed by 5-hmC immunoprecipitation and methylation sensitive restriction enzyme digest. Chromatin-immunoprecipitation (ChIP) assay showed that T3 induced TET3 recruitment to these loci. Furthermore, the messenger ribonucleic acid for several genes encoding DNA demethylation enzymes were induced by T3 in a time-dependent manner in tadpole brain. A TR ChIP-sequencing experiment identified putative TR binding sites at several of these genes, and we provide multiple lines of evidence to support that tet2 contains a bona fide T3 response element. Our findings show that T3 can promote DNA demethylation in developing tadpole brain, in part by promoting TET3 recruitment to discrete genomic regions, and by inducing genes that encode DNA demethylation enzymes.

DNA methylation dynamics underlie metamorphic gene regulation programs in Xenopus tadpole brain.

Methylation of cytosine residues in DNA influences chromatin structure and gene transcription, and its regulation is crucial for brain development. There is mounting evidence that DNA methylation can be modulated by hormone signaling. We analyzed genome-wide changes in DNA methylation and their relationship to gene regulation in the brain of Xenopus tadpoles during metamorphosis, a thyroid hormone-dependent developmental process. We studied the region of the tadpole brain containing neurosecretory neurons that control pituitary hormone secretion, a region that is highly responsive to thyroid hormone action. Using Methylated DNA Capture sequencing (MethylCap-seq) we discovered a diverse landscape of DNA methylation across the tadpole neural cell genome, and pairwise stage comparisons identified several thousand differentially methylated regions (DMRs). During the pre-to pro-metamorphic period, the number of DMRs was lowest (1,163), with demethylation predominating. From pre-metamorphosis to metamorphic climax DMRs nearly doubled (2,204), with methylation predominating. The largest changes in DNA methylation were seen from metamorphic climax to the completion of metamorphosis (2960 DMRs), with 80% of the DMRs representing demethylation. Using RNA sequencing, we found negative correlations between differentially expressed genes and DMRs localized to gene bodies and regions upstream of transcription start sites. DNA demethylation at metamorphosis revealed by MethylCap-seq was corroborated by increased immunoreactivity for the DNA demethylation intermediates 5-hydroxymethylcytosine and 5-carboxymethylcytosine, and the methylcytosine dioxygenase ten eleven translocation 3 that catalyzes DNA demethylation. Our findings show that the genome of tadpole neural cells undergoes significant changes in DNA methylation during metamorphosis, and these changes likely influence chromatin architecture, and gene regulation programs occurring during this developmental period.

Quantification, Dynamic Visualization, and Validation of Bias in ATAC-Seq Data with ataqv.

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) has become the preferred method for mapping chromatin accessibility due to its time and input material efficiency. However, it can be difficult to evaluate data quality and identify sources of technical bias across samples. Here, we present ataqv, a computational toolkit for efficiently measuring, visualizing, and comparing quality control (QC) results across samples and experiments. We use ataqv to analyze 2,009 public ATAC-seq datasets; their QC metrics display a 10-fold range. Tn5 dosage experiments and statistical modeling show that technical variation in the ratio of Tn5 transposase to nuclei and sequencing flowcell density induces systematic bias in ATAC-seq data by changing the enrichment of reads across functional genomic annotations including promoters, enhancers, and transcription-factor-bound regions, with the notable exception of CTCF. ataqv can be integrated into existing computational pipelines and is freely available at https://github.com/ParkerLab/ataqv/.

Genomic annotation of disease-associated variants reveals shared functional contexts.

Variation in non-coding DNA, encompassing gene regulatory regions such as enhancers and promoters, contributes to risk for complex disorders, including type 2 diabetes. While genome-wide association studies have successfully identified hundreds of type 2 diabetes loci throughout the genome, the vast majority of these reside in non-coding DNA, which complicates the process of determining their functional significance and level of priority for further study. Here we review the methods used to experimentally annotate these non-coding variants, to nominate causal variants and to link them to diabetes pathophysiology. In recent years, chromatin profiling, massively parallel sequencing, high-throughput reporter assays and CRISPR gene editing technologies have rapidly become indispensable tools. Rather than treating individual variants in isolation, we discuss the importance of accounting for context, both genetic (such as flanking DNA sequence) and environmental (such as cellular state or environmental exposure). Incorporating these features shows promise in terms of revealing biologically convergent molecular signatures across distant and seemingly unrelated loci. Studying regulatory elements in the proper context will be crucial for interpreting the functional significance of disease-associated variants and applying the resulting knowledge to improve patient care.

Global dynamics of stage-specific transcription factor binding during thymocyte development.

In vertebrates, multiple transcription factors (TFs) bind to gene regulatory elements (promoters, enhancers, and silencers) to execute developmental expression changes. ChIP experiments are often used to identify where TFs bind to regulatory elements in the genome, but the requirement of TF-specific antibodies hampers analyses of tens of TFs at multiple loci. Here we tested whether TF binding predictions using ATAC-seq can be used to infer the identity of TFs that bind to functionally validated enhancers of the Cd4, Cd8, and Gata3 genes in thymocytes. We performed ATAC-seq at four distinct stages of development in mouse thymus, probing the chromatin accessibility landscape in double negative (DN), double positive (DP), CD4 single positive (SP4) and CD8 SP (SP8) thymocytes. Integration of chromatin accessibility with TF motifs genome-wide allowed us to infer stage-specific occupied TF binding sites within known and potentially novel regulatory elements. Our results provide genome-wide stage-specific T cell open chromatin profiles, and allow the identification of candidate TFs that drive thymocyte differentiation at each developmental stage.

Developmental and Thyroid Hormone Regulation of the DNA Methyltransferase 3a Gene in Xenopus Tadpoles.

Thyroid hormone is essential for normal development in vertebrates. In amphibians, T3 controls metamorphosis by inducing tissue-specific gene regulation programs. A hallmark of T3 action is the modification of chromatin structure, which underlies changes in gene transcription. We found that mRNA for the de novo DNA methyltransferase (DNMT) dnmt3a, but not dnmt1, increased in the brain of Xenopus tadpoles during metamorphosis in parallel with plasma [T3]. Addition of T3 to the rearing water caused a time-dependent increase in dnmt3a mRNA in tadpole brain, tail, and hind limb. By analyzing data from a genome-wide analysis of T3 receptor (TR) binding in tadpole tail, we identified several putative T3 response elements (TREs) within the dnmt3a locus. Using in vitro DNA binding, transient transfection-reporter, and chromatin immunoprecipitation assays for TRs, we identified two functional TREs at -7.1 kb and +5.1 kb relative to the dnmt3a transcription start site. Sequence alignment showed that these TREs are conserved between two related frog species, X. laevis and X. tropicalis, but not with amniotes. Our previous findings showed that this gene is directly regulated by liganded TRs in mouse brain, and whereas the two mouse TREs are conserved among Eutherian mammals, they are not conserved in Xenopus species. Thus, although T3 regulation of dnmt3a may be an ancient pathway in vertebrates, the genomic sites responsible for hormone regulation may have diverged or arisen by convergent evolution. We hypothesize that direct T3 regulation of dnmt3a may be an important mechanism for modulating global changes in DNA methylation.

Liganded Thyroid Hormone Receptors Transactivate the DNA Methyltransferase 3a Gene in Mouse Neuronal Cells.

Thyroid hormone (T3) is essential for proper neurological development. The hormone, bound to its receptors, regulates gene transcription in part by modulating posttranslational modifications of histones. Methylation of DNA, which is established by the de novo DNA methyltransferase (DNMT)3a and DNMT3b, and maintained by DNMT1 is another epigenetic modification influencing gene transcription. The expression of Dnmt3a, but not other Dnmt genes, increases in mouse brain in parallel with the postnatal rise in plasma [T3]. We found that treatment of the mouse neuroblastoma cell line Neuro2a[TRß1] with T3 caused rapid induction of Dnmt3a mRNA, which was resistant to protein synthesis inhibition, supporting that it is a direct T3-response gene. Injection of T3 into postnatal day 6 mice increased Dnmt3a mRNA in the brain by 1 hour. Analysis of two chromatin immunoprecipitation-sequencing datasets, and targeted analyses using chromatin immunoprecipitation, transfection-reporter assays, and in vitro DNA binding identified 2 functional T3-response elements (TREs) at the mouse Dnmt3a locus located +30.3 and +49.3 kb from the transcription start site. Thyroid hormone receptors associated with both of these regions in mouse brain chromatin, but with only 1 (+30.3 kb) in Neuro2a[TRß1] cells. Deletion of the +30.3-kb TRE using CRISPR/Cas9 genome editing eliminated or strongly reduced the Dnmt3a mRNA response to T3. Bioinformatics analysis showed that both TREs are highly conserved among eutherian mammals. Thyroid regulation of Dnmt3a may be an evolutionarily conserved mechanism for modulating global changes in DNA methylation during postnatal neurological development.

Deciphering the regulatory logic of an ancient, ultraconserved nuclear receptor enhancer module.

Cooperative, synergistic gene regulation by nuclear hormone receptors can increase sensitivity and amplify cellular responses to hormones. We investigated thyroid hormone (TH) and glucocorticoid (GC) synergy on the Krüppel-like factor 9 (Klf9) gene, which codes for a zinc finger transcription factor involved in development and homeostasis of diverse tissues. We identified regions of the Xenopus and mouse Klf9 genes 5-6 kb upstream of the transcription start sites that supported synergistic transactivation by TH plus GC. Within these regions, we found an orthologous sequence of approximately 180 bp that is highly conserved among tetrapods, but absent in other chordates, and possesses chromatin marks characteristic of an enhancer element. The Xenopus and mouse approximately 180-bp DNA element conferred synergistic transactivation by hormones in transient transfection assays, so we designate this the Klf9 synergy module (KSM). We identified binding sites within the mouse KSM for TH receptor, GC receptor, and nuclear factor κB. TH strongly increased recruitment of liganded GC receptor and serine 5 phosphorylated (initiating) RNA polymerase II to chromatin at the KSM, suggesting a mechanism for transcriptional synergy. The KSM is transcribed to generate long noncoding RNAs, which are also synergistically induced by combined hormone treatment, and the KSM interacts with the Klf9 promoter and a far upstream region through chromosomal looping. Our findings support that the KSM plays a central role in hormone regulation of vertebrate Klf9 genes, it evolved in the tetrapod lineage, and has been maintained by strong stabilizing selection.

Visualizing stimulus convergence in amygdala neurons during associative learning.

A central feature of models of associative memory formation is the reliance on information convergence from pathways responsive to the conditioned stimulus (CS) and unconditioned stimulus (US). In particular, cells receiving coincident input are held to be critical for subsequent plasticity. Yet identification of neurons in the mammalian brain that respond to such coincident inputs during a learning event remains elusive. Here we use Arc cellular compartmental analysis of temporal gene transcription by fluorescence in situ hybridization (catFISH) to locate populations of neurons in the mammalian brain that respond to both the CS and US during training in a one-trial learning task, conditioned taste aversion (CTA). Individual neurons in the basolateral nucleus of the amygdala (BLA) responded to both the CS taste and US drug during conditioning. Coincident activation was not evident, however, when stimulus exposure was altered so as to be ineffective in promoting learning (backward conditioning, latent inhibition). Together, these data provide clear visualization of neurons in the mammalian brain receiving convergent information about the CS and US during acquisition of a learned association.