High-Throughput Genotyping of Green Algal Mutants Reveals Random Distribution of Mutagenic Insertion Sites and Endonucleolytic Cleavage of Transforming DNA

A high-throughput genetic screening platform in a single-celled photosynthetic eukaryote would be a transformative addition to the plant biology toolbox. Here, we present ChlaMmeSeq (Chlamydomonas MmeI-based insertion site Sequencing), a tool for simultaneous mapping of tens of thousands of mutagenic insertion sites in the eukaryotic unicellular green alga Chlamydomonas reinhardtii. We first validated ChlaMmeSeq by in-depth characterization of individual insertion sites. We then applied ChlaMmeSeq to a mutant pool and mapped 11,478 insertions, covering 39% of annotated protein coding genes. We observe that insertions are distributed in a manner largely indistinguishable from random, indicating that mutants in nearly all genes can be obtained efficiently. The data reveal that sequence-specific endonucleolytic activities cleave the transforming DNA and allow us to propose a simple model to explain the origin of the poorly understood exogenous sequences that sometimes surround insertion sites. ChlaMmeSeq is quantitatively reproducible, enabling its use for pooled enrichment screens and for the generation of indexed mutant libraries. Additionally, ChlaMmeSeq allows genotyping of hits from Chlamydomonas screens on an unprecedented scale, opening the door to comprehensive identification of genes with roles in photosynthesis, algal lipid metabolism, the algal carbon-concentrating mechanism, phototaxis, the biogenesis and function of cilia, and other processes for which C. reinhardtii is a leading model system.     

Time- and Cost-Efficient Identification of T-DNA Insertion Sites through Targeted Genomic Sequencing

Forward genetic screens enable the unbiased identification of genes involved in biological processes. In Arabidopsis, several mutant collections are publicly available, which greatly facilitates such practice. Most of these collections were generated by agrotransformation of a T-DNA at random sites in the plant genome. However, precise mapping of T-DNA insertion sites in mutants isolated from such screens is a laborious and time-consuming task. Here we report a simple, low-cost and time efficient approach to precisely map T-DNA insertions simultaneously in many different mutants. By combining sequence capture, next-generation sequencing and 2D-PCR pooling, we developed a new method that allowed the rapid localization of T-DNA insertion sites in 55 out of 64 mutant plants isolated in a screen for gyrase inhibition hypersensitivity.     

An Inverse Pcr Screen for the Detection of P Element Insertions in Cloned Genomic Intervals in Drosophila Melanogaster

We developed a screening approach that utilizes an inverse polymerase chain reaction (PCR) to detect P element insertions in or near previously cloned genes in Drosophila melanogaster. We used this approach in a large scale genetic screen in which P elements were mobilized from sites on the X chromosome to new autosomal locations. Mutagenized flies were combined in pools, and our screening approach was used to generate probes corresponding to the sequences flanking each site of insertion. These probes then were used for hybridization to cloned genomic intervals, allowing individuals carrying insertions in them to be detected. We used the same approach to perform repeated rounds of sib-selection to generate stable insertion lines. We screened 16,100 insert bearing individuals and recovered 11 insertions in five intervals containing genes encoding members of the kinesin superfamily in Drosophila melanogaster. In addition, we recovered an insertion in the region including the Larval Serum Protein-2 gene. Examination by Southern hybridization confirms that the lines we recovered represent genuine insertions in the corresponding genomic intervals. Our data indicates that this approach will be very efficient both for P element mutagenesis of new genomic regions and for detection and recovery of ``local' P element transposition events. In addition, our data constitutes a survey of preferred P element insertion sites in the Drosophila genome and suggests that insertion sites that are mutable at a rate of ~10(-4) are distributed every 40-50 kb.     

Non-random distribution of T-DNA insertions at various levels of the genome hierarchy as revealed by analyzing 13 804 T-DNA flanking sequences from an enhancer-trap mutant library

We isolated 13 804 T-DNA flanking sequence tags (FSTs) from a T-DNA insertion library of rice. A comprehensive analysis of the 13 804 FSTs revealed a number of features demonstrating a highly non-random distribution of the T-DNA insertions in the rice genome: T-DNA insertions were biased towards large chromosomes, not only in the absolute number of insertions but also in the relative density; within chromosomes the insertions occurred more densely in the distal ends, and less densely in the centromeric regions; the distribution of the T-DNA insertions was highly correlated with that of full-length cDNAs, but the correlations were highly heterogeneous among the chromosomes; T-DNA insertions strongly disfavored transposable element (TE)-related sequences, but favored genic sequences with a strong bias toward the 5' upstream and 3' downstream regions of the genes; T-DNA insertions preferentially occurred among the various classes of functional genes, such that the numbers of insertions were in excess in certain functional categories but were deficient in other categories. The analysis of DNA sequence compositions around the T-DNA insertion sites also revealed several prominent features, including an elevated bendability from -200 to 200 bp relative to the insertion sites, an inverse relationship between the GC and TA skews, and reversed GC and TA skews in sequences upstream and downstream of the insertion sites, with both GC and TA skews equal to zero at the insertion sites. It was estimated that 365 380 insertions are needed to saturate the genome with P = 0.95, and that the 45 441 FSTs that have been isolated so far by various groups tagged 14 287 of the 42 653 non-TE related genes.     

