Fanconi anemia-isogenic head and neck cancer cell line pairs: A basic and translational science resource.

Fanconi anemia (FA) is a heritable malformation, bone marrow failure and cancer predisposition syndrome that confers an exceptionally high risk of squamous carcinomas. These carcinomas originate in epithelia lining the mouth, proximal esophagus, vulva and anus: their origins are not understood, and no effective ways have been identified to prevent or delay their appearance. Many FA-associated carcinomas are also therapeutically challenging: they may be multi-focal and stage-advanced at diagnosis, and most individuals with FA cannot tolerate standard-of-care systemic therapies such as DNA cross-linking drugs or ionizing radiation due to constitutional DNA damage hypersensitivity. We developed the Fanconi Anemia Cancer Cell Line Resource (FA-CCLR) to foster new work on the origins, treatment and prevention of FA-associated carcinomas. The FA-CCLR consists of Fanconi-isogenic head and neck squamous cell carcinoma (HNSCC) cell line pairs generated from five individuals with FA-associated HNSCC, and five individuals with sporadic HNSCC. Sporadic, isogenic HNSCC cell line pairs were generated in parallel with FA patient-derived isogenic cell line pairs to provide comparable experimental material to use to identify cell and molecular phenotypes driven by germline or somatic loss of Fanconi pathway function, and the subset of these FA-dependent phenotypes that can be modified, complemented or suppressed. All 10 FANC-isogenic cell line pairs are available to academic, non-profit and industry investigators via the "Fanconi Anemia Research Materials" Resource and Repository at Oregon Health & Sciences University, Portland OR.

Unsupervised discovery of dynamic cell phenotypic states from transmitted light movies.

Identification of cell phenotypic states within heterogeneous populations, along with elucidation of their switching dynamics, is a central challenge in modern biology. Conventional single-cell analysis methods typically provide only indirect, static phenotypic readouts. Transmitted light images, on the other hand, provide direct morphological readouts and can be acquired over time to provide a rich data source for dynamic cell phenotypic state identification. Here, we describe an end-to-end deep learning platform, UPSIDE (Unsupervised Phenotypic State IDEntification), for discovering cell states and their dynamics from transmitted light movies. UPSIDE uses the variational auto-encoder architecture to learn latent cell representations, which are then clustered for state identification, decoded for feature interpretation, and linked across movie frames for transition rate inference. Using UPSIDE, we identified distinct blood cell types in a heterogeneous dataset. We then analyzed movies of patient-derived acute myeloid leukemia cells, from which we identified stem-cell associated morphological states as well as the transition rates to and from these states. UPSIDE opens up the use of transmitted light movies for systematic exploration of cell state heterogeneity and dynamics in biology and medicine.

High-throughput, microscope-based sorting to dissect cellular heterogeneity.

Microscopy is a powerful tool for characterizing complex cellular phenotypes, but linking these phenotypes to genotype or RNA expression at scale remains challenging. Here, we present Visual Cell Sorting, a method that physically separates hundreds of thousands of live cells based on their visual phenotype. Automated imaging and phenotypic analysis directs selective illumination of Dendra2, a photoconvertible fluorescent protein expressed in live cells; these photoactivated cells are then isolated using fluorescence-activated cell sorting. First, we use Visual Cell Sorting to assess hundreds of nuclear localization sequence variants in a pooled format, identifying variants that improve nuclear localization and enabling annotation of nuclear localization sequences in thousands of human proteins. Second, we recover cells that retain normal nuclear morphologies after paclitaxel treatment, and then derive their single-cell transcriptomes to identify pathways associated with paclitaxel resistance in cancers. Unlike alternative methods, Visual Cell Sorting depends on inexpensive reagents and commercially available hardware. As such, it can be readily deployed to uncover the relationships between visual cellular phenotypes and internal states, including genotypes and gene expression programs.

The Werner syndrome RECQ helicase targets G4 DNA in human cells to modulate transcription.

