RNA-Seq mapping and detection of gene fusions with a suffix array algorithm.

High-throughput RNA sequencing enables quantification of transcripts (both known and novel), exon/exon junctions and fusions of exons from different genes. Discovery of gene fusions-particularly those expressed with low abundance- is a challenge with short- and medium-length sequencing reads. To address this challenge, we implemented an RNA-Seq mapping pipeline within the LifeScope software. We introduced new features including filter and junction mapping, annotation-aided pairing rescue and accurate mapping quality values. We combined this pipeline with a Suffix Array Spliced Read (SASR) aligner to detect chimeric transcripts. Performing paired-end RNA-Seq of the breast cancer cell line MCF-7 using the SOLiD system, we called 40 gene fusions among over 120,000 splicing junctions. We validated 36 of these 40 fusions with TaqMan assays, of which 25 were expressed in MCF-7 but not the Human Brain Reference. An intra-chromosomal gene fusion involving the estrogen receptor alpha gene ESR1, and another involving the RPS6KB1 (Ribosomal protein S6 kinase beta-1) were recurrently expressed in a number of breast tumor cell lines and a clinical tumor sample.

Sequence and structural variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding.

We describe the genome sequencing of an anonymous individual of African origin using a novel ligation-based sequencing assay that enables a unique form of error correction that improves the raw accuracy of the aligned reads to >99.9%, allowing us to accurately call SNPs with as few as two reads per allele. We collected several billion mate-paired reads yielding approximately 18x haploid coverage of aligned sequence and close to 300x clone coverage. Over 98% of the reference genome is covered with at least one uniquely placed read, and 99.65% is spanned by at least one uniquely placed mate-paired clone. We identify over 3.8 million SNPs, 19% of which are novel. Mate-paired data are used to physically resolve haplotype phases of nearly two-thirds of the genotypes obtained and produce phased segments of up to 215 kb. We detect 226,529 intra-read indels, 5590 indels between mate-paired reads, 91 inversions, and four gene fusions. We use a novel approach for detecting indels between mate-paired reads that are smaller than the standard deviation of the insert size of the library and discover deletions in common with those detected with our intra-read approach. Dozens of mutations previously described in OMIM and hundreds of nonsynonymous single-nucleotide and structural variants in genes previously implicated in disease are identified in this individual. There is more genetic variation in the human genome still to be uncovered, and we provide guidance for future surveys in populations and cancer biopsies.

Genome-wide analysis of host responses to the Pseudomonas aeruginosa type III secretion system yields synergistic effects.

The type III secretion system (TTSS) is a dedicated bacterial pathogen protein targeting system that directly affects host cell signalling and response pathways. Our goal was to identify host responses to the Pseudomonas aeruginosa effectors, introduced into target cells utilizing the TTSS. We carried out expression profiling of a human lung pneumocyte cell line A549 exposed to isogenic mutants of P. aeruginosa PAK lacking individual or a combination of TTSS components. We then devised a data analysis method to isolate the key responses to specific secreted bacterial effector proteins as well as components of the TTSS machinery. Individually, the effector proteins elicited host responses consistent with their known functions, many of which were cell cycle-related. However, our analysis has shown that the effector proteins elicit a distinct host transcriptional response when present in combination, suggesting a synergistic effect. Furthermore, the pattern of host transcriptional responses is consistent with the pore forming ability of the TTSS needle complex. This study shows that the individual components of the TTSS define an integrated system and that a systems biology approach is required to fully understand the complex interplay between pathogen and host.

Influence of cystic fibrosis transmembrane conductance regulator on gene expression in response to Pseudomonas aeruginosa infection of human bronchial epithelial cells.

