Autophagy is required during high MUC2 mucin biosynthesis in colonic goblet cells to contend metabolic stress.

Goblet cells are specialized for the production and secretion of MUC2 glycoproteins that forms a thick layer covering the mucosal epithelium as a protective barrier against noxious substances and invading microbes. High MUC2 mucin biosynthesis induces endoplasmic reticulum (ER) stress and apoptosis in goblet cells during inflammatory and infectious diseases. Autophagy is an intracellular degradation process required for maintenance of intestinal homeostasis. In this study, we hypothesized that autophagy was triggered during high MUC2 mucin biosynthesis from colonic goblet cells to cope with metabolic stress. To interrogate this, we analyzed the autophagy process in high MUC2-producing human HT29-H and a clone HT29-L silenced for MUC2 expression by lentivirus-mediated shRNA, and WT and CRISPR/Cas9 MUC2 KO LS174T cells. Autophagy was constitutively increased in high MUC2-producing cells characterized by elevated pULK1S555 expression and increased numbers of autophagosomes as compared with MUC2 silenced or gene edited cells. Similarly, colonoids from Muc2+/+ but not Muc2-/- littermates differentiated into goblet cells showed increased autophagy. IL-22 treatment corrected misfolded MUC2 protein and alleviated the autophagy process in LS174T cells. This study highlights that autophagy plays an essential role in goblet cells to survive during high mucin biosynthesis by regulating cellular homeostasis.NEW & NOTEWORTHY It is unclear how colonic goblet cells survive by producing high output MUC2 mucin that triggers endoplasmic stress by misfolded MUC2 proteins. To cope with metabolic stress, we interrogated if autophagy played an essential role in regulating cellular homeostasis. Indeed, high MUC2 mucin biosynthesis dysregulated autophagy processes that was regulated by IL-22 to maintain gut barrier innate host defenses.

Entamoeba histolytica exploits the autophagy pathway in macrophages to trigger inflammation in disease pathogenesis.

The mechanism whereby Entamoeba histolytica (Eh) binding with macrophages at the intercellular junction triggers aggressive pro-inflammatory responses in disease pathogenesis is not well understood. The host intracellular protein degradation process autophagy and its regulatory proteins are involved in maintenance of cellular homeostasis and excessive inflammatory responses. In this study we unraveled how Eh hijacks the autophagy process in macrophages to dysregulate pro-inflammatory responses. Direct contact of live Eh with macrophages activated caspase-6 that induced rapid proteolytic degradation of the autophagy ATG16L1 protein complex independent of NLRP3 inflammasome and caspase-3/8 activation. Crohn's disease susceptible ATG16L1 T300A variant was highly susceptible to Eh-mediated degradation that augmented pro-inflammatory cytokines in mice. Quantitative proteomics revealed downregulation of autophagy and vesicle-mediated transport and upregulation of cysteine-type endopeptidase pathways in response to Eh. We conclude during Eh-macrophage outside-in signaling, ATG16L1 protein complex plays an overlooked regulatory role in shaping the pro-inflammatory landscape in amebiasis.

Role of inflammasomes in innate host defense against Entamoeba histolytica.

Intestinal amebiasis is the disease caused by the extracellular protozoan parasite Entamoeba histolytica (Eh) that induces a dynamic and heterogeneous interaction profile with the host immune system during disease pathogenesis. In 90% of asymptomatic infection, Eh resides with indigenous microbiota in the outer mucus layer of the colon without prompting an immune response. However, for reasons that remain unclear, in a minority of the Eh-infected individuals, this fine tolerated relationship is switched to a pathogenic phenotype and advanced to an increasingly complex host-parasite interaction. Eh disease susceptibility depends on parasite virulence factors and their interactions with indigenous bacteria, disruption of the mucus bilayers, and adherence to the epithelium provoking host immune cells to evoke a robust pro-inflammatory response mediated by inflammatory caspases and inflammasome activation. To understand Eh pathogenicity and innate host immune responses, this review highlights recent advances in our understanding of how Eh induces outside-in signaling via Mϕs to activate inflammatory caspases and inflammasome to regulate pro-inflammatory responses.

