Whole-genome analysis of Fusarium graminearum insertional mutants identifies virulence associated genes and unmasks untagged chromosomal deletions

Background

Identifying pathogen virulence genes required to cause disease is crucial to understand the mechanisms underlying the pathogenic process. Plasmid insertion mutagenesis of fungal protoplasts is frequently used for this purpose in filamentous ascomycetes. Post transformation, the mutant population is screened for loss of virulence to a specific plant or animal host. Identifying the insertion event has previously met with varying degrees of success, from a cleanly disrupted gene with minimal deletion of nucleotides at the insertion point to multiple-copy insertion events and large deletions of chromosomal regions. Currently, extensive mutant collections exist in laboratories globally where it was hitherto impossible to identify all the affected genes.

Results

We used a whole-genome sequencing (WGS) approach using Illumina HiSeq 2000 technology to investigate DNA tag insertion points and chromosomal deletion events in mutagenised, reduced virulence F. graminearum isolates identified in disease tests on wheat (Triticum aestivum). We developed the FindInsertSeq workflow to localise the DNA tag insertions to the nucleotide level. The workflow was tested using four mutants showing evidence of single and multi-copy insertions in DNA blot analysis. FindInsertSeq was able to identify both single and multi-copy concatenation insertion sites. By comparing sequencing coverage, unexpected molecular recombination events such as large tagged and untagged chromosomal deletions, and DNA amplification were observed in three of the analysed mutants. A random data sampling approach revealed the minimum genome coverage required to survey the F. graminearum genome for alterations.

Conclusions

This study demonstrates that whole-genome re-sequencing to 22x fold genome coverage is an efficient tool to characterise single and multi-copy insertion mutants in the filamentous ascomycete Fusarium graminearum. In some cases insertion events are accompanied with large untagged chromosomal deletions while in other cases a straight-forward insertion event could be confirmed. The FindInsertSeq analysis workflow presented in this study enables researchers to efficiently characterise insertion and deletion mutants.

Electronic supplementary material

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

Genetic and molecular characterization of embryonic mutants identified following seed transformation in Arabidopsis

Over 5000 transgenic families of Arabidopsis thaliana produced following seed transformation with Agrobacterium tumefaciens were screened for embryonic lethals, defectives, and pattern mutants. One hundred and seventy-eight mutants with a wide range of developmental abnormalities were identified. Forty-one mutants appear from genetic studies to be tagged (36% of the 115 mutants examined in detail). Mapping with visible markers demonstrated that mutant genes were randomly distributed throughout the genome. Seven mutant families appeared to contain chromosomal translocations because the mutant genes exhibited linkage to visible markers on two different chromosomes. Chromosomal rearrangements may therefore be widespread following seed transformation. DNA gel blot hybridizations with 34 tagged mutants and three T-DNA probes revealed a wide range of insertion patterns. Models of T-DNA structure at each mutant locus were constructed to facilitate gene isolation. The value of such models was demonstrated by using plasmid rescue to clone flanking plant DNA from four tagged mutants. Further analysis of genes isolated from these insertional mutants should help to elucidate the relationship between gene function and plant embryogenesis.     

Isolation of the structural genes for the Rieske Fe-S protein, cytochrome b and cytochrome c1 all components of the ubiquinol: cytochrome c2 oxidoreductase complex of Rhodopseudomonas capsulata

The structural genes for the Rieske Fe-S protein (petA), cytochrome b (petB) and cytochrome c1 (petC) subunits of the ubiquinol:cytochrome c2 oxidoreductase (bc1 complex) of Rhodopseudomonas capsulata have been cloned by complementation, using a mutant defective in this complex. The location of these genes on the obtained plasmid, pR14A, was determined using synthetic mixed oligonucleotide probes corresponding to highly conserved amino acid sequences of these proteins from various organisms. Their correct identity was established by partial sequencing. The petA, petB and petC genes were found to lie close to each other in this order, spanning two adjacent EcoRI fragments of 2.7 X 10(3) and 1.3 X 10(3) base-pairs, respectively. An insertion-deletion mutation, covering most of petB and all of petC and an insertion mutation, located in petB were constructed in vitro and were introduced into the chromosome of an otherwise wild-type strain by gene transfer agent-mediated genetic crosses. The bc-1 mutants obtained were defective in photosynthesis but, as expected, they could grow by respiration because of a branched respiratory pathway. Therefore, in R. capsulata a functional bc1 complex is essential in vivo for photosynthesis but not for respiration. Further, in the respiratory pathway the branch point must be before the bc1 complex, most likely at the quinone pool. These mutants were also proficient in anaerobic growth in the presence of dimethylsulfoxide, indicating that a functional bc1 complex is not required for this pathway. Several other insertions and deletions, located outside of the pet gene cluster, were also constructed. The ability of these latter mutants to grow photosynthetically suggested that no other gene essential for photosynthesis is located in the proximity of the pet cluster. The plasmid pR14A was shown to complement in trans the bc-1 insertion or insertion-deletion mutants, indicating that the pet genes were expressed in R. capsulata. Cross-hybridization experiments showed that the pet cluster was quite distinct from other known genes involved in photosynthesis.     

