The spatial organization of microbial communities during range expansion.

Microbes in nature often live in dense and diverse communities exhibiting a variety of spatial structures. Microbial range expansion is a universal ecological process that enables populations to form spatial patterns. It can be driven by both passive and active processes, for example, mechanical forces from cell growth and bacterial motility. In this review, we provide a taste of recent creative and sophisticated efforts being made to address basic questions in spatial ecology and pattern formation during range expansion. We especially highlight the role of motility to shape community structures, and discuss the research challenges and future directions.

General quantitative relations linking cell growth and the cell cycle in Escherichia coli.

Growth laws emerging from studies of cell populations provide essential constraints on the global mechanisms that coordinate cell growth1-3. The foundation of bacterial cell cycle studies relies on two interconnected dogmas that were proposed more than 50 years ago-the Schaechter-Maaloe-Kjeldgaard growth law that relates cell mass to growth rate1 and Donachie's hypothesis of a growth-rate-independent initiation mass4. These dogmas spurred many efforts to understand their molecular bases and physiological consequences5-14. Although they are generally accepted in the fast-growth regime, that is, for doubling times below 1 h, extension of these dogmas to the slow-growth regime has not been consistently achieved. Here, through a quantitative physiological study of Escherichia coli cell cycles over an extensive range of growth rates, we report that neither dogma holds in either the slow- or fast-growth regime. In their stead, linear relations between the cell mass and the rate of chromosome replication-segregation were found across the range of growth rates. These relations led us to propose an integral-threshold model in which the cell cycle is controlled by a licensing process, the rate of which is related in a simple way to chromosomal dynamics. These results provide a quantitative basis for predictive understanding of cell growth-cell cycle relationships.

Microfluidic Synchronizer Using a Synthetic Nanoparticle-Capped Bacterium.

Conventional techniques to synchronize bacterial cells often require manual manipulations and lengthy incubation lacking precise temporal control. An automated microfluidic device was recently developed to overcome these limitations. However, it exploits the stalk property of Caulobacter crescentus that undergoes asymmetric stalked and swarmer cell cycle stages and is therefore restricted to this species. To address this shortcoming, we have engineered Escherichia coli cells to adhere to microchannel walls via a synthetic and inducible "stalk". The pole of E. coli is capped by magnetic fluorescent nanoparticles via a polar-localized outer membrane protein. A mass of cells is immobilized in a microfluidic chamber by an externally applied magnetic field. Daughter cells are formed without the induced stalk and hence are flushed out, yielding a synchronous population of "baby" cells. The stalks can be tracked by GFP and nanoparticle fluorescence; no fluorescence signal is detected in the eluted cell population, indicating that it consists solely of daughters. The collected daughter cells display superb synchrony. The results demonstrate a new on-chip method to synchronize the model bacterium E. coli and likely other bacterial species, and also foster the application of synthetic biology to the study of the bacterial cell cycle.

No wisdom in the crowd: genome annotation in the era of big data - current status and future prospects.

Science and engineering rely on the accumulation and dissemination of knowledge to make discoveries and create new designs. Discovery-driven genome research rests on knowledge passed on via gene annotations. In response to the deluge of sequencing big data, standard annotation practice employs automated procedures that rely on majority rules. We argue this hinders progress through the generation and propagation of errors, leading investigators into blind alleys. More subtly, this inductive process discourages the discovery of novelty, which remains essential in biological research and reflects the nature of biology itself. Annotation systems, rather than being repositories of facts, should be tools that support multiple modes of inference. By combining deduction, induction and abduction, investigators can generate hypotheses when accurate knowledge is extracted from model databases. A key stance is to depart from 'the sequence tells the structure tells the function' fallacy, placing function first. We illustrate our approach with examples of critical or unexpected pathways, using MicroScope to demonstrate how tools can be implemented following the principles we advocate. We end with a challenge to the reader.

A human gut phage catalog correlates the gut phageome with type 2 diabetes.

