SAMBLASTER: fast duplicate marking and structural variant read extraction

Global analysis of genomic instability caused by DNA replication stress in Saccharomyces cerevisiae

DNA replication stress (DRS)-induced genomic instability is an important factor driving cancer development. To understand the mechanisms of DRS-associated genomic instability, we measured the rates of genomic alterations throughout the genome in a yeast strain with lowered expression of the replicative DNA polymerase δ. By a genetic test, we showed that most recombinogenic DNA lesions were introduced during S or G₂ phase, presumably as a consequence of broken replication forks. We observed a high rate of chromosome loss, likely reflecting a reduced capacity of the low-polymerase strains to repair double-stranded DNA breaks (DSBs). We also observed a high frequency of deletion events within tandemly repeated genes such as the ribosomal RNA genes. By whole-genome sequencing, we found that low levels of DNA polymerase δ elevated mutation rates, both single-base mutations and small insertions/deletions. Finally, we showed that cells with low levels of DNA polymerase δ tended to accumulate small promoter mutations that increased the expression of this polymerase. These deletions conferred a selective growth advantage to cells, demonstrating that DRS can be one factor driving phenotypic evolution.

Spatial intratumoral heterogeneity and temporal clonal evolution in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is among the most common malignancies, but little is known about its spatial intratumor heterogeneity (ITH) and temporal clonal evolutionary processes. To address this, we performed multiregion whole-exome sequencing on 51 tumor regions from 13 ESCCs, and multiregion global methylation profiling on three of these 13 cases. We found an average of 35.8% heterogeneous somatic mutations with strong evidence of ITH. Half of driver mutations located on the branches targeted oncogenes, including PIK3CA, NFE2L2, MTOR, etc. By contrast, the majority of truncal and clonal driver mutations occurred in tumor suppressor genes, including TP53, KMT2D, ZNF750, etc. Interestingly, the phyloepigenetic trees robustly recapitulated the topologic structures of the phylogenetic ones, indicating the possible relationship between genetic and epigenetic alterations. Our integrated investigations of the spatial ITH and clonal evolution provide an important molecular foundation for enhanced understanding of the tumorigenesis and progression of ESCC.

The Evolutionary Fates of a Large Segmental Duplication in Mouse

Gene duplication and loss are major sources of genetic polymorphism in populations, and are important forces shaping the evolution of genome content and organization. We have reconstructed the origin and history of a 127-kbp segmental duplication, R2d, in the house mouse (Mus musculus). R2d contains a single protein-coding gene, Cwc22 De novo assembly of both the ancestral (R2d1) and the derived (R2d2) copies reveals that they have been subject to nonallelic gene conversion events spanning tens of kilobases. R2d2 is also a hotspot for structural variation: its diploid copy number ranges from zero in the mouse reference genome to >80 in wild mice sampled from around the globe. Hemizygosity for high copy-number alleles of R2d2 is associated in cis with meiotic drive; suppression of meiotic crossovers; and copy-number instability, with a mutation rate in excess of 1 per 100 transmissions in some laboratory populations. Our results provide a striking example of allelic diversity generated by duplication and demonstrate the value of de novo assembly in a phylogenetic context for understanding the mutational processes affecting duplicate genes.

Loss of BRCA1 or BRCA2 markedly increases the rate of base substitution mutagenesis and has distinct effects on genomic deletions

Loss-of-function mutations in the BRCA1 and BRCA2 genes increase the risk of cancer. Owing to their function in homologous recombination repair, much research has focused on the unstable genomic phenotype of BRCA1/2 mutant cells manifest mainly as large-scale rearrangements. We used whole-genome sequencing of multiple isogenic chicken DT40 cell clones to precisely determine the consequences of BRCA1/2 loss on all types of genomic mutagenesis. Spontaneous base substitution mutation rates increased sevenfold upon the disruption of either BRCA1 or BRCA2, and the arising mutation spectra showed strong and specific correlation with a mutation signature associated with BRCA1/2 mutant tumours. To model endogenous alkylating damage, we determined the mutation spectrum caused by methyl methanesulfonate (MMS), and showed that MMS also induces more base substitution mutations in BRCA1/2-deficient cells. Spontaneously arising and MMS-induced insertion/deletion mutations and large rearrangements were also more common in BRCA1/2 mutant cells compared with the wild-type control. A difference in the short deletion phenotypes of BRCA1 and BRCA2 suggested distinct roles for the two proteins in the processing of DNA lesions, as BRCA2 mutants contained more short deletions, with a wider size distribution, which frequently showed microhomology near the breakpoints resembling repair by non-homologous end joining. An increased and prolonged gamma-H2AX signal in MMS-treated BRCA1/2 cells suggested an aberrant processing of stalled replication forks as the cause of increased mutagenesis. The high rate of base substitution mutagenesis demonstrated by our experiments is likely to significantly contribute to the oncogenic effect of the inactivation of BRCA1 or BRCA2.