Identification of Yersinia enterocolitica genes affecting survival in an animal host using signature-tagged transposon mutagenesis

Pathogenic Yersinia species are associated with both localized and systemic infections in mammalian hosts. In this study, signature-tagged transposon mutagenesis was used to identify Yersinia enterocolitica genes required for survival in a mouse model of infection. Approximately 2000 transposon insertion mutants were screened for attenuation. This led to the identification of 55 mutants defective for survival in the animal host, as judged by their ability to compete with the wild-type strain in mixed infections. A total of 28 mutants had transposon insertions in the virulence plasmid, validating the screen. Two of the plasmid mutants with severe virulence defects had insertions in an uncharacterized region. Several of the chromosomal insertions were in a gene cluster involved in O-antigen biosynthesis. Other chromosomal insertions identified genes not previously demonstrated as being required for in vivo survival of Y. enterocolitica. These include genes involved in the synthesis of outer membrane components, stress response and nutrient acquisition. One severely attenuated mutant had an insertion in a homologue of the pspC gene (phage shock protein C) of Escherichia coli. The phage shock protein operon has no known biochemical or physiological function in E. coli, but is apparently essential for the survival of Y. enterocolitica during infection.     

Random Mutagenesis Identifies Novel Genes Involved in the Secretion of Antimicrobial, Cell Wall-Lytic Enzymes by Lactococcus lactis

Lactococcus lactis is a gram-positive bacterium that is widely used in the food industry and is therefore desirable as a candidate for the production and secretion of recombinant proteins. Previously, we generated a L. lactis strain that expressed and secreted the antimicrobial cell wall-lytic enzyme lysostaphin. To identify lactococcal gene products that affect the production of lysostaphin, we isolated and characterized mutants generated by random transposon mutagenesis that had altered lysostaphin activity. Out of 35,000 mutants screened, only one with no lysostaphin activity was identified, and it was found to contain an insertion in the lysostaphin expression cassette. Ten mutants with higher lysostaphin activity contained insertions in only four different genes, which encode an uncharacterized putative transmembrane protein (llmg_0609) (three mutants), an enzyme catalyzing the first step in peptidoglycan biosynthesis (murA2) (five mutants), a putative regulator of peptidoglycan modification (trmA) (one mutant), and an uncharacterized enzyme possibly involved in ubiquinone biosynthesis (llmg_2148) (one mutant). These mutants were found to secrete larger amounts of lysostaphin than the control strain (MG1363[lss]), and the greatest increase in secretion was 9.8- to 16.1-fold, for the llmg_0609 mutants. The lysostaphin-oversecreting llmg_0609, murA2, and trmA mutants were also found to secrete larger amounts of another cell wall-lytic enzyme (the Listeria monocytogenes bacteriophage endolysin Ply511) than the control strain, indicating that the phenotype is not limited to lysostaphin.     

An Assessment of Heavy Ion Irradiation Mutagenesis for Reverse Genetics in Wheat (Triticum aestivum L.)

Reverse genetic techniques harnessing mutational approaches are powerful tools that can provide substantial insight into gene function in plants. However, as compared to diploid species, reverse genetic analyses in polyploid plants such as bread wheat can present substantial challenges associated with high levels of sequence and functional similarity amongst homoeologous loci. We previously developed a high-throughput method to identify deletions of genes within a physically mutagenized wheat population. Here we describe our efforts to combine multiple homoeologous deletions of three candidate disease susceptibility genes (TaWRKY11, TaPFT1 and TaPLDß1). We were able to produce lines featuring homozygous deletions at two of the three homoeoloci for all genes, but this was dependent on the individual mutants used in crossing. Intriguingly, despite extensive efforts, viable lines possessing homozygous deletions at all three homoeoloci could not be produced for any of the candidate genes. To investigate deletion size as a possible reason for this phenomenon, we developed an amplicon sequencing approach based on synteny to Brachypodium distachyon to assess the size of the deletions removing one candidate gene (TaPFT1) in our mutants. These analyses revealed that genomic deletions removing the locus are relatively large, resulting in the loss of multiple additional genes. The implications of this work for the use of heavy ion mutagenesis for reverse genetic analyses in wheat are discussed.     