Construction of Signature-tagged Mutant Library in Mesorhizobium loti as a Powerful Tool for Functional Genomics

Rhizobia are nitrogen-fixing soil bacteria that establish endosymbiosis with some leguminous plants. The completion of several rhizobial genome sequences provides opportunities for genome-wide functional studies of the physiological roles of many rhizobial genes. In order to carry out genome-wide phenotypic screenings, we have constructed a large mutant library of the nitrogen-fixing symbiotic bacterium, Mesorhizobium loti, by transposon mutagenesis. Transposon insertion mutants were generated using the signature-tagged mutagenesis (STM) technique and a total of 29 330 independent mutants were obtained. Along with the collection of transposon mutants, we have determined the transposon insertion sites for 7892 clones, and confirmed insertions in 3680 non-redundant M. loti genes (50.5% of the total number of M. loti genes). Transposon insertions were randomly distributed throughout the M. loti genome without any bias toward G+C contents of insertion target sites and transposon plasmids used for the mutagenesis. We also show the utility of STM mutants by examining the specificity of signature tags and test screenings for growth- and nodulation-deficient mutants. This defined mutant library allows for genome-wide forward- and reverse-genetic functional studies of M. loti and will serve as an invaluable resource for researchers to further our understanding of rhizobial biology.Key words: Mesorhizobium loti, signature-tagged mutagenesis, mutant library, reverse genetics     

Technical advance: a high throughput system for transposon tagging and promoter trapping in tomato

We describe new tools for functional analysis of the tomato genome based on insertional mutagenesis with the maize Ac/Ds transposable elements in the background of the miniature cultivar Micro-Tom. 2932 F3 families, in which Ds elements transposed and were stabilized, were screened for phenotypic mutations. Out of 10 families that had a clear mutant phenotype, only one mutant was Ds-tagged. In addition, we developed promoter trapping using the firefly luciferase reporter gene and enhancer trapping, using beta-glucuronidase (GUS). We show that luciferase can be used as a non-invasive reporter to identify, isolate and regenerate somatic sectors, to study the time course of mutant expression, and to identify inducible genes. Out of 108 families screened for luciferase activity 55% showed expression in the flower, 11% in the fruit and 4% in seedlings, suggesting a high rate of Ds insertion into genes. Preferential insertion into genes was supported by the analysis of Ds flanking sequences: 28 out of 50 sequenced Ds insertion sites were similar to known genes or to ESTs. In summary, the 2932 lines described here contain 2-3 Ds inserts per line, representing a collection of approximately 7500 Ds insertions. This collection has potential for use in high-throughput functional analysis of genes and promoter isolation in tomato.     

Clonal Expansion Analysis of Transposon Insertions by High-Throughput Sequencing Identifies Candidate Cancer Genes in a PiggyBac Mutagenesis Screen

Somatic transposon mutagenesis in mice is an efficient strategy to investigate the genetic mechanisms of tumorigenesis. The identification of tumor driving transposon insertions traditionally requires the generation of large tumor cohorts to obtain information about common insertion sites. Tumor driving insertions are also characterized by their clonal expansion in tumor tissue, a phenomenon that is facilitated by the slow and evolving transformation process of transposon mutagenesis. We describe here an improved approach for the detection of tumor driving insertions that assesses the clonal expansion of insertions by quantifying the relative proportion of sequence reads obtained in individual tumors. To this end, we have developed a protocol for insertion site sequencing that utilizes acoustic shearing of tumor DNA and Illumina sequencing. We analyzed various solid tumors generated by PiggyBac mutagenesis and for each tumor >10⁶ reads corresponding to >10⁴ insertion sites were obtained. In each tumor, 9 to 25 insertions stood out by their enriched sequence read frequencies when compared to frequencies obtained from tail DNA controls. These enriched insertions are potential clonally expanded tumor driving insertions, and thus identify candidate cancer genes. The candidate cancer genes of our study comprised many established cancer genes, but also novel candidate genes such as Mastermind-like1 (Mamld1) and Diacylglycerolkinase delta (Dgkd). We show that clonal expansion analysis by high-throughput sequencing is a robust approach for the identification of candidate cancer genes in insertional mutagenesis screens on the level of individual tumors.     