Specific amino acids in HIV-1 Vpr are significantly associated with differences in patient neurocognitive status

Even in the era of combination antiretroviral therapies used to combat human immunodeficiency virus type 1 (HIV-1) infection, up to 50 % of well-suppressed HIV-1-infected patients are still diagnosed with mild neurological deficits referred to as HIV-associated neurocognitive disorders (HAND). The multifactorial nature of HAND likely involves the HIV-1 accessory protein viral protein R (Vpr) as an agent of neuropathogenesis. To investigate the effect of naturally occurring variations in Vpr on HAND in well-suppressed HIV-1-infected patients, bioinformatic analyses were used to correlate peripheral blood-derived Vpr sequences with patient neurocognitive performance, as measured by comprehensive neuropsychological assessment and the resulting Global Deficit Score (GDS). Our studies revealed unique associations between GDS and the presence of specific amino acid changes in peripheral blood-derived Vpr sequences [neuropsychological impairment Vpr (niVpr) variants]. Amino acids N41 and A55 in the Vpr sequence were associated with more pronounced neurocognitive deficits (higher GDS). In contrast, amino acids I37 and S41 were connected to measurably lower GDS. All niVpr variants were also detected in DNA isolated from HIV-1-infected brain tissues. The implication of these results is that niVpr variants alter the genesis and/or progression of HAND through differences in Vpr-mediated effects in the peripheral blood and/or the brain.

Third Chromosome Balancer Inversions Disrupt Protein-Coding Genes and Influence Distal Recombination Events in Drosophila melanogaster

Balancer chromosomes are multiply inverted chromosomes that suppress meiotic crossing over and prevent the recovery of crossover products. Balancers are commonly used in Drosophila melanogaster to maintain deleterious alleles and in stock construction. They exist for all three major chromosomes, yet the molecular location of the breakpoints and the exact nature of many of the mutations carried by the second and third chromosome balancers has not been available. Here, we precisely locate eight of 10 of the breakpoints on the third chromosome balancer TM3, six of eight on TM6, and nine of 11 breakpoints on TM6B. We find that one of the inversion breakpoints on TM3 bisects the highly conserved tumor suppressor gene p53—a finding that may have important consequences for a wide range of studies in Drosophila. We also identify evidence of single and double crossovers between several TM3 and TM6B balancers and their normal-sequence homologs that have created genetic diversity among these chromosomes. Overall, this work demonstrates the practical importance of precisely identifying the position of inversion breakpoints of balancer chromosomes and characterizing the mutant alleles carried by them.

Extensive sequencing of seven human genomes to characterize benchmark reference materials

The Genome in a Bottle Consortium, hosted by the National Institute of Standards and Technology (NIST) is creating reference materials and data for human genome sequencing, as well as methods for genome comparison and benchmarking. Here, we describe a large, diverse set of sequencing data for seven human genomes; five are current or candidate NIST Reference Materials. The pilot genome, NA12878, has been released as NIST RM 8398. We also describe data from two Personal Genome Project trios, one of Ashkenazim Jewish ancestry and one of Chinese ancestry. The data come from 12 technologies: BioNano Genomics, Complete Genomics paired-end and LFR, Ion Proton exome, Oxford Nanopore, Pacific Biosciences, SOLiD, 10X Genomics GemCode WGS, and Illumina exome and WGS paired-end, mate-pair, and synthetic long reads. Cell lines, DNA, and data from these individuals are publicly available. Therefore, we expect these data to be useful for revealing novel information about the human genome and improving sequencing technologies, SNP, indel, and structural variant calling, and de novo assembly.

Identification of ETV6-RUNX1-like and DUX4-rearranged subtypes in paediatric B-cell precursor acute lymphoblastic leukaemia

