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.     

Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-602) contains supplementary material, which is available to authorized users.     

Mono-allelic retrotransposon insertion addresses epigenetic transcriptional repression in human genome

Background

Retrotransposons have been extensively studied in plants and animals and have been shown to have an impact on human genome dynamics and evolution. Their ability to move within genomes gives retrotransposons to affect genome instability.

Methods

we examined the polymorphic inserted AluYa5, evolutionary young Alu, in the progesterone receptor gene to determine the effects of Alu insertion on molecular environment. We used mono-allelic inserted cell lines which carry both Alu-present and Alu-absent alleles. To determine the epigenetic change and gene expression, we performed restriction enzyme digestion, Pyrosequencing, and Chromatin Immunoprecipitation.

Results

We observed that the polymorphic insertion of evolutionally young Alu causes increasing levels of DNA methylation in the surrounding genomic area and generates inactive histone tail modifications. Consequently the Alu insertion deleteriously inactivates the neighboring gene expression.

Conclusion

The mono-allelic Alu insertion cell line clearly showed that polymorphic inserted repetitive elements cause the inactivation of neighboring gene expression, bringing aberrant epigenetic changes.     

Structural and sequence diversity of the transposon Galileo in the Drosophila willistoni genome

Background

Galileo is one of three members of the P superfamily of DNA transposons. It was originally discovered in Drosophila buzzatii, in which three segregating chromosomal inversions were shown to have been generated by ectopic recombination between Galileo copies. Subsequently, Galileo was identified in six of 12 sequenced Drosophila genomes, indicating its widespread distribution within this genus. Galileo is strikingly abundant in Drosophila willistoni, a neotropical species that is highly polymorphic for chromosomal inversions, suggesting a role for this transposon in the evolution of its genome.

Results

We carried out a detailed characterization of all Galileo copies present in the D. willistoni genome. A total of 191 copies, including 133 with two terminal inverted repeats (TIRs), were classified according to structure in six groups. The TIRs exhibited remarkable variation in their length and structure compared to the most complete copy. Three copies showed extended TIRs due to internal tandem repeats, the insertion of other transposable elements (TEs), or the incorporation of non-TIR sequences into the TIRs. Phylogenetic analyses of the transposase (TPase)-encoding and TIR segments yielded two divergent clades, which we termed Galileo subfamilies V and W. Target-site duplications (TSDs) in D. willistoni Galileo copies were 7- or 8-bp in length, with the consensus sequence GTATTAC. Analysis of the region around the TSDs revealed a target site motif (TSM) with a 15-bp palindrome that may give rise to a stem-loop secondary structure.

Conclusions

There is a remarkable abundance and diversity of Galileo copies in the D. willistoni genome, although no functional copies were found. The TIRs in particular have a dynamic structure and extend in different ways, but their ends (required for transposition) are more conserved than the rest of the element. The D. willistoni genome harbors two Galileo subfamilies (V and W) that diverged ~9 million years ago and may have descended from an ancestral element in the genome. Galileo shows a significant insertion preference for a 15-bp palindromic TSM.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-792) contains supplementary material, which is available to authorized users.     

A new platform for ultra-high density Staphylococcus aureus transposon libraries

Background

Staphylococcus aureus readily develops resistance to antibiotics and achieving effective therapies to overcome resistance requires in-depth understanding of S. aureus biology. High throughput, parallel-sequencing methods for analyzing transposon mutant libraries have the potential to revolutionize studies of S. aureus, but the genetic tools to take advantage of the power of next generation sequencing have not been fully developed.

Results

Here we report a phage-based transposition system to make ultra-high density transposon libraries for genome-wide analysis of mutant fitness in any Φ11-transducible S. aureus strain. The high efficiency of the delivery system has made it possible to multiplex transposon cassettes containing different regulatory elements in order to make libraries in which genes are over- or under-expressed as well as deleted. By incorporating transposon-specific barcodes into the cassettes, we can evaluate how null mutations and changes in gene expression levels affect fitness in a single sequencing data set. Demonstrating the power of the system, we have prepared a library containing more than 690,000 unique insertions. Because one unique feature of the phage-based approach is that temperature-sensitive mutants are retained, we have carried out a genome-wide study of S. aureus genes involved in withstanding temperature stress. We find that many genes previously identified as essential are temperature sensitive and also identify a number of genes that, when disrupted, confer a growth advantage at elevated temperatures.

