3D genomic mapping reveals multifocality of human pancreatic precancers.

Pancreatic intraepithelial neoplasias (PanINs) are the most common precursors of pancreatic cancer, but their small size and inaccessibility in humans make them challenging to study1. Critically, the number, dimensions and connectivity of human PanINs remain largely unknown, precluding important insights into early cancer development. Here, we provide a microanatomical survey of human PanINs by analysing 46 large samples of grossly normal human pancreas with a machine-learning pipeline for quantitative 3D histological reconstruction at single-cell resolution. To elucidate genetic relationships between and within PanINs, we developed a workflow in which 3D modelling guides multi-region microdissection and targeted and whole-exome sequencing. From these samples, we calculated a mean burden of 13 PanINs per cm3 and extrapolated that the normal intact adult pancreas harbours hundreds of PanINs, almost all with oncogenic KRAS hotspot mutations. We found that most PanINs originate as independent clones with distinct somatic mutation profiles. Some spatially continuous PanINs were found to contain multiple KRAS mutations; computational and in situ analyses demonstrated that different KRAS mutations localize to distinct cell subpopulations within these neoplasms, indicating their polyclonal origins. The extensive multifocality and genetic heterogeneity of PanINs raises important questions about mechanisms that drive precancer initiation and confer differential progression risk in the human pancreas. This detailed 3D genomic mapping of molecular alterations in human PanINs provides an empirical foundation for early detection and rational interception of pancreatic cancer.

Cell-type-resolved mosaicism reveals clonal dynamics of the human forebrain.

Debate remains around the anatomical origins of specific brain cell subtypes and lineage relationships within the human forebrain1-7. Thus, direct observation in the mature human brain is critical for a complete understanding of its structural organization and cellular origins. Here we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific mosaic variant barcode analysis. From four hemispheres and two different human neurotypical donors, we identified 287 and 780 mosaic variants, respectively, that were used to deconvolve clonal dynamics. Clonal spread and allele fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted than resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome transcriptome analysis at both a cell-type-specific and a single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of mosaic variants across 17 locations within one parietal lobe reveals that restriction of clonal spread in the anterior-posterior axis precedes restriction in the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus, cell-type-resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Extraction of the CDRH3 sequence of the mouse antibody repertoire selected upon influenza virus infection by subtraction of the background antibody repertoire.

Historically, antibody reactivity to pathogens and vaccine antigens has been evaluated using serological measurements of antigen-specific antibodies. However, it is difficult to evaluate all antibodies that contribute to various functions in a single assay, such as the measurement of the neutralizing antibody titer. Bulk antibody repertoire analysis using next-generation sequencing is a comprehensive method for analyzing the overall antibody response; however, it is unreliable for estimating antigen-specific antibodies due to individual variation. To address this issue, we propose a method to subtract the background signal from the repertoire of data of interest. In this study, we analyzed changes in antibody diversity and inferred the heavy-chain complementarity-determining region 3 (CDRH3) sequences of antibody clones that were selected upon influenza virus infection in a mouse model using bulk repertoire analysis. A decrease in the diversity of the antibody repertoire was observed upon viral infection, along with an increase in neutralizing antibody titers. Using kernel density estimation of sequences in a high-dimensional sequence space with background signal subtraction, we identified several clusters of CDRH3 sequences induced upon influenza virus infection. Most of these repertoires were detected more frequently in infected mice than in uninfected control mice, suggesting that infection-specific antibody sequences can be extracted using this method. Such an accurate extraction of antigen- or infection-specific repertoire information will be a useful tool for vaccine evaluation in the future. IMPORTANCE: As specific interactions between antigens and cell-surface antibodies trigger the proliferation of B-cell clones, the frequency of each antibody sequence in the samples reflects the size of each clonal population. Nevertheless, it is extremely difficult to extract antigen-specific antibody sequences from the comprehensive bulk antibody sequences obtained from blood samples due to repertoire bias influenced by exposure to dietary antigens and other infectious agents. This issue can be addressed by subtracting the background noise from the post-immunization or post-infection repertoire data. In the present study, we propose a method to quantify repertoire data from comprehensive repertoire data. This method allowed subtraction of the background repertoire, resulting in more accurate extraction of expanded antibody repertoires upon influenza virus infection. This accurate extraction of antigen- or infection-specific repertoire information is a useful tool for vaccine evaluation.

Deciphering cell states and genealogies of human haematopoiesis.

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Chemically defined cytokine-free expansion of human haematopoietic stem cells.