The Werner syndrome (WS) is a prototypic adult Mendelian progeroid syndrome in which signs of premature aging are associated with genomic instability and an elevated risk of cancer. The WRN RECQ helicase protein binds and unwinds G-quadruplex (G4) DNA substrates in vitro, and we identified significant enrichment in G4 sequence motifs at the transcription start site and 5' ends of first introns (false discovery rate < 0.001) of genes down-regulated in WS patient fibroblasts. This finding provides strong evidence that WRN binds G4 DNA structures at many chromosomal sites to modulate gene expression. WRN appears to bind a distinct subpopulation of G4 motifs in human cells, when compared with the related Bloom syndrome RECQ helicase protein. Functional annotation of the genes and miRNAs altered in WS provided new insight into WS disease pathogenesis. WS patient fibroblasts displayed altered expression of multiple, mechanistically distinct, senescence-associated gene expression programs, with altered expression of disease-associated miRNAs, and dysregulation of canonical pathways that regulate cell signaling, genome stability and tumorigenesis. WS fibroblasts also displayed a highly statistically significant and distinct gene expression signature, with coordinate overexpression of nearly all of the cytoplasmic tRNA synthetases and associated ARS-interacting multifunctional protein genes. The 'non-canonical' functions of many of these upregulated tRNA charging proteins may together promote WS disease pathogenesis. Our results identify the human WRN RECQ protein as a G4 helicase that modulates gene expression in G4-dependent fashion at many chromosomal sites and provide several new and unexpected mechanistic insights into WS disease pathogenesis.

Metformin improves defective hematopoiesis and delays tumor formation in Fanconi anemia mice.

Fanconi anemia (FA) is an inherited bone marrow failure disorder associated with a high incidence of leukemia and solid tumors. Bone marrow transplantation is currently the only curative therapy for the hematopoietic complications of this disorder. However, long-term morbidity and mortality remain very high, and new therapeutics are badly needed. Here we show that the widely used diabetes drug metformin improves hematopoiesis and delays tumor formation in Fancd2-/- mice. Metformin is the first compound reported to improve both of these FA phenotypes. Importantly, the beneficial effects are specific to FA mice and are not seen in the wild-type controls. In this preclinical model of FA, metformin outperformed the current standard of care, oxymetholone, by improving peripheral blood counts in Fancd2-/- mice significantly faster. Metformin increased the size of the hematopoietic stem cell compartment and enhanced quiescence in hematopoietic stem and progenitor cells. In tumor-prone Fancd2-/-Trp53+/- mice, metformin delayed the onset of tumors and significantly extended the tumor-free survival time. In addition, we found that metformin and the structurally related compound aminoguanidine reduced DNA damage and ameliorated spontaneous chromosome breakage and radials in human FA patient-derived cells. Our results also indicate that aldehyde detoxification might be one of the mechanisms by which metformin reduces DNA damage in FA cells.

Genome Sequence and Transcriptome Analyses of Chrysochromulina tobin: Metabolic Tools for Enhanced Algal Fitness in the Prominent Order Prymnesiales (Haptophyceae).

Haptophytes are recognized as seminal players in aquatic ecosystem function. These algae are important in global carbon sequestration, form destructive harmful blooms, and given their rich fatty acid content, serve as a highly nutritive food source to a broad range of eco-cohorts. Haptophyte dominance in both fresh and marine waters is supported by the mixotrophic nature of many taxa. Despite their importance the nuclear genome sequence of only one haptophyte, Emiliania huxleyi (Isochrysidales), is available. Here we report the draft genome sequence of Chrysochromulina tobin (Prymnesiales), and transcriptome data collected at seven time points over a 24-hour light/dark cycle. The nuclear genome of C. tobin is small (59 Mb), compact (∼ 40% of the genome is protein coding) and encodes approximately 16,777 genes. Genes important to fatty acid synthesis, modification, and catabolism show distinct patterns of expression when monitored over the circadian photoperiod. The C. tobin genome harbors the first hybrid polyketide synthase/non-ribosomal peptide synthase gene complex reported for an algal species, and encodes potential anti-microbial peptides and proteins involved in multidrug and toxic compound extrusion. A new haptophyte xanthorhodopsin was also identified, together with two "red" RuBisCO activases that are shared across many algal lineages. The Chrysochromulina tobin genome sequence provides new information on the evolutionary history, ecology and economic importance of haptophytes.

