InMut-finder: a software tool for insertion identification in mutagenesis using Nanopore long reads

Background Biological mutagens (such as transposon) with sequences inserted, play a crucial role to link observed phenotype and genotype in reverse genetic studies. For this reason, accurate and efficient software tools for identifying insertion sites based on the analysis of sequencing reads are desired. Results We developed a bioinformatics tool, a Finder, to identify genome-wide Insertions in Mutagenesis (named as “InMut-Finder”), based on target sequences and flanking sequences from long reads, such as Oxford Nanopore Sequencing. InMut-Finder succeeded in identify > 100 insertion sites in Medicago truncatula and soybean mutants based on sequencing reads of whole-genome DNA or enriched insertion-site DNA fragments. Insertion sites discovered by InMut-Finder were validated by PCR experiments. Conclusion InMut-Finder is a comprehensive and powerful tool for automated insertion detection from Nanopore long reads. The simplicity, efficiency, and flexibility of InMut-Finder make it a valuable tool for functional genomics and forward and reverse genetics. InMut-Finder was implemented with Perl, R, and Shell scripts, which are independent of the OS. The source code and instructions can be accessed at https://github.com/jsg200830/InMut-Finder.

Correlation between Morphological Features of the Anterior Cruciate Ligament: A Quantitative Study Using a Porcine Model

Background: Inadequate restoration of the rotational knee stability is often reported after ACL reconstruction, mainly due to insufficient restoration of the morphology of the ACL. This study aimed to explore the quantitative correlation between morphological features of the ACL, thus to provide useful information for improving anatomical reconstruction techniques and designing artificial ligaments. Methods: Porcine model was used in this study mainly because the morphology and function of porcine ACLs have been reported to be similar to those of the human ACL. 19 porcine knees (1 year, male) were fixed at full extension using 10% formalin and were dissected to expose the ACL. ACL lengths were measured using a caliper. Mid-substances of the ACL were cut and scanned using X-ray microscopy, and the cross-sectional area (CSA) was measured at the isthmus. Margins of direct and indirect bone insertion sites were distinguished and marked. Measurements were performed on digital photographs to obtain the areas of bone insertions. Statistical analysis using nonlinear regression was used to identify potential correlations among the measurements. Results: The results showed large individual differences on all of the measurements (e.g. the CSA at ACL isthmus ranged from 44.7 to 87.2 mm2, when the area of its femoral and tibial insertion sites ranged from 188.1 and 233.6 mm2 to 258.4 and 412.0 mm2, respectively). The CSA at the isthmus was significantly correlated with the total area of the bone insertion sites and the area of tibial insertion. The area of the tibial insertion was significantly correlated with the area of its direct insertion site. In contrast, the area of the femoral insertion was significantly correlated with the area of its indirect insertion site. The area of the indirect tibial insertion showed a weak correlation with the length of ACL, whereas the length of the ACL was not able to predict or be predicted by any other parameters.Conclusions: The CSA at the ACL isthmus might be more suitable for assessing the morphology of the ACL. However, ACL length has little correlation with the CSA of the isthmus or bone insertion sites, thus should be evaluated independently before surgery.

Characterization of an Acute Myeloid Leukemia Murine Model Driven By MLL/AF9: Effect of Retroviral Insertion Sites and Somatic Mutations on Gene Expression