Haematopoietic stem cells (HSCs) are a rare cell type that reconstitute the entire blood and immune systems after transplantation and can be used as a curative cell therapy for a variety of haematological diseases1,2. However, the low number of HSCs in the body makes both biological analyses and clinical application difficult, and the limited extent to which human HSCs can be expanded ex vivo remains a substantial barrier to the wider and safer therapeutic use of HSC transplantation3. Although various reagents have been tested in attempts to stimulate the expansion of human HSCs, cytokines have long been thought to be essential for supporting HSCs ex vivo4. Here we report the establishment of a culture system that allows the long-term ex vivo expansion of human HSCs, achieved through the complete replacement of exogenous cytokines and albumin with chemical agonists and a caprolactam-based polymer. A phosphoinositide 3-kinase activator, in combination with a thrombopoietin-receptor agonist and the pyrimidoindole derivative UM171, were sufficient to stimulate the expansion of umbilical cord blood HSCs that are capable of serial engraftment in xenotransplantation assays. Ex vivo HSC expansion was further supported by split-clone transplantation assays and single-cell RNA-sequencing analysis. Our chemically defined expansion culture system will help to advance clinical HSC therapies.

Long-primed germinal centres with enduring affinity maturation and clonal migration.

Germinal centres are the engines of antibody evolution. Here, using human immunodeficiency virus (HIV) Env protein immunogen priming in rhesus monkeys followed by a long period without further immunization, we demonstrate germinal centre B (BGC) cells that last for at least 6 months. A 186-fold increase in BGC cells was present by week 10 compared with conventional immunization. Single-cell transcriptional profiling showed that both light- and dark-zone germinal centre states were sustained. Antibody somatic hypermutation of BGC cells continued to accumulate throughout the 29-week priming period, with evidence of selective pressure. Env-binding BGC cells were still 49-fold above baseline at 29 weeks, which suggests that they could remain active for even longer periods of time. High titres of HIV-neutralizing antibodies were generated after a single booster immunization. Fully glycosylated HIV trimer protein is a complex antigen, posing considerable immunodominance challenges for B cells1,2. Memory B cells generated under these long priming conditions had higher levels of antibody somatic hypermutation, and both memory B cells and antibodies were more likely to recognize non-immunodominant epitopes. Numerous BGC cell lineage phylogenies spanning more than the 6-month germinal centre period were identified, demonstrating continuous germinal centre activity and selection for at least 191 days with no further antigen exposure. A long-prime, slow-delivery (12 days) immunization approach holds promise for difficult vaccine targets and suggests that patience can have great value for tuning of germinal centres to maximize antibody responses.

Clonal dynamics of haematopoiesis across the human lifespan.

Age-related change in human haematopoiesis causes reduced regenerative capacity1, cytopenias2, immune dysfunction3 and increased risk of blood cancer4-6, but the reason for such abrupt functional decline after 70 years of age remains unclear. Here we sequenced 3,579 genomes from single cell-derived colonies of haematopoietic cells across 10 human subjects from 0 to 81 years of age. Haematopoietic stem cells or multipotent progenitors (HSC/MPPs) accumulated a mean of 17 mutations per year after birth and lost 30 base pairs per year of telomere length. Haematopoiesis in adults less than 65 years of age was massively polyclonal, with high clonal diversity and a stable population of 20,000-200,000 HSC/MPPs contributing evenly to blood production. By contrast, haematopoiesis in individuals aged over 75 showed profoundly decreased clonal diversity. In each of the older subjects, 30-60% of haematopoiesis was accounted for by 12-18 independent clones, each contributing 1-34% of blood production. Most clones had begun their expansion before the subject was 40 years old, but only 22% had known driver mutations. Genome-wide selection analysis estimated that between 1 in 34 and 1 in 12 non-synonymous mutations were drivers, accruing at constant rates throughout life, affecting more genes than identified in blood cancers. Loss of the Y chromosome conferred selective benefits in males. Simulations of haematopoiesis, with constant stem cell population size and constant acquisition of driver mutations conferring moderate fitness benefits, entirely explained the abrupt change in clonal structure in the elderly. Rapidly decreasing clonal diversity is a universal feature of haematopoiesis in aged humans, underpinned by pervasive positive selection acting on many more genes than currently identified.

The longitudinal dynamics and natural history of clonal haematopoiesis.