Human RECQ Helicase Pathogenic Variants, Population Variation and "Missing" Diseases.

Heritable loss of function mutations in the human RECQ helicase genes BLM, WRN, and RECQL4 cause Bloom, Werner, and Rothmund-Thomson syndromes, cancer predispositions with additional developmental or progeroid features. In order to better understand RECQ pathogenic and population variation, we systematically analyzed genetic variation in all five human RECQ helicase genes. A total of 3,741 unique base pair-level variants were identified, across 17,605 potential mutation sites. Direct counting of BLM, RECQL4, and WRN pathogenic variants was used to determine aggregate and disease-specific carrier frequencies. The use of biochemical and model organism data, together with computational prediction, identified over 300 potentially pathogenic population variants in RECQL and RECQL5, the two RECQ helicases that are not yet linked to a heritable deficiency syndrome. Despite the presence of these predicted pathogenic variants in the human population, we identified no individuals homozygous for any biochemically verified or predicted pathogenic RECQL or RECQL5 variant. Nor did we find any individual heterozygous for known pathogenic variants in two or more of the disease-associated RECQ helicase genes BLM, RECQL4, or WRN. Several postulated RECQ helicase deficiency syndromes-RECQL or RECQL5 loss of function, or compound haploinsufficiency for the disease-associated RECQ helicases-may remain missing, as they likely incompatible with life.

A synthetic homing endonuclease-based gene drive system in the human malaria mosquito.

Genetic methods of manipulating or eradicating disease vector populations have long been discussed as an attractive alternative to existing control measures because of their potential advantages in terms of effectiveness and species specificity. The development of genetically engineered malaria-resistant mosquitoes has shown, as a proof of principle, the possibility of targeting the mosquito's ability to serve as a disease vector. The translation of these achievements into control measures requires an effective technology to spread a genetic modification from laboratory mosquitoes to field populations. We have suggested previously that homing endonuclease genes (HEGs), a class of simple selfish genetic elements, could be exploited for this purpose. Here we demonstrate that a synthetic genetic element, consisting of mosquito regulatory regions and the homing endonuclease gene I-SceI, can substantially increase its transmission to the progeny in transgenic mosquitoes of the human malaria vector Anopheles gambiae. We show that the I-SceI element is able to invade receptive mosquito cage populations rapidly, validating mathematical models for the transmission dynamics of HEGs. Molecular analyses confirm that expression of I-SceI in the male germline induces high rates of site-specific chromosomal cleavage and gene conversion, which results in the gain of the I-SceI gene, and underlies the observed genetic drive. These findings demonstrate a new mechanism by which genetic control measures can be implemented. Our results also show in principle how sequence-specific genetic drive elements like HEGs could be used to take the step from the genetic engineering of individuals to the genetic engineering of populations.

Regulation of gene expression by the BLM helicase correlates with the presence of G-quadruplex DNA motifs.

Bloom syndrome is a rare autosomal recessive disorder characterized by genetic instability and cancer predisposition, and caused by mutations in the gene encoding the Bloom syndrome, RecQ helicase-like (BLM) protein. To determine whether altered gene expression might be responsible for pathological features of Bloom syndrome, we analyzed mRNA and microRNA (miRNA) expression in fibroblasts from individuals with Bloom syndrome and in BLM-depleted control fibroblasts. We identified mRNA and miRNA expression differences in Bloom syndrome patient and BLM-depleted cells. Differentially expressed mRNAs are connected with cell proliferation, survival, and molecular mechanisms of cancer, and differentially expressed miRNAs target genes involved in cancer and in immune function. These and additional altered functions or pathways may contribute to the proportional dwarfism, elevated cancer risk, immune dysfunction, and other features observed in Bloom syndrome individuals. BLM binds to G-quadruplex (G4) DNA, and G4 motifs were enriched at transcription start sites (TSS) and especially within first introns (false discovery rate ≤ 0.001) of differentially expressed mRNAs in Bloom syndrome compared with normal cells, suggesting that G-quadruplex structures formed at these motifs are physiologic targets for BLM. These results identify a network of mRNAs and miRNAs that may drive the pathogenesis of Bloom syndrome.