BACKGROUND: Substantial efforts have been made to link the gut bacterial community to many complex human diseases. Nevertheless, the gut phages are often neglected. RESULTS: In this study, we used multiple bioinformatic methods to catalog gut phages from whole-community metagenomic sequencing data of fecal samples collected from both type II diabetes (T2D) patients (n = 71) and normal Chinese adults (n = 74). The definition of phage operational taxonomic units (pOTUs) and identification of large phage scaffolds (n = 2567, ≥ 10 k) revealed a comprehensive human gut phageome with a substantial number of novel sequences encoding genes that were unrelated to those in known phages. Interestingly, we observed a significant increase in the number of gut phages in the T2D group and, in particular, identified 7 pOTUs specific to T2D. This finding was further validated in an independent dataset of 116 T2D and 109 control samples. Co-occurrence/exclusion analysis of the bacterial genera and pOTUs identified a complex core interaction between bacteria and phages in the human gut ecosystem, suggesting that the significant alterations of the gut phageome cannot be explained simply by co-variation with the altered bacterial hosts. CONCLUSIONS: Alterations in the gut bacterial community have been linked to the chronic disease T2D, but the role of gut phages therein is not well understood. This is the first study to identify a T2D-specific gut phageome, indicating the existence of other mechanisms that might govern the gut phageome in T2D patients. These findings suggest the importance of the phageome in T2D risk, which warrants further investigation.

Transcription and methylation analyses of preleukemic promyelocytes indicate a dual role for PML/RARA in leukemia initiation.

Acute promyelocytic leukemia is an aggressive malignancy characterized by the accumulation of promyelocytes in the bone marrow. PML/RARA is the primary abnormality implicated in this pathology, but the mechanisms by which this chimeric fusion protein initiates disease are incompletely understood. Identifying PML/RARA targets in vivo is critical for comprehending the road to pathogenesis. Utilizing a novel sorting strategy, we isolated highly purified promyelocyte populations from normal and young preleukemic animals, carried out microarray and methylation profiling analyses, and compared the results from the two groups of animals. Surprisingly, in the absence of secondary lesions, PML/RARA had an overall limited impact on both the transcriptome and methylome. Of interest, we did identify down-regulation of secondary and tertiary granule genes as the first step engaging the myeloid maturation block. Although initially not sufficient to arrest terminal granulopoiesis in vivo, such alterations set the stage for the later, complete differentiation block seen in leukemia. Further, gene set enrichment analysis revealed that PML/RARA promyelocytes exhibit a subtle increase in expression of cell cycle genes, and we show that this leads to both increased proliferation of these cells and expansion of the promyelocyte compartment. Importantly, this proliferation signature was absent from the poorly leukemogenic p50/RARA fusion model, implying a critical role for PML in the altered cell-cycle kinetics and ability to initiate leukemia. Thus, our findings challenge the predominant model in the field and we propose that PML/RARA initiates leukemia by subtly shifting cell fate decisions within the promyelocyte compartment.

Differential Requirements for eIF4E Dose in Normal Development and Cancer.

eIF4E, the major cap-binding protein, has long been considered limiting for translating the mammalian genome. However, the eIF4E dose requirement at an organismal level remains unexplored. By generating an Eif4e haploinsufficient mouse, we found that a 50% reduction in eIF4E expression, while compatible with normal development and global protein synthesis, significantly impeded cellular transformation. Genome-wide translational profiling uncovered a translational program induced by oncogenic transformation and revealed a critical role for the dose of eIF4E, specifically in translating a network of mRNAs enriched for a unique 5' UTR signature. In particular, we demonstrate that the dose of eIF4E is essential for translating mRNAs that regulate reactive oxygen species, fueling transformation and cancer cell survival in vivo. Our findings indicate eIF4E is maintained at levels in excess for normal development that are hijacked by cancer cells to drive a translational program supporting tumorigenesis.

Gene expression profile identifies tyrosine kinase c-Met as a targetable mediator of antiangiogenic therapy resistance.