The MLL/AF9 fusion is found in approximately 30% of MLL-rearranged leukemias and has an intermediate prognosis. Genomically well-characterized murine leukemia models enable us to understand leukemogenesis. We generated a retroviral transduction murine bone marrow transplantation leukemia model (MBMTLM) using the MLL/AF9 fusion gene. Fifteen of 20 mice transplanted with syngeneic bone marrow transduced with a MLL/AF9 carrying retrovirus developed leukemia after a median latency of 149 days. Half a million leukemic bone marrow (LBM) cells from two of these primary leukemias, MA03-P and MA86-P, were transplanted into irradiated recipient mice to establish secondary leukemias, MA03-S (n=3) and MA86-S (n=4). Half a million LBM cells from these secondary leukemias were further transplanted into irradiated recipient mice to generate tertiary leukemias, MA03-T (n=3) and MA86-T (n=4). The latency of the leukemias shortened from 141 days in MA03-P to 18 and 22 days in MA03-S and MA03-T, respectively. Similarly, MA86-P had a latency of 98 days, and the latency was reduced to about 28 days in MA986-S and MA986-T. We used retroviral insertion sites (RISs) to track leukemia clones during serial transplantation. We identified 5 RISs in MA03-P. One RIS, RIS#1-03 at chromosome 7:4602500-4609499 accounted for 52.5% of the total RIS-related reads in MA03-P, while the other four RISs were each represented by fewer than 5% of the reads. Only RIS#1-03 was detected in all of the MA03 secondary and tertiary leukemias , indicating that the cells with RIS#1-03 were the dominant clone in MA03 leukemias. Two RISs were detected in MA86-P: RIS#1-86 at chromosome 19:41338500-41341999 and RIS#2-86 at chromosome 10:127106000-127109499 at 46.7% and 2.5%, respectively . RIS#1-986 was contained in the dominant clone as only this RIS was subsequently detected in the secondary and tertiary MA86 leukemias. The relatively long latency to leukemia development in our MLL/AF9 model was most likely due to the requirement of cooperating somatic mutations. We performed whole exome sequencing on DNA from LBM (n=15) and DNA from their corresponding germline (n=2). An average of 4.5 of single nucleotide variants (SNVs) and 11.4 indels affecting protein coding sequences were found in the MA03 family of leukemias (n=7) which, among others, mutated genes involved in tyrosine kinase pathways such as Epha5 and Pik3r1. We identified an average of 14.8 (SNVs) and 0.5 indels per exome in the MA86 leukemias (n=8). Transcription regulator (Brd1) and tumor suppressor genes (Stk11 and Trp53) were affected by somatic changes in the MA86 family. RNA sequencing was performed on LBM (n=15) and healthy bone marrow (HBM) (n=8). Principal component analysis (PCA) on the expression profiles showed that LBM samples clustered together. Differential gene expression analysis identified genes such as Six1, Eya1 and Bcor which had been reported in previous studies to be essential for leukemogenesis in MLL/AF9 murine model. We also observed downregulation of genes such as Gata2, Btg1, Ifitm1, which had been implicated in other types of leukemias. We next investigated the effect of the RISs and somatic mutations on gene expression. RIS#1-903 was in intron 1 of Ppp6r1. A reduction in fragments per kilobase of transcript per Million mapped reads (FPKM) of Ppp6r1 was observed in MA03 family leukemias compared to leukemias of the MA86 family which did not have RIS#1-03 and showed no difference to HBM samples (MA03: 87.71±1.5; MA86: 132.1±5.1; HBM: 77.56±1.7, p< 0.001). We then determined the expression of Tm9sf3 as it is located 600bp away from RIS#1-986. The FPKM of Tm9sf3 was significantly higher in LBM (both of MA903 and MA986 leukemias) than in HBM (LBM: 146.0±12.7; HBM: 64.66±2.8, p<0.001). In MA86 leukemias which all have RIS#1-86, the FPKM of Tm9sf was two fold higher than in MA03 leukemias without RIS#1-86 (MA86: 189.3±4.4; MA03:97.59±1.7, p < 0.001). In contrast, none of the somatic mutations had a significant effect on the expression of any of the mutated genes. In conclusion, we have established a MBMTLM driven by the MLL/AF9 fusion gene. This well-characterized model provides insights to further understand leukemia development and drug testing. Moreover, we demonstrated that RISs can have an impact on gene expression. Future work on whether Ppp6r1 and Tm9sf3 identified by our RIS analysis are drivers in MLL/AF9 leukemias is warranted. Disclosures Browett: MSD: Membership on an entity's Board of Directors or advisory committees; Janssen: Membership on an entity's Board of Directors or advisory committees; AbbVie: Honoraria.

Mapping Transgene Insertion Sites Reveals the α-Cre Transgene Expression in Both Developing Retina and Olfactory Neurons

The Tg(Pax6-cre,GFP)2Pgr (α-Cre) mouse (MGI:3052661) is a commonly used Cre line thought to be retinal-specific. Using targeted locus amplification (TLA), we mapped the insertion site of the transgene, and defined primers useful to deduce zygosity. Further analyses revealed four tandem copies of the transgene. The insertion site mapped to clusters of vomeronasal and olfactory receptor genes. Using R26R and Ai14 Cre reporter mice, we confirmed retinal Cre activity, but also detected expression in olfactory neurons, implicating the influence of neighbouring regulatory elements. RT-PCR and the buried food pellet (BFP) test did not detect any effects of the transgene on flanking genes in the nasal mucosa and retina. Together, these data precisely map α-Cre, show that it does not affect surrounding loci, but reveal previously unanticipated transgene expression in olfactory neurons. The α-Cre mouse can be a valuable tool in both retinal and olfactory research.

Genomic Analysis of Prophages from Klebsiella pneumoniae Clinical Isolates