Global and disease-associated genetic variation in the human Fanconi anemia gene family.

Fanconi anemia (FA) is a human recessive genetic disease resulting from inactivating mutations in any of 16 FANC (Fanconi) genes. Individuals with FA are at high risk of developmental abnormalities, early bone marrow failure and leukemia. These are followed in the second and subsequent decades by a very high risk of carcinomas of the head and neck and anogenital region, and a small continuing risk of leukemia. In order to characterize base pair-level disease-associated (DA) and population genetic variation in FANC genes and the segregation of this variation in the human population, we identified 2948 unique FANC gene variants including 493 FA DA variants across 57,240 potential base pair variation sites in the 16 FANC genes. We then analyzed the segregation of this variation in the 7578 subjects included in the Exome Sequencing Project (ESP) and the 1000 Genomes Project (1KGP). There was a remarkably high frequency of FA DA variants in ESP/1KGP subjects: at least 1 FA DA variant was identified in 78.5% (5950 of 7578) individuals included in these two studies. Six widely used functional prediction algorithms correctly identified only a third of the known, DA FANC missense variants. We also identified FA DA variants that may be good candidates for different types of mutation-specific therapies. Our results demonstrate the power of direct DNA sequencing to detect, estimate the frequency of and follow the segregation of deleterious genetic variation in human populations.

Essential role for Cdk2 inhibitory phosphorylation during replication stress revealed by a human Cdk2 knockin mutation.

Cyclin-dependent kinases (Cdks) coordinate cell division, and their activities are tightly controlled. Phosphorylation of threonine 14 (T14) and tyrosine 15 (Y15) inhibits Cdks and regulates their activities in numerous physiologic contexts. Although the roles of Cdk1 inhibitory phosphorylation during mitosis are well described, studies of Cdk2 inhibitory phosphorylation during S phrase have largely been indirect. To specifically study the functions of Cdk2 inhibitory phosphorylation, we used gene targeting to make an endogenous Cdk2 knockin allele in human cells, termed Cdk2AF, which prevents Cdk2 T14 and Y15 phosphorylation. Cdk2AF caused premature S-phase entry, rapid cyclin E degradation, abnormal DNA replication, and genome instability. Cdk2AF cells also exhibited strikingly abnormal responses to replication stress, accumulated irreparable DNA damage, and permanently exited the cell cycle after transient exposure to S-phase inhibitors. Our results reveal the specific and essential roles of Cdk2 inhibitory phosphorylation in the successful execution of the replication stress checkpoint response and in maintaining genome integrity.

CpG island methylator phenotype is associated with response to adjuvant irinotecan-based therapy for stage III colon cancer.

BACKGROUND & AIMS: The CpG island methylator phenotype (CIMP), defined by a high frequency of aberrantly methylated genes, is a characteristic of a subclass of colon tumors with distinct clinical and molecular features. Cohort studies have produced conflicting results on responses of CIMP-positive tumors to chemotherapy. We assessed the association between tumor CIMP status and survival of patients receiving adjuvant fluorouracil and leucovorin alone or with irinotecan (IFL). METHODS: We analyzed data from patients with stage III colon adenocarcinoma randomly assigned to groups given fluorouracil and leucovorin or IFL after surgery, from April 1999 through April 2001. The primary end point of the trial was overall survival and the secondary end point was disease-free survival. DNA isolated from available tumor samples (n = 615) was used to determine CIMP status based on methylation patterns at the CACNA1G, IGF2, NEUROG1, RUNX3, and SOCS1 loci. The effects of CIMP on survival were modeled using Kaplan-Meier and Cox proportional hazards; interactions with treatment and BRAF, KRAS, and mismatch repair (MMR) status were also investigated. RESULTS: Of the tumor samples characterized for CIMP status, 145 were CIMP positive (23%). Patients with CIMP-positive tumors had shorter overall survival times than patients with CIMP-negative tumors (hazard ratio = 1.36; 95% confidence interval: 1.01-1.84). Treatment with IFL showed a trend toward increased overall survival for patients with CIMP-positive tumors, compared with treatment with fluorouracil and leucovorin (hazard ratio = 0.62; 95% CI: 0.37-1.05; P = .07), but not for patients with CIMP-negative tumors (hazard ratio = 1.38; 95% CI: 1.00-1.89; P = .049). In a 3-way interaction analysis, patients with CIMP-positive, MMR-intact tumors benefited most from the addition of irinotecan to fluorouracil and leucovorin therapy (for the interaction, P = .01). CIMP was more strongly associated with response to IFL than MMR status. Results for disease-free survival times were comparable among all analyses. CONCLUSIONS: Patients with stage III, CIMP-positive, MMR-intact colon tumors have longer survival times when irinotecan is added to combination therapy with fluorouracil and leucovorin.

