Is amplification bias consequential in transposon sequencing (TnSeq) assays? A case study with a Staphylococcus aureus TnSeq library subjected to PCR-based and amplification-free enrichment methods

As transposon sequencing (TnSeq) assays have become prolific in the microbiology field, it is of interest to scrutinize their potential drawbacks. TnSeq data consist of millions of nucleotide sequence reads that are generated by PCR amplification of transposon-genomic junctions. Reads mapping to the junctions are enumerated thus providing information on the number of transposon insertion mutations in each individual gene. Here we explore the possibility that PCR amplification of transposon insertions in a TnSeq library skews the results by introducing bias into the detection and/or enumeration of insertions. We compared the detection and frequency of mapped insertions when altering the number of PCR cycles, and when including a nested PCR, in the enrichment step. Additionally, we present nCATRAs - a novel, amplification-free TnSeq method where the insertions are enriched via CRISPR/Cas9-targeted transposon cleavage and subsequent Oxford Nanopore MinION sequencing. nCATRAs achieved 54 and 23% enrichment of the transposons and transposon-genomic junctions, respectively, over background genomic DNA. These PCR-based and PCR-free experiments demonstrate that, overall, PCR amplification does not significantly bias the results of TnSeq insofar as insertions in the majority of genes represented in our library were similarly detected regardless of PCR cycle number and whether or not PCR amplification was employed. However, the detection of a small subset of genes which had been previously described as essential is sensitive to the number of PCR cycles. We conclude that PCR-based enrichment of transposon insertions in a TnSeq assay is reliable, but researchers interested in profiling putative essential genes should carefully weigh the number of amplification cycles employed in their library preparation protocols. In addition, nCATRAs is comparable to traditional PCR-based methods (Kendall’s correlation=0.896–0.897) although the latter remain superior owing to their accessibility and high sequencing depth.

Mapping the genetic landscape of biomineralization in Magnetospirillum magneticum AMB-1 with RB-Tnseq

Magnetotactic bacteria (MTB) are a phylogenetically diverse group of bacteria remarkable for their ability to biomineralize magnetite (Fe3O4) or greigite (Fe3S4) in organelles called magnetosomes. The majority of genes required for magnetosome formation are encoded by a magnetosome gene island (MAI). Here, we conducted random barcoded transposon mutagenesis (RB-TnSeq) in Magnetospirillum magneticum AMB-1 to identify the global genetic requirements for magnetosome formation under different growth conditions. We generated a library of 184,710 unique strains in a wild-type background, generating ~34 mutant strains for each gene. RB-TnSeq also allowed us to determine the essential gene set of AMB-1 under standard laboratory growth conditions. To pinpoint novel genes that are important for magnetosome formation, we subjected the library to magnetic selection screens in varied growth conditions. We compared biomineralization in standard growth conditions to biomineralization in high iron and anaerobic conditions, respectively. Strains with transposon insertions in the MAI gene mamT had an exacerbated biomineralization defect under both high iron and anerobic conditions compared to standard conditions, adding to our knowledge of the role of MamT in magnetosome formation. Mutants in amb4151, a gene outside of the MAI, are more magnetic than wild-type cells under anaerobic conditions. All three of these phenotypes were validated by creating a markerless deletion strain of the gene and evaluating with TEM imaging. Overall, our results indicate that growth conditions affect which genes are required for biomineralization and that some MAI genes may have more nuanced functions than was previously understood.

SETBP1 overexpression acts in the place of class-defining mutations to drive FLT3-ITD–mutant AML

Advances in cancer genomics have revealed genomic classes of acute myeloid leukemia (AML) characterized by class-defining mutations, such as chimeric fusion genes or in genes such as NPM1, MLL, and CEBPA. These class-defining mutations frequently synergize with internal tandem duplications in FLT3 (FLT3-ITDs) to drive leukemogenesis. However, ∼20% of FLT3-ITD–positive AMLs bare no class-defining mutations, and mechanisms of leukemic transformation in these cases are unknown. To identify pathways that drive FLT3-ITD mutant AML in the absence of class-defining mutations, we performed an insertional mutagenesis (IM) screening in Flt3-ITD mice, using Sleeping Beauty transposons. All mice developed acute leukemia (predominantly AML) after a median of 73 days. Analysis of transposon insertions in 38 samples from Flt3-ITD/IM leukemic mice identified recurrent integrations at 22 loci, including Setbp1 (20/38), Ets1 (11/38), Ash1l (8/38), Notch1 (8/38), Erg (7/38), and Runx1 (5/38). Insertions at Setbp1 led exclusively to AML and activated a transcriptional program similar, but not identical, to those of NPM1-mutant and MLL-rearranged AMLs. Guide RNA targeting of Setbp1 was highly detrimental to Flt3ITD/+/Setbp1IM+, but not to Flt3ITD/+/Npm1cA/+, AMLs. Also, analysis of RNA-sequencing data from hundreds of human AMLs revealed that SETBP1 expression is significantly higher in FLT3-ITD AMLs lacking class-defining mutations. These findings propose that SETBP1 overexpression collaborates with FLT3-ITD to drive a subtype of human AML. To identify genetic vulnerabilities of these AMLs, we performed genome-wide CRISPR-Cas9 screening in Flt3ITD/+/Setbp1IM+ AMLs and identified potential therapeutic targets, including Kdm1a, Brd3, Ezh2, and Hmgcr. Our study gives new insights into epigenetic pathways that can drive AMLs lacking class-defining mutations and proposes therapeutic approaches against such cases.