Coupling the CRISPR/Cas9 System with Lambda Red Recombineering Enables Simplified Chromosomal Gene Replacement in Escherichia coli

To date, most genetic engineering approaches coupling the type II Streptococcus pyogenes clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system to lambda Red recombineering have involved minor single nucleotide mutations. Here we show that procedures for carrying out more complex chromosomal gene replacements in Escherichia coli can be substantially enhanced through implementation of CRISPR/Cas9 genome editing. We developed a three-plasmid approach that allows not only highly efficient recombination of short single-stranded oligonucleotides but also replacement of multigene chromosomal stretches of DNA with large PCR products. By systematically challenging the proposed system with respect to the magnitude of chromosomal deletion and size of DNA insertion, we demonstrated DNA deletions of up to 19.4 kb, encompassing 19 nonessential chromosomal genes, and insertion of up to 3 kb of heterologous DNA with recombination efficiencies permitting mutant detection by colony PCR screening. Since CRISPR/Cas9-coupled recombineering does not rely on the use of chromosome-encoded antibiotic resistance, or flippase recombination for antibiotic marker recycling, our approach is simpler, less labor-intensive, and allows efficient production of gene replacement mutants that are both markerless and “scar”-less.     

Roles of the Structural Symbiosis Polysaccharide (syp) Genes in Host Colonization,Biofilm Formation,and Polysaccharide Biosynthesis in Vibrio fischeri

The symbiosis polysaccharide locus, syp, is required for Vibrio fischeri to form a symbiotic association with the squid Euprymna scolopes. It is also required for biofilm formation induced by the unlinked regulator RscS. The syp locus includes 18 genes that can be classified into four groups based on putative function: 4 genes encode putative regulators, 6 encode glycosyltransferases, 2 encode export proteins, and the remaining 6 encode proteins with other functions, including polysaccharide modification. To understand the roles of each of the 14 structural syp genes in colonization and biofilm formation, we generated nonpolar in-frame deletions of each gene. All of the deletion mutants exhibited defects in their ability to colonize juvenile squid, although the impact of the loss of SypB or SypI was modest. Consistent with their requirement for colonization, most of the structural genes were also required for RscS-induced biofilm formation. In particular, the production of wrinkled colonies, pellicles, and the matrix on the colony surface was eliminated or severely decreased in all mutants except for the sypB and sypI mutants; in contrast, only a subset of genes appeared to play a role in attachment to glass. Finally, immunoblotting data suggested that the structural Syp proteins are involved in polysaccharide production and/or export. These results provide important insights into the requirements for the syp genes under different environmental conditions and thus lay the groundwork for a more complete understanding of the matrix produced by V. fischeri to enhance cell-cell interactions and promote symbiotic colonization.     

Nucleotide sequences of Bacillus subtilis flagellar biosynthetic genes fliP and fliQ and identification of a novel flagellar gene, fliZ.

Three genes from the Bacillus subtilis major che-fla operon have been cloned and sequenced. Two of the genes encode proteins that are homologous to the Escherichia coli and Salmonella typhimurium flagellar biosynthetic proteins FliP and FliQ. The third gene, designated fliZ, encodes a 219-amino-acid protein with a predicted molecular mass of 24,872 Da. FliZ is not significantly homologous to any known proteins. Null mutants in fliP and fliZ do not have flagella; however, motility can be restored to the fliZ null mutant by expression of fliZ from a plasmid. FliZ has a conventional N-terminal signal sequence that does not direct secretion of the protein but appears to target the protein to the membrane. Two possible models of insertion of FliZ into the membrane are described.     

Genome-wide mutagenesis of <Emphasis Type="Italic">Zea mays</Emphasis> L. using <Emphasis Type="Italic">RescueMu</Emphasis> transposons

Derived from the maize Mu1 transposon, RescueMu provides strategies for maize gene discovery and mutant phenotypic analysis. 9.92 Mb of gene-enriched sequences next to RescueMu insertion sites were co-assembled with expressed sequence tags and analyzed. Multiple plasmid recoveries identified probable germinal insertions and screening of RescueMu plasmid libraries identified plants containing probable germinal insertions. Although frequently recovered parental insertions and insertion hotspots reduce the efficiency of gene discovery per plasmid, RescueMu targets a large variety of genes and produces knockout mutants.     

Genetic and molecular characterization of embryonic mutants identified following seed transformation in Arabidopsis