Fusion genes are potent driver mutations in cancer. In this study, we delineate the fusion gene landscape in a consecutive series of 195 paediatric B-cell precursor acute lymphoblastic leukaemia (BCP ALL). Using RNA sequencing, we find in-frame fusion genes in 127 (65%) cases, including 27 novel fusions. We describe a subtype characterized by recurrent IGH-DUX4 or ERG-DUX4 fusions, representing 4% of cases, leading to overexpression of DUX4 and frequently co-occurring with intragenic ERG deletions. Furthermore, we identify a subtype characterized by an ETV6-RUNX1-like gene-expression profile and coexisting ETV6 and IKZF1 alterations. Thus, this study provides a detailed overview of fusion genes in paediatric BCP ALL and adds new pathogenetic insights, which may improve risk stratification and provide therapeutic options for this disease.

DeepSNVMiner: a sequence analysis tool to detect emergent, rare mutations in subsets of cell populations

Background. Massively parallel sequencing technology is being used to sequence highly diverse populations of DNA such as that derived from heterogeneous cell mixtures containing both wild-type and disease-related states. At the core of such molecule tagging techniques is the tagging and identification of sequence reads derived from individual input DNA molecules, which must be first computationally disambiguated to generate read groups sharing common sequence tags, with each read group representing a single input DNA molecule. This disambiguation typically generates huge numbers of reads groups, each of which requires additional variant detection analysis steps to be run specific to each read group, thus representing a significant computational challenge. While sequencing technologies for producing these data are approaching maturity, the lack of available computational tools for analysing such heterogeneous sequence data represents an obstacle to the widespread adoption of this technology. Results. Using synthetic data we successfully detect unique variants at dilution levels of 1 in a 1,000,000 molecules, and find DeeepSNVMiner obtains significantly lower false positive and false negative rates compared to popular variant callers GATK, SAMTools, FreeBayes and LoFreq, particularly as the variant concentration levels decrease. In a dilution series with genomic DNA from two cells lines, we find DeepSNVMiner identifies a known somatic variant when present at concentrations of only 1 in 1,000 molecules in the input material, the lowest concentration amongst all variant callers tested. Conclusions. Here we present DeepSNVMiner; a tool to disambiguate tagged sequence groups and robustly identify sequence variants specific to subsets of starting DNA molecules that may indicate the presence of a disease. DeepSNVMiner is an automated workflow of custom sequence analysis utilities and open source tools able to differentiate somatic DNA variants from artefactual sequence variants that likely arose during DNA amplification. The workflow remains flexible such that it may be customised to variants of the data production protocol used, and supports reproducible analysis through detailed logging and reporting of results. DeepSNVMiner is available for academic non-commercial research purposes at https://github.com/mattmattmattmatt/DeepSNVMiner.

Genotyping by Sequencing Using Specific Allelic Capture to Build a High-Density Genetic Map of Durum Wheat

Targeted sequence capture is a promising technology which helps reduce costs for sequencing and genotyping numerous genomic regions in large sets of individuals. Bait sequences are designed to capture specific alleles previously discovered in parents or reference populations. We studied a set of 135 RILs originating from a cross between an emmer cultivar (Dic2) and a recent durum elite cultivar (Silur). Six thousand sequence baits were designed to target Dic2 vs. Silur polymorphisms discovered in a previous RNAseq study. These baits were exposed to genomic DNA of the RIL population. Eighty percent of the targeted SNPs were recovered, 65% of which were of high quality and coverage. The final high density genetic map consisted of more than 3,000 markers, whose genetic and physical mapping were consistent with those obtained with large arrays.

A comprehensive survey of the mutagenic impact of common cancer cytotoxics

BACKGROUND

Genomic mutations caused by cytotoxic agents used in cancer chemotherapy may cause secondary malignancies as well as contribute to the evolution of treatment-resistant tumour cells. The stable diploid genome of the chicken DT40 lymphoblast cell line, an established DNA repair model system, is well suited to accurately assay genomic mutations.

RESULTS