Conclusions

The platform described here reliably provides mutant collections of unparalleled genotypic diversity and will enable a wide range of functional genomic studies in S. aureus.

Electronic supplementary material

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

Disruption of Mycobacterium avium subsp. paratuberculosis-specific genes impairs in vivo fitness

Background

Mycobacterium avium subsp. paratuberculosis (MAP) is an obligate intracellular pathogen that infects many ruminant species. The acquisition of foreign genes via horizontal gene transfer has been postulated to contribute to its pathogenesis, as these genetic elements are absent from its putative ancestor, M. avium subsp. hominissuis (MAH), an environmental organism with lesser pathogenicity. In this study, high-throughput sequencing of MAP transposon libraries were analyzed to qualitatively and quantitatively determine the contribution of individual genes to bacterial survival during infection.

Results

Out of 52384 TA dinucleotides present in the MAP K-10 genome, 12607 had a MycoMarT7 transposon in the input pool, interrupting 2443 of the 4350 genes in the MAP genome (56%). Of 96 genes situated in MAP-specific genomic islands, 82 were disrupted in the input pool, indicating that MAP-specific genomic regions are dispensable for in vitro growth (odds ratio = 0.21). Following 5 independent in vivo infections with this pool of mutants, the correlation between output pools was high for 4 of 5 (R = 0.49 to 0.61) enabling us to define genes whose disruption reproducibly reduced bacterial fitness in vivo. At three different thresholds for reduced fitness in vivo, MAP-specific genes were over-represented in the list of predicted essential genes. We also identified additional genes that were severely depleted after infection, and several of them have orthologues that are essential genes in M. tuberculosis.

Conclusions

This work indicates that the genetic elements required for the in vivo survival of MAP represent a combination of conserved mycobacterial virulence genes and MAP-specific genes acquired via horizontal gene transfer. In addition, the in vitro and in vivo essential genes identified in this study may be further characterized to offer a better understanding of MAP pathogenesis, and potentially contribute to the discovery of novel therapeutic and vaccine targets.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-415) contains supplementary material, which is available to authorized users.     

Simulating the impact of genetic diversity of Medicago truncatula on germination and emergence using a crop emergence model for ideotype breeding

Background and Aims

Germination and heterotrophic growth are crucial steps for stand establishment. Numerical experiments based on the modelling of these early stages in relation to major environmental factors at sowing were used as a powerful tool to browse the effects of the genetic diversity of Medicago truncatula, one of the model legume species, under a range of agronomic scenarios, and to highlight the most important plant parameters for emergence. To this end, the emergence of several genotypes of M. truncatula was simulated under a range of sowing conditions with a germination and emergence simulation model.

Methods

After testing the predictive quality of the model by comparing simulations to field observations of several genotypes of M. truncatula, numerical experiments were performed under a wide range of environmental conditions (sowing dates × years × seedbed structure). Germination and emergence was simulated for a set of five genotypes previously parameterized and for two virtual genotypes engineered to maximize the potential effects of genetic diversity.

Key Results

The simulation results gave an average value of 5–10 % difference in final emergence between genotypes, which was low, but the analysis underlined considerable inter-annual variation. The effects of parameters describing germination and emergence processes were quantified and ranked according to their contribution to the variation in emergence. Seedling non-emergence was mainly related to mechanical obstacles (40–50 %). More generally, plant parameters that accelerated the emergence time course significantly contributed to limiting the risk of soil surface crusting occurring before seedling emergence.

Conclusions

The model-assisted analysis of the effects of genetic diversity demonstrated its usefulness in helping to identify the parameters which have most influence that could be improved by breeding programmes. These results should also enable a deeper analysis of the genetic determinism of the main plant parameters influencing emergence, using the genomic tools available for this model plant.     

