Introgression browser: high‐throughput whole‐genome SNP visualization

Breeding by introgressive hybridization is a pivotal strategy to broaden the genetic basis of crops. Usually, the desired traits are monitored in consecutive crossing generations by marker‐assisted selection, but their analyses fail in chromosome regions where crossover recombinants are rare or not viable. Here, we present the Introgression Browser (iBrowser ), a bioinformatics tool aimed at visualizing introgressions at nucleotide or SNP (Single Nucleotide Polymorphisms) accuracy. The software selects homozygous SNPs from Variant Call Format (VCF) information and filters out heterozygous SNPs, multi‐nucleotide polymorphisms (MNPs) and insertion–deletions (InDels). For data analysis iBrowser makes use of sliding windows, but if needed it can generate any desired fragmentation pattern through General Feature Format (GFF) information. In an example of tomato (Solanum lycopersicum) accessions we visualize SNP patterns and elucidate both position and boundaries of the introgressions. We also show that our tool is capable of identifying alien DNA in a panel of the closely related S. pimpinellifolium by examining phylogenetic relationships of the introgressed segments in tomato. In a third example, we demonstrate the power of the iBrowser in a panel of 597 Arabidopsis accessions, detecting the boundaries of a SNP‐free region around a polymorphic 1.17 Mbp inverted segment on the short arm of chromosome 4. The architecture and functionality of iBrowser makes the software appropriate for a broad set of analyses including SNP mining, genome structure analysis, and pedigree analysis. Its functionality, together with the capability to process large data sets and efficient visualization of sequence variation, makes iBrowser a valuable breeding tool.     

Systematic identification of cancer driving signaling pathways based on mutual exclusivity of genomic alterations

We present a novel method for the identification of sets of mutually exclusive gene alterations in a given set of genomic profiles. We scan the groups of genes with a common downstream effect on the signaling network, using a mutual exclusivity criterion that ensures that each gene in the group significantly contributes to the mutual exclusivity pattern. We test the method on all available TCGA cancer genomics datasets, and detect multiple previously unreported alterations that show significant mutual exclusivity and are likely to be driver events.

Electronic supplementary material

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

Genomic landscape of rat strain and substrain variation

Background

Since the completion of the rat reference genome in 2003, whole-genome sequencing data from more than 40 rat strains have become available. These data represent the broad range of strains that are used in rat research including commonly used substrains. Currently, this wealth of information cannot be used to its full extent, because the variety of different variant calling algorithms employed by different groups impairs comparison between strains. In addition, all rat whole genome sequencing studies to date used an outdated reference genome for analysis (RGSC3.4 released in 2004).

Results

Here we present a comprehensive, multi-sample and uniformly called set of genetic variants in 40 rat strains, including 19 substrains. We reanalyzed all primary data using a recent version of the rat reference assembly (RGSC5.0 released in 2012) and identified over 12 million genomic variants (SNVs, indels and structural variants) among the 40 strains. 28,318 SNVs are specific to individual substrains, which may be explained by introgression from other unsequenced strains and ongoing evolution by genetic drift. Substrain SNVs may have a larger predicted functional impact compared to older shared SNVs.

Conclusions

In summary we present a comprehensive catalog of uniformly analyzed genetic variants among 40 widely used rat inbred strains based on the RGSC5.0 assembly. This represents a valuable resource, which will facilitate rat functional genomic research. In line with previous observations, our genome-wide analyses do not show evidence for contribution of multiple ancestral founder rat subspecies to the currently used rat inbred strains, as is the case for mouse. In addition, we find that the degree of substrain variation is highly variable between strains, which is of importance for the correct interpretation of experimental data from different labs.

Electronic supplementary material

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

Dissecting the complete lipoprotein biogenesis pathway in Streptomyces scabies

Following translocation, bacterial lipoproteins are lipidated by lipoprotein diacylglycerol transferase (Lgt) and cleaved of their signal peptides by lipoprotein signal peptidase (Lsp). In Gram-negative bacteria and mycobacteria, lipoproteins are further lipidated by lipoprotein N-acyl transferase (Lnt), to give triacylated lipoproteins. Streptomyces are unusual amongst Gram-positive bacteria because they export large numbers of lipoproteins via the twin arginine protein transport (Tat) pathway. Furthermore, some Streptomyces species encode two Lgt homologues and all Streptomyces species encode two homologues of Lnt. Here we characterize lipoprotein biogenesis in the plant pathogen Streptomyces scabies and report that lgt and lsp mutants are defective in growth and development while only moderately affected in virulence. Lipoproteins are lost from the membrane in an S. scabies lgt mutant but restored by expression of Streptomyces coelicolor lgt1 or lgt2 confirming that both encode functional Lgt enzymes. Furthermore, lipoproteins are N-acylated in Streptomyces with efficient N-acylation dependent on Lnt1 and Lnt2. However, deletion of lnt1 and lnt2 has no effect on growth, development or virulence. We thus present a detailed study of lipoprotein biogenesis in Streptomyces, the first study of Lnt function in a monoderm bacterium and the first study of bacterial lipoproteins as virulence factors in a plant pathogen.     

Major remodelling of the murine stem cell kinome following differentiation in the hematopoietic compartment

The changes in signal transduction associated with the acquisition of specific cell fates remain poorly understood. We performed massive parallel assessment of kinase signatures of the radiations of the hematopoietic system, including long-term repopulating hematopoietic stem cells (LT-HSC), short-term repopulating HSC (ST-HSC), immature natural killer (iNK) cells, NK cells, B cells, T cells, and myeloid cells. The LT-HSC kinome is characterized by noncanonical Wnt, Ca(2+) and classical protein kinase C (PKC)-driven signaling, which is lost upon the transition to ST-HSC, whose kinome signature prominently features receptor tyrosine kinase (RTK) activation of the Ras/MAPK signaling cassette. Further differentiation to iNK maintains signaling through this cassette but simultaneously leads to activation of a PI3K/PKB/Rac signaling, which becomes the dominant trait in the kinase signature following full differentiation toward NK cells. Differentiation along the myeloid and B cell lineages is accompanied by hyperactivation of both the Ras/MAPK and PI3K/PKB/Rac signaling cassette. T cells, however, deactivate signaling and only display residual G protein-coupled pathways. Thus, differentiation along the hematopoietic lineage is associated with major remodelling of cellular kinase signature.