Entamoeba histolytica stimulates the alarmin molecule HMGB1 from macrophages to amplify innate host defenses.

Even though Entamoeba histolytica (Eh)-induced host pro-inflammatory responses play a critical role in disease, we know very little about the host factors that regulate this response. Direct contact between host cell and Eh signify the highest level of danger, and to eliminate this threat, the host immune system elicits an augmented immune response. To understand the mechanisms of this response, we investigated the induction and release of the endogenous alarmin molecule high-mobility group box 1 (HMGB1) that act as a pro-inflammatory cytokine and chemoattractant during Eh infection. Eh in contact with macrophage induced a dose- and time-dependent secretion of HMGB1 in the absence of cell death. Secretion of HMGB1 was facilitated by Eh surface Gal-lectin-activated phosphoinositide 3-kinase and nuclear factor-κB signaling and up-regulation of histone acetyltransferase activity to trigger acetylated HMGB1 translocation from the nucleus. Unlike lipopolysaccharide, Eh-induced HMGB1 release was independent of caspase-1-mediated inflammasome and gasdermin D pores. In vivo, Eh inoculation in specific pathogen-free but not germ-free mice was associated with high levels of pro-inflammatory cytokines such as tumor necrosis factor-α, interleukin-1ß, and keratinocyte-derived chemokine, which was suppressed with HMGB1 neutralization. This study reveals that Eh-induced active secretion of the HMGB1 plays a key role in shaping the pro-inflammatory landscape critical in innate host defense against amebiasis.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

The macrophage cytoskeleton acts as a contact sensor upon interaction with Entamoeba histolytica to trigger IL-1ß secretion.

Entamoeba histolytica (Eh) is the causative agent of amebiasis, one of the major causes of dysentery-related morbidity worldwide. Recent studies have underlined the importance of the intercellular junction between Eh and host cells as a determinant in the pathogenesis of amebiasis. Despite the fact that direct contact and ligation between Eh surface Gal-lectin and EhCP-A5 with macrophage α5ß1 integrin are absolute requirements for NLRP3 inflammasome activation and IL-1ß release, many other undefined molecular events and downstream signaling occur at the interface of Eh and macrophage. In this study, we investigated the molecular events at the intercellular junction that lead to recognition of Eh through modulation of the macrophage cytoskeleton. Upon Eh contact with macrophages key cytoskeletal-associated proteins were rapidly post-translationally modified only with live Eh but not with soluble Eh proteins or fragments. Eh ligation with macrophages rapidly activated caspase-6 dependent cleavage of the cytoskeletal proteins talin, Pyk2 and paxillin and caused robust release of the pro-inflammatory cytokine, IL-1ß. Macrophage cytoskeletal cleavages were dependent on Eh cysteine proteinases EhCP-A1 and EhCP-A4 but not EhCP-A5 based on pharmacological blockade of Eh enzyme inhibitors and EhCP-A5 deficient parasites. These results unravel a model where the intercellular junction between macrophages and Eh form an area of highly interacting proteins that implicate the macrophage cytoskeleton as a sensor for Eh contact that leads downstream to subsequent inflammatory immune responses.

BCL9L Dysfunction Impairs Caspase-2 Expression Permitting Aneuploidy Tolerance in Colorectal Cancer.

Chromosomal instability (CIN) contributes to cancer evolution, intratumor heterogeneity, and drug resistance. CIN is driven by chromosome segregation errors and a tolerance phenotype that permits the propagation of aneuploid genomes. Through genomic analysis of colorectal cancers and cell lines, we find frequent loss of heterozygosity and mutations in BCL9L in aneuploid tumors. BCL9L deficiency promoted tolerance of chromosome missegregation events, propagation of aneuploidy, and genetic heterogeneity in xenograft models likely through modulation of Wnt signaling. We find that BCL9L dysfunction contributes to aneuploidy tolerance in both TP53-WT and mutant cells by reducing basal caspase-2 levels and preventing cleavage of MDM2 and BID. Efforts to exploit aneuploidy tolerance mechanisms and the BCL9L/caspase-2/BID axis may limit cancer diversity and evolution.

Optimization of Quantitative PCR Methods for Enteropathogen Detection.