PURPOSE: To identify mediators of glioblastoma antiangiogenic therapy resistance and target these mediators in xenografts. EXPERIMENTAL DESIGN: We conducted microarray analysis comparing bevacizumab-resistant glioblastomas (BRG) with pretreatment tumors from the same patients. We established novel xenograft models of antiangiogenic therapy resistance to target candidate resistance mediator(s). RESULTS: BRG microarray analysis revealed upregulation versus pretreatment of receptor tyrosine kinase c-Met, which underwent further investigation because of its prior biologic plausibility as a bevacizumab resistance mediator. BRGs exhibited increased hypoxia versus pretreatment in a manner correlating with their c-Met upregulation, increased c-Met phosphorylation, and increased phosphorylation of c-Met-activated focal adhesion kinase and STAT3. We developed 2 novel xenograft models of antiangiogenic therapy resistance. In the first model, serial bevacizumab treatment of an initially responsive xenograft generated a xenograft with acquired bevacizumab resistance, which exhibited upregulated c-Met expression versus pretreatment. In the second model, a BRG-derived xenograft maintained refractoriness to the MRI tumor vasculature alterations and survival-promoting effects of bevacizumab. Growth of this BRG-derived xenograft was inhibited by a c-Met inhibitor. Transducing these xenograft cells with c-Met short hairpin RNA inhibited their invasion and survival in hypoxia, disrupted their mesenchymal morphology, and converted them from bevacizumab-resistant to bevacizumab-responsive. Engineering bevacizumab-responsive cells to express constitutively active c-Met caused these cells to form bevacizumab-resistant xenografts. CONCLUSION: These findings support the role of c-Met in survival in hypoxia and invasion, features associated with antiangiogenic therapy resistance, and growth and therapeutic resistance of xenografts resistant to antiangiogenic therapy. Therapeutically targeting c-Met could prevent or overcome antiangiogenic therapy resistance.

Epistatic interactions between Tgfb1 and genetic loci, Tgfbm2 and Tgfbm3, determine susceptibility to an asthmatic stimulus.

TGFß activation and signaling have been extensively studied in experimental models of allergen-induced asthma as potential therapeutic targets during chronic or acute phases of the disease. Outcomes of experimental manipulation of TGFß activity have been variable, in part due to use of different model systems. Using an ovalbumin (OVA)-induced mouse model of asthma, we here show that innate variation within TGFß1 genetic modifier loci, Tgfbm2 and Tgfbm3, alters disease susceptibility. Specifically, Tgfbm2(129) and Tgfbm3(C57) synergize to reverse accentuated airway hyperresponsiveness (AHR) caused by low TGFß1 levels in Tgfb1(+/-) mice of the NIH/OlaHsd strain. Moreover, epistatic interaction between Tgfbm2(129) and Tgfbm3(C57) uncouples the inflammatory response to ovalbumin from those of airway remodeling and airway hyperresponsiveness, illustrating independent genetic control of these responses. We conclude that differential inheritance of genetic variants of Tgfbm genes alters biological responses to reduced TGFß1 signaling in an experimental asthma model. TGFß antagonists for treatment of lung diseases might therefore give diverse outcomes, dependent on genetic variation.

Microarray analysis verifies two distinct phenotypes of glioblastomas resistant to antiangiogenic therapy.

PURPOSE: To identify mechanisms and mediators of resistance to antiangiogenic therapy in human glioblastoma. EXPERIMENTAL DESIGN: We carried out microarray gene expression analysis and immunohistochemistry comparing 21 recurrent glioblastomas progressing during antiangiogenic treatment with VEGF neutralizing antibody bevacizumab to paired pretreatment tumors from the same patients. RESULTS: Microarray analysis revealed that bevacizumab-resistant glioblastomas (BRG) had two clustering patterns defining subtypes that reflect radiographic growth patterns. Enhancing BRGs (EBRG) exhibited MRI enhancement, a long-established criterion for glioblastoma progression, and expressed mitogen-activated protein kinases, neural cell adhesion molecule-1 (NCAM-1), and aquaporin 4. Compared with their paired pretreatment tumors, EBRGs had unchanged vascularity and hypoxia, with increased proliferation. Nonenhancing BRGs (NBRG) exhibited minimal MRI enhancement but had FLAIR-bright expansion, a newer criterion for glioblastoma recurrence since the advent of antiangiogenic therapy, and expressed integrin α5, laminin, fibronectin1, and PDGFRß. NBRGs had less vascularity, more hypoxia, and unchanged proliferation than their paired pretreatment tumors. Primary NBRG cells exhibited more stellate morphology with a 3-fold increased shape factor and were nearly 4-fold more invasive in Matrigel chambers than primary cells from EBRGs or bevacizumab-naive glioblastomas (P < 0.05). CONCLUSION: Using microarray analysis, we found two resistance patterns during antiangiogenic therapy with distinct molecular profiles and radiographic growth patterns. These studies provide valuable biologic insight into the resistance that has limited antiangiogenic therapy to date.