Three common symbiotic ABC-B transporters in Medicago truncatula are regulated by a NIN-independent branch of the symbiosis signalling pathway.

Several ATP-Binding Cassette (ABC) transporters involved in the arbuscular mycorrhizal symbiosis and nodulation have been identified. We describe three previously-unreported ABC subfamily-B transporters, named ABCB for Mycorrhization and Nodulation (AMN1, AMN2, and AMN3), that are expressed early during infection by rhizobia and arbuscular mycorrhizal fungi. These ABCB transporters are strongly expressed in symbiotically infected tissues, including in root hair cells with rhizobial infection threads and arbusculated cells. During nodulation, the expression of these genes is highly induced by rhizobia and purified Nod factors, and was dependent on DMI3, but is not dependent on other known major regulators of infection such as NIN, NSP1, or NSP2. During mycorrhization their expression is dependent on DMI3 and RAM1, but not on NSP1 and NSP2. Therefore, they may be commonly regulated through a distinct branch of the common symbiotic pathway. Mutants with exonic Tnt1-transposon insertions were isolated for all three genes. None of the single or double mutants showed any differences in colonization by either rhizobia or mycorrhizal fungi, but the triple amn1 amn2 amn3 mutant showed an increase in nodule number. Further studies are needed to identify potential substrates of these transporters and understand their roles in these beneficial symbioses.

Isolating an active and inactive CACTA transposon from lettuce color mutants and characterizing their family

Dietary flavonoids play an important role in human nutrition and health. Flavonoid biosynthesis genes have recently been identified in lettuce (Lactuca sativa); however, few mutants have been characterized. We now report the causative mutations in Green Super Lettuce (GSL), a natural light green mutant derived from red cultivar NAR; and GSL-Dark Green (GSL-DG), an olive-green natural derivative of GSL. GSL harbors CACTA 1 (LsC1), a 3.9-kb active nonautonomous CACTA superfamily transposon inserted in the 5′ untranslated region of anthocyanidin synthase (ANS), a gene coding for a key enzyme in anthocyanin biosynthesis. Both terminal inverted repeats (TIRs) of this transposon were intact, enabling somatic excision of the mobile element, which led to the restoration of ANS expression and the accumulation of red anthocyanins in sectors on otherwise green leaves. GSL-DG harbors CACTA 2 (LsC2), a 1.1-kb truncated copy of LsC1 that lacks one of the TIRs, rendering the transposon inactive. RNA-sequencing and reverse transcription quantitative PCR of NAR, GSL, and GSL-DG indicated the relative expression level of ANS was strongly influenced by the transposon insertions. Analysis of flavonoid content indicated leaf cyanidin levels correlated positively with ANS expression. Bioinformatic analysis of the cv Salinas lettuce reference genome led to the discovery and characterization of an LsC1 transposon family with a putative transposon copy number greater than 1,700. Homologs of tnpA and tnpD, the genes encoding two proteins necessary for activation of transposition of CACTA elements, were also identified in the lettuce genome.

Genetic Causes of Non-pathogenic Pseudomonas syringae pv. actinidiae Isolates in Kiwifruit Orchards