Detection and quantification of enteropathogens in stool specimens is useful for diagnosing the cause of diarrhea but is technically challenging. Here we evaluate several important determinants of quantification: specimen collection, nucleic acid extraction, and extraction and amplification efficiency. First, we evaluate the molecular detection and quantification of pathogens in rectal swabs versus stool, using paired flocked rectal swabs and whole stool collected from 129 children hospitalized with diarrhea in Tanzania. Swabs generally yielded a higher quantification cycle (Cq) (average 29.7, standard deviation 3.5 vs. 25.3 ± 2.9 from stool, P<0.001) but were still able to detect 80% of pathogens with a Cq < 30 in stool. Second, a simplified total nucleic acid (TNA) extraction procedure was compared to separate DNA and RNA extractions and showed 92% (318/344) sensitivity and 98% (951/968) specificity, with no difference in Cq value for the positive results (ΔCq(DNA+RNA-TNA) = -0.01 ± 1.17, P = 0.972, N = 318). Third, we devised a quantification scheme that adjusts pathogen quantity to the specimen's extraction and amplification efficiency, and show that this better estimates the quantity of spiked specimens than the raw target Cq. In sum, these methods for enteropathogen quantification, stool sample collection, and nucleic acid extraction will be useful for laboratories studying enteric disease.

Impact of enterovirus and other enteric pathogens on oral polio and rotavirus vaccine performance in Bangladeshi infants.

BACKGROUND: Oral polio vaccine (OPV) and rotavirus vaccine (RV) exhibit poorer performance in low-income settings compared to high-income settings. Prior studies have suggested an inhibitory effect of concurrent non-polio enterovirus (NPEV) infection, but the impact of other enteric infections has not been comprehensively evaluated. METHODS: In urban Bangladesh, we tested stools for a broad range of enteric viruses, bacteria, parasites, and fungi by quantitative PCR from infants at weeks 6 and 10 of life, coincident with the first OPV and RV administration respectively, and examined the association between enteropathogen quantity and subsequent OPV serum neutralizing titers, serum rotavirus IgA, and rotavirus diarrhea. RESULTS: Campylobacter and enterovirus (EV) quantity at the time of administration of the first dose of OPV was associated with lower OPV1-2 serum neutralizing titers, while enterovirus quantity was also associated with diminished rotavirus IgA (-0.08 change in log titer per tenfold increase in quantity; P=0.037), failure to seroconvert (OR 0.78, 95% CI: 0.64-0.96; P=0.022), and breakthrough rotavirus diarrhea (OR 1.34, 95% CI: 1.05-1.71; P=0.020) after adjusting for potential confounders. These associations were not observed for Sabin strain poliovirus quantity. CONCLUSION: In this broad survey of enteropathogens and oral vaccine performance we find a particular association between EV carriage, particularly NPEV, and OPV immunogenicity and RV protection. Strategies to reduce EV infections may improve oral vaccine responses. ClinicalTrials.gov Identifier: NCT01375647.

Role of the Gut Microbiota of Children in Diarrhea Due to the Protozoan Parasite Entamoeba histolytica.

BACKGROUND: An estimated 1 million children die each year before their fifth birthday from diarrhea. Previous population-based surveys of pediatric diarrheal diseases have identified the protozoan parasite Entamoeba histolytica, the etiological agent of amebiasis, as one of the causes of moderate-to-severe diarrhea in sub-Saharan Africa and South Asia. METHODS: We prospectively studied the natural history of E. histolytica colonization and diarrhea among infants in an urban slum of Dhaka, Bangladesh. RESULTS: Approximately 80% of children were infected with E. histolytica by the age of 2 years. Fecal anti-galactose/N-acetylgalactosamine lectin immunoglobulin A was associated with protection from reinfection, while a high parasite burden and expansion of the Prevotella copri level was associated with diarrhea. CONCLUSIONS: E. histolytica infection was prevalent in this population, with most infections asymptomatic and diarrhea associated with both the amount of parasite and the composition of the microbiota.

Immune Evasion Mechanisms of Entamoeba histolytica: Progression to Disease.