We use whole genome sequencing of multiple DT40 clones to determine the mutagenic effect of eight common cytotoxics used for the treatment of millions of patients worldwide. We determine the spontaneous mutagenesis rate at 2.3 × 10(-10) per base per cell division and find that cisplatin, cyclophosphamide and etoposide induce extra base substitutions with distinct spectra. After four cycles of exposure, cisplatin induces 0.8 mutations per Mb, equivalent to the median mutational burden in common leukaemias. Cisplatin-induced mutations, including short insertions and deletions, are mainly located at sites of putative intrastrand crosslinks. We find two of the newly defined cisplatin-specific mutation types as causes of the reversion of BRCA2 mutations in emerging cisplatin-resistant tumours or cell clones. Gemcitabine, 5-fluorouracil, hydroxyurea, doxorubicin and paclitaxel have no measurable mutagenic effect. The cisplatin-induced mutation spectrum shows good correlation with cancer mutation signatures attributed to smoking and other sources of guanine-directed base damage.

CONCLUSION

This study provides support for the use of cell line mutagenesis assays to validate or predict the mutagenic effect of environmental and iatrogenic exposures. Our results suggest genetic reversion due to cisplatin-induced mutations as a distinct mechanism for developing resistance.

Transformed Recombinant Enrichment Profiling Rapidly Identifies HMW1 as an Intracellular Invasion Locus in Haemophilus influenza

Many bacterial species actively take up and recombine homologous DNA into their genomes, called natural competence, a trait that offers a means to identify the genetic basis of naturally occurring phenotypic variation. Here, we describe “transformed recombinant enrichment profiling” (TREP), in which natural transformation is used to generate complex pools of recombinants, phenotypic selection is used to enrich for specific recombinants, and deep sequencing is used to survey for the genetic variation responsible. We applied TREP to investigate the genetic architecture of intracellular invasion by the human pathogen Haemophilus influenzae, a trait implicated in persistence during chronic infection. TREP identified the HMW1 adhesin as a crucial factor. Natural transformation of the hmw1 operon from a clinical isolate (86-028NP) into a laboratory isolate that lacks it (Rd KW20) resulted in ~1,000-fold increased invasion into airway epithelial cells. When a distinct recipient (Hi375, already possessing hmw1 and its paralog hmw2) was transformed by the same donor, allelic replacement of hmw2A_Hi375 by hmw1A_86-028NP resulted in a ~100-fold increased intracellular invasion rate. The specific role of hmw1A_86-028NP was confirmed by mutant and western blot analyses. Bacterial self-aggregation and adherence to airway cells were also increased in recombinants, suggesting that the high invasiveness induced by hmw1A_86-028NP might be a consequence of these phenotypes. However, immunofluorescence results found that intracellular hmw1A_86-028NP bacteria likely invaded as groups, instead of as individual bacterial cells, indicating an emergent invasion-specific consequence of hmw1A-mediated self-aggregation.

Many bacteria are naturally competent, actively taking up DNA from their surroundings and incorporating it into their genomes by homologous recombination. This cellular process has had a large impact on the evolution of these species, for example by enabling pathogens to acquire virulence factors and antibiotic resistances from their relatives. But natural competence can also be exploited by researchers to identify the underlying genetic variation responsible for naturally varying phenotypic traits, similar to how eukaryotic geneticists use meiotic recombination during sexual reproduction to create genetically admixed populations. Here we exploited natural competence, phenotypic selection, and deep sequencing to rapidly identify the hmw1 locus as a major contributor to intracellular invasion of airway epithelial cells by the human pathogen Haemophilus influenzae, a trait that likely allows bacterial cells to evade the immune system and therapeutic interventions during chronic infections. Genetic variation in this locus can strongly modulate bacterial intracellular invasion rates, and possession of a certain allele favors adhesion and self-aggregation, which appear to prompt bacteria to invade airway cells as groups, rather than as individuals. Overall, our findings indicate that targeting HMW1 could block the ability of H. influenzae to invade airway cells, which would make antibiotic therapy to treat chronic lung infections more effective. Furthermore, our new approach to identifying the genetic basis of natural phenotypic variation is applicable to a wide-range of phenotypically selectable traits within the widely distributed naturally competent bacterial species, including pathogenesis traits in many human pathogens.

Comprehensive mutational analysis of primary and relapse acute promyelocytic leukemia

