Patrick Flaherty

Functional profiling of the Saccharomyces cerevisiae genome

Determining the effect of gene deletion is a fundamental approach to understanding gene function. Conventional genetic screens exhibit biases, and genes contributing to a phenotype are often missed. We systematically constructed a nearly complete collection of gene-deletion mutants (96% of annotated open reading frames, or ORFs) of the yeast Saccharomyces cerevisiae. DNA sequences dubbed ‘molecular bar codes’ uniquely identify each strain, enabling their growth to be analysed in parallel and the fitness contribution of each gene to be quantitatively assessed by hybridization to high-density oligonucleotide arrays. We show that previously known and new genes are necessary for optimal growth under six well-studied conditions: high salt, sorbitol, galactose, pH 8, minimal medium and nystatin treatment. Less than 7% of genes that exhibit a significant increase in messenger RNA expression are also required for optimal growth in four of the tested conditions. Our results validate the yeast gene-deletion collection as a valuable resource for functional genomics.

Oncogenic BRAF mutation with CDKN2A inactivation is characteristic of a subset of pediatric malignant astrocytomas.

Malignant astrocytomas are a deadly solid tumor in children. Limited understanding of their underlying genetic basis has contributed to modest progress in developing more effective therapies. In an effort to identify such alterations, we performed a genome-wide search for DNA copy number aberrations (CNA) in a panel of 33 tumors encompassing grade 1 through grade 4 tumors. Genomic amplifications of 10-fold or greater were restricted to grade 3 and 4 astrocytomas and included the MDM4 (1q32), PDGFRA (4q12), MET (7q21), CMYC (8q24), PVT1 (8q24), WNT5B (12p13), and IGF1R (15q26) genes. Homozygous deletions of CDKN2A (9p21), PTEN (10q26), and TP53 (17p3.1) were evident among grade 2 to 4 tumors. BRAF gene rearrangements that were indicated in three tumors prompted the discovery of KIAA1549-BRAF fusion transcripts expressed in 10 of 10 grade 1 astrocytomas and in none of the grade 2 to 4 tumors. In contrast, an oncogenic missense BRAF mutation (BRAF(V600E)) was detected in 7 of 31 grade 2 to 4 tumors but in none of the grade 1 tumors. BRAF(V600E) mutation seems to define a subset of malignant astrocytomas in children, in which there is frequent concomitant homozygous deletion of CDKN2A (five of seven cases). Taken together, these findings highlight BRAF as a frequent mutation target in pediatric astrocytomas, with distinct types of BRAF alteration occurring in grade 1 versus grade 2 to 4 tumors.

Genome-Wide Requirements for Resistance to Functionally Distinct DNA-Damaging Agents

The mechanistic and therapeutic differences in the cellular response to DNA-damaging compounds are not completely understood, despite intense study. To expand our knowledge of DNA damage, we assayed the effects of 12 closely related DNA-damaging agents on the complete pool of ~4,700 barcoded homozygous deletion strains of Saccharomyces cerevisiae. In our protocol, deletion strains are pooled together and grown competitively in the presence of compound. Relative strain sensitivity is determined by hybridization of PCR-amplified barcodes to an oligonucleotide array carrying the barcode complements. These screens identified genes in well-characterized DNA-damage-response pathways as well as genes whose role in the DNA-damage response had not been previously established. High-throughput individual growth analysis was used to independently confirm microarray results. Each compound produced a unique genome-wide profile. Analysis of these data allowed us to determine the relative importance of DNA-repair modules for resistance to each of the 12 profiled compounds. Clustering the data for 12 distinct compounds uncovered both known and novel functional interactions that comprise the DNA-damage response and allowed us to define the genetic determinants required for repair of interstrand cross-links. Further genetic analysis allowed determination of epistasis for one of these functional groups.

Ultrasensitive detection of rare mutations using next-generation targeted resequencing