Conversion of a laboratory-based test for phenylalanine detection to a simple paper-based format and implications for PKU screening in low-resource settings.

Laboratory-based testing does not reach many individuals in lower-resource settings who could benefit from access to appropriate tests for diagnosis and therapy. A critical issue is laboratory-based testing often requires an environment with a high level of resources and supporting infrastructure that is not available in many areas of the world. The current report describes the conversion of a laboratory-based test for phenylalanine detection to a simple paper-based test appropriate for use in low-resource settings. The paper-based test is easy to operate, with all reagents stored dry on the card, is compatible with visible detection for clinically relevant concentrations of phenylalanine, and has a time to result of 10 minutes. Next steps for test development are discussed in the context of the potential for the paper-based Phe test to be used as a newborn PKU screening test in settings that are not well served by existing screening approaches.

The mitochondrial and chloroplast genomes of the haptophyte Chrysochromulina tobin contain unique repeat structures and gene profiles.

BACKGROUND: Haptophytes are widely and abundantly distributed in both marine and freshwater ecosystems. Few genomic analyses of representatives within this taxon have been reported, despite their early evolutionary origins and their prominent role in global carbon fixation. RESULTS: The complete mitochondrial and chloroplast genome sequences of the haptophyte Chrysochromulina tobin (Prymnesiales) provide insight into the architecture and gene content of haptophyte organellar genomes. The mitochondrial genome (~34 kb) encodes 21 protein coding genes and contains a complex, 9 kb tandem repeat region. Similar to other haptophytes and rhodophytes, but not cryptophytes or stramenopiles, the mitochondrial genome has lost the nad7, nad9 and nad11 genes. The ~105 kb chloroplast genome encodes 112 protein coding genes, including ycf39 which has strong structural homology to NADP-binding nitrate transcriptional regulators; a divergent 'CheY-like' two-component response regulator (ycf55) and Tic/Toc (ycf60 and ycf80) membrane transporters. Notably, a zinc finger domain has been identified in the rpl36 ribosomal protein gene of all chloroplasts sequenced to date with the exception of haptophytes and cryptophytes--algae that have gained (via lateral gene transfer) an alternative rpl36 lacking the zinc finger motif. The two C. tobin chloroplast ribosomal RNA operon spacer regions differ in tRNA content. Additionally, each ribosomal operon contains multiple single nucleotide polymorphisms (SNPs)--a pattern observed in rhodophytes and cryptophytes, but few stramenopiles. Analysis of small (<200 bp) chloroplast encoded tandem and inverted repeats in C. tobin and 78 other algal chloroplast genomes show that repeat type, size and location are correlated with gene identity and taxonomic clade. CONCLUSION: The Chrysochromulina tobin organellar genomes provide new insight into organellar function and evolution. These are the first organellar genomes to be determined for the prymnesiales, a taxon that is present in both oceanic and freshwater systems and represents major primary photosynthetic producers and contributors to global ecosystem stability.

Roles of ATM and NBS1 in chromatin structure modulation and DNA double-strand break repair.