Acute promyelocytic leukemia (APL) is a subtype of myeloid leukemia characterized by differentiation block at the promyelocyte stage. Besides the presence of chromosomal rearrangement t(15;17), leading to the formation of PML-RARA (promyelocytic leukemia-retinoic acid receptor alpha) fusion, other genetic alterations have also been implicated in APL. Here, we performed comprehensive mutational analysis of primary and relapse APL to identify somatic alterations, which cooperate with PML-RARA in the pathogenesis of APL. We explored the mutational landscape using whole-exome (n=12) and subsequent targeted sequencing of 398 genes in 153 primary and 69 relapse APL. Both primary and relapse APL harbored an average of eight non-silent somatic mutations per exome. We observed recurrent alterations of FLT3, WT1, NRAS and KRAS in the newly diagnosed APL, whereas mutations in other genes commonly mutated in myeloid leukemia were rarely detected. The molecular signature of APL relapse was characterized by emergence of frequent mutations in PML and RARA genes. Our sequencing data also demonstrates incidence of loss-of-function mutations in previously unidentified genes, ARID1B and ARID1A, both of which encode for key components of the SWI/SNF complex. We show that knockdown of ARID1B in APL cell line, NB4, results in large-scale activation of gene expression and reduced in vitro differentiation potential.

Whole genome sequencing in cats, identifies new models for blindness in AIPL1 and somite segmentation in HES7

BACKGROUND

The reduced cost and improved efficiency of whole genome sequencing (WGS) is drastically improving the development of cats as biomedical models. Persian cats are models for Leber's congenital amaurosis (LCA), the most severe and earliest onset form of visual impairment in humans. Cats with innocuous breed-defining traits, such as a bobbed tail, can also be models for somite segmentation and vertebral column development.

METHODS

The first WGS in cats was conducted on a trio segregating for LCA and the bobbed tail abnormality. Variants were identified using FreeBayes and effects predicted using SnpEff. Variants within a known haplotype block for cat LCA and specific candidate genes for both phenotypes were prioritized by the predicted variant effect on the proteins and concordant segregation within the trio. The efficiency of WGS of a single trio of domestic cats was evaluated.

RESULTS

A stop gain was identified at position c.577C > T in cat AIPL1, a predicted p.Arg193*. A c.5A > G variant causing a p.V2A was identified in HES7. The variants segregated concordantly in a Persian - Japanese bobtail pedigree. Over 1700 cats from 40 different breeds and populations were genotyped for the AIPL1 variant, defining an allelic frequency in only Persian -related breeds of 1.15%. A sub-set of cats was genotyped for the HES7 variant, supporting the variant as private to the Japanese bobtail breed. Approximately 18 million SNPs were identified for application in cat research. The cat AIPL1 variant would have been considered a high priority variant for evaluation, regardless of a priori knowledge from previous genetic studies.

CONCLUSIONS

This study represents the first effort of the 99 Lives Cat Genome Sequencing Initiative to identify disease--causing variants in the domestic cat using WGS. The current cat reference assembly is efficient for gene and variant identification. However, as the feline variant database improves, development of cats as biomedical models for human disease will be more efficient, providing an alternative, large animal model for drug and gene therapy trials. Undiagnosed human patients with early-onset blindness should be screened for this AIPL1 variant. The HES7 variant should further calibrate the somite segmentation clock.

Rare recombination events generate sequence diversity among balancer chromosomes in Drosophila melanogaster

Multiply inverted balancer chromosomes that suppress exchange with their homologs are an essential part of the Drosophila melanogaster genetic toolkit. Despite their widespread use, the organization of balancer chromosomes has not been characterized at the molecular level, and the degree of sequence variation among copies of balancer chromosomes is unknown. To map inversion breakpoints and study potential diversity in descendants of a structurally identical balancer chromosome, we sequenced a panel of laboratory stocks containing the most widely used X chromosome balancer, First Multiple 7 (FM7). We mapped the locations of FM7 breakpoints to precise euchromatic coordinates and identified the flanking sequence of breakpoints in heterochromatic regions. Analysis of SNP variation revealed megabase-scale blocks of sequence divergence among currently used FM7 stocks. We present evidence that this divergence arose through rare double-crossover events that replaced a female-sterile allele of the singed gene (sn(X2)) on FM7c with a sequence from balanced chromosomes. We propose that although double-crossover events are rare in individual crosses, many FM7c chromosomes in the Bloomington Drosophila Stock Center have lost sn(X2) by this mechanism on a historical timescale. Finally, we characterize the original allele of the Bar gene (B(1)) that is carried on FM7, and validate the hypothesis that the origin and subsequent reversion of the B(1) duplication are mediated by unequal exchange. Our results reject a simple nonrecombining, clonal mode for the laboratory evolution of balancer chromosomes and have implications for how balancer chromosomes should be used in the design and interpretation of genetic experiments in Drosophila.