Over 5000 transgenic families of Arabidopsis thaliana produced following seed transformation with Agrobacterium tumefaciens were screened for embryonic lethals, defectives, and pattern mutants. One hundred and seventy-eight mutants with a wide range of developmental abnormalities were identified. Forty-one mutants appear from genetic studies to be tagged (36% of the 115 mutants examined in detail). Mapping with visible markers demonstrated that mutant genes were randomly distributed throughout the genome. Seven mutant families appeared to contain chromosomal translocations because the mutant genes exhibited linkage to visible markers on two different chromosomes. Chromosomal rearrangements may therefore be widespread following seed transformation. DNA gel blot hybridizations with 34 tagged mutants and three T-DNA probes revealed a wide range of insertion patterns. Models of T-DNA structure at each mutant locus were constructed to facilitate gene isolation. The value of such models was demonstrated by using plasmid rescue to clone flanking plant DNA from four tagged mutants. Further analysis of genes isolated from these insertional mutants should help to elucidate the relationship between gene function and plant embryogenesis.     

Directed mutagenesis of chromosomal genes of Escherichia coli in vivo: construction of RecA- and HtpR- mutants

The deletions in Escherichia coli chromosomal genes recA and htpR were constructed using the site-directed mutagenesis techniques. The obtained RecA- mutants are UV-sensitive and have a phenotype defective for the homologous DNA recombination. HtpR- mutant is temperature sensitive for growth and deficient in intracellular proteolysis. As a result a HtpR- mutant seems to be a preferable candidate for attempting to synthesize efficiently any alien protein in Escherichia coli cells.     

Genome-Wide Screening with Hydroxyurea Reveals a Link between Nonessential Ribosomal Proteins and Reactive Oxygen Species Production

We performed a screening of hydroxyurea (HU)-sensitive mutants using a single-gene-deletion mutant collection in Escherichia coli. HU inhibits ribonucleotide reductase (RNR), which leads to arrest of the replication fork. Surprisingly, the wild-type was less resistant to HU than the average for the Keio Collection. Respiration-defective mutants were significantly more resistant to HU, suggesting that the generation of reactive oxygen species (ROS) contributes to cell death. High-throughput screening revealed that 15 mutants were completely sensitive on plates containing 7.5 mM HU. Unexpectedly, translation-related mutants based on COG categorization were the most enriched, and three of them were deletion mutants of nonessential ribosomal proteins (L1, L32, and L36). We found that, in these mutants, an increased membrane stress response was provoked, resulting in increased ROS generation. The addition of OH radical scavenger thiourea rescued the HU sensitivity of these mutants, suggesting that ROS generation is the direct cause of cell death. Conversely, both the deletion of rpsF and the deletion of rimK, which encode S6 and S6 modification enzymes, respectively, showed an HU-resistant phenotype. These mutants increased the copy number of the p15A-based plasmid and exhibited reduced basal levels of SOS response. The data suggest that nonessential proteins indirectly affect the DNA-damaging process.     

Establishment of Tn5096-Based Transposon Mutagenesis in Gordonia polyisoprenivorans

The transposons Tn5, Tn10, Tn611, and Tn5096 were characterized regarding transposition in Gordonia polyisoprenivorans strain VH2. No insertional mutants were obtained employing Tn5 or Tn10. The thermosensitive plasmid pCG79 harboring Tn611 integrated into the chromosome of G. polyisoprenivorans; however, the insertional mutants were fairly unstable und reverted frequently to the wild-type phenotype. In contrast, various stable mutants were obtained employing Tn5096-mediated transposon mutagenesis. Auxotrophic mutants, mutants defective or deregulated in carotenoid biosynthesis, and mutants defective in utilization of rubber and/or highly branched isoprenoid hydrocarbons were obtained by integration of plasmid pMA5096 harboring Tn5096 as a whole into the genome. From about 25,000 isolated mutants, the insertion loci of pMA5096 were subsequently mapped in 20 independent mutants in genes which could be related to the above-mentioned metabolic pathways or to putative regulation proteins. Analyses of the genotypes of pMA5096-mediated mutants defective in biodegradation of poly(cis-1,4-isoprene) did not reveal homologues to recently identified genes coding for enzymes catalyzing the initial cleavage of poly(cis-1,4-isoprene). One rubber-negative mutant was disrupted in mcr, encoding an α-methylacyl-coenzyme A racemase. This mutant was defective in degradation of poly(cis-1,4-isoprene) and also of highly branched isoprenoid hydrocarbons.     

Genome structure of Halobacterium halobium: plasmid dynamics in gas vacuole deficient mutants