Chronic lung infection by Pseudomonas aeruginosa causes significant morbidity in cystic fibrosis patients initiated by the failure of innate immune responses. We used microarray analysis and real-time PCR to detect transcriptional changes associated with cytokine production in isogenic bronchial epithelial cell lines with either wild-type (WT) or mutant cystic fibrosis transmembrane conductance regulator (CFTR) in response to P. aeruginosa infection. The transcription of four NF-kappaB-regulated cytokine genes was maximal in the presence of WT CFTR: the interleukin-8 (IL-8), IL-6, CXCL1, and intracellular adhesion molecule 1 (ICAM-1) genes. Analysis of protein expression in two cell lines paired for wild-type and mutant CFTR with three P. aeruginosa strains showed IL-6 and IL-8 expressions were consistently enhanced by the presence of WT CFTR in both cell lines with all three strains of P. aeruginosa, although some strains gave small IL-8 increases in cells with mutant CFTR. CXCL1 production showed consistent enhancement in cells with WT CFTR using all three bacterial strains in one cell line, whereas in the other cell line, CXCL1 showed a significant increase in cells with either WT or mutant CFTR. ICAM-1 was unchanged at the protein level in one of the cell lines but did show mild enhancement with WT CFTR in the other cell pair. Inhibitions of NF-kappaB prior to infection indicated differing degrees of dependence on NF-kappaB for production of the cytokines, contingent on the cell line. Cytokine effectors of innate immunity to P. aeruginosa were found to be positively influenced by the presence of WT CFTR, indicating a role in resistance to P. aeruginosa infection.

Cystic fibrosis airway epithelial cell polarity and bacterial flagellin determine host response to Pseudomonas aeruginosa.

The role of epithelial polarity and bacterial factors in the control of the innate immune response of airway epithelial cells to Pseudomonas aeruginosa PAK was investigated using a human, nasal cystic fibrosis (DeltaF508/DeltaF508) epithelial cell line CF15 grown as confluent layers on permeable supports. Addition of PAK to the basal surface of CF15 layers caused significant expression changes in 1525 different genes (out of 12 625 examined), including the cytokines IL-6, IL-8, IL-1beta and TNF-alpha, as well as genes associated with leucocyte adhesion, antibacterial factors, and NF-kappaB signalling. Confocal microscopy showed that nuclear migration of NF-kappaB in all CF15 cells was preceded by PAK binding to the basal and lateral surfaces of some cells. Addition of PAK to the apical surface of CF15 monolayers elicited changes in expression of only 602 genes, including 256 not affected during basolateral PAK exposure. Over time, cytokine expression during apical PAK was similar to that exhibited by basal PAK, but the magnitudes during apical treatment were much smaller with little/no nuclear migration of NF-kappaB in CF15 cells. Furthermore, these responses depended on the presence of flagellin, but not pili on the bacteria. Thus, P. aeruginosa triggered a strong innate immune response that depended on the apical versus basolateral polarity of CF15 cells and the presence of flagellin on the bacteria.

Genome mosaicism is conserved but not unique in Pseudomonas aeruginosa isolates from the airways of young children with cystic fibrosis.

Pseudomonas aeruginosa strains from the chronic lung infections of cystic fibrosis (CF) patients are phenotypically and genotypically diverse. Using strain PAO1 whole genome DNA microarrays, we assessed the genomic variation in P. aeruginosa strains isolated from young children with CF (6 months to 8 years of age) as well as from the environment. Eighty-nine to 97% of the PAO1 open reading frames were detected in 20 strains by microarray analysis, while subsets of 38 gene islands were absent or divergent. No specific pattern of genome mosaicism defined strains associated with CF. Many mosaic regions were distinguished by their low G + C content; their inclusion of phage related or pyocin genes; or by their linkage to a vgr gene or a tRNA gene. Microarray and phenotypic analysis of sequential isolates from individual patients revealed two deletions of greater than 100 kbp formed during evolution in the lung. The gene loss in these sequential isolates raises the possibility that acquisition of pyomelanin production and loss of pyoverdin uptake each may be of adaptive significance. Further characterization of P. aeruginosa diversity within the airways of individual CF patients may reveal common adaptations, perhaps mediated by gene loss, that suggest new opportunities for therapy.