Characterization of T-DNA insertion sites in Arabidopsis thaliana and the implications for saturation mutagenesis

A key component of a sound functional genomics infrastructure is the availability of a knockout mutant for every gene in the genome. A fruitful approach to systematically knockingout genes in the plant Arabidopsis thaliana has been the use of transferred-DNA (T-DNA) from Agrobacterium tumefaciens as an insertional mutagen. One of the assumptions underlying the use of T-DNA as a mutagen is that the insertion of these DNA elements into the Arabidopsis genome occurs at randomly selected locations. We have directly investigated the distribution of T-DNA insertions sites in populations of transformed Arabidopsis using two different approaches. To begin with, we utilized a polymerase chain reaction (PCR) procedure to systematically catalog the precise locations of all the T-DNA elements inserted within a 65 kb segment of chromosome IV. Of the 47 T-DNA insertions identified, 30% were found within the coding regions of genes. We also documented the insertion of T-DNA elements within the centromeric region of chromosome IV. In addition to these targeted T-DNA screens, we also mapped the genomic locations of 583 randomly chosen T-DNA elements by sequencing the genomic DNA flanking the insertion sites from individual T-DNA-transformed lines. 35% of these randomly chosen T-DNA insertions were located within the coding regions of genes. For comparison, coding sequences account for 44% of the Arabidopsis genome. Our results demonstrate that there is a small bias towards recovering T-DNA insertions within intergenic regions. However, this bias does not limit the utility of T-DNA as an effective insertional mutagen for use in reverse-genetic strategies.     

An Arabidopsis thaliana T-DNA mutagenized population (GABI-Kat) for flanking sequence tag-based reverse genetics

The GABI-Kat population of T-DNA mutagenized Arabidopsis thaliana lines with sequence-characterized insertion sites is used extensively for efficient progress in plant functional genomics. Here we provide details about the establishment of the material, demonstrate the population's functionality and discuss results from quality control studies. T-DNA insertion mutants of the accession Columbia (Col-0) were created by Agrobacterium tumefaciens-mediated transformation. To allow selection of transformed plants under greenhouse conditions, a sulfadiazine resistance marker was employed. DNA from leaves of T1 plants was extracted and used as a template for PCR-based amplification of DNA fragments spanning insertion site borders. After sequencing, the data were placed in a flanking sequence tag (FST) database describing which mutant allele was present in which line. Analysis of the distribution of T-DNA insertions revealed a clear bias towards intergenic regions. Insertion sites appeared more frequent in regions in front of the ATG and after STOP codons of predicted genes. Segregation analysis for sulfadiazine resistance showed that 62% of the transformants contain an insertion at only one genetic locus. In quality control studies with gene-specific primers in combination with T-DNA primers, 76% of insertions could be confirmed. Finally, the functionality of the GABI-Kat population was demonstrated by exemplary confirmation of several new transparent testa alleles, as well as a number of other mutants, which were identified on the basis of the FST data.     

Genome-wide LORE1 retrotransposon mutagenesis and high-throughput insertion detection in Lotus japonicus

Use of insertion mutants facilitates functional analysis of genes, but it has been difficult to identify a suitable mutagen and to establish large populations for reverse genetics in most plant species. The main challenge is developing efficient high-throughput procedures for both mutagenesis and identification of insertion sites. To date, only floral-dip T-DNA transformation of Arabidopsis has produced independent germinal insertions, thereby allowing generation of mutant populations from seeds of single plants. In addition, advances in insertion detection have been hampered by a lack of protocols, including software for automated data analysis, that take full advantage of high-throughput next-generation sequencing. We have addressed these challenges by developing the FSTpoolit protocol and software package, and here we demonstrate its efficacy by detecting 8935 LORE1 insertions in 3744 Lotus japonicus plants. The identified insertions show that the endogenous LORE1 retrotransposon is well suited for insertion mutagenesis due to homogenous gene targeting and exonic insertion preference. As LORE1 transposition occurs in the germline, harvesting seeds from a single founder line and cultivating progeny generates a complete mutant population. This ease of LORE1 mutagenesis, combined with the efficient FSTpoolit protocol, which exploits 2D pooling, Illumina sequencing and automated data analysis, allows highly cost-efficient development of a comprehensive reverse genetic resource.     