With next-generation DNA sequencing technologies, one can interrogate a specific genomic region of interest at very high depth of coverage and identify less prevalent, rare mutations in heterogeneous clinical samples. However, the mutation detection levels are limited by the error rate of the sequencing technology as well as by the availability of variant-calling algorithms with high statistical power and low false positive rates. We demonstrate that we can robustly detect mutations at 0.1% fractional representation. This represents accurate detection of one mutant per every 1000 wild-type alleles. To achieve this sensitive level of mutation detection, we integrate a high accuracy indexing strategy and reference replication for estimating sequencing error variance. We employ a statistical model to estimate the error rate at each position of the reference and to quantify the fraction of variant base in the sample. Our method is highly specific (99%) and sensitive (100%) when applied to a known 0.1% sample fraction admixture of two synthetic DNA samples to validate our method. As a clinical application of this method, we analyzed nine clinical samples of H1N1 influenza A and detected an oseltamivir (antiviral therapy) resistance mutation in the H1N1 neuraminidase gene at a sample fraction of 0.18%.

Phase II Study of Gemcitabine, Carboplatin, and Iniparib As Neoadjuvant Therapy for Triple-Negative and BRCA1/2 Mutation–Associated Breast Cancer With Assessment of a Tumor-Based Measure of Genomic Instability: PrECOG 0105

PURPOSE This study was designed to assess efficacy, safety, and predictors of response to iniparib in combination with gemcitabine and carboplatin in early-stage triple-negative and BRCA1/2 mutation-associated breast cancer. PATIENTS AND METHODS This single-arm phase II study enrolled patients with stage I to IIIA (T ≥ 1 cm) estrogen receptor-negative (≤ 5%), progesterone receptor-negative (≤ 5%), and human epidermal growth factor receptor 2-negative or BRCA1/2 mutation-associated breast cancer. Neoadjuvant gemcitabine (1,000 mg/m(2) intravenously [IV] on days 1 and 8), carboplatin (area under curve of 2 IV on days 1 and 8), and iniparib (5.6 mg/kg IV on days 1, 4, 8, and 11) were administered every 21 days for four cycles, until the protocol was amended to six cycles. The primary end point was pathologic complete response (no invasive carcinoma in breast or axilla). All patients underwent comprehensive BRCA1/2 genotyping, and homologous recombination deficiency was assessed by loss of heterozygosity (HRD-LOH) in pretreatment core breast biopsies. RESULTS Among 80 patients, median age was 48 years; 19 patients (24%) had germline BRCA1 or BRCA2 mutations; clinical stage was I (13%), IIA (36%), IIB (36%), and IIIA (15%). Overall pathologic complete response rate in the intent-to-treat population (n = 80) was 36% (90% CI, 27 to 46). Mean HRD-LOH scores were higher in responders compared with nonresponders (P = .02) and remained significant when BRCA1/2 germline mutations carriers were excluded (P = .021). CONCLUSION Preoperative combination of gemcitabine, carboplatin, and iniparib is active in the treatment of early-stage triple-negative and BRCA1/2 mutation-associated breast cancer. The HRD-LOH assay was able to identify patients with sporadic triple-negative breast cancer lacking a BRCA1/2 mutation, but with an elevated HRD-LOH score, who achieved a favorable pathologic response. Confirmatory controlled trials are warranted.

Metastatic tumor evolution and organoid modeling implicate TGFBR2 as a cancer driver in diffuse gastric cancer

