Combinatorial Mutations Research Articles

Quivers for silting mutation

We give a combinatorial mutation rule for Aihara and Iyama's silting mutation.As an application, we reprove Keller–Yang's mutation rule for Ginzburg algebras, and obtain an analog of that rule for Ginzburg algebras of arbitrary Calabi–Yau dimension.

Highly sensitive detection of the PIK3CA (H1047R) mutation in colorectal cancer using a novel PCR-RFLP method.

BackgroundThe PIK3CAH1047R mutation is considered to be a potential predictive biomarker for EGFR-targeted therapies. In this study, we developed a novel PCR-PFLP approach to detect the PIK3CAH1047R mutation in high effectiveness.MethodsA 126-bp fragment of PIK3CA exon-20 was amplified by PCR, digested with FspI restriction endonuclease and separated by 3 % agarose gel electrophoresis for the PCR-RFLP analysis. The mutant sequence of the PIK3CAH1047R was spiked into the corresponding wild-type sequence in decreasing ratios for sensitivity analysis. Eight-six cases of formalin-fixed paraffin-embedded colorectal cancer (CRC) specimens were subjected to PCR-RFLP to evaluate the applicability of the method.ResultsThe PCR-RFLP method had a capability to detect as litter as 0.4 % of mutation, and revealed 16.3 % of the PIK3CAH1047R mutation in 86 CRC tissues, which was significantly higher than that discovered by DNA sequencing (9.3 %). A positive association between the PIK3CAH1047R mutation and the patients’ age was first found, except for the negative relationship with the degree of tumor differentiation. In addition, the highly sensitive detection of a combinatorial mutation of PIK3CA, KRAS and BRAF was achieved using individual PCR-RFLP methods.ConclusionsWe developed a sensitive, simple and rapid approach to detect the low-abundance PIK3CAH1047R mutation in real CRC specimens, providing an effective tool for guiding cancer targeted therapy.

Hierarchy within the mammary STAT5-driven Wap super-enhancer.

Super-enhancers comprise of dense transcription factor platforms highly enriched for active chromatin marks. A paucity of functional data led us to investigate their role in the mammary gland, an organ characterized by exceptional gene regulatory dynamics during pregnancy. ChIP-Seq for the master regulator STAT5, the glucocorticoid receptor, H3K27ac and MED1, identified 440 mammary-specific super-enhancers, half of which were associated with genes activated during pregnancy. We interrogated the Wap super-enhancer, generating mice carrying mutations in STAT5 binding sites within its three constituent enhancers. Individually, only the most distal site displayed significant enhancer activity. However, combinatorial mutations showed that the 1,000-fold gene induction relied on all enhancers. Disabling the binding sites of STAT5, NFIB and ELF5 in the proximal enhancer incapacitated the entire super-enhancer, suggesting an enhancer hierarchy. The identification of mammary-specific super-enhancers and the mechanistic exploration of the Wap locus provide insight into the complexity of cell-specific and hormone-regulated genes.

Exploring phenotype patterns of breast cancer within somatic mutations: a modicum in the intrinsic code.

Triple-negative (TN) breast cancer (BC) patients have limited treatment options and poor prognosis even after extant treatments and standard chemotherapeutic regimens. Linking TN patients to clinically known phenotypes with appropriate treatments is vital. Location-specific sequence variants are expected to be useful for this purpose by identifying subgroups within a disease population. Single gene mutational signatures have been widely reported, with related phenotypes in literature. We thoroughly survey currently available mutations (and mutated genes), linked to BC phenotypes, to demonstrate their limited performance as sole predictors/biomarkers to assign phenotypes to patients. We then explore mutational combinations, as a pilot study, using The Cancer Genome Atlas Research Network mutational data of BC and three machine learning methods: association rules (limitless arity multiple procedure), decision tree and hierarchical disjoint clustering. The study results in a patient classification scheme through combinatorial mutations in Phosphatidylinositol-4,5-Bisphosphate 3-Kinase and tumor protein 53, being consistent with all three methods, implying its validity from a diverse viewpoint. However, it would warrant further research to select multi-gene signatures to identify phenotypes specifically and be clinically used routinely.

Thermostabilization of the Human Endothelin Type B Receptor