TEAseq-based identification of 35,696 Dissociation insertional mutations facilitates functional genomic studies in maize

In plants, transposable element (TE)-triggered mutants are important resources for functional genomic studies. However, conventional approaches for genome-wide identification of TE insertion sites are costly and laborious. This study developed a novel, rapid, and high-throughput TE insertion site identification workflow based on next-generation sequencing and named it Transposable Element Amplicon Sequencing (TEAseq). Using TEAseq, we systemically profiled the Dissociation (Ds) insertion sites in 1606 independent Ds insertional mutants in advanced backcross generation using K17 as background. The Ac-containing individuals were excluded for getting rid of the potential somatic insertions. We characterized 35,696 germinal Ds insertions tagging 10,323 genes, representing approximately 23.3% of the total genes in the maize genome. The insertion sites were presented in chromosomal hotspots around the ancestral Ds loci, and insertions occurred preferentially in gene body regions. Furthermore, we mapped a loss-of-function AGL2 gene using bulked segregant RNA-sequencing assay and proved that AGL2 is essential for seed development. We additionally established an open-access database named MEILAM for easy access to Ds insertional mutations. Overall, our results have provided an efficient workflow for TE insertion identification and rich sequence-indexed mutant resources for maize functional genomic studies.     

Insertion mutagenesis of the gene encoding the ferrichrome-iron receptor of Escherichia coli K-12.

The ferrichrome-iron receptor of Escherichia coli K-12 encoded by the fhuA gene is a multifunctional outer membrane receptor with an Mr of 78,000. It is required for the binding and uptake of ferrichrome and is the receptor for bacteriophages T5, T1, phi 80, and UC-1 as well as for colicin M. The fhuA gene was cloned into pBR322, and the recombinant plasmid pGC01 was mutagenized by the insertion of 6-base-pair TAB (two amino acid Barany) linkers into CfoI and HpaII restriction sites distributed throughout the coding region. A library of 18 TAB linker insertions in fhuA was generated; 8 of the mutations were at CfoI sites and 10 were at HpaII sites. All mutations inserted a hexamer that encoded a unique SacI site. A large deletion in fhuA was also isolated by TAB linker mutagenesis. Except for the deletion mutant, all of the linker insertion mutant FhuA proteins were found in the outer membrane in amounts similar to those found in the wild type. Five of the linker insertion mutants were susceptible to cleavage by endogenous proteolytic activity: a second FhuA-related band that migrated at approximately 72 kilodaltons could be detected on Coomassie blue-stained gels and on Western blots (immunoblots) by using a carboxy terminus-specific anti-peptide antibody. Receptor functions were measured with the mutated genes present in a single copy on the chromosome. Some of the receptors conferred wild-type phenotypes: they demonstrated growth promotion by ferrichrome and the same efficiency of plating as that of wild-type FhuA; killing by colicin M was also unaffected. Several mutants were altered in their sensitivities to the lethal agents. TAB linker insertions after amino acids 69 and 128 abolished all receptor functions. Phage T5 id not bind to these mutant FhuA proteins in detergent extracts. The deletion mutant was also defective in all FhuA functions. Sensitivity to the lethal agents of cellsl that expressed mutant FhuAs with insertions after amino acids 59 and 135 was reduced by several orders of magnitude. Insertion at other selected sites decreased some or all receptor functions only slightly. An insertion after amino acid 321 selectively eliminated ferrichrome growth promotion. Finally, a strain carrying a mutant fhuA gene on the chromosome in which the linker insertion occurred after amino acid 82 showed a tonB phenotype. These subtle perturbations that were introduced into the FhuA protein resulted in changes in its stability and in the binding and uptake of its cognate ligands.     

ITIS,a bioinformatics tool for accurate identification of transposon insertion sites using next-generation sequencing data

Background

Transposable elements constitute an important part of the genome and are essential in adaptive mechanisms. Transposition events associated with phenotypic changes occur naturally or are induced in insertional mutant populations. Transposon mutagenesis results in multiple random insertions and recovery of most/all the insertions is critical for forward genetics study. Using genome next-generation sequencing data and appropriate bioinformatics tool, it is plausible to accurately identify transposon insertion sites, which could provide candidate causal mutations for desired phenotypes for further functional validation.

Results

We developed a novel bioinformatics tool, ITIS (Identification of Transposon Insertion Sites), for localizing transposon insertion sites within a genome. It takes next-generation genome re-sequencing data (NGS data), transposon sequence, and reference genome sequence as input, and generates a list of highly reliable candidate insertion sites as well as zygosity information of each insertion. Using a simulated dataset and a case study based on an insertional mutant line from Medicago truncatula, we showed that ITIS performed better in terms of sensitivity and specificity than other similar algorithms such as RelocaTE, RetroSeq, TEMP and TIF. With the case study data, we demonstrated the efficiency of ITIS by validating the presence and zygosity of predicted insertion sites of the Tnt1 transposon within a complex plant system, M. truncatula.