BackgroundGastric cancer is the second-leading cause of global cancer deaths, with metastatic disease representing the primary cause of mortality. To identify candidate drivers involved in oncogenesis and tumor evolution, we conduct an extensive genome sequencing analysis of metastatic progression in a diffuse gastric cancer. This involves a comparison between a primary tumor from a hereditary diffuse gastric cancer syndrome proband and its recurrence as an ovarian metastasis.ResultsBoth the primary tumor and ovarian metastasis have common biallelic loss-of-function of both the CDH1 and TP53 tumor suppressors, indicating a common genetic origin. While the primary tumor exhibits amplification of the Fibroblast growth factor receptor 2 (FGFR2) gene, the metastasis notably lacks FGFR2 amplification but rather possesses unique biallelic alterations of Transforming growth factor-beta receptor 2 (TGFBR2), indicating the divergent in vivo evolution of a TGFBR2-mutant metastatic clonal population in this patient. As TGFBR2 mutations have not previously been functionally validated in gastric cancer, we modeled the metastatic potential of TGFBR2 loss in a murine three-dimensional primary gastric organoid culture. The Tgfbr2 shRNA knockdown within Cdh1-/-; Tp53-/- organoids generates invasion in vitro and robust metastatic tumorigenicity in vivo, confirming Tgfbr2 metastasis suppressor activity.ConclusionsWe document the metastatic differentiation and genetic heterogeneity of diffuse gastric cancer and reveal the potential metastatic role of TGFBR2 loss-of-function. In support of this study, we apply a murine primary organoid culture method capable of recapitulating in vivo metastatic gastric cancer. Overall, we describe an integrated approach to identify and functionally validate putative cancer drivers involved in metastasis.

Identification of a New Class of Antifungals Targeting the Synthesis of Fungal Sphingolipids

ABSTRACT Recent estimates suggest that >300 million people are afflicted by serious fungal infections worldwide. Current antifungal drugs are static and toxic and/or have a narrow spectrum of activity. Thus, there is an urgent need for the development of new antifungal drugs. The fungal sphingolipid glucosylceramide (GlcCer) is critical in promoting virulence of a variety of human-pathogenic fungi. In this study, we screened a synthetic drug library for compounds that target the synthesis of fungal, but not mammalian, GlcCer and found two compounds [N′-(3-bromo-4-hydroxybenzylidene)-2-methylbenzohydrazide (BHBM) and its derivative, 3-bromo-N′-(3-bromo-4-hydroxybenzylidene) benzohydrazide (D0)] that were highly effective in vitro and in vivo against several pathogenic fungi. BHBM and D0 were well tolerated in animals and are highly synergistic or additive to current antifungals. BHBM and D0 significantly affected fungal cell morphology and resulted in the accumulation of intracellular vesicles. Deep-sequencing analysis of drug-resistant mutants revealed that four protein products, encoded by genes APL5, COS111, MKK1, and STE2, which are involved in vesicular transport and cell cycle progression, are targeted by BHBM. IMPORTANCE Fungal infections are a significant cause of morbidity and mortality worldwide. Current antifungal drugs suffer from various drawbacks, including toxicity, drug resistance, and narrow spectrum of activity. In this study, we have demonstrated that pharmaceutical inhibition of fungal glucosylceramide presents a new opportunity to treat cryptococcosis and various other fungal infections. In addition to being effective against pathogenic fungi, the compounds discovered in this study were well tolerated by animals and additive to current antifungals. These findings suggest that these drugs might pave the way for the development of a new class of antifungals. Fungal infections are a significant cause of morbidity and mortality worldwide. Current antifungal drugs suffer from various drawbacks, including toxicity, drug resistance, and narrow spectrum of activity. In this study, we have demonstrated that pharmaceutical inhibition of fungal glucosylceramide presents a new opportunity to treat cryptococcosis and various other fungal infections. In addition to being effective against pathogenic fungi, the compounds discovered in this study were well tolerated by animals and additive to current antifungals. These findings suggest that these drugs might pave the way for the development of a new class of antifungals.

A dual receptor crosstalk model of G-protein-coupled signal transduction.