The peptide hormone endothelin, produced by the vascular endothelium, is involved in several physiological functions, including maintenance of vascular tone and humoral homeostasis. Endothelin transmits signals through the endothelin receptor, a G-protein-coupled receptor. Structural studies of the endothelin type B receptor (ETBR) have been unsuccessful due to its structural flexibility and instability in detergent-solubilized solution. To overcome these problems, we explored thermostabilization of human ETBR by establishing an ETBR expression system in Escherichia coli, followed by systematic alanine scanning mutagenesis. Among 297 point mutations, 11 thermostabilizing residues were selected and further mutated to other amino acids. The thermostability indices of these residues, represented by the ratios of endothelin-1 (ET-1) binding activities with or without heat treatment at 27°C for 30min in a ligand-free form, were compared. The ligand affinity and apparent melting temperature (Tm) of the five most thermostable mutants, R124Y, D154A, K270A, S342A, and I381A, were then examined. The apparent Tm of three single mutants, R124Y, D154A, and K270A, was approximately 7°C higher than that of the wild type. The apparent Tm value of a combination of the five residues, named the Y5 ETBR mutant, was 17°C higher than that of the wild type. The Y5 ETBR mutant exhibited an affinity for ET-1 and activated Gq similar to the wild type. Further investigation of the pharmacological properties affected by combinatorial mutations of ET-1, ET-3, TxET-1, and K8794 suggested that Y5 ETBR is highly suitable for representing a ligand-free form of ETBR and is potentially applicable for studying an ET-1-bound form.

Combinatorial Evolution of Enzymes and Synthetic Pathways Using One-Step PCR.

DNA engineering is the fundamental motive driving the rapid development of modern biotechnology. Here, we present a versatile evolution method termed "rapidly efficient combinatorial oligonucleotides for directed evolution" (RECODE) for rapidly introducing multiple combinatorial mutations to the target DNA by combined action of a thermostable high-fidelity DNA polymerase and a thermostable DNA Ligase in one reaction system. By applying this method, we rapidly constructed a variant library of the rpoS promoters (with activity of 8-460%), generated a novel heparinase from the highly specific leech hyaluronidase (with more than 30 mutant residues) and optimized the heme biosynthetic pathway by combinatorial evolution of regulatory elements and pathway enzymes (2500 ± 120 mg L(-1) with 20-fold increase). The simple RECODE method enabled researchers the unparalleled ability to efficiently create diverse mutant libraries for rapid evolution and optimization of enzymes and synthetic pathways.

Functional Hallmarks of a Catalytic DNA that Makes Lariat RNA.

Catalytic DNAs, also known as deoxyribozymes, are of practical value for the synthesis of structurally or topologically complex RNAs, but little is known about the molecular details of DNA catalysis. We have investigated a deoxyribozyme that catalyzes the formation of a specific intramolecular 2',5'-phosphodiester bond to produce lariat RNA, which is an important biological intermediate in eukaryotic mRNA splicing. The results of combinatorial mutation interference analysis (CoMA) allowed us to shrink the catalytic core to 70 % of its original length and revealed that the essential part of the deoxyribozyme sequence contained more than 50 % guanosines. Nucleotide analogue interference mapping (dNAIM) and dimethyl sulfate interference (DMSi) experiments provided atomic details of individual guanosine functional groups. Additional spectroscopic experiments and structural probing data identified conformational changes upon metal-ion binding and catalysis. Overall, this comprehensive analysis of the DNA-catalyzed reaction has provided specific insights into the synthesis of 2',5'-branched RNA, and suggested the general features of deoxyribozymes that catalyze nucleic acid ligation reactions.

Assessment of the In Vivo Activity of PI3K and MEK Inhibitors in Genetically Defined Models of Colorectal Cancer.

The objective of tailoring medicines for cancer patients according to the molecular profile of their disease holds great promise for the improvement of cancer therapy. Nevertheless, this approach has been limited, in part, due to the lack of predictive and informative preclinical studies. Herein, we describe an assessment of the therapeutic potential of targeting PI3K/mTOR and MAPK signaling in genetically defined mouse models of colorectal cancer mirroring disease subtypes targeted for novel therapy in the FOCUS4 trial. Our studies demonstrate that dual PI3K/mTOR inhibition is highly effective in invasive adenocarcinoma models characterized by combinatorial mutations in Apc and Pten; Apc and Kras; and Apc, Pten and Kras. MEK inhibition was effective in the combinatorial Apc and Kras setting, but had no impact in either Apc Pten mutants or in Apc Pten Kras triple mutants. Furthermore, we describe the importance of scheduling for combination studies and show that although no additional benefit is gained in Apc Pten mice, combination of PI3K/mTOR and MAPK inhibition leads to an additive benefit in survival in Apc Kras mice and a synergistic increase in survival in Apc Pten Kras mice. This is the first study using robust colorectal cancer genetically engineered mouse models to support the validity of PI3K/mTOR and MEK inhibitors as tailored therapies for colorectal cancer and highlight the potential importance of drug scheduling in the clinic.