Conclusion

This study showed that ITIS is a robust and powerful tool for forward genetic studies in identifying transposable element insertions causing phenotypes. ITIS is suitable in various systems such as cell culture, bacteria, yeast, insect, mammal and plant.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0507-2) contains supplementary material, which is available to authorized users.     

Structural and energetic differences between insertions and substitutions in staphylococcal nuclease.

In a previous study, the small protein staphylococcal nuclease was shown to readily accommodate single alanine and glycine insertions, with average losses in stability comparable to substitutions at the same sites (PROT. 7:299-305, 1990). To more fully explore this unexpected adaptability to changes in residue spacing, 2 double amino acid insertions (alanyl-glycine, glycyl-glycine) and 3 additional single amino acid insertions with dissimilar side chains (proline, leucine, and glutamine) were constructed at 10 of the sites previously studied. At 8 of these sites, the type of amino acid side chain on the inserted residue significantly influenced the stability of the mutant protein. However, at 9 of the 10 sites, the double insertions were found to be no more destabilizing than the single alanine or glycine insertions. In contrast, double substitution mutations of staphylococcal nuclease, which replace two adjacent residues with alanine, do not show this striking degree of non-additivity. A comparison of the effects of single glutamine and single glycine insertions with alanyl-glycine insertions indicates that insertion of alanine into the peptide backbone is, on average, less destabilizing than appending the equivalent atoms onto the side chain of a glycine insertion. To explain their very different energetic effects, we propose that, unlike most substitutions, the inserted residue(s) must induce lateral displacements of the polypeptide chain, forcing the folded conformation away from that of wild type. The resulting obligatory shifts in the positioning of residues flanking the insertion generate a large number of degrees of freedom around which the mutant structure can relax.(ABSTRACT TRUNCATED AT 250 WORDS)     

High-throughput linkage analysis of Mutator insertion sites in maize

Insertional mutagenesis is a cornerstone of functional genomics. High-copy transposable element systems such as Mutator ( Mu ) in maize ( Zea mays ) afford the advantage of high forward mutation rates but pose a challenge for identifying the particular element responsible for a given mutation. Several large mutant collections have been generated in Mu -active genetic stocks, but current methods limit the ability to rapidly identify the causal Mu insertions. Here we present a method to rapidly assay Mu insertions that are genetically linked to a mutation of interest. The method combines elements of MuTAIL (thermal asymmetrically interlaced) and amplification of insertion mutagenized sites (AIMS) protocols and is applicable to the analysis of single mutants or to high-throughput analyses of mutant collections. Briefly, genomic DNA is digested with a restriction enzyme and adapters are ligated. Polymerase chain reaction is performed with TAIL cycling parameters, using a fluorescently labeled Mu primer, which results in the preferential amplification and labeling of Mu -containing genomic fragments. Products from a segregating line are analyzed on a capillary sequencer. To recover a fragment of interest, PCR products are cloned and sequenced. Sequences with lengths matching the size of a band that co-segregates with the mutant phenotype represent candidate linked insertion sites, which are then confirmed by PCR. We demonstrate the utility of the method by identifying Mu insertion sites linked to seed-lethal mutations with a preliminary success rate of nearly 50%.     

Integration Profiling of Gene Function With Dense Maps of Transposon Integration

Understanding how complex networks of genes integrate to produce dividing cells is an important goal that is limited by the difficulty in defining the function of individual genes. Current resources for the systematic identification of gene function such as siRNA libraries and collections of deletion strains are costly and organism specific. We describe here integration profiling, a novel approach to identify the function of eukaryotic genes based upon dense maps of transposon integration. As a proof of concept, we used the transposon Hermes to generate a library of 360,513 insertions in the genome of Schizosaccharomyces pombe. On average, we obtained one insertion for every 29 bp of the genome. Hermes integrated more often into nucleosome free sites and 33% of the insertions occurred in ORFs. We found that ORFs with low integration densities successfully identified the genes that are essential for cell division. Importantly, the nonessential ORFs with intermediate levels of insertion correlated with the nonessential genes that have functions required for colonies to reach full size. This finding indicates that integration profiles can measure the contribution of nonessential genes to cell division. While integration profiling succeeded in identifying genes necessary for propagation, it also has the potential to identify genes important for many other functions such as DNA repair, stress response, and meiosis.