Entamoeba histolytica (Eh) is a protozoan parasite that infects 10% of the world's population and results in 100,000 deaths/year from amebic dysentery and/or liver abscess. In most cases, this extracellular parasite colonizes the colon by high affinity binding to MUC2 mucin without disease symptoms, whereas in some cases, Eh triggers an aggressive inflammatory response upon invasion of the colonic mucosa. The specific host-parasite factors critical for disease pathogenesis are still not well characterized. From the parasite, the signature events that lead to disease progression are cysteine protease cleavage of the C-terminus of MUC2 that dissolves the mucus layer followed by Eh binding and cytotoxicity of the mucosal epithelium. The host mounts an ineffective excessive host pro-inflammatory response following contact with host cells that causes tissue damage and participates in disease pathogenesis as Eh escapes host immune clearance by mechanisms that are not completely understood. Ameba can modulate or destroy effector immune cells by inducing neutrophil apoptosis and suppressing respiratory burst or nitric oxide (NO) production from macrophages. Eh adherence to the host cells also induce multiple cytotoxic effects that can promote cell death through phagocytosis, apoptosis or by trogocytosis (ingestion of living cells) that might play critical roles in immune evasion. This review focuses on the immune evasion mechanisms that Eh uses to survive and induce disease manifestation in the host.

Tracking the genomic evolution of esophageal adenocarcinoma through neoadjuvant chemotherapy.

UNLABELLED: Esophageal adenocarcinomas are associated with a dismal prognosis. Deciphering the evolutionary history of this disease may shed light on therapeutically tractable targets and reveal dynamic mutational processes during the disease course and following neoadjuvant chemotherapy (NAC). We exome sequenced 40 tumor regions from 8 patients with operable esophageal adenocarcinomas, before and after platinum-containing NAC. This revealed the evolutionary genomic landscape of esophageal adenocarcinomas with the presence of heterogeneous driver mutations, parallel evolution, early genome-doubling events, and an association between high intratumor heterogeneity and poor response to NAC. Multiregion sequencing demonstrated a significant reduction in thymine to guanine mutations within a CpTpT context when comparing early and late mutational processes and the presence of a platinum signature with enrichment of cytosine to adenine mutations within a CpC context following NAC. Esophageal adenocarcinomas are characterized by early chromosomal instability leading to amplifications containing targetable oncogenes persisting through chemotherapy, providing a rationale for future therapeutic approaches. SIGNIFICANCE: This work illustrates dynamic mutational processes occurring during esophageal adenocarcinoma evolution and following selective pressures of platinum exposure, emphasizing the iatrogenic impact of therapy on cancer evolution. Identification of amplifications encoding targetable oncogenes maintained through NAC suggests the presence of stable vulnerabilities, unimpeded by cytotoxics, suitable for therapeutic intervention.

Recurrent chromosomal gains and heterogeneous driver mutations characterise papillary renal cancer evolution.

Papillary renal cell carcinoma (pRCC) is an important subtype of kidney cancer with a problematic pathological classification and highly variable clinical behaviour. Here we sequence the genomes or exomes of 31 pRCCs, and in four tumours, multi-region sequencing is undertaken. We identify BAP1, SETD2, ARID2 and Nrf2 pathway genes (KEAP1, NHE2L2 and CUL3) as probable drivers, together with at least eight other possible drivers. However, only ~10% of tumours harbour detectable pathogenic changes in any one driver gene, and where present, the mutations are often predicted to be present within cancer sub-clones. We specifically detect parallel evolution of multiple SETD2 mutations within different sub-regions of the same tumour. By contrast, large copy number gains of chromosomes 7, 12, 16 and 17 are usually early, monoclonal changes in pRCC evolution. The predominance of large copy number variants as the major drivers for pRCC highlights an unusual mode of tumorigenesis that may challenge precision medicine approaches.

Kinetics of poliovirus shedding following oral vaccination as measured by quantitative reverse transcription-PCR versus culture.