Two distinct routes to oral cancer differing in genome instability and risk for cervical node metastasis.

PURPOSE: Problems in management of oral cancers or precancers include identification of patients at risk for metastasis, tumor recurrence, and second primary tumors or risk for progression of precancers (dysplasia) to cancer. Thus, the objective of this study was to clarify the role of genomic aberrations in oral cancer progression and metastasis. EXPERIMENTAL DESIGN: The spectrum of copy number alterations in oral dysplasia and squamous cell carcinomas (SCC) was determined by array comparative genomic hybridization. Associations with clinical characteristics were studied and results confirmed in an independent cohort. RESULTS: The presence of one or more of the chromosomal aberrations +3q24-qter, -8pter-p23.1, +8q12-q24.2, and +20 distinguishes a major subgroup (70%-80% of lesions, termed 3q8pq20 subtype) from the remainder (20%-30% of lesions, non-3q8pq20). The 3q8pq20 subtype is associated with chromosomal instability and differential methylation in the most chromosomally unstable tumors. The two subtypes differ significantly in clinical outcome with risk for cervical (neck) lymph node metastasis almost exclusively associated with the 3q8pq20 subtype in two independent oral SCC cohorts. CONCLUSIONS: Two subtypes of oral lesions indicative of at least two pathways for oral cancer development were distinguished that differ in chromosomal instability and risk for metastasis, suggesting that +3q,-8p, +8q, and +20 constitute a biomarker with clinical utility for identifying patients at risk for metastasis. Moreover, although increased numbers of genomic alterations can be harbingers of progression to cancer, dysplastic lesions lacking copy number changes cannot be considered benign as they are potential precursors to non-3q8pq20 locally invasive, yet not metastatic oral SCC.

EGAN: exploratory gene association networks.

UNLABELLED: Exploratory Gene Association Networks (EGAN) is a Java desktop application that provides a point-and-click environment for contextual graph visualization of high-throughput assay results. By loading the entire network of genes, pathways, interactions, annotation terms and literature references directly into memory, EGAN allows a biologist to repeatedly query and interpret multiple experimental results without incurring additional delays for data download/integration. Other compelling features of EGAN include: support for diverse -omics technologies, a simple and interactive graph display, sortable/searchable data tables, links to external web resources including > or = 240 000 articles at PubMed, hypergeometric and GSEA-like enrichment statistics, pipeline-compatible automation via scripting and the ability to completely customize and/or supplement the network with new/proprietary data. AVAILABILITY: Runs on most operating systems via Java; downloadable from http://akt.ucsf.edu/EGAN/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetic and phenotypic characteristics of pleomorphic lobular carcinoma in situ of the breast.

The clinical, pathologic, and molecular features of pleomorphic lobular carcinoma in situ (PLCIS) and the relationship of PLCIS to classic LCIS (CLCIS) are poorly defined. In this study, we analyzed 31 cases of PLCIS (13 apocrine and 18 nonapocrine subtypes) and compared the clinical, pathologic, immunophenotypic, and genetic characteristics of these cases with those of 24 cases of CLCIS. Biomarker expression was examined using immunostaining for E-cadherin, gross cystic disease fluid protein-15, estrogen, progesterone, androgen receptor, human epidermal growth factor receptor2, CK5/6, and Ki67. Array-based comparative genomic hybridization to assess the genomic alterations was performed using microdissected formalin-fixed paraffin-embedded samples. Patients with PLCIS presented with mammographic abnormalities. Histologically, the tumor cells were dyshesive and showed pleomorphic nuclei, and there was often associated necrosis and microcalcifications. All lesions were E-cadherin negative. Compared with CLCIS, PLCIS showed significantly higher Ki67 index, lower estrogen receptor and progesterone receptor expression, and higher incidence of HER2 gene amplification. The majority of PLCIS and CLCIS demonstrated loss of 16q and gain of 1q. Apocrine PLCIS had significantly more genomic alterations than CLCIS and nonapocrine PLCIS. Although lack of E-cadherin expression and the 16q loss and 1q gain-array-based comparative genomic hybridization pattern support a relationship to CLCIS, PLCIS has clinical, mammographic, histologic, immunophenotypic, and genetic features that distinguish it from CLCIS. The histologic features, biomarker profile, and genomic instability observed in PLCIS suggest a more aggressive phenotype than CLCIS. However, clinical follow-up studies will be required to define the natural history and most appropriate management of these lesions.