Thermostability enhancement of cellobiose 2-epimerase from Caldicellulosiruptor saccharolyticus by site-directed mutagenesis

Cellobiose 2-epimerase from the thermophile Caldicellulosiruptor saccharolyticus (CsCE) catalyzes the isomerization of lactose into lactulose, a non-digestible disaccharide widely used in food and pharmaceutical industries. Semi-rational approaches were applied to enhance the thermostability of CsCE. A total of eight single-site mutants were designed, and five of them showed prolonged half-life of inactivation at 80°C. Combinatorial mutations were subsequently introduced, and the superior mutant was double mutant E161D/N365P. The half-life was approximately 4-fold higher than that of the wild type enzyme. In addition, the reaction temperature for maximum activity increased from 80°C to 87.5°C, and catalytic efficiency (kcat/Km) for lactulose production was increased 29%. Moreover, this mutant E161D/N365P was more stable against chemical denaturation and showed also a broader pH profile. The second most stable variant were mutant E161D/S180P/S351G with a 3.3-fold increase in half-life. These results provided new insights into the thermostability of CsCE and suggested further potential industrial applications.

Parkin-mediated responses against infection and wound involve TSPO-VDAC complex in Drosophila

Parkin, an E3 ubuquitin ligase associated with Parkinson's disease (PD), has recently been implicated in mediating innate immunity. However, molecular details regarding parkin-mediated immune response remain to be elucidated. Here, we identified mitochondrial TSPO-VDAC complex to genetically interact with parkin in mediating responses against infection and wound in Drosophila. The loss-of-function mutation in parkin results in defective immune response against bacterial infection. Additionally, parkin mutant larvae showed hypersensitivity against wound regardless of bacterial infection. Interestingly, the combinatorial trans-heterozygotic mutations in parkin and TSPO, or parkin and VDAC showed similar lethal tendency with parkin homozygous mutants. Furthermore, knockdown of TSPO alone also resulted in defective responses to infection and wound analogously to parkin mutants. Taken together, we propose that parkin cooperates with TSPO-VDAC complex to mediate responses against infection and wound.

Multiplexed tracking of combinatorial genomic mutations in engineered cell populations.

Multiplexed genome engineering approaches can be used to generate targeted genetic diversity in cell populations on laboratory timescales, but methods to track mutations and link them to phenotypes have been lacking. We present an approach for tracking combinatorial engineered libraries (TRACE) through the simultaneous mapping of millions of combinatorially engineered genomes at single-cell resolution. Distal genomic sites are assembled into individual DNA constructs that are compatible with next-generation sequencing strategies. We used TRACE to map growth selection dynamics for Escherichia coli combinatorial libraries created by recursive multiplex recombineering at a depth 10(4)-fold greater than before. TRACE was used to identify genotype-to-phenotype correlations and to map the evolutionary trajectory of two individual combinatorial mutants in E. coli. Combinatorial mutations in the human ES2 ovarian carcinoma cell line were also assessed with TRACE. TRACE completes the combinatorial engineering cycle and enables more sophisticated approaches to genome engineering in both bacteria and eukaryotic cells than are currently possible.

Protein engineering in the 21st century.

Protein engineering, the process by which novel proteins with desired properties are developed, has grown by leaps and bounds since the first examples of protein mutagenesis were described over three decades ago. In the last 30 years, protein engineers have successfully created a wide range of proteins tailored to specific industrial, medical, and research applications. Yet challenges remain that prevent the engineering of complex protein functions on demand. With the aim of creating a forum for the protein engineering community, the inaugural Protein Engineering Canada Conference was held on June 20–22, 2014, in Ottawa, Canada. Over the course of the 2-day conference, 115 protein scientists from over 30 institutions representing 5 countries shared ideas, networked with colleagues, and learned about exciting new research covering a broad array of topics such as enzyme engineering, computational protein design, X-ray crystallography, protein NMR, molecular modeling, and protein evolution. This special issue highlights the breadth of research that was covered during the Protein Engineering Canada Conference by including previously unpublished research articles submitted by speakers as well as some selections that were submitted directly to Protein Science.