Macrophage cells that are stimulated by two different ligands that bind to G-protein-coupled receptors (GPCRs) usually respond as if the stimulus effects are additive, but for a minority of ligand combinations the response is synergistic. The G-protein-coupled receptor system integrates signaling cues from the environment to actuate cell morphology, gene expression, ion homeostasis, and other physiological states. We analyze the effects of the two signaling molecules complement factors 5a (C5a) and uridine diphosphate (UDP) on the intracellular second messenger calcium to elucidate the principles that govern the processing of multiple signals by GPCRs. We have developed a formal hypothesis, in the form of a kinetic model, for the mechanism of action of this GPCR signal transduction system using data obtained from RAW264.7 macrophage cells. Bayesian statistical methods are employed to represent uncertainty in both data and model parameters and formally tie the model to experimental data. When the model is also used as a tool in the design of experiments, it predicts a synergistic region in the calcium peak height dose response that results when cells are simultaneously stimulated by C5a and UDP. An analysis of the model reveals a potential mechanism for crosstalk between the Gαi-coupled C5a receptor and the Gαq-coupled UDP receptor signaling systems that results in synergistic calcium release.

Systematic genomic identification of colorectal cancer genes delineating advanced from early clinical stage and metastasis

BackgroundColorectal cancer is the third leading cause of cancer deaths in the United States. The initial assessment of colorectal cancer involves clinical staging that takes into account the extent of primary tumor invasion, determining the number of lymph nodes with metastatic cancer and the identification of metastatic sites in other organs. Advanced clinical stage indicates metastatic cancer, either in regional lymph nodes or in distant organs. While the genomic and genetic basis of colorectal cancer has been elucidated to some degree, less is known about the identity of specific cancer genes that are associated with advanced clinical stage and metastasis.MethodsWe compiled multiple genomic data types (mutations, copy number alterations, gene expression and methylation status) as well as clinical meta-data from The Cancer Genome Atlas (TCGA). We used an elastic-net regularized regression method on the combined genomic data to identify genetic aberrations and their associated cancer genes that are indicators of clinical stage. We ranked candidate genes by their regression coefficient and level of support from multiple assay modalities.ResultsA fit of the elastic-net regularized regression to 197 samples and integrated analysis of four genomic platforms identified the set of top gene predictors of advanced clinical stage, including: WRN, SYK, DDX5 and ADRA2C. These genetic features were identified robustly in bootstrap resampling analysis.ConclusionsWe conducted an analysis integrating multiple genomic features including mutations, copy number alterations, gene expression and methylation. This integrated approach in which one considers all of these genomic features performs better than any individual genomic assay. We identified multiple genes that robustly delineate advanced clinical stage, suggesting their possible role in colorectal cancer metastatic progression.

Prediction of multiple infections after severe burn trauma: a prospective cohort study.

OBJECTIVE To develop predictive models for early triage of burn patients based on hypersusceptibility to repeated infections. BACKGROUND Infection remains a major cause of mortality and morbidity after severe trauma, demanding new strategies to combat infections. Models for infection prediction are lacking. METHODS Secondary analysis of 459 burn patients (≥16 years old) with 20% or more total body surface area burns recruited from 6 US burn centers. We compared blood transcriptomes with a 180-hour cutoff on the injury-to-transcriptome interval of 47 patients (≤1 infection episode) to those of 66 hypersusceptible patients [multiple (≥2) infection episodes (MIE)]. We used LASSO regression to select biomarkers and multivariate logistic regression to built models, accuracy of which were assessed by area under receiver operating characteristic curve (AUROC) and cross-validation. RESULTS Three predictive models were developed using covariates of (1) clinical characteristics; (2) expression profiles of 14 genomic probes; (3) combining (1) and (2). The genomic and clinical models were highly predictive of MIE status [AUROCGenomic = 0.946 (95% CI: 0.906-0.986); AUROCClinical = 0.864 (CI: 0.794-0.933); AUROCGenomic/AUROCClinical P = 0.044]. Combined model has an increased AUROCCombined of 0.967 (CI: 0.940-0.993) compared with the individual models (AUROCCombined/AUROCClinical P = 0.0069). Hypersusceptible patients show early alterations in immune-related signaling pathways, epigenetic modulation, and chromatin remodeling. CONCLUSIONS Early triage of burn patients more susceptible to infections can be made using clinical characteristics and/or genomic signatures. Genomic signature suggests new insights into the pathophysiology of hypersusceptibility to infection may lead to novel potential therapeutic or prophylactic targets.