Acquired genomic aberrations associated with methotrexate resistance vary with background genomic instability.

Tumors vary widely in chromosomal level genome instability. To gain a better understanding of the underlying defects which foster specific types of aberrations, we investigated the response of cells of related genetic backgrounds to challenge with methotrexate. We studied mismatch repair deficient HCT116 cells, two derivatives also deficient in XRCC5 (HCT116 Ku86+/-) or BLM (HCT116 BLM-/-), and mismatch repair competent HCT116+chr3 cells. We show that colony formation occurred at a significantly higher frequency in HCT116 cells and HCT116 Ku86+/- cells compared to HCT116 BLM-/- and HCT116+chr3 cells. Visible colonies arose most rapidly in HCT116 Ku86+/- cells, whereas they formed most slowly in HCT116+chr3 cells. Copy number changes acquired by the methotrexate resistant HCT116 and HCT116 BLM-/- cells most often included whole chromosome gains or losses or no acquired copy number changes, whereas resistance in HCT116+chr3 and HCT116 Ku86+/- cells was associated with amplification of DHFR and copy number transitions leading to increased copy number of DHFR, respectively. The additional copies of DHFR were present on unstable chromosomes and organized as inverted repeats in HCT116+chr3 cells, while they were most often present as direct repeats in HCT116 Ku86+/- cells. These observations suggest that different mutational mechanisms promote drug resistance in these genetic backgrounds; mismatch repair deficiency in HCT116, high rates of chromosomal instability in HCT116 Ku86+/-, and low rates of chromosomal instability in HCT116+chr3. On the other hand, it appears that loss of BLM function suppresses the mismatch repair mutator mechanism in mismatch repair and BLM deficient HCT116 BLM-/- cells.

Genome position and gene amplification.

BACKGROUND: Amplifications, regions of focal high-level copy number change, lead to overexpression of oncogenes or drug resistance genes in tumors. Their presence is often associated with poor prognosis; however, the use of amplification as a mechanism for overexpression of a particular gene in tumors varies. To investigate the influence of genome position on propensity to amplify, we integrated a mutant form of the gene encoding dihydrofolate reductase into different positions in the human genome, challenged cells with methotrexate and then studied the genomic alterations arising in drug resistant cells. RESULTS: We observed site-specific differences in methotrexate sensitivity, amplicon organization and amplification frequency. One site was uniquely associated with a significantly enhanced propensity to amplify and recurrent amplicon boundaries, possibly implicating a rare folate-sensitive fragile site in initiating amplification. Hierarchical clustering of gene expression patterns and subsequent gene enrichment analysis revealed two clusters differing significantly in expression of MYC target genes independent of integration site. CONCLUSION: These studies suggest that genome context together with the particular challenges to genome stability experienced during the progression to cancer contribute to the propensity to amplify a specific oncogene or drug resistance gene, whereas the overall functional response to drug (or other) challenge may be independent of the genomic location of an oncogene.

Detection of single clone deletions using array CGH: identification of submicroscopic deletions in the 22q11.2 deletion syndrome as a model system.

Constitutional submicroscopic DNA copy number alterations have been shown to cause numerous medical genetic syndromes, and are suspected to occur in a portion of cases for which the causal events remain undiscovered. Array comparative genomic hybridization (array CGH) allows high-throughput, high-resolution genome scanning for DNA dosage aberrations and thus offers an attractive approach for both clinical diagnosis and discovery efforts. Here we assess this capability by applying array CGH to the analysis of copy number alterations in 44 patients with a phenotype of the 22q11.2 deletion syndrome. Twenty-five patients had the deletion on chromosome 22 characteristic of this syndrome as determined by fluorescence in situ hybridization (FISH). The array measurements were in complete concordance with the FISH analysis, supporting their diagnostic utility. These data show that a genome-scanning microarray has the level of sensitivity and specificity required to prospectively interrogate and identify single copy number aberrations in a clinical setting. We demonstrate that such technology is ideally suited for microdeletion syndromes such as 22q11.2.

