RAF1 rearrangements are common in pancreatic acinar cell carcinomas.

There is now evidence that gene fusions activating the MAPK pathway are relatively common in pancreatic acinar cell carcinoma with potentially actionable BRAF or RET fusions being found in ~30%. We sought to investigate the incidence of RAF1 fusions in pancreatic malignancies with acinar cell differentiation. FISH testing for RAF1 was undertaken on 30 tumors comprising 25 'pure' acinar cell carcinomas, 2 mixed pancreatic acinar-neuroendocrine carcinomas, 1 mixed acinar cell-low grade neuroendocrine tumor and 2 pancreatoblastomas. RAF1 rearrangements were identified in 5 cases and confirmed by DNA and RNA sequencing to represent oncogenic fusions (GATM-RAF1, GOLGA4-RAF1, PDZRN3-RAF1, HERPUD1-RAF1 and TRIM33-RAF1) and to be mutually exclusive with BRAF and RET fusions, as well as KRAS mutations. Large genome-wide copy number changes were common and included 1q gain and/or 1p loss in all five RAF1 FISH-positive acinar cell carcinomas. RAF1 expression by immunohistochemistry was found in 3 of 5 (60%) of fusion-positive cases and no FISH-negative cases. Phospho-ERK1/2 expression was found in 4 of 5 RAF1-fusion-positive cases. Expression of both RAF1 and phospho-ERK1/2 was heterogeneous and often only detected at the tumor-stroma interface, thus limiting their clinical utility. We conclude that RAF1 gene rearrangements are relatively common in pancreatic acinar cell carcinomas (14.3% to 18.5% of cases) and can be effectively identified by FISH with follow up molecular testing. The combined results of several studies now indicate that BRAF, RET or RAF1 fusions occur in between one third and one-half of these tumors but are extremely rare in other pancreatic malignancies. As these fusions are potentially actionable with currently available therapies, a strong argument can be made to perform FISH or molecular testing on all pancreatic acinar cell carcinomas.

Decreased vancomycin susceptibility in Staphylococcus aureus caused by IS256 tempering of WalKR expression.

Vancomycin-intermediate Staphylococcus aureus (VISA) strains often arise by mutations in the essential two-component regulator walKR; however their impact on walKR function has not been definitively established. Here, we investigated 10 MRSA strains recovered serially after exposure of vancomycin-susceptible S. aureus (VSSA) JKD6009 to simulated human vancomycin dosing regimens (500 mg to 4,000 mg every 12 h) using a 10-day hollow fiber infection model. After continued exposure to the vancomycin regimens, two isolates displayed reduced susceptibility to both vancomycin and daptomycin, developing independent IS256 insertions in the walKR 5' untranslated region (5' UTR). Quantitative reverse transcription-PCR (RT-PCR) revealed a 50% reduction in walKR gene expression in the IS256 mutants compared to the VSSA parent. Green fluorescent protein (GFP) reporter analysis, promoter mapping, and site-directed mutagenesis confirmed these findings and showed that the IS256 insertions had replaced two SigA-like walKR promoters with weaker, hybrid promoters. Removal of IS256 reverted the phenotype to VSSA, showing that reduced expression of WalKR did induce the VISA phenotype. Analysis of selected WalKR-regulated autolysins revealed upregulation of ssaA but no change in expression of sak and sceD in both IS256 mutants. Whole-genome sequencing of the two mutants revealed an additional IS256 insertion within agrC for one mutant, and we confirmed that this mutation abolished agr function. These data provide the first substantial analysis of walKR promoter function and show that prolonged vancomycin exposure can result in VISA through an IS256-mediated reduction in walKR expression; however, the mechanisms by which this occurs remain to be determined.

Evolution of multidrug resistance during Staphylococcus aureus infection involves mutation of the essential two component regulator WalKR.