Inter- and intra-tumor profiling of multi-regional colon cancer and metastasis

Intra- and inter-tumor heterogeneity may hinder personalized molecular-target treatment that depends on the somatic mutation profiles. We performed mutation profiling of formalin-fixed paraffin embedded tumors of multi-regional colon cancer and characterized the consequences of intra- and inter-tumor heterogeneity and metastasis using targeted re-sequencing.We performed targeted re-sequencing on multiple spatially separated samples obtained from multi-regional primary colon carcinoma and associated metastatic sites in two patients using next-generation sequencing. In Patient 1 with four primary tumors (P1-1, P1-2, P1-3, and P1-4) and one liver metastasis (H1), mutually exclusive pattern of mutations was observed in four primary tumors. Mutations in primary tumors were identified in three regions; KARS (G13D) and APC (R876*) in P1-2, TP53 (A161S) in P1-3, and KRAS (G12D), PIK3CA (Q546R), and ERBB4 (T272A) in P1-4. Similar combinatorial mutations were observed between P1-4 and H1. The ERBB4 (T272A) mutation observed in P1-4, however, disappeared in H1. In Patient 2 with two primary tumors (P2-1 and P2-2) and one liver metastasis (H2), mutually exclusive pattern of mutations were observed in two primary tumors. We identified mutations; KRAS (G12V), SMAD4 (N129K, R445*, and G508D), TP53 (R175H), and FGFR3 (R805W) in P2-1, and NRAS (Q61K) and FBXW7 (R425C) in P2-2. Similar combinatorial mutations were observed between P2-1 and H2. The SMAD4 (N129K and G508D) mutations observed in P2-1, however, were nor detected in H2. These results suggested that different clones existed in primary tumors and metastatic tumor in Patient 1 and 2 likely originated from P1-4 and P2-1, respectively.In conclusion, we detected the muti-clonalities between intra- and inter-tumors based on mutational profiling in multi-regional colon cancer using next-generation sequencing. Primary region from which metastasis originated could be speculated by mutation profile. Characterization of inter- and inter-tumor heterogeneity can lead to underestimation of the tumor genomics landscape and treatment strategy of personal medicine.

Prognostic significance of novel katG mutations in Mycobacterium tuberculosis

BackgroundBy using whole genome sequencing (WGS), researchers are beginning to understand the genetic diversity of Mycobacterium tuberculosis (MTB) and its consequences for the diagnosis of multidrug-resistant tuberculosis (MDR–TB) on a genomic scale. The Global Consortium for Drug-resistant TB Diagnostics (GCDD) conducted a genome scale variant analyses of 366 clinical MTB genomes (mostly MDR/XDR [extensively drug resistant]) from four countries in order to inform the development of rapid molecular diagnostics. This project has been extended by performing an evolutionary analysis of isoniazid (INH)-resistant isolates for prognostic purposes. Methods151 (130 INHR, 21 INHS) clinical MTB isolates from India (19: 17 INHR, 2 INHS), Moldova (48: 42 INHR, 6 INHS), the Philippines (26: 20 INHR, 6 INHS), and South Africa (58: 51 INHR, 7 INHS) were included in this study. INH drug susceptibility was determined by using MGIT 960 and WHO (World Health Organization)-recommended critical concentration of 0.1mg/L. Isolates were sequenced using PacBio RS WGS platform. A genome-wide variant analysis was conducted using a proprietary pipeline (PacDAP) developed at San Diego State University. To infer the amino acid changes in katG that confer resistance, PAML was utilized to detect sites in silico that are under positive selection. The dN/dS method was used in combination with Bayes empirical Bayes to determine sites under positive selection and Chi-Squared analysis to determine the significance of the selected sites. ResultsPacDAP variant analysis revealed 22 novel catalase-peroxidase (katG product) mutations. Of these, 14 were single nucleotide polymorphisms, while 8 novel mutations appeared in combination with katG S315T and/or with inhA promoter C-15T. These SNPs have not been previously reported. Additionally, 11 previously observed, but uncommon, katG mutations were also observed in these clinical isolates. These results suggest that 17 amino acids in the enzyme are under positive selective pressure; most significantly in South Africa and the Philippines. No selective pressure on codons other than 315 was observed in isolates from Moldova. Due to the low number of isolates from India, the significance of the sites under positive selection was low and no prediction for India could be made based on this study. ConclusionsEleven of the 14 SNPs are resistance conferring, and it is believed that the remaining 8 combinatorial mutations are either compensatory in nature or, in combination with known SNPs, could increase resistance levels. Positive selection results indicate a diversifying evolutionary path to resistance more in line with long tail statistics and therefore indicate a departure from the traditional point mutation (or “hotspot”) model that current molecular diagnostics are based on. Positive selection pressures indicate a future with elevated diagnostic and prognostic significance of the “long tail” (i.e., alternative mechanisms of resistance) and potentially diminishing significance of the canonical mutations (especially in South Africa and the Philippines), which could have significant future implications on narrowly targeting molecular diagnostics.