Amid polio eradication efforts, detection of oral polio vaccine (OPV) virus in stool samples can provide information about rates of mucosal immunity and allow estimation of the poliovirus reservoir. We developed a multiplex one-step quantitative reverse transcription-PCR (qRT-PCR) assay for detection of OPV Sabin strains 1, 2, and 3 directly in stool samples with an external control to normalize samples for viral quantity and compared its performance with that of viral culture. We applied the assay to samples from infants in Dhaka, Bangladesh, after the administration of trivalent OPV (tOPV) at weeks 14 and 52 of life (on days 0 [pre-OPV], +4, +11, +18, and +25 relative to vaccination). When 1,350 stool samples were tested, the sensitivity and specificity of the quantitative PCR (qPCR) assay were 89 and 91% compared with culture. A quantitative relationship between culture(+)/qPCR(+) and culture(-)/qPCR(+) stool samples was observed. The kinetics of shedding revealed by qPCR and culture were similar. qPCR quantitative cutoffs based on the day +11 or +18 stool samples could be used to identify the culture-positive shedders, as well as the long-duration or high-frequency shedders. Interestingly, qPCR revealed that a small minority (7%) of infants contributed the vast majority (93 to 100%) of the total estimated viral excretion across all subtypes at each time point. This qPCR assay for OPV can simply and quantitatively detect all three Sabin strains directly in stool samples to approximate shedding both qualitatively and quantitatively.

Molecular genotyping and quantitation assay for rotavirus surveillance.

Rotavirus genotyping is useful for surveillance purposes especially in areas where rotavirus vaccination has been or will be implemented. RT-PCR based molecular methods have been applied widely, but quantitative assays targeting a broad spectrum of genotypes have not been developed. Three real time RT-PCR panels were designed to identify G1, G2, G9, G12 (panel GI), G3, G4, G8, G10 (panel GII), and P[4], P[6], P[8], P[10], P[11] (panel P), respectively. An assay targeting NSP3 was included in both G panels as an internal control. The cognate assays were also formulated as one RT-PCR-Luminex panel for simultaneous detection of all the genotypes listed above plus P[9]. The assays were evaluated with various rotavirus isolates and 89 clinical samples from Virginia, Bangladesh and Tanzania, and exhibited 95% (81/85) sensitivity compared with the conventional RT-PCR-Gel-electrophoresis method, and 100% concordance with sequencing. Real time assays identified a significantly higher rate of mixed genotypes in Bangladeshi samples than the conventional gel-electrophoresis-based RT-PCR assay (32.5% versus 12.5%, P<0.05). In these mixed infections, the relative abundance of the rotavirus types could be estimated by Cq values. These typing assays detect and discriminate a broad range of G/P types circulating in different geographic regions with high sensitivity and specificity and can be used for rotavirus surveillance.

Spatial and temporal diversity in genomic instability processes defines lung cancer evolution.

Spatial and temporal dissection of the genomic changes occurring during the evolution of human non-small cell lung cancer (NSCLC) may help elucidate the basis for its dismal prognosis. We sequenced 25 spatially distinct regions from seven operable NSCLCs and found evidence of branched evolution, with driver mutations arising before and after subclonal diversification. There was pronounced intratumor heterogeneity in copy number alterations, translocations, and mutations associated with APOBEC cytidine deaminase activity. Despite maintained carcinogen exposure, tumors from smokers showed a relative decrease in smoking-related mutations over time, accompanied by an increase in APOBEC-associated mutations. In tumors from former smokers, genome-doubling occurred within a smoking-signature context before subclonal diversification, which suggested that a long period of tumor latency had preceded clinical detection. The regionally separated driver mutations, coupled with the relentless and heterogeneous nature of the genome instability processes, are likely to confound treatment success in NSCLC.

Development of synchronous VHL syndrome tumors reveals contingencies and constraints to tumor evolution.