Propionibacterium acnes: Disease-Causing Agent or Common Contaminant? Detection in Diverse Patient Samples by Next-Generation Sequencing

Propionibacterium acnes is the most abundant bacterium on human skin, particularly in sebaceous areas. P. acnes is suggested to be an opportunistic pathogen involved in the development of diverse medical conditions but is also a proven contaminant of human clinical samples and surgical wounds. Its significance as a pathogen is consequently a matter of debate. In the present study, we investigated the presence of P. acnes DNA in 250 next-generation sequencing data sets generated from 180 samples of 20 different sample types, mostly of cancerous origin. The samples were subjected to either microbial enrichment, involving nuclease treatment to reduce the amount of host nucleic acids, or shotgun sequencing. We detected high proportions of P. acnes DNA in enriched samples, particularly skin tissue-derived and other tissue samples, with the levels being higher in enriched samples than in shotgun-sequenced samples. P. acnes reads were detected in most samples analyzed, though the proportions in most shotgun-sequenced samples were low. Our results show that P. acnes can be detected in practically all sample types when molecular methods, such as next-generation sequencing, are employed. The possibility of contamination from the patient or other sources, including laboratory reagents or environment, should therefore always be considered carefully when P. acnes is detected in clinical samples. We advocate that detection of P. acnes always be accompanied by experiments validating the association between this bacterium and any clinical condition.

SpeedSeq: ultra-fast personal genome analysis and interpretation

SpeedSeq is an open-source genome analysis platform that accomplishes alignment, variant detection and functional annotation of a 50× human genome in 13 h on a low-cost server and alleviates a bioinformatics bottleneck that typically demands weeks of computation with extensive hands-on expert involvement. SpeedSeq offers performance competitive with or superior to current methods for detecting germline and somatic single-nucleotide variants, structural variants, insertions and deletions, and it includes novel functionality for streamlined interpretation.

COLQ variant associated with Devon Rex and Sphynx feline hereditary myopathy

Some Devon Rex and Sphynx cats have a variably progressive myopathy characterized by appendicular and axial muscle weakness, megaesophagus, pharyngeal weakness and fatigability with exercise. Muscle biopsies from affected cats demonstrated variable pathological changes ranging from dystrophic features to minimal abnormalities. Affected cats have exacerbation of weakness following anticholinesterase dosing, a clue that there is an underlying congenital myasthenic syndrome (CMS). A genome-wide association study and whole-genome sequencing suggested a causal variant for this entity was a c.1190G>A variant causing a cysteine to tyrosine substitution (p.Cys397Tyr) within the C-terminal domain of collagen-like tail subunit (single strand of homotrimer) of asymmetric acetylcholinesterase (COLQ). Alpha-dystroglycan expression, which is associated with COLQ anchorage at the motor end-plate, has been shown to be deficient in affected cats. Eighteen affected cats were identified by genotyping, including cats from the original clinical descriptions in 1993 and subsequent publications. Eight Devon Rex and one Sphynx not associated with the study were identified as carriers, suggesting an allele frequency of ~2.0% in Devon Rex. Over 350 tested cats from other breeds did not have the variant. Characteristic clinical features and variant presence in all affected cats suggest a model for COLQ CMS. The association between the COLQ variant and this CMS affords clinicians the opportunity to confirm diagnosis via genetic testing and permits owners and breeders to identify carriers in the population. Moreover, accurate diagnosis increases available therapeutic options for affected cats based on an understanding of the pathophysiology and experience from human CMS associated with COLQ variants.

ISMapper: identifying transposase insertion sites in bacterial genomes from short read sequence data

Background

Insertion sequences (IS) are small transposable elements, commonly found in bacterial genomes. Identifying the location of IS in bacterial genomes can be useful for a variety of purposes including epidemiological tracking and predicting antibiotic resistance. However IS are commonly present in multiple copies in a single genome, which complicates genome assembly and the identification of IS insertion sites. Here we present ISMapper, a mapping-based tool for identification of the site and orientation of IS insertions in bacterial genomes, directly from paired-end short read data.

Results