Clonal expansions driven by somatic mutations become pervasive across human tissues with age, including in the haematopoietic system, where the phenomenon is termed clonal haematopoiesis1-4. The understanding of how and when clonal haematopoiesis develops, the factors that govern its behaviour, how it interacts with ageing and how these variables relate to malignant progression remains limited5,6. Here we track 697 clonal haematopoiesis clones from 385 individuals 55 years of age or older over a median of 13 years. We find that 92.4% of clones expanded at a stable exponential rate over the study period, with different mutations driving substantially different growth rates, ranging from 5% (DNMT3A and TP53) to more than 50% per year (SRSF2P95H). Growth rates of clones with the same mutation differed by approximately ±5% per year, proportionately affecting slow drivers more substantially. By combining our time-series data with phylogenetic analysis of 1,731 whole-genome sequences of haematopoietic colonies from 7 individuals from an older age group, we reveal distinct patterns of lifelong clonal behaviour. DNMT3A-mutant clones preferentially expanded early in life and displayed slower growth in old age, in the context of an increasingly competitive oligoclonal landscape. By contrast, splicing gene mutations drove expansion only later in life, whereas TET2-mutant clones emerged across all ages. Finally, we show that mutations driving faster clonal growth carry a higher risk of malignant progression. Our findings characterize the lifelong natural history of clonal haematopoiesis and give fundamental insights into the interactions between somatic mutation, ageing and clonal selection.

Multicolor strategies for investigating clonal expansion and tissue plasticity.

Understanding the generation of complexity in living organisms requires the use of lineage tracing tools at a multicellular scale. In this review, we describe the different multicolor strategies focusing on mouse models expressing several fluorescent reporter proteins, generated by classical (MADM, Brainbow and its multiple derivatives) or acute (StarTrack, CLoNe, MAGIC Markers, iOn, viral vectors) transgenesis. After detailing the multi-reporter genetic strategies that serve as a basis for the establishment of these multicolor mouse models, we briefly mention other animal and cellular models (zebrafish, chicken, drosophila, iPSC) that also rely on these constructs. Then, we highlight practical applications of multicolor mouse models to better understand organogenesis at single progenitor scale (clonal analyses) in the brain and briefly in several other tissues (intestine, skin, vascular, hematopoietic and immune systems). In addition, we detail the critical contribution of multicolor fate mapping strategies in apprehending the fine cellular choreography underlying tissue morphogenesis in several models with a particular focus on brain cytoarchitecture in health and diseases. Finally, we present the latest technological advances in multichannel and in-depth imaging, and automated analyses that enable to better exploit the large amount of data generated from multicolored tissues.

Signature of long-lived memory CD8+ T cells in acute SARS-CoV-2 infection.

Immunological memory is a hallmark of adaptive immunity and facilitates an accelerated and enhanced immune response upon re-infection with the same pathogen1,2. Since the outbreak of the ongoing COVID-19 pandemic, a key question has focused on which SARS-CoV-2-specific T cells stimulated during acute infection give rise to long-lived memory T cells3. Here, using spectral flow cytometry combined with cellular indexing of transcriptomes and T cell receptor sequencing, we longitudinally characterized individual SARS-CoV-2-specific CD8+ T cells of patients with COVID-19 from acute infection to 1 year into recovery and found a distinct signature identifying long-lived memory CD8+ T cells. SARS-CoV-2-specific memory CD8+ T cells persisting 1 year after acute infection express CD45RA, IL-7 receptor-α and T cell factor 1, but they maintain low expression of CCR7, thus resembling CD45RA+ effector memory T cells. Tracking individual clones of SARS-CoV-2-specific CD8+ T cells, we reveal that an interferon signature marks clones that give rise to long-lived cells, whereas prolonged proliferation and mechanistic target of rapamycin signalling are associated with clonal disappearance from the blood. Collectively, we describe a transcriptional signature that marks long-lived, circulating human memory CD8+ T cells following an acute viral infection.

Butyrate and Metformin Affect Energy Metabolism Independently of the Metabolic Phenotype in the Tumor Therapy Model.

The BALB/c cell transformation assay (BALB-CTA) considers inter- and intra-tumor heterogeneities and affords the possibility of a direct comparison between untransformed and malignant cells. In the present study, we established monoclonal cell lines that originate from the BALB-CTA and mimic heterogeneous tumor cell populations, in order to investigate phenotype-specific effects of the anti-diabetic drug metformin and the short-chain fatty acid butyrate. Growth inhibitory effects were measured with a ViCell XR cell counter. The BALB/c tumor therapy model (BALB-TTM) was performed, and the extracellular glucose level was measured in the medium supernatant. Using a Seahorse Analyzer, the metabolic phenotypes of four selected clones were characterized, and effects on energy metabolism were investigated. Anti-carcinogenic effects and reduced glucose uptake after butyrate application were observed in the BALB-TTM. Metabolic characterization of the cell clones revealed three different phenotypes. Surprisingly, treatment with metformin or butyrate induced opposite metabolic shifts with similar patterns in all cell clones tested. In conclusion, the BALB-TTM is a relevant model for mechanistic cancer research, and the generation of monoclonal cell lines offers a novel possibility to investigate specific drug effects in a heterogeneous tumor cell population. The results indicate that induced alterations in energy metabolism seem to be independent of the original metabolic phenotype.

