Gabriel Piedrafita

Human somatic tissues accumulate mutations during normal aging. Some of these affect cancer-associated driver genes and confer mutant progenitor cells a competitive advantage that leads to clonal expansions. The human esophageal epithelium exemplifies this phenomenon, becoming a dense mosaic of competing mutant clones by adulthood. However, the phenotypic consequence of those mutations and their possible role in carcinogenesis remains unknown. Novel bioinformatic tools for de novo mutant detection from single-cell transcriptomics (scRNA-seq) could potentially leverage on the wealth of publicly available data to help draw mutant cell phenotypes in vivo . In this study we test SComatic algorithm’s ability to identify somatic mutations in the normal, polyclonal esophageal epithelium. We analyze a public scRNA-seq dataset from a human cohort with multiple esophageal samples per donor, and an independent study in mice subjected to experimental mutagenesis where samples have been re-sequenced for validation. These unconventional experimental designs allow us to control unspecificity. We observe scRNA-seq variant calling output is heavily affected by undesired technical artifacts and germline variants, which we are able to reduce following a customized series of rational filters that enrich in somatic mutations. Final candidate mutations are then used to reconstitute clonal lineages and map them to differentiation trajectories in the UMAP embeddings. We find low read depth and sparse cellular sampling favor detection of passenger mutations and hinder driver mutant phenotypic inferences. Altogether, we showcase current limitations of scRNA-seq-derived mutation calling, while we offer methodological indications that should be considered for future studies aimed at investigating mutant clone behavior in normal polyclonal tissues from single-cell transcriptomics.

Aging epithelia are colonized by somatic mutations, which are subjected to selection influenced by intrinsic and extrinsic factors. The lack of suitable culture systems has slowed the study of this and other long-term biological processes. Here, we describe epithelioids, a facile, cost-effective method of culturing multiple mouse and human epithelia. Esophageal epithelioids self-maintain without passaging for at least 1 year, maintaining a three-dimensional structure with proliferative basal cells that differentiate into suprabasal cells, which eventually shed and retain genomic stability. Live imaging over 5 months showed that epithelioids replicate in vivo cell dynamics. Epithelioids support genetic manipulation and enable the study of mutant cell competition and selection in three-dimensional epithelia, and show how anti-cancer treatments modulate competition between transformed and wild-type cells. Finally, a targeted CRISPR–Cas9 screen shows that epithelioids recapitulate mutant gene selection in aging human esophagus and identifies additional drivers of clonal expansion, resolving the genetic networks underpinning competitive fitness.

Oncogenic PIK3CA mutations generate large clones in aging human esophagus. Here we investigate the behavior of Pik3ca mutant clones in the normal esophageal epithelium of transgenic mice. Expression of a heterozygous Pik3ca^H1047R mutation drives clonal expansion by tilting cell fate toward proliferation. CRISPR screening and inhibitor treatment of primary esophageal keratinocytes confirmed the PI3K–mTOR pathway increased mutant cell competitive fitness. The antidiabetic drug metformin reduced mutant cell advantage in vivo and in vitro. Conversely, metabolic conditions such as type 1 diabetes or diet-induced obesity enhanced the competitive fitness of Pik3ca^H1047R cells. Consistently, we found a higher density of PIK3CA gain-of-function mutations in the esophagus of individuals with high body mass index compared with those with normal weight. We conclude that the metabolic environment selectively influences the evolution of the normal epithelial mutational landscape. Clinically feasible interventions to even out signaling imbalances between wild-type and mutant cells may limit the expansion of oncogenic mutants in normal tissues.

NOTCH1 mutant clones occupy the majority of normal human esophagus by middle age but are comparatively rare in esophageal cancers, suggesting NOTCH1 mutations drive clonal expansion but impede carcinogenesis. Here we test this hypothesis. Sequencing NOTCH1 mutant clones in aging human esophagus reveals frequent biallelic mutations that block NOTCH1 signaling. In mouse esophagus, heterozygous Notch1 mutation confers a competitive advantage over wild-type cells, an effect enhanced by loss of the second allele. Widespread Notch1 loss alters transcription but has minimal effects on the epithelial structure and cell dynamics. In a carcinogenesis model, Notch1 mutations were less prevalent in tumors than normal epithelium. Deletion of Notch1 reduced tumor growth, an effect recapitulated by anti-NOTCH1 antibody treatment. Notch1 null tumors showed reduced proliferation. We conclude that Notch1 mutations in normal epithelium are beneficial as wild-type Notch1 favors tumor expansion. NOTCH1 blockade may have therapeutic potential in preventing esophageal squamous cancer.

Studying long-term biological processes such as the colonization of aging epithelia by somatic mutant clones has been slowed by the lack of suitable culture systems. Here we describe epithelioids, a facile, cost-effective method of culturing multiple mouse and human epithelia. Esophageal epithelioids self-maintain without passaging for at least a year, recapitulating the 3D structure, cell dynamics, transcriptome, and genomic stability of the esophagus. Live imaging over 5 months showed epithelioids replicate in vivo cell dynamics. Epithelioids enable the study of cell competition and mutant selection in 3D epithelia, and how anti-cancer treatments modulate the competition between transformed and wild type cells. Epithelioids are a novel method with a wide range of applications in epithelial tissues, particularly the study of long term processes, that cannot be accessed using other culture models.