Antimicrobial resistance in Staphylococcus aureus is a major public health threat, compounded by emergence of strains with resistance to vancomycin and daptomycin, both last line antimicrobials. Here we have performed high throughput DNA sequencing and comparative genomics for five clinical pairs of vancomycin-susceptible (VSSA) and vancomycin-intermediate ST239 S. aureus (VISA); each pair isolated before and after vancomycin treatment failure. These comparisons revealed a frequent pattern of mutation among the VISA strains within the essential walKR two-component regulatory locus involved in control of cell wall metabolism. We then conducted bi-directional allelic exchange experiments in our clinical VSSA and VISA strains and showed that single nucleotide substitutions within either walK or walR lead to co-resistance to vancomycin and daptomycin, and caused the typical cell wall thickening observed in resistant clinical isolates. Ion Torrent genome sequencing confirmed no additional regulatory mutations had been introduced into either the walR or walK VISA mutants during the allelic exchange process. However, two potential compensatory mutations were detected within putative transport genes for the walK mutant. The minimal genetic changes in either walK or walR also attenuated virulence, reduced biofilm formation, and led to consistent transcriptional changes that suggest an important role for this regulator in control of central metabolism. This study highlights the dramatic impacts of single mutations that arise during persistent S. aureus infections and demonstrates the role played by walKR to increase drug resistance, control metabolism and alter the virulence potential of this pathogen.

Proteogenomic analysis of polymorphisms and gene annotation divergences in prokaryotes using a clustered mass spectrometry-friendly database.

Precise annotation of genes or open reading frames is still a difficult task that results in divergence even for data generated from the same genomic sequence. This has an impact in further proteomic studies, and also compromises the characterization of clinical isolates with many specific genetic variations that may not be represented in the selected database. We recently developed software called multistrain mass spectrometry prokaryotic database builder (MSMSpdbb) that can merge protein databases from several sources and be applied on any prokaryotic organism, in a proteomic-friendly approach. We generated a database for the Mycobacterium tuberculosis complex (using three strains of Mycobacterium bovis and five of M. tuberculosis), and analyzed data collected from two laboratory strains and two clinical isolates of M. tuberculosis. We identified 2561 proteins, of which 24 were present in M. tuberculosis H37Rv samples, but not annotated in the M. tuberculosis H37Rv genome. We were also able to identify 280 nonsynonymous single amino acid polymorphisms and confirm 367 translational start sites. As a proof of concept we applied the database to whole-genome DNA sequencing data of one of the clinical isolates, which allowed the validation of 116 predicted single amino acid polymorphisms and the annotation of 131 N-terminal start sites. Moreover we identified regions not present in the original M. tuberculosis H37Rv sequence, indicating strain divergence or errors in the reference sequence. In conclusion, we demonstrated the potential of using a merged database to better characterize laboratory or clinical bacterial strains.

Using a label-free proteomics method to identify differentially abundant proteins in closely related hypo- and hypervirulent clinical Mycobacterium tuberculosis Beijing isolates.

Although the genome of the Mycobacterium tuberculosis H37Rv laboratory strain has been available for over 10 years, it is only recently that genomic information from clinical isolates has been used to generate the hypothesis of virulence differences between different strains. In addition, the relationship between strains displaying differing virulence in an epidemiological setting and their behavior in animal models has received little attention. The potential causes for variation in virulence between strains, as determined by differential protein expression, have similarly been a neglected area of investigation. In this study, we used a label-free quantitative proteomics approach to estimate differences in protein abundance between two closely related Beijing genotypes that have been shown to be hyper- and hypovirulent on the basis of both epidemiological and mouse model studies. We were able to identify a total of 1668 proteins from both samples, and protein abundance calculations revealed that 48 proteins were over-represented in the hypovirulent isolate, whereas 53 were over-represented in the hypervirulent. Functional classification of these results shows that molecules of cell wall organization and DNA transcription regulatory proteins may have a critical influence in defining the level of virulence. The reduction in the presence of ESAT-6, other Esx-like proteins, and FbpD (MPT51) in the hypervirulent strain indicates that changes in the repertoire of highly immunogenic proteins can be a defensive process undertaken by the virulent cell. In addition, most of the previously well characterized gene targets related to virulence were found to be similarly expressed in our model. Our data support the use of proteomics as a complementary tool for genomic comparisons to understand the biology of M. tuberculosis virulence.

Rifampicin reduces susceptibility to ofloxacin in rifampicin-resistant Mycobacterium tuberculosis through efflux.