Unusual genome complexity in Lactobacillus salivarius JCM1046

Background

Lactobacillus salivarius strains are increasingly being exploited for their probiotic properties in humans and animals. Dissemination of antibiotic resistance genes among species with food or probiotic-association is undesirable and is often mediated by plasmids or integrative and conjugative elements. L. salivarius strains typically have multireplicon genomes including circular megaplasmids that encode strain-specific traits for intestinal survival and probiotic activity. Linear plasmids are less common in lactobacilli and show a very limited distribution in L. salivarius. Here we present experimental evidence that supports an unusually complex multireplicon genome structure in the porcine isolate L. salivarius JCM1046.

Results

JCM1046 harbours a 1.83 Mb chromosome, and four plasmids which constitute 20% of the genome. In addition to the known 219 kb repA-type megaplasmid pMP1046A, we identified and experimentally validated the topology of three additional replicons, the circular pMP1046B (129 kb), a linear plasmid pLMP1046 (101 kb) and pCTN1046 (33 kb) harbouring a conjugative transposon. pMP1046B harbours both plasmid-associated replication genes and paralogues of chromosomally encoded housekeeping and information-processing related genes, thus qualifying it as a putative chromid. pLMP1046 shares limited sequence homology or gene synteny with other L. salivarius plasmids, and its putative replication-associated protein is homologous to the RepA/E proteins found in the large circular megaplasmids of L. salivarius. Plasmid pCTN1046 harbours a single copy of an integrated conjugative transposon (Tn6224) which appears to be functionally intact and includes the tetracycline resistance gene tetM.

Conclusion

Experimental validation of sequence assemblies and plasmid topology resolved the complex genome architecture of L. salivarius JCM1046. A high-coverage draft genome sequence would not have elucidated the genome complexity in this strain. Given the expanding use of L. salivarius as a probiotic, it is important to determine the genotypic and phenotypic organization of L. salivarius strains. The identification of Tn6224-like elements in this species has implications for strain selection for probiotic applications.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-771) contains supplementary material, which is available to authorized users.     

Sequencing methods and datasets to improve functional interpretation of sleeping beauty mutagenesis screens

Background

Animal models of cancer are useful to generate complementary datasets for comparison to human tumor data. Insertional mutagenesis screens, such as those utilizing the Sleeping Beauty (SB) transposon system, provide a model that recapitulates the spontaneous development and progression of human disease. This approach has been widely used to model a variety of cancers in mice. Comprehensive mutation profiles are generated for individual tumors through amplification of transposon insertion sites followed by high-throughput sequencing. Subsequent statistical analyses identify common insertion sites (CISs), which are predicted to be functionally involved in tumorigenesis. Current methods utilized for SB insertion site analysis have some significant limitations. For one, they do not account for transposon footprints – a class of mutation generated following transposon remobilization. Existing methods also discard quantitative sequence data due to uncertainty regarding the extent to which it accurately reflects mutation abundance within a heterogeneous tumor. Additionally, computational analyses generally assume that all potential insertion sites have an equal probability of being detected under non-selective conditions, an assumption without sufficient relevant data. The goal of our study was to address these potential confounding factors in order to enhance functional interpretation of insertion site data from tumors.

Results

We describe here a novel method to detect footprints generated by transposon remobilization, which revealed minimal evidence of positive selection in tumors. We also present extensive characterization data demonstrating an ability to reproducibly assign semi-quantitative information to individual insertion sites within a tumor sample. Finally, we identify apparent biases for detection of inserted transposons in several genomic regions that may lead to the identification of false positive CISs.

Conclusion

The information we provide can be used to refine analyses of data from insertional mutagenesis screens, improving functional interpretation of results and facilitating the identification of genes important in cancer development and progression.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1150) contains supplementary material, which is available to authorized users.     

Genomic factors related to tissue tropism in Chlamydia pneumoniae infection

Background