Identification of key residues in Debaryomyces hansenii carbonyl reductase for highly productive preparation of (S)-aryl halohydrins.

Key residues of Debaryomyces hansenii carbonyl reductase in the determination of the reducing activity towards aryl haloketones were identified through combinatorial mutation of conserved residues. This study provides a green and efficient biocatalyst for the synthesis of (S)-aryl halohydrins.

CONSTRICTOR: constraint modification provides insight into design of biochemical networks.

Advances in computational methods that allow for exploration of the combinatorial mutation space are needed to realize the potential of synthetic biology based strain engineering efforts. Here, we present Constrictor, a computational framework that uses flux balance analysis (FBA) to analyze inhibitory effects of genetic mutations on the performance of biochemical networks. Constrictor identifies engineering interventions by classifying the reactions in the metabolic model depending on the extent to which their flux must be decreased to achieve the overproduction target. The optimal inhibition of various reaction pathways is determined by restricting the flux through targeted reactions below the steady state levels of a baseline strain. Constrictor generates unique in silico strains, each representing an “expression state”, or a combination of gene expression levels required to achieve the overproduction target. The Constrictor framework is demonstrated by studying overproduction of ethylene in Escherichia coli network models iAF1260 and iJO1366 through the addition of the heterologous ethylene-forming enzyme from Pseudomonas syringae. Targeting individual reactions as well as combinations of reactions reveals in silico mutants that are predicted to have as high as 25% greater theoretical ethylene yields than the baseline strain during simulated exponential growth. Altering the degree of restriction reveals a large distribution of ethylene yields, while analysis of the expression states that return lower yields provides insight into system bottlenecks. Finally, we demonstrate the ability of Constrictor to scan networks and provide targets for a range of possible products. Constrictor is an adaptable technique that can be used to generate and analyze disparate populations of in silico mutants, select gene expression levels and provide non-intuitive strategies for metabolic engineering.

Design and construction of small perturbation mutagenesis libraries for antibody affinity maturation using massive microchip-synthesized oligonucleotides

We report a rational strategy to design and construct multiple small perturbation mutagenesis (SPM) libraries using massively parallel synthesis of oligonucleotides on a microchip for affinity maturation of an engineered anti-ErbB2 antibody chA21. On the basis of a comprehensive analysis of the sequence and structural relationships of six complementary determination regions (CDRs) in the Kabatman database, a computational algorithm was developed to introduce single-site and double-site mutations into variable CDR positions using ambiguous nucleotides. The six SPM libraries were composed of 419 degenerate oligonucleotides that can be expanded into 161,832 unique CDR sequences with a high coverage ratio of 95% natural amino acid diversity. We used Illumina next-generation sequencing to demonstrate that the synthetic CDR library sequences, as well as relative quantities per sequence, can be controlled precisely by adjusting reaction chamber assignment and input nucleoside composition. The microchip-synthesized oligonucleotides were used for construction of single-chain antibody fragment (scFv) phage libraries through one-step mutagenic PCR of double-stranded plasmids with >106E. coli transformants. A variant with combinatorial mutations from four individual CDRs achieved more than 19-fold affinity increase. The strategy described herein should be broadly applicable to affinity and selectivity studies of antibodies and other proteins.

A systematic analysis of the in vitro and in vivo functions of the HD-GYP domain proteins of Vibrio cholerae.