BACKGROUND: Genomic analysis of multi-focal renal cell carcinomas from an individual with a germline VHL mutation offers a unique opportunity to study tumor evolution. RESULTS: We perform whole exome sequencing on four clear cell renal cell carcinomas removed from both kidneys of a patient with a germline VHL mutation. We report that tumors arising in this context are clonally independent and harbour distinct secondary events exemplified by loss of chromosome 3p, despite an identical genetic background and tissue microenvironment. We propose that divergent mutational and copy number anomalies are contingent upon the nature of 3p loss of heterozygosity occurring early in tumorigenesis. However, despite distinct 3p events, genomic, proteomic and immunohistochemical analyses reveal evidence for convergence upon the PI3K-AKT-mTOR signaling pathway. Four germline tumors in this young patient, and in a second, older patient with VHL syndrome demonstrate minimal intra-tumor heterogeneity and mutational burden, and evaluable tumors appear to follow a linear evolutionary route, compared to tumors from patients with sporadic clear cell renal cell carcinoma. CONCLUSIONS: In tumors developing from a germline VHL mutation, the evolutionary principles of contingency and convergence in tumor development are complementary. In this small set of patients with early stage VHL-associated tumors, there is reduced mutation burden and limited evidence of intra-tumor heterogeneity.

Development and assessment of molecular diagnostic tests for 15 enteropathogens causing childhood diarrhoea: a multicentre study.

BACKGROUND: Childhood diarrhoea can be caused by many pathogens that are difficult to assay in the laboratory. Molecular diagnostic techniques provide a uniform method to detect and quantify candidate enteropathogens. We aimed to develop and assess molecular tests for identification of enteropathogens and their association with disease. METHODS: We developed and assessed molecular diagnostic tests for 15 enteropathogens across three platforms-PCR-Luminex, multiplex real-time PCR, and TaqMan array card-at five laboratories worldwide. We judged the analytical and clinical performance of these molecular techniques against comparator methods (bacterial culture, ELISA, and PCR) using 867 diarrhoeal and 619 non-diarrhoeal stool specimens. We also measured molecular quantities of pathogens to predict the association with diarrhoea, by univariate logistic regression analysis. FINDINGS: The molecular tests showed very good analytical and clinical performance at all five laboratories. Comparator methods had limited sensitivity compared with the molecular techniques (20-85% depending on the target) but good specificity (median 97·3%, IQR 96·5-98·9; mean 95·2%, SD 9·1). Positive samples by comparator methods usually had higher molecular quantities of pathogens than did negative samples, across almost all platforms and for most pathogens (p<0·05). The odds ratio for diarrhoea at a given quantity (measured by quantification cycle, Cq) showed that for most pathogens associated with diarrhoea-including Campylobacter jejuni and Campylobacter coli, Cryptosporidium spp, enteropathogenic Escherichia coli, heat-stable enterotoxigenic E coli, rotavirus, Shigella spp and enteroinvasive E coli, and Vibrio cholerae-the strength of association with diarrhoea increased at higher pathogen loads. For example, Shigella spp at a Cq range of 15-20 had an odds ratio of 8·0 (p<0·0001), but at a Cq range of 25-30 the odds ratio fell to 1·7 (p=0·043). INTERPRETATION: Molecular diagnostic tests can be implemented successfully and with fidelity across laboratories around the world. In the case of diarrhoea, these techniques can detect pathogens with high sensitivity and ascribe diarrhoeal associations based on quantification, including in mixed infections, providing rich and unprecedented measurements of infectious causes. FUNDING: Bill & Melinda Gates Foundation Next Generation Molecular Diagnostics Project.

Identification of an LGP2-associated MDA5 agonist in picornavirus-infected cells.

The RIG-I-like receptors RIG-I, LGP2, and MDA5 initiate an antiviral response that includes production of type I interferons (IFNs). The nature of the RNAs that trigger MDA5 activation in infected cells remains unclear. Here, we purify and characterise LGP2/RNA complexes from cells infected with encephalomyocarditis virus (EMCV), a picornavirus detected by MDA5 and LGP2 but not RIG-I. We show that those complexes contain RNA that is highly enriched for MDA5-stimulatory activity and for a specific sequence corresponding to the L region of the EMCV antisense RNA. Synthesis of this sequence by in vitro transcription is sufficient to generate an MDA5 stimulatory RNA. Conversely, genomic deletion of the L region in EMCV generates viruses that are less potent at stimulating MDA5-dependent IFN production. Thus, the L region antisense RNA of EMCV is a key determinant of innate immunity to the virus and represents an RNA that activates MDA5 in virally-infected cells. DOI: http://dx.doi.org/10.7554/eLife.01535.001.