Genomic and transcriptional aberrations linked to breast cancer pathophysiologies.

This study explores the roles of genome copy number abnormalities (CNAs) in breast cancer pathophysiology by identifying associations between recurrent CNAs, gene expression, and clinical outcome in a set of aggressively treated early-stage breast tumors. It shows that the recurrent CNAs differ between tumor subtypes defined by expression pattern and that stratification of patients according to outcome can be improved by measuring both expression and copy number, especially high-level amplification. Sixty-six genes deregulated by the high-level amplifications are potential therapeutic targets. Nine of these (FGFR1, IKBKB, ERBB2, PROCC, ADAM9, FNTA, ACACA, PNMT, and NR1D1) are considered druggable. Low-level CNAs appear to contribute to cancer progression by altering RNA and cellular metabolism.

Breast tumor copy number aberration phenotypes and genomic instability.

BACKGROUND: Genomic DNA copy number aberrations are frequent in solid tumors, although the underlying causes of chromosomal instability in tumors remain obscure. Genes likely to have genomic instability phenotypes when mutated (e.g. those involved in mitosis, replication, repair, and telomeres) are rarely mutated in chromosomally unstable sporadic tumors, even though such mutations are associated with some heritable cancer prone syndromes. METHODS: We applied array comparative genomic hybridization (CGH) to the analysis of breast tumors. The variation in the levels of genomic instability amongst tumors prompted us to investigate whether alterations in processes/genes involved in maintenance and/or manipulation of the genome were associated with particular types of genomic instability. RESULTS: We discriminated three breast tumor subtypes based on genomic DNA copy number alterations. The subtypes varied with respect to level of genomic instability. We find that shorter telomeres and altered telomere related gene expression are associated with amplification, implicating telomere attrition as a promoter of this type of aberration in breast cancer. On the other hand, the numbers of chromosomal alterations, particularly low level changes, are associated with altered expression of genes in other functional classes (mitosis, cell cycle, DNA replication and repair). Further, although loss of function instability phenotypes have been demonstrated for many of the genes in model systems, we observed enhanced expression of most genes in tumors, indicating that over expression, rather than deficiency underlies instability. CONCLUSION: Many of the genes associated with higher frequency of copy number aberrations are direct targets of E2F, supporting the hypothesis that deregulation of the Rb pathway is a major contributor to chromosomal instability in breast tumors. These observations are consistent with failure to find mutations in sporadic tumors in genes that have roles in maintenance or manipulation of the genome.

Fractional genomic alteration detected by array-based comparative genomic hybridization independently predicts survival after hepatic resection for metastatic colorectal cancer.

PURPOSE: Although liver resection is the primary curative therapy for patients with colorectal hepatic metastases, most patients have a recurrence. Identification of molecular markers that predict patients at highest risk for recurrence may help to target further therapy. EXPERIMENTAL DESIGN: Array-based comparative genomic hybridization was used to investigate the association of DNA copy number alterations with outcome in patients with colorectal liver metastasis resected with curative intent. DNA from 50 liver metastases was labeled and hybridized onto an array consisting of 2,463 bacterial artificial chromosome clones covering the entire genome. The total fraction of genome altered (FGA) in the metastases and the patient's clinical risk score (CRS) were calculated to identify independent prognostic factors for survival. RESULTS: An average of 30 +/- 14% of the genome was altered in the liver metastases (14% gained and 16% lost). As expected, a lower CRS was an independent predictor of overall survival (P = 0.03). In addition, a high FGA also was an independent predictor of survival (P = 0.01). The median survival time in patients with a low CRS (score 0-2) and a high (> or =20%) FGA was 38 months compared with 18 months in patients with a low CRS and a low FGA. Supervised analyses, using Prediction Analysis of Microarrays and Significance Analysis of Microarrays, identified a set of clones, predominantly located on chromosomes 7 and 20, which best predicted survival. CONCLUSIONS: Both FGA and CRS are independent predictors of survival in patients with resected hepatic colorectal cancer metastases. The greater the FGA, the more likely the patient is to survive.