BackgroundThe second messenger cyclic diguanylate (c-di-GMP) plays a central role in bacterial adaptation to extracellular stimuli, controlling processes such as motility, biofilm development, cell development and, in some pathogens, virulence. The intracellular level of c-di-GMP is controlled by the complementary activities of diguanylate cyclases containing a GGDEF domain and two classes of c-di-GMP phosphodiesterases containing an EAL or HD-GYP hydrolytic domain. Compared to the GGDEF and EAL domains, the functions of HD-GYP domain family proteins are poorly characterized. The human diarrheal pathogen Vibrio cholerae encodes nine putative HD-GYP domain proteins. To determine the contributions of HD-GYP domain proteins to c-di-GMP signaling in V. cholerae, we systematically analyzed the enzymatic functionality of each protein and their involvement in processes known to be regulated by c-di-GMP: motility, biofilm development and virulence.ResultsComplementary in vitro and in vivo experiments showed that four HD-GYP domain proteins are active c-di-GMP phosphodiesterases: VC1295, VC1348, VCA0210 and VCA0681. Mutation of individual HD-GYP domain genes, as well as combinatorial mutations of multiple HD-GYP domain genes, had no effect on motility or biofilm formation of V. cholerae under the conditions tested. Furthermore, no single HD-GYP domain gene affected intestinal colonization by V. cholerae in an infant mouse model. However, inactivation of multiple HD-GYP domain genes, including the four encoding functional phosphodiesterases, significantly attenuated colonization.ConclusionsThese results indicate that the HD-GYP family of c-di-GMP phosphodiesterases impacts signaling by this second messenger during infection. Altogether, this work greatly furthers the understanding of this important family of c-di-GMP metabolic enzymes and demonstrates a role for HD-GYP domain proteins in the virulence of V. cholerae.Electronic supplementary materialThe online version of this article (doi:10.1186/s12866-014-0272-9) contains supplementary material, which is available to authorized users.

Directed evolution of human heavy chain variable domain (VH) using in vivo protein fitness filter.

Human immunoglobulin heavy chain variable domains (VH) are promising scaffolds for antigen binding. However, VH is an unstable and aggregation-prone protein, hindering its use for therapeutic purposes. To evolve the VH domain, we performed in vivo protein solubility selection that linked antibiotic resistance to the protein folding quality control mechanism of the twin-arginine translocation pathway of E. coli. After screening a human germ-line VH library, 95% of the VH proteins obtained were identified as VH3 family members; one VH protein, MG2x1, stood out among separate clones expressing individual VH variants. With further screening of combinatorial framework mutation library of MG2x1, we found a consistent bias toward substitution with tryptophan at the position of 50 and 58 in VH. Comparison of the crystal structures of the VH variants revealed that those substitutions with bulky side chain amino acids filled the cavity in the VH interface between heavy and light chains of the Fab arrangement along with the increased number of hydrogen bonds, decreased solvation energy, and increased negative charge. Accordingly, the engineered VH acquires an increased level of thermodynamic stability, reversible folding, and soluble expression. The library built with the VH variant as a scaffold was qualified as most of VH clones selected randomly were expressed as soluble form in E. coli regardless length of the combinatorial CDR. Furthermore, a non-aggregation feature of the selected VH conferred a free of humoral response in mice, even when administered together with adjuvant. As a result, this selection provides an alternative directed evolution pathway for unstable proteins, which are distinct from conventional methods based on the phage display.

Oncogenic transformation of diverse gastrointestinal tissues in primary organoid culture.

The application of primary organoid cultures containing epithelial and mesenchymal elements to cancer modeling holds promise for combining the accurate multilineage differentiation and physiology of in vivo systems with the facile in vitro manipulation of transformed cell lines. Here, a single air-liquid interface culture method was used without modification to engineer oncogenic mutations into primary epithelial/mesenchymal organoids from mouse colon, stomach and pancreas. Pancreatic and gastric organoids exhibited dysplasia upon KrasG12D expression and/or p53 loss, and readily generated adenocarcinoma upon in vivo transplantation. In contrast, primary colon organoids required combinatorial Apc, p53, KrasG12D and Smad4 mutations for progressive transformation to invasive adenocarcinoma-like histology in vitro and tumorigenicity in vivo, recapitulating multi-hit models of colorectal cancer (CRC), and versus more promiscuous transformation of small intestinal organoids. Colon organoid culture functionally validated the microRNA miR-483 as a dominant driver oncogene at the Insulin-like growth factor-2 (IGF2) 11p15.5 CRC amplicon, inducing dysplasia in vitro and tumorigenicity in vivo. These studies demonstrate the general utility of a highly tractable primary organoid system for cancer modeling and driver oncogene validation in diverse gastrointestinal tissues.