Chlamydia pneumoniae (Cpn) are obligate intracellular bacteria that cause acute infections of the upper and lower respiratory tract and have been implicated in chronic inflammatory diseases. Although of significant clinical relevance, complete genome sequences of only four clinical Cpn strains have been obtained. All of them were isolated from the respiratory tract and shared more than 99% sequence identity. Here we investigate genetic differences on the whole-genome level that are related to Cpn tissue tropism and pathogenicity.

Results

We have sequenced the genomes of 18 clinical isolates from different anatomical sites (e.g. lung, blood, coronary arteries) of diseased patients, and one animal isolate. In total 1,363 SNP loci and 184 InDels have been identified in the genomes of all clinical Cpn isolates. These are distributed throughout the whole chlamydial genome and enriched in highly variable regions. The genomes show clear evidence of recombination in at least one potential region but no phage insertions. The tyrP gene was always encoded as single copy in all vascular isolates. Phylogenetic reconstruction revealed distinct evolutionary lineages containing primarily non-respiratory Cpn isolates. In one of these, clinical isolates from coronary arteries and blood monocytes were closely grouped together. They could be distinguished from all other isolates by characteristic nsSNPs in genes involved in RB to EB transition, inclusion membrane formation, bacterial stress response and metabolism.

Conclusions

This study substantially expands the genomic data of Cpn and elucidates its evolutionary history. The translation of the observed Cpn genetic differences into biological functions and the prediction of novel pathogen-oriented diagnostic strategies have to be further explored.

Electronic supplementary material

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

Tn5/7-lux: a versatile tool for the identification and capture of promoters in Gram-negative bacteria

Background

The combination of imaging technologies and luciferase-based bioluminescent bacterial reporter strains provide a sensitive and simple non-invasive detection method (photonic bioimaging) for the study of diverse biological processes, as well as efficacy of therapeutic interventions, in live animal models of disease. The engineering of bioluminescent bacteria required for photonic bioimaging is frequently hampered by lack of promoters suitable for strong, yet stable luciferase gene expression.

Results

We devised a novel method for identification of constitutive native promoters in Gram-negative bacteria. The method is based on a Tn5/7 transposon that exploits the unique features of Tn5 (random transposition) and Tn7 (site-specific transposition). The transposons are designed such that Tn5 transposition will allow insertion of a promoter-less bacterial luxCDABE operon downstream of a bacterial gene promoter. Cloning of DNA fragments from luminescent isolates results in a plasmid that replicates in pir⁺ hosts. Sequencing of the lux-chromosomal DNA junctions on the plasmid reveals transposon insertion sites within genes or operons. The plasmid is also a mini-Tn7-lux delivery vector that can be used to introduce the promoter-lux operon fusion into other derivatives of the bacterium of interest in an isogenic fashion. Alternatively, promoter-containing sequences can be PCR-amplified from plasmid or chromosomal DNA and cloned into a series of accompanying mini-Tn7-lux vectors. The mini-Tn5/7-lux and mini-Tn7-lux vectors are equipped with diverse selection markers and thus applicable in numerous Gram-negative bacteria. Various mini-Tn5/7-lux vectors were successfully tested for transposition and promoter identification by imaging in Acinetobacter baumannii, Escherichia coli, and Burkholderia pseudomallei. Strong promoters were captured for lux expression in E. coli and A. baumannii. Some mini-Tn7-lux vectors are also equipped with attB sites for swapping of the lux operon with other reporter genes using Gateway technology.

Conclusions

Although mini-Tn5-lux and mini-Tn7-lux elements have previously been developed and used for bacterial promoter identification and chromosomal insertion of promoter-lux gene fusions, respectively, the newly developed mini-Tn5/7-lux and accompanying accessory plasmids streamline and accelerate the promoter discovery and bioluminescent strain engineering processes. Availability of vectors with diverse selection markers greatly extend the host-range of promoter probe and lux gene fusion vectors.

Electronic supplementary material

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

Divergences in gene repertoire among the reference Prevotella genomes derived from distinct body sites of human

Background