ISMapper was validated using three types of short read data: (i) simulated reads from a variety of species, (ii) Illumina reads from 5 isolates for which finished genome sequences were available for comparison, and (iii) Illumina reads from 7 Acinetobacter baumannii isolates for which predicted IS locations were tested using PCR. A total of 20 genomes, including 13 species and 32 distinct IS, were used for validation. ISMapper correctly identified 97 % of known IS insertions in the analysis of simulated reads, and 98 % in real Illumina reads. Subsampling of real Illumina reads to lower depths indicated ISMapper was able to correctly detect insertions for average genome-wide read depths >20x, although read depths >50x were required to obtain confident calls that were highly-supported by evidence from reads. All ISAba1 insertions identified by ISMapper in the A. baumannii genomes were confirmed by PCR. In each A. baumannii genome, ISMapper successfully identified an IS insertion upstream of the ampC beta-lactamase that could explain phenotypic resistance to third-generation cephalosporins. The utility of ISMapper was further demonstrated by profiling genome-wide IS6110 insertions in 138 publicly available Mycobacterium tuberculosis genomes, revealing lineage-specific insertions and multiple insertion hotspots.

Conclusions

ISMapper provides a rapid and robust method for identifying IS insertion sites directly from short read data, with a high degree of accuracy demonstrated across a wide range of bacteria.

FermiKit: assembly-based variant calling for Illumina resequencing data

FermiKit is a variant calling pipeline for Illumina whole-genome germline data. It de novo assembles short reads and then maps the assembly against a reference genome to call SNPs, short insertions/deletions and structural variations. FermiKit takes about one day to assemble 30-fold human whole-genome data on a modern 16-core server with 85 GB RAM at the peak, and calls variants in half an hour to an accuracy comparable to the current practice. FermiKit assembly is a reduced representation of raw data while retaining most of the original information.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lh3/fermikit

CONTACT

hengli@broadinstitute.org.

elPrep: High-Performance Preparation of Sequence Alignment/Map Files for Variant Calling

elPrep is a high-performance tool for preparing sequence alignment/map files for variant calling in sequencing pipelines. It can be used as a replacement for SAMtools and Picard for preparation steps such as filtering, sorting, marking duplicates, reordering contigs, and so on, while producing identical results. What sets elPrep apart is its software architecture that allows executing preparation pipelines by making only a single pass through the data, no matter how many preparation steps are used in the pipeline. elPrep is designed as a multithreaded application that runs entirely in memory, avoids repeated file I/O, and merges the computation of several preparation steps to significantly speed up the execution time. For example, for a preparation pipeline of five steps on a whole-exome BAM file (NA12878), we reduce the execution time from about 1:40 hours, when using a combination of SAMtools and Picard, to about 15 minutes when using elPrep, while utilising the same server resources, here 48 threads and 23GB of RAM. For the same pipeline on whole-genome data (NA12878), elPrep reduces the runtime from 24 hours to less than 5 hours. As a typical clinical study may contain sequencing data for hundreds of patients, elPrep can remove several hundreds of hours of computing time, and thus substantially reduce analysis time and cost.

Sambamba: fast processing of NGS alignment formats

Summary: Sambamba is a high-performance robust tool and library for working with SAM, BAM and CRAM sequence alignment files; the most common file formats for aligned next generation sequencing data. Sambamba is a faster alternative to samtools that exploits multi-core processing and dramatically reduces processing time. Sambamba is being adopted at sequencing centers, not only because of its speed, but also because of additional functionality, including coverage analysis and powerful filtering capability.

Availability and implementation: Sambamba is free and open source software, available under a GPLv2 license. Sambamba can be downloaded and installed from http://www.open-bio.org/wiki/Sambamba.

Sambamba v0.5.0 was released with doi:10.5281/zenodo.13200.

Contact: j.c.p.prins@umcutrecht.nl

Churchill: an ultra-fast, deterministic, highly scalable and balanced parallelization strategy for the discovery of human genetic variation in clinical and population-scale genomics

While advances in genome sequencing technology make population-scale genomics a possibility, current approaches for analysis of these data rely upon parallelization strategies that have limited scalability, complex implementation and lack reproducibility. Churchill, a balanced regional parallelization strategy, overcomes these challenges, fully automating the multiple steps required to go from raw sequencing reads to variant discovery. Through implementation of novel deterministic parallelization techniques, Churchill allows computationally efficient analysis of a high-depth whole genome sample in less than two hours. The method is highly scalable, enabling full analysis of the 1000 Genomes raw sequence dataset in a week using cloud resources. http://churchill.nchri.org/.