Klebsiella pneumoniae is an increasing threat to public health and represents one of the most concerning pathogens involved in life-threatening infections. The resistant and virulence determinants are coded by mobile genetic elements which can easily spread between bacteria populations and co-evolve with its genomic host. In this study, we present the full genomic sequences, insertion sites and phylogenetic analysis of 150 prophages found in 40 K. pneumoniae clinical isolates obtained from an outbreak in a Portuguese hospital. All strains harbored at least one prophage and we identified 104 intact prophages (69.3%). The prophage size ranges from 29.7 to 50.6 kbp, coding between 32 and 78 putative genes. The prophage GC content is 51.2%, lower than the average GC content of 57.1% in K. pneumoniae. Complete prophages were classified into three families in the order Caudolovirales: Myoviridae (59.6%), Siphoviridae (38.5%) and Podoviridae (1.9%). In addition, an alignment and phylogenetic analysis revealed nine distinct clusters. Evidence of recombination was detected within the genome of some prophages but, in most cases, proteins involved in viral structure, transcription, replication and regulation (lysogenic/lysis) were maintained. These results support the knowledge that prophages are diverse and widely disseminated in K. pneumoniae genomes, contributing to the evolution of this species and conferring additional phenotypes. Moreover, we identified K. pneumoniae prophages in a set of endolysin genes, which were found to code for proteins with lysozyme activity, cleaving the β-1,4 linkages between N-acetylmuramic acid and N-acetyl-D-glucosamine residues in the peptidoglycan network and thus representing genes with the potential for lysin phage therapy.

Insertion of Exogenous Genes within the ORF1a Coding Region of Porcine Astrovirus

A tagged or reporter astrovirus can be a valuable tool for the analysis of various aspects of the virus life cycle, and to aid in the development of genetically engineered astroviruses as vectors. Here, transposon-mediated insertion mutagenesis was used to insert a 15-nucleotide (nt) sequence into random sites of open reading frame 1a (ORF1a) based on an infectious full-length cDNA clone of porcine astrovirus (PAstV). Five sites in the predicted coiled-coil structures (CC), genome-linked protein (VPg), and hypervariable region (HVR) in ORF1a of the PAstV genome were identified that could tolerate random 15 nt insertions. Incorporation of the commonly used epitope tags, His, Flag, and HA, into four of the five insertion sites permitted the production of infectious viruses and allowed recognition by specifically tagged monoclonal antibodies. The results of immuno-fluorescent assays showed that Flag-tagged ORF1a protein overlapped partially with capsid and ORF2b proteins in the cytoplasm. Improved light-oxygen-voltage (iLOV) gene was also introduced at the insertion sites of CC, VPg, and HVR. Only one viable recombinant reporter PAstV expressing iLOV inserted in HVR was recovered. Biological analysis of the reporter virus showed that it displayed similar growth characteristics, and yet produced less infectious virus particles, when compared with the parental virus. The recombinant virus carrying the iLOV fused with the HVR of ORF1a protein maintained its stability and showed green fluorescence after 15 passages in cell cultures. The resultant fluorescently tagged virus could provide a promising tool for the rapid screening of antiviral drugs as well as allowing the visualization of PAstV infection and replication in living cells.

Effect of Subinhibitory Concentrations of Antibiotics and Disinfectants on ISAba-Mediated Inactivation of Lipooligosaccharide Biosynthesis Genes in Acinetobacter baumannii

Inactivation of the lipooligosaccharide (LOS) biosynthesis genes lpxA, lpxC and lpxD by ISAba insertion elements results in high-level resistance to colistin in A. baumannii. In the present study, we quantify the rate of spontaneous insertional inactivation of LOS biosynthesis genes by ISAba elements in the ATCC 19606-type strain and two multidrug clinical isolates. Using insertional inactivation of lpxC by ISAba11 in the ATCC 19606 strain as a model, we determine the effect of several subinhibitory concentrations of the antibiotics, namely tetracycline, ciprofloxacin, meropenem, kanamycin and rifampicin, as well as the disinfectants ethanol and chlorhexidine on ISAba11 insertion frequencies. Notably, subinhibitory concentrations of tetracycline significantly increased ISAba11 insertion, and rifampicin completely inhibited the emergence of colistin resistance due to ISAba11 inactivation of lpxC. Sequencing of ISAba11 insertion sites within the lpxC gene demonstrated that insertions clustered between nucleotides 382 and 618 (58.3% of unique insertions detected), indicating that this may be a hotspot for ISAba11 insertion. The alignment of insertion sites revealed a semi-conserved AT-rich consensus sequence upstream of the ISAba11 insertion site, suggesting that ISAba11 insertion sites may be sequence-dependent. This study explores previously uncharacterized aspects regarding the acquisition of colistin resistance through insertional activation in LOS biosynthesis genes in A. baumannii.