We developed a novel system to create DNA double-strand breaks (DSBs) at defined endogenous sites in the human genome, and used this system to detect protein recruitment and loss at and around these breaks by chromatin immunoprecipitation (ChIP). The detection of human ATM protein at site-specific DSBs required functional NBS1 protein, ATM kinase activity and ATM autophosphorylation on Ser 1981. DSB formation led to the localized disruption of nucleosomes, a process that depended on both functional NBS1 and ATM. These two proteins were also required for efficient recruitment of the repair cofactor XRCC4 to DSBs, and for efficient DSB repair. These results demonstrate the functional importance of ATM kinase activity and phosphorylation in the response to DSBs, and support a model in which ordered chromatin structure changes that occur after DNA breakage depend on functional NBS1 and ATM, and facilitate DNA DSB repair.

DNA polymerases and human disease.

The human genome encodes at least 14 DNA-dependent DNA polymerases--a surprisingly large number. These include the more abundant, high-fidelity enzymes that replicate the bulk of genomic DNA, together with eight or more specialized DNA polymerases that have been discovered in the past decade. Although the roles of the newly recognized polymerases are still being defined, one of their crucial functions is to allow synthesis past DNA damage that blocks replication-fork progression. We explore the reasons that might justify the need for so many DNA polymerases, describe their function and mode of regulation, and finally consider links between mutations in DNA polymerases and human disease.

Comprehensive homing endonuclease target site specificity profiling reveals evolutionary constraints and enables genome engineering applications.

Homing endonucleases (HEs) promote the evolutionary persistence of selfish DNA elements by catalyzing element lateral transfer into new host organisms. The high site specificity of this lateral transfer reaction, termed homing, reflects both the length (14-40 bp) and the limited tolerance of target or homing sites for base pair changes. In order to better understand molecular determinants of homing, we systematically determined the binding and cleavage properties of all single base pair variant target sites of the canonical LAGLIDADG homing endonucleases I-CreI and I-MsoI. These Chlorophyta algal HEs have very similar three-dimensional folds and recognize nearly identical 22 bp target sites, but use substantially different sets of DNA-protein contacts to mediate site-specific recognition and cleavage. The site specificity differences between I-CreI and I-MsoI suggest different evolutionary strategies for HE persistence. These differences also provide practical guidance in target site finding, and in the generation of HE variants with high site specificity and cleavage activity, to enable genome engineering applications.

James German and the Quest to Understand Human RECQ Helicase Deficiencies.

James German's work to establish the natural history and cancer risk associated with Bloom syndrome (BS) has had a strong influence on the generation of scientists and clinicians working to understand other RECQ deficiencies and heritable cancer predisposition syndromes. I summarize work by us and others below, inspired by James German's precedents with BS, to understand and compare BS with the other heritable RECQ deficiency syndromes with a focus on Werner syndrome (WS). What we know, unanswered questions and new opportunities are discussed, as are potential ways to treat or modify WS-associated disease mechanisms and pathways.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Mutation Patterns Predict Drug Sensitivity in Acute Myeloid Leukemia.

PURPOSE: The inherent genetic heterogeneity of acute myeloid leukemia (AML) has challenged the development of precise and effective therapies. The objective of this study was to elucidate the genomic basis of drug resistance or sensitivity, identify signatures for drug response prediction, and provide resources to the research community. EXPERIMENTAL DESIGN: We performed targeted sequencing, high-throughput drug screening, and single-cell genomic profiling on leukemia cell samples derived from patients with AML. Statistical approaches and machine learning models were applied to identify signatures for drug response prediction. We also integrated large public datasets to understand the co-occurring mutation patterns and further investigated the mutation profiles in the single cells. The features revealed in the co-occurring or mutual exclusivity pattern were further subjected to machine learning models. RESULTS: We detected genetic signatures associated with sensitivity or resistance to specific agents, and identified five co-occurring mutation groups. The application of single-cell genomic sequencing unveiled the co-occurrence of variants at the individual cell level, highlighting the presence of distinct subclones within patients with AML. Using the mutation pattern for drug response prediction demonstrates high accuracy in predicting sensitivity to some drug classes, such as MEK inhibitors for RAS-mutated leukemia. CONCLUSIONS: Our study highlights the importance of considering the gene mutation patterns for the prediction of drug response in AML. It provides a framework for categorizing patients with AML by mutations that enable drug sensitivity prediction.