RATIONALE: Central dogma suggests that rifampicin resistance in Mycobacterium tuberculosis develops solely through rpoB gene mutations. OBJECTIVE: To determine whether rifampicin induces efflux pumps activation in rifampicin resistant M. tuberculosis strains thereby defining rifampicin resistance levels and reducing ofloxacin susceptibility. METHODS: Rifampicin and/or ofloxacin minimum inhibitory concentrations (MICs) were determined in rifampicin resistant strains by culture in BACTEC 12B medium. Verapamil and reserpine were included to determine their effect on rifampicin and ofloxacin susceptibility. RT-qPCR was applied to assess expression of efflux pump/transporter genes after rifampicin exposure. To determine whether verapamil could restore susceptibility to first-line drugs, BALB/c mice were infected with a MDR-TB strain and treated with first-line drugs with/without verapamil. MEASUREMENTS AND MAIN FINDINGS: Rifampicin MICs varied independently of rpoB mutation and genetic background. Addition reserpine and verapamil significantly restored rifampicin susceptibility (p = 0.0000). RT-qPCR demonstrated that rifampicin induced differential expression of efflux/transporter genes in MDR-TB isolates. Incubation of rifampicin mono-resistant strains in rifampicin (2 µg/ml) for 7 days induced ofloxacin resistance (MIC > 2 µg/ml) in strains with an rpoB531 mutation. Ofloxacin susceptibility was restored by exposure to efflux pump inhibitors. Studies in BALB/c mice showed that verapamil in combination with first-line drugs significantly reduced pulmonary CFUs after 1 and 2 months treatment (p < 0.05). CONCLUSION: Exposure of rifampicin resistant M. tuberculosis strains to rifampicin can potentially compromise the efficacy of the second-line treatment regimens containing ofloxacin, thereby emphasising the need for rapid diagnostics to guide treatment. Efflux pump inhibitors have the potential to improve the efficacy of anti-tuberculosis drug treatment.

The Saccharomyces cerevisiae linker histone Hho1p is essential for chromatin compaction in stationary phase and is displaced by transcription.

The importance of core histones in the regulation of DNA function by chromatin is clear. However, little is known about the role of the linker histone. We investigated the role of H1 in Saccharomyces cerevisiae during extensive transcriptional reprogramming in stationary phase. Although the levels of linker histone Hho1p remained constant during growth to semiquiescence, there was a genome-wide increase in binding to chromatin. Hho1p was essential for compaction of chromatin in stationary phase, but not for general transcriptional repression. A clear, genome-wide anticorrelation was seen between the level of bound Hho1p and gene expression. Surprisingly, the rank order of gene activity was maintained even in the absence of Hho1p. Based on these findings, we suggest that linker histone Hho1p has a limited role in transcriptional regulation and that the dynamically exchanging linker histone may be evicted from chromatin by transcriptional activity.

Evidence for a rapid rate of molecular evolution at the hypervariable and immunogenic Mycobacterium tuberculosis PPE38 gene region.

BACKGROUND: PPE38 (Rv2352c) is a member of the large PPE gene family of Mycobacterium tuberculosis and related mycobacteria. The function of PPE proteins is unknown but evidence suggests that many are cell-surface associated and recognised by the host immune system. Previous studies targeting other PPE gene members suggest that some display high levels of polymorphism and it is thought that this might represent a means of providing antigenic variation. We have analysed the genetic variability of the PPE38 genomic region on a cohort of M. tuberculosis clinical isolates representing all of the major phylogenetic lineages, along with the ancestral M. tuberculosis complex (MTBC) member M. canettii, and supplemented this with analysis of publicly available whole genome sequences representing additional M. tuberculosis clinical isolates, other MTBC members and non tuberculous mycobacteria (NTM). Where possible we have extended this analysis to include the adjacent plcABC and PPE39/40 genomic regions. RESULTS: We show that the ancestral MTBC PPE38 region comprises 2 homologous PPE genes (PPE38 and PPE71), separated by 2 esat-6 (esx)-like genes and that this structure derives from an esx/esx/PPE duplication in the common ancestor of M. tuberculosis and M. marinum. We also demonstrate that this region of the genome is hypervariable due to frequent IS6110 integration, IS6110-associated recombination, and homologous recombination and gene conversion events between PPE38 and PPE71. These mutations result in combinations of gene deletion, gene truncation and gene disruption in the majority of clinical isolates. These mutations were generally found to be IS6110 strain lineage-specific, although examples of additional within-lineage and even within-cluster mutations were observed. Furthermore, we provide evidence that the published M. tuberculosis H37Rv whole genome sequence is inaccurate regarding this region. CONCLUSION: Our results show that this antigen-encoding region of the M. tuberculosis genome is hypervariable. The observation that numerous different mutations have become fixed within specific lineages demonstrates that this genomic region is undergoing rapid molecular evolution and that further lineage-specific evolutionary expansion and diversification has occurred subsequent to the lineage-defining mutational events. We predict that functional loss of these genes could aid immune evasion. Finally, we also show that the PPE38 region of the published M. tuberculosis H37Rv whole genome sequence is not representative of the ATCC H37Rv reference strain.