Bacterial canker disease has become the largest threat to kiwifruit cultivation and production. A monomorphic subpopulation of Pseudomonas syringae pv. actinidiae biovar 3 (Psa3) is responsible for the pandemic worldwide. Diversity in pathogenicity has been found in the pandemic subpopulation and in other Psa3 subpopulations causing epidemics in China. However, the genetic bases have not yet been elucidated. In this study, 117 Psa3 isolates were identified by Psa- and Psa3-specific primers, and evaluated for pathogenicity. Three isolates G4, G40, and S2 are not pathogenic to kiwifruit and do not elicit hypersensitivity responses (HRs) in non-host Nicotiana benthamiana leaves. Two isolates, G25 and G35, exhibited attenuated HR-eliciting activity in non-host N. benthamiana, but they exhibited greatly and slightly reduced pathogenicity in host plants, respectively. The genomes of the five isolates were sequenced and compared with closely related isolates revealed by MLVA and whole-genome typing methods. The candidate genetic loci responsible for the changes in pathogenicity and HR elicitation, were further evaluated by allele replacement experiments. We found that the three non-pathogenic isolates were formed due to the independent, identical insertion events of ISPsy36 transposon in the hrpR gene, encoding a key regulator of type III secretion system (T3SS) and type III effectors (T3Es). In the symptomatic sample from which G4 was isolated, 27% HR negative isolates were detected. In isolate G25, transposon insertion of ISPsy32 at the non-coding sequence upstream of the hrpR gene was detected, similar to a previously reported low-virulent Psa3 strain M227. In isolate G35, we detected disruptions of T3Es hopBB1-1 and hopBB1-2, which induce HR in N. benthamiana leaves revealed by Agrobacterium tumefaciens infiltration. These phenotype-changed isolates were formed at low frequencies during the course of pathogen infection in host plants, supported by the binding assay of ISPsy32 and the non-coding DNA sequences upstream of the hrpR gene, the co-isolation of the virulent isolates belonging to the same MLVA clade, and the low levels of transcription of the transposon genes. Taken together, in terms of short-term field evolution, transposon insertions in the T3SS-related genes resulted in the formation of non-pathogenic and low-virulent Psa3 isolates.

Assessing the Involvement of Selected Phenotypes of Pseudomonas simiae PICF7 in Olive Root Colonization and Biological Control of Verticillium dahliae

Pseudomonas simiae PICF7 is an indigenous inhabitant of the olive (Olea europaea L.) rhizosphere/root endosphere and an effective biocontrol agent against Verticillium wilt of olive (VWO), caused by the soil-borne fungus Verticillium dahliae. This study aimed to evaluate the potential involvement of selected phenotypes of strain PICF7 in root colonization ability and VWO biocontrol. Therefore, a random transposon-insertion mutant bank of P. simiae PICF7 was screened for the loss of phenotypes likely involved in rhizosphere/soil persistence (copper resistance), root colonization (biofilm formation) and plant growth promotion (phytase activity). Transposon insertions in genes putatively coding for the transcriptional regulator CusR or the chemotaxis protein CheV were found to affect copper resistance, whereas an insertion in fleQ gene putatively encoding a flagellar regulatory protein hampered the ability to form a biofilm. However, these mutants displayed the same antagonistic effect against V. dahliae as the parental strain. Remarkably, two mutants impaired in biofilm formation were never found inside olive roots, whereas their ability to colonize the root exterior and to control VWO remained unaffected. Endophytic colonization of olive roots was unaltered in mutants impaired in copper resistance and phytase production. Results demonstrated that the phenotypes studied were irrelevant for VWO biocontrol.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5′ splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5′ splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

Dimerisation of the PICTS complex via LC8/Cut-up drives co-transcriptional transposon silencing in Drosophila

In animal gonads, the PIWI-interacting RNA (piRNA) pathway guards genome integrity in part through the co-transcriptional gene silencing of transposon insertions. In Drosophila ovaries, piRNA-loaded Piwi detects nascent transposon transcripts and instructs heterochromatin formation through the Panoramix-induced co-transcriptional silencing (PICTS) complex, containing Panoramix, Nxf2 and Nxt1. Here, we report that the highly conserved dynein light chain LC8/Cut-up (Ctp) is an essential component of the PICTS complex. Loss of Ctp results in transposon de-repression and a reduction in repressive chromatin marks specifically at transposon loci. In turn, Ctp can enforce transcriptional silencing when artificially recruited to RNA and DNA reporters. We show that Ctp drives dimerisation of the PICTS complex through its interaction with conserved motifs within Panoramix. Artificial dimerisation of Panoramix bypasses the necessity for its interaction with Ctp, demonstrating that conscription of a protein from a ubiquitous cellular machinery has fulfilled a fundamental requirement for a transposon silencing complex.

Dimerisation of the PICTS complex via LC8/Cut-up drives co-transcriptional transposon silencing in Drosophila

In animal gonads, the PIWI-interacting RNA (piRNA) pathway guards genome integrity in part through the co-transcriptional gene silencing of transposon insertions. In Drosophila ovaries, piRNA-loaded Piwi detects nascent transposon transcripts and instructs heterochromatin formation through the Panoramix-induced co-transcriptional silencing (PICTS) complex, containing Panoramix, Nxf2 and Nxt1. Here, we report that the highly conserved dynein light chain LC8/Cut-up (Ctp) is an essential component of the PICTS complex. Loss of Ctp results in transposon de-repression and a reduction in repressive chromatin marks specifically at transposon loci. In turn, Ctp can enforce transcriptional silencing when artificially recruited to RNA and DNA reporters. We show that Ctp drives dimerisation of the PICTS complex through its interaction with conserved motifs within Panoramix. Artificial dimerisation of Panoramix bypasses the necessity for its interaction with Ctp, demonstrating that conscription of a protein from a ubiquitous cellular machinery has fulfilled a fundamental requirement for a transposon silencing complex.