Halobacterium halobium contains two gas vacuole protein genes that are located in plasmid pHH1 (p-vac) and in the chromosomal DNA (c-vac). The mutation frequency for these genes is different: the constitutively expressed p-vac gene is mutated with a frequency of 10(-2), while the chromosomal gene expressed in the stationary phase of growth is mutated with a frequency of 10(-5). The difference in the mutation susceptibility is due to the dynamics of plasmid pHH1. p-vac gene mutations are caused (i) by the integration of an insertion element or (ii) by a deletion event encompassing the p-vac gene region. In contrast, c-vac mutants analyzed to date incurred neither insertion elements nor deletions. Deletion events within pHH1 occur at high frequencies during the development of a H. halobium culture. The investigation of the fusion regions resulting from deletion events indicates that insertion elements are involved. The analysis of pHH1 deletion variants led to a 4 kilobase pair DNA region containing the origin of replication of the pHH1 plasmid.     

Comprehensive genetic analysis of adhesin proteins and their role in virulence of Candida albicans

Candida albicans is a microbial fungus that exists as a commensal member of the human microbiome and an opportunistic pathogen. Cell surface-associated adhesin proteins play a crucial role in C. albicans’ ability to undergo cellular morphogenesis, develop robust biofilms, colonize, and cause infection in a host. However, a comprehensive analysis of the role and relationships between these adhesins has not been explored. We previously established a CRISPR-based platform for efficient generation of single- and double-gene deletions in C. albicans, which was used to construct a library of 144 mutants, comprising 12 unique adhesin genes deleted singly, and every possible combination of double deletions. Here, we exploit this adhesin mutant library to explore the role of adhesin proteins in C. albicans virulence. We perform a comprehensive, high-throughput screen of this library, using Caenorhabditis elegans as a simplified model host system, which identified mutants critical for virulence and significant genetic interactions. We perform follow-up analysis to assess the ability of high- and low-virulence strains to undergo cellular morphogenesis and form biofilms in vitro, as well as to colonize the C. elegans host. We further perform genetic interaction analysis to identify novel significant negative genetic interactions between adhesin mutants, whereby combinatorial perturbation of these genes significantly impairs virulence, more than expected based on virulence of the single mutant constituent strains. Together, this study yields important new insight into the role of adhesins, singly and in combinations, in mediating diverse facets of virulence of this critical fungal pathogen.     

Characterization of DNA Repair Deficient Strains of Chlamydomonas reinhardtii Generated by Insertional Mutagenesis

While the mechanisms governing DNA damage response and repair are fundamentally conserved, cross-kingdom comparisons indicate that they differ in many aspects due to differences in life-styles and developmental strategies. In photosynthetic organisms these differences have not been fully explored because gene-discovery approaches are mainly based on homology searches with known DDR/DNA repair proteins. Here we performed a forward genetic screen in the green algae Chlamydomonas reinhardtii to identify genes deficient in DDR/DNA repair. We isolated five insertional mutants that were sensitive to various genotoxic insults and two of them exhibited altered efficiency of transgene integration. To identify genomic regions disrupted in these mutants, we established a novel adaptor-ligation strategy for the efficient recovery of the insertion flanking sites. Four mutants harbored deletions that involved known DNA repair factors, DNA Pol zeta, DNA Pol theta, SAE2/COM1, and two neighbouring genes encoding ERCC1 and RAD17. Deletion in the last mutant spanned two Chlamydomonas-specific genes with unknown function, demonstrating the utility of this approach for discovering novel factors involved in genome maintenance.     

Mutational analysis of the Streptomyces lividans recA gene suggests that only mutants with residual activity remain viable

Temperature-sensitive integration plasmids carrying internal fragments of the Streptomyces lividans TK24 recA gene were constructed and used to inactivate the chromosomal recA gene of S. lividans by gene disruption and gene replacement. Integration of these plasmids resulted in recA mutants expressing C-terminally truncated RecA proteins, as deduced from Southern hybridization experiments. Mutants FRECD2 in which the last 42 amino acids, comprising the variable part of bacterial RecA proteins, had been deleted retained the wild-type phenotype. The S. lividans recA mutant FRECD3 produced a RecA protein lacking 87 amino acids probably including the interfilament contact site. FRECD3 was more sensitive to UV and MMS than the wild-type. Its ability to undergo homologous recombination was impaired, but not completely abolished. Integration of the disruption plasmid pFRECD3 in S. coelicolor“Müller” caused the same mutant phenotype as S. lividans FRECD3. In spite of many attempts no S. lividans recA mutants with deletions of 165 C-terminal amino acids or more were isolated. Furthermore, the recA gene could not be replaced by a kanamycin resistance cassette. These experiments indicate a crucial role of the recA gene in ensuring viability of Streptomyces.