The community composition of the human microbiome is known to vary at distinct anatomical niches. But little is known about the nature of variations, if any, at the genome/sub-genome levels of a specific microbial community across different niches. The present report aims to explore, as a case study, the variations in gene repertoire of 28 Prevotella reference genomes derived from different body-sites of human, as reported earlier by the Human Microbiome Consortium.

Results

The pan-genome for Prevotella remains “open”. On an average, 17% of predicted protein-coding genes of any particular Prevotella genome represent the conserved core genes, while the remaining 83% contribute to the flexible and singletons. The study reveals exclusive presence of 11798, 3673, 3348 and 934 gene families and exclusive absence of 17, 221, 115 and 645 gene families in Prevotella genomes derived from human oral cavity, gastro-intestinal tracts (GIT), urogenital tract (UGT) and skin, respectively. Distribution of various functional COG categories differs significantly among the habitat-specific genes. No niche-specific variations could be observed in distribution of KEGG pathways.

Conclusions

Prevotella genomes derived from different body sites differ appreciably in gene repertoire, suggesting that these microbiome components might have developed distinct genetic strategies for niche adaptation within the host. Each individual microbe might also have a component of its own genetic machinery for host adaptation, as appeared from the huge number of singletons.

Electronic supplementary material

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

Ontogenetic contingency of tolerance mechanisms in response to apical damage

Background and Aims

Plants are able to tolerate tissue loss through vigorous branching which is often triggered by release from apical dominance and activation of lateral meristems. However, damage-induced branching might not be a mere physiological outcome of released apical dominance, but an adaptive response to environmental signals, such as damage timing and intensity. Here, branching responses to both factors were examined in the annual plant Medicago truncatula.

Methods

Branching patterns and allocation to reproductive traits were examined in response to variable clipping intensities and timings in M. truncatula plants from two populations that vary in the onset of reproduction. Phenotypic selection analysis was used to evaluate the strength and direction of selection on branching under the damage treatments.

Key Results

Plants of both populations exhibited an ontogenetic shift in tolerance mechanisms: while early damage induced greater meristem activation, late damage elicited investment in late-determined traits, including mean pod and seed biomass, and supported greater germination rates. Severe damage mostly elicited simultaneous development of multiple-order lateral branches, but this response was limited to early damage. Selection analyses revealed positive directional selection on branching in plants under early- compared with late- or no-damage treatments.

Conclusions

The results demonstrate that damage-induced meristem activation is an adaptive response that could be modified according to the plant''s developmental stage, severity of tissue loss and their interaction, stressing the importance of considering these effects when studying plastic responses to apical damage.     

Cytogeography and genome size variation in the Claytonia perfoliata (Portulacaceae) polyploid complex

Background and Aims

Genome duplication is a central process in plant evolution and contributes to patterns of variation in genome size within and among lineages. Studies that combine cytogeography with genome size measurements contribute to our basic knowledge of cytotype distributions and their associations with variation in genome size.

Methods

Ploidy and genome size were assessed with direct chromosome counts and flow cytometry for 78 populations within the Claytonia perfoliata complex, comprised of three diploid taxa with numerous polyploids that range to the decaploid level. The relationship between genome size and temperature and precipitation was investigated within and across cytotypes to test for associations between environmental factors and nuclear DNA content.

Key Results

A euploid series (n = 6) of diploids to octoploids was documented through chromosome counts, and decaploids were suggested by flow cytometry. Increased variation in genome size among populations was found at higher ploidy levels, potentially associated with differential contributions of diploid parental genomes, variation in rates of genomic loss or gain, or undetected hybridization. Several accessions were detected with atypical genome sizes, including a diploid population of C. parviflora ssp. grandiflora with an 18 % smaller genome than typical, and hexaploids of C. perfoliata and C. parviflora with genomes 30 % larger than typical. There was a slight but significant association of larger genome sizes with colder winter temperature across the C. perfoliata complex as a whole, and a strong association between lower winter temperatures and large genome size for tetraploid C. parviflora.

Conclusions

The C. perfoliata complex is characterized by polyploids ranging from tetraploid to decaploid, with large magnitude variation in genome size at higher ploidy levels, associated in part with environmental variation in temperature.