The role of IS6110 in the evolution of Mycobacterium tuberculosis.

Members of the Mycobacterium tuberculosis complex contain the transposable element IS6110 which, due to its high numerical and positional polymorphism, has become a widely used marker in epidemiological studies. Here, we review the evidence that IS6110 is not simply a passive or 'junk' DNA sequence, but that, through its transposable activity, it is able to generate genotypic variation that translates into strain-specific phenotypic variation. We also speculate on the role that this variation has played in the evolution of M. tuberculosis and conclude that the presence of a moderate IS6110 copy number within the genome may provide the pathogen with a selective advantage that has aided its virulence.

Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints.

Mycobacterium tuberculosis complex (MTBC) genomes contain 2 large gene families termed pe and ppe. The function of pe/ppe proteins remains enigmatic but studies suggest that they are secreted or cell surface associated and are involved in bacterial virulence. Previous studies have also shown that some pe/ppe genes are polymorphic, a finding that suggests involvement in antigenic variation. Using comparative sequence analysis of 18 publicly available MTBC whole genome sequences, we have performed alignments of 33 pe (excluding pe_pgrs) and 66 ppe genes in order to detect the frequency and nature of genetic variation. This work has been supplemented by whole gene sequencing of 14 pe/ppe (including 5 pe_pgrs) genes in a cohort of 40 diverse and well defined clinical isolates covering all the main lineages of the M. tuberculosis phylogenetic tree. We show that nsSNP's in pe (excluding pgrs) and ppe genes are 3.0 and 3.3 times higher than in non-pe/ppe genes respectively and that numerous other mutation types are also present at a high frequency. It has previously been shown that non-pe/ppe M. tuberculosis genes display a remarkably low level of purifying selection. Here, we also show that compared to these genes those of the pe/ppe families show a further reduction of selection pressure that suggests neutral evolution. This is inconsistent with the positive selection pressure of "classical" antigenic variation. Finally, by analyzing such a large number of genes we were able to detect large differences in mutation type and frequency between both individual genes and gene sub-families. The high variation rates and absence of selective constraints provides valuable insights into potential pe/ppe function. Since pe/ppe proteins are highly antigenic and have been studied as potential vaccine components these results should also prove informative for aspects of M. tuberculosis vaccine design.

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes.

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.

Multiple, independent, identical IS6110 insertions in Mycobacterium tuberculosis PPE genes.

IS6110 is a transposable element found in Mycobacterium tuberculosis complex members. Regions of preferential IS6110 integration occur within the M. tuberculosis genome but the element has not previously been shown to exhibit any sequence-specific integration preferences. Here we provide evidence for multiple independent IS6110 insertions into identical, or near-identical, positions within the highly homologous 5' region of 3 PPE genes.

Evidence for whole chromosome 6 loss and duplication of the remaining chromosome in acute lymphoblastic leukemia.

HLA class I molecules serve the essential immunological function of presenting antigen to CD8+ T lymphocytes. Tumor cells may present tumor-specific antigen to T cells via these molecules, but many tumors show a loss or down-regulation of HLA class I expression and this may serve as an immune escape mechanism. Using a microsatellite marker-based method, we have searched for loss of heterozygosity (LOH) mutations at 3 genomic regions implicated in HLA class I expression in a cohort of 56 acute lymphoblastic leukemia (ALL) samples. The regions analyzed consisted of the HLA class I heavy chain genes located within the MHC genomic region on chromosome arm 6p, the HLA class I light chain (beta-2-microglobulin, B2M) gene on chromosome arm 15q, and the putative HLA modifier of methylation gene (MEMO1) located on chromosome arm 1q. Results revealed low frequencies of B2M (2/55) and MEMO1 (5/42) LOH but a high frequency of MHC LOH (19/56) that was usually associated with whole chromosome 6 loss (13/19). Cytogenetic data were available for 30 samples, including nine of those that exhibited apparent whole chromosome 6 loss. No cases of chromosome 6 monosomy were observed. We propose that whole chromosome 6 loss with reduplication of the remaining chromosome is common in ALL and that it is driven by the presence of tumor-inhibiting factors on chromosome arm 6p (the HLA loci) along with previously localized tumor-suppressor genes on chromosome arm 6q.