Transposable elements (TE) have played a major role in configuring the structures of mammalian genomes through evolution. In normal conditions, the expression of these elements is repressed by different epigenetic regulation mechanisms such as DNA methylation, histone modification and regulation by small RNAs. TE re-activation is associated with stemness potential acquisition, regulation of innate immunity and disease, such as cancer. However, the vast majority of current knowledge in the field is based on bulk expression studies, and very little is known on cell-type- or state-specific expression of TE-derived transcripts. Therefore, cost-efficient single-cell-resolution TE expression analytical approaches are needed. We have implemented an analytical approach based on pseudoalignment to consensus sequences to incorporate TE expression information to scRNAseq data. All the data and code implemented are available as Supplementary data and in: https://github.com/jmzvillarreal/kallisto_TE_scRNAseq. Supplementary data are available at Bioinformatics online.

The structure of the metabolic network is highly conserved, but we know little about its evolutionary origins. Key for explaining the early evolution of metabolism is solving a chicken–egg dilemma, which describes that enzymes are made from the very same molecules they produce. The recent discovery of several nonenzymatic reaction sequences that topologically resemble central metabolism has provided experimental support for a “metabolism first” theory, in which at least part of the extant metabolic network emerged on the basis of nonenzymatic reactions. But how could evolution kick-start on the basis of a metal catalyzed reaction sequence, and how could the structure of nonenzymatic reaction sequences be imprinted on the metabolic network to remain conserved for billions of years? We performed an in vitro screening where we add the simplest components of metabolic enzymes, proteinogenic amino acids, to a nonenzymatic, iron-driven reaction network that resembles glycolysis and the pentose phosphate pathway (PPP). We observe that the presence of the amino acids enhanced several of the nonenzymatic reactions. Particular attention was triggered by a reaction that resembles a rate-limiting step in the oxidative PPP. A prebiotically available, proteinogenic amino acid cysteine accelerated the formation of RNA nucleoside precursor ribose-5-phosphate from 6-phosphogluconate. We report that iron and cysteine interact and have additive effects on the reaction rate so that ribose-5-phosphate forms at high specificity under mild, metabolism typical temperature and environmental conditions. We speculate that accelerating effects of amino acids on rate-limiting nonenzymatic reactions could have facilitated a stepwise enzymatization of nonenzymatic reaction sequences, imprinting their structure on the evolving metabolic network.

Human epithelial tissues accumulate cancer-driver mutations with age^{1,2,3,4,5,6,7,8,9}, yet tumour formation remains rare. The positive selection of these mutations suggests that they alter the behaviour and fitness of proliferating cells^10,11,12. Thus, normal adult tissues become a patchwork of mutant clones competing for space and survival, with the fittest clones expanding by eliminating their less competitive neighbours^11,12,13,14. However, little is known about how such dynamic competition in normal epithelia influences early tumorigenesis. Here we show that the majority of newly formed oesophageal tumours are eliminated through competition with mutant clones in the adjacent normal epithelium. We followed the fate of nascent, microscopic, pre-malignant tumours in a mouse model of oesophageal carcinogenesis and found that most were rapidly lost with no indication of tumour cell death, decreased proliferation or an anti-tumour immune response. However, deep sequencing of ten-day-old and one-year-old tumours showed evidence of selection on the surviving neoplasms. Induction of highly competitive clones in transgenic mice increased early tumour removal, whereas pharmacological inhibition of clonal competition reduced tumour loss. These results support a model in which survival of early neoplasms depends on their competitive fitness relative to that of mutant clones in the surrounding normal tissue. Mutant clones in normal epithelium have an unexpected anti-tumorigenic role in purging early tumours through cell competition, thereby preserving tissue integrity.

How metabolic pathways emerged in early evolution remains largely unknown. Recently discovered chemical networks driven by iron and sulfur resemble reaction sequences found within glycolysis, gluconeogenesis, the oxidative and reductive Krebs cycle, the Wood Ljungdahl as well as the S-adenosylmethionine pathways, components of the core cellular metabolic network. These findings suggest that the evolution of central metabolism was primed by environmental chemical reactions, implying that non-enzymatic reaction networks served as a “template” in the evolution of enzymatic activities. We speculated that the turning point for this transition would depend on the catalytic properties of the simplest structural components of proteins, single amino acids. Here, we systematically combine constituents of Fe(II)-driven non-enzymatic reactions resembling glycolysis and pentose phosphate pathway (PPP), with single proteinogenic amino acids. Multiple reaction rates are enhanced by amino acids. In particular, cysteine is able to replace (and/or complement) the metal ion Fe(II) in driving the non-enzymatic formation of the RNA-backbone metabolite ribose 5-phosphate from 6-phosphogluconate, a rate-limiting reaction of the oxidative PPP. In the presence of both Fe(II) and cysteine, a complex is formed, enabling the non-enzymatic reaction to proceed at a wide range of temperatures. At mundane temperatures, this ‘minimal enzyme-like complex’ achieves a much higher specificity in the formation of ribose 5-phosphate than the Fe(II)-driven reaction at high temperatures. Hence, simple amino acids can accelerate key steps within metal-promoted metabolism-like chemical networks. Our results imply a stepwise scenario, in which environmental chemical networks served as primers in the early evolution of the metabolic network structure. Significance Statement: The evolutionary roots of metabolic pathways are barely understood. Here we show results consistent with a stepwise scenario during the evolution of (enzymatic) metabolism, starting from non-enzymatic chemical networks. By systematic screening of metabolic-like reactivities in vitro , and using high-throughput analytical techniques, we identify an iron/cysteine complex to act as a ‘minimal enzymelike complex’, which consists of a metal ion, an amino acid, and a sugar phosphate ligand. Integrated in a metal-driven, non-enzymatic pentose phosphate pathway, it promotes the formation of the RNA-backbone precursor ribose 5-phosphate at ambient temperature.