Mutant clones in normal epithelium outcompete and eliminate emerging tumours.

Human epithelial tissues accumulate cancer-driver mutations with age1-9, yet tumour formation remains rare. The positive selection of these mutations suggests that they alter the behaviour and fitness of proliferating cells10-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 neighbours11-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.

Differential Marker Expression between Keratinocyte Stem Cells and Their Progeny Generated from a Single Colony.

The stemness in keratinocyte stem cells (KSCs) is determined by their gene expression patterns. KSCs are crucial in maintaining epidermal homeostasis and wound repair and are widely used candidates for therapeutic applications. Although several studies have reported their positive identifiers, unique biomarkers for KSCs remain elusive. Here, we aim to identify potential candidate stem cell markers. Human epidermal keratinocytes (HEKs) from neonatal foreskin tissues were isolated and cultured. Single-cell clonal analysis identified and characterized three types of cells: KSCs (holoclones), transient amplifying cells (TACs; meroclones), and differentiated cells (DSCs; paraclones). The clonogenic potential of KSCs demonstrated the highest proliferation potential of KSCs, followed by TACs and DSCs, respectively. Whole-transcriptome analysis using microarray technology unraveled the molecular signatures of these cells. These results were validated by quantitative real-time polymerase chain reaction and flow cytometry analysis. A total of 301 signature upregulated and 149 downregulated differentially expressed genes (DEGs) were identified in the KSCs, compared to TACs and DSCs. Furthermore, DEG analyses revealed new sets of genes related to cell proliferation, cell adhesion, surface makers, and regulatory factors. In conclusion, this study provides a useful source of information for the identification of potential SC-specific candidate markers.

A body map of somatic mutagenesis in morphologically normal human tissues.

Somatic mutations that accumulate in normal tissues are associated with ageing and disease1,2. Here we performed a comprehensive genomic analysis of 1,737 morphologically normal tissue biopsies of 9 organs from 5 donors. We found that somatic mutation accumulations and clonal expansions were widespread, although to variable extents, in morphologically normal human tissues. Somatic copy number alterations were rarely detected, except for in tissues from the oesophagus and cardia. Endogenous mutational processes with the SBS1 and SBS5 mutational signatures are ubiquitous among normal tissues, although they exhibit different relative activities. Exogenous mutational processes operate in multiple tissues from the same donor. We reconstructed the spatial somatic clonal architecture with sub-millimetre resolution. In the oesophagus and cardia, macroscopic somatic clones that expanded to hundreds of micrometres were frequently seen, whereas in tissues such as the colon, rectum and duodenum, somatic clones were microscopic in size and evolved independently, possibly restricted by local tissue microstructures. Our study depicts a body map of somatic mutations and clonal expansions from the same individual.

Genome Editing in Human Induced Pluripotent Stem Cells (hiPSCs).

Cardiomyocytes differentiated from human induced pluripotent stem cells (hiPSCs) are powerful tools for elucidating the pathology behind inherited cardiomyopathies. Genome editing technologies enable targeted genome replacement and the generation of isogenic hiPSCs, allowing investigators to precisely determine the roles of identified mutations. Here, we describe a protocol to obtain isogenic hiPSCs with the corrected allele via homology-directed repair (HDR) using CRISPR/Cas9 genome editing under feeder-free conditions. Seeding hiPSCs in a 24-well plate and conducting the initial evaluation using direct genomic sequencing after 1 week is cost- and time-effective. Following optimization of the protocol, sequence confirmation of the corrected HDR clone is completed within 21 days.

High-risk clones of extended-spectrum ß-lactamase-producing Klebsiella pneumoniae isolated from the University Hospital Establishment of Oran, Algeria (2011-2012).

The purpose of the study was to characterize the resistome, virulome, mobilome and Clustered Regularly Interspaced Short Palindromic Repeats-associated (CRISPR-Cas) system of extended-spectrum ß-lactamase producing Klebsiella pneumoniae (ESBL-KP) clinical isolates and to determine their phylogenetic relatedness. The isolates were from Algeria, isolated at the University Hospital Establishment of Oran, between 2011 and 2012. ESBL-KP isolates (n = 193) were screened for several antibiotic resistance genes (ARGs) using qPCR followed by Pulsed-Field Gel Electrophoresis (PFGE). Representative isolates were selected from PFGE clusters and subjected to whole-genome sequencing (WGS). Genomic characterization of the WGS data by studying prophages, CRISPR-Cas systems, Multi-Locus Sequence Typing (MLST), serotype, ARGs, virulence genes, plasmid replicons, and their pMLST. Phylogenetic and comparative genomic were done using core genome MLST and SNP-Based analysis. Generally, the ESBL-KP isolates were polyclonal. The whole genome sequences of nineteen isolates were taken of main PFGE clusters. Sixteen sequence types (ST) were found including high-risk clones ST14, ST23, ST37, and ST147. Serotypes K1 (n = 1), K2 (n = 2), K3 (n = 1), K31 (n = 1), K62 (n = 1), and K151 (n = 1) are associated with hyper-virulence. CRISPR-Cas system was found in 47.4%, typed I-E and I-E*. About ARGs, from 193 ESBL-KP, the majority of strains were multidrug-resistant, the CTX-M-1 enzyme was predominant (99%) and the prevalence of plasmid-mediated quinolone resistance (PMQR) genes was high with aac(6')-lb-cr (72.5%) and qnr's (65.8%). From 19 sequenced isolates we identified ESBL, AmpC, and carbapenemase genes: blaCTX-M-15 (n = 19), blaOXA-48 (n = 1), blaCMY-2 (n = 2), and blaCMY-16 (n = 2), as well as non-ESBL genes: qnrB1 (n = 12), qnrS1 (n = 1) and armA (n = 2). We found IncF, IncN, IncL/M, IncA/C2, and Col replicon types, at least once per isolate. This study is the first to report qnrS in ESBL-KP in Algeria. Our analysis shows the concerning co-existence of virulence and resistance genes and would support that genomic surveillance should be a high priority in the hospital environment.

SARS-CoV-2 mRNA vaccines induce persistent human germinal centre responses.

SARS-CoV-2 mRNA-based vaccines are about 95% effective in preventing COVID-191-5. The dynamics of antibody-secreting plasmablasts and germinal centre B cells induced by these vaccines in humans remain unclear. Here we examined antigen-specific B cell responses in peripheral blood (n = 41) and draining lymph nodes in 14 individuals who had received 2 doses of BNT162b2, an mRNA-based vaccine that encodes the full-length SARS-CoV-2 spike (S) gene1. Circulating IgG- and IgA-secreting plasmablasts that target the S protein peaked one week after the second immunization and then declined, becoming undetectable three weeks later. These plasmablast responses preceded maximal levels of serum anti-S binding and neutralizing antibodies to an early circulating SARS-CoV-2 strain as well as emerging variants, especially in individuals who had previously been infected with SARS-CoV-2 (who produced the most robust serological responses). By examining fine needle aspirates of draining axillary lymph nodes, we identified germinal centre B cells that bound S protein in all participants who were sampled after primary immunization. High frequencies of S-binding germinal centre B cells and plasmablasts were sustained in these draining lymph nodes for at least 12 weeks after the booster immunization. S-binding monoclonal antibodies derived from germinal centre B cells predominantly targeted the receptor-binding domain of the S protein, and fewer clones bound to the N-terminal domain or to epitopes shared with the S proteins of the human betacoronaviruses OC43 and HKU1. These latter cross-reactive B cell clones had higher levels of somatic hypermutation as compared to those that recognized only the SARS-CoV-2 S protein, which suggests a memory B cell origin. Our studies demonstrate that SARS-CoV-2 mRNA-based vaccination of humans induces a persistent germinal centre B cell response, which enables the generation of robust humoral immunity.

Lineage tracing of human development through somatic mutations.

The ontogeny of the human haematopoietic system during fetal development has previously been characterized mainly through careful microscopic observations1. Here we reconstruct a phylogenetic tree of blood development using whole-genome sequencing of 511 single-cell-derived haematopoietic colonies from healthy human fetuses at 8 and 18 weeks after conception, coupled with deep targeted sequencing of tissues of known embryonic origin. We found that, in healthy fetuses, individual haematopoietic progenitors acquire tens of somatic mutations by 18 weeks after conception. We used these mutations as barcodes and timed the divergence of embryonic and extra-embryonic tissues during development, and estimated the number of blood antecedents at different stages of embryonic development. Our data support a hypoblast origin of the extra-embryonic mesoderm and primitive blood in humans.