Engineered mammalian cells are key for biotechnology by enabling broad applications ranging from in vitro model systems to therapeutic biofactories. Engineered cell lines exist as a population containing sub-lineages of cell clones that exhibit substantial genetic and phenotypic heterogeneity. There is still a limited understanding of the source of this inter-clonal heterogeneity as well as its implications for biotechnological applications. Here, we developed a genomic barcoding strategy for a targeted integration (TI)-based CHO antibody producer cell line development process. This technology provided novel insights about clone diversity during stable cell line selection on pool level, enabled an imaging-independent monoclonality assessment after single cell cloning, and eventually improved hit-picking of antibody producer clones by monitoring of cellular lineages during the cell line development (CLD) process. Specifically, we observed that CHO producer pools generated by TI of two plasmids at a single genomic site displayed a low diversity (< 0.1% RMCE efficiency), which further depends on the expressed molecules, and underwent rapid population skewing towards dominant clones during routine cultivation. Clonal cell lines from one individual TI event demonstrated a significantly lower variance regarding production-relevant and phenotypic parameters as compared to cell lines from distinct TI events. This implies that the observed cellular diversity lies within pre-existing cell-intrinsic factors and that the majority of clonal variation did not develop during the CLD process, especially during single cell cloning. Using cellular barcodes as a proxy for cellular diversity, we improved our CLD screening workflow and enriched diversity of production-relevant parameters substantially. This work, by enabling clonal diversity monitoring and control, paves the way for an economically valuable and data-driven CLD process.

Recurrent exposures to a pathogenic antigen remodel the CD8+ T cell compartment and generate a functional memory repertoire that is polyclonal and complex. At the clonotype level, the response to the conserved influenza antigen, M158-66 has been well characterized in healthy individuals, but not in patients receiving immunosuppressive therapy or with aberrant immunity, such as those with juvenile idiopathic arthritis (JIA). Here we show that patients with JIA have a reduced number of M158-66 specific RS/RA clonotypes, indicating decreased clonal richness and, as a result, have lower repertoire diversity. By using a rank-frequency approach to analyze the distribution of the repertoire, we found several characteristics of the JIA T cell repertoire to be akin to repertoires seen in healthy adults, including an amplified RS/RA-specific antigen response, representing greater clonal unevenness. Unlike mature repertoires, however, there is more fluctuation in clonotype distribution, less clonotype stability, and more variable IFNy response of the M158-66 specific RS/RA clonotypes in JIA. This indicates that functional clonal expansion is altered in patients with JIA on immunosuppressive therapies. We propose that the response to the influenza M158-66 epitope described here is a general phenomenon for JIA patients receiving immunosuppressive therapy, and that the changes in clonal richness and unevenness indicate a retarded and uneven generation of a mature immune response.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Mutations are abundantly present in tissues of healthy individuals, including the breast epithelium. Yet it remains unknown whether mutant cells directly induce lesion formation or first spread, leading to a field of mutant cells that is predisposed towards lesion formation. To study the clonal and spatial relationships between morphologically normal breast epithelium adjacent to pre-cancerous lesions, we developed a three-dimensional (3D) imaging pipeline combined with spatially resolved genomics on archival, formalin-fixed breast tissue with the non-obligate breast cancer precursor ductal carcinoma in situ (DCIS). Using this 3D image-guided characterization method, we built high-resolution spatial maps of DNA copy number aberration (CNA) profiles within the DCIS lesion and the surrounding normal mammary ducts. We show that the local heterogeneity within a DCIS lesion is limited. However, by mapping the CNA profiles back onto the 3D reconstructed ductal subtree, we find that in eight out of 16 cases the healthy epithelium adjacent to the DCIS lesions has overlapping structural variations with the CNA profile of the DCIS. Together, our study indicates that pre-malignant breast transformations frequently develop within mutant clonal fields of morphologically normal-looking ducts. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

In somatic tissue differentiation, chromatin accessibility changes govern priming and precursor commitment towards cellular fates1-3. Therefore, somatic mutations are likely to alter chromatin accessibility patterns, as they disrupt differentiation topologies leading to abnormal clonal outgrowth. However, defining the impact of somatic mutations on the epigenome in human samples is challenging due to admixed mutated and wild-type cells. Here, to chart how somatic mutations disrupt epigenetic landscapes in human clonal outgrowths, we developed genotyping of targeted loci with single-cell chromatin accessibility (GoT-ChA). This high-throughput platform links genotypes to chromatin accessibility at single-cell resolution across thousands of cells within a single assay. We applied GoT-ChA to CD34+ cells from patients with myeloproliferative neoplasms with JAK2V617F-mutated haematopoiesis. Differential accessibility analysis between wild-type and JAK2V617F-mutant progenitors revealed both cell-intrinsic and cell-state-specific shifts within mutant haematopoietic precursors, including cell-intrinsic pro-inflammatory signatures in haematopoietic stem cells, and a distinct profibrotic inflammatory chromatin landscape in megakaryocytic progenitors. Integration of mitochondrial genome profiling and cell-surface protein expression measurement allowed expansion of genotyping onto DOGMA-seq through imputation, enabling single-cell capture of genotypes, chromatin accessibility, RNA expression and cell-surface protein expression. Collectively, we show that the JAK2V617F mutation leads to epigenetic rewiring in a cell-intrinsic and cell type-specific manner, influencing inflammation states and differentiation trajectories. We envision that GoT-ChA will empower broad future investigations of the critical link between somatic mutations and epigenetic alterations across clonal populations in malignant and non-malignant contexts.

T cells are required for protective immunity against Mycobacterium tuberculosis. We recently described a cohort of Ugandan household contacts of tuberculosis cases who appear to \"resist\" M. tuberculosis infection (resisters; RSTRs) and showed that these individuals harbor IFN-γ-independent T cell responses to M. tuberculosis-specific peptide antigens. However, T cells also recognize nonprotein antigens via antigen-presenting systems that are independent of genetic background, known as donor-unrestricted T cells (DURTs). We used tetramer staining and flow cytometry to characterize the association between DURTs and \"resistance\" to M. tuberculosis infection. Peripheral blood frequencies of most DURT subsets were comparable between RSTRs and latently infected controls (LTBIs). However, we observed a 1.65-fold increase in frequency of MR1-restricted T (MR1T) cells among RSTRs in comparison with LTBIs. Single-cell RNA sequencing of 18,251 MR1T cells sorted from 8 donors revealed 5,150 clonotypes that expressed a common transcriptional program, the majority of which were private. Sequencing of the T cell receptor α/T cell receptor δ (TCRα/δ) repertoire revealed several DURT clonotypes were expanded among RSTRs, including 2 MR1T clonotypes that recognized mycobacteria-infected cells in a TCR-dependent manner. Overall, our data reveal unexpected donor-specific diversity in the TCR repertoire of human MR1T cells as well as associations between mycobacteria-reactive MR1T clonotypes and resistance to M. tuberculosis infection.

During brain development, neural progenitors expand through symmetric divisions before giving rise to differentiating cell types via asymmetric divisions. Transition between those modes varies among individual neural stem cells, resulting in clones of different sizes. Imaging-based lineage tracing allows for lineage analysis at high cellular resolution but systematic approaches to analyse clonal behaviour of entire tissues are currently lacking. Here we implement whole-tissue lineage tracing by genomic DNA barcoding in 3D human cerebral organoids, to show that individual stem cell clones produce progeny on a vastly variable scale. By using stochastic modelling we find that variable lineage sizes arise because a subpopulation of lineages retains symmetrically dividing cells. We show that lineage sizes can adjust to tissue demands after growth perturbation via chemical ablation or genetic restriction of a subset of cells in chimeric organoids. Our data suggest that adaptive plasticity of stem cell populations ensures robustness of development in human brain organoids.

Insulin is a central autoantigen in the pathogenesis of T1D, and thymic epithelial cell expression of insulin under the control of the Autoimmune Regulator (Aire) is thought to be a key component of maintaining tolerance to insulin. In spite of this general working model, direct detection of this thymic selection on insulin-specific T cells has been somewhat elusive. Here, we used a combination of highly sensitive T cell receptor transgenic models for detecting thymic selection and sorting and sequencing of Insulin-specific CD4+ T cells from Aire-deficient mice as a strategy to further define their selection. This analysis revealed a number of unique t cell receptor (TCR) clones in Aire-deficient hosts with high affinity for insulin/major histocompatibility complex (MHC) ligands. We then modeled the thymic selection of one of these clones in Aire-deficient versus wild-type hosts and found that this model clone could escape thymic negative selection in the absence of thymic Aire. Together, these results suggest that thymic expression of insulin plays a key role in trimming and removing high-affinity insulin-specific T cells from the repertoire to help promote tolerance.

B cells and T cells are important components of the adaptive immune system and mediate anticancer immunity. The T cell landscape in cancer is well characterized, but the contribution of B cells to anticancer immunosurveillance is less well explored. Here we show an integrative analysis of the B cell and T cell receptor repertoire from individuals with metastatic breast cancer and individuals with early breast cancer during neoadjuvant therapy. Using immune receptor, RNA and whole-exome sequencing, we show that both B cell and T cell responses seem to coevolve with the metastatic cancer genomes and mirror tumor mutational and neoantigen architecture. B cell clones associated with metastatic immunosurveillance and temporal persistence were more expanded and distinct from site-specific clones. B cell clonal immunosurveillance and temporal persistence are predictable from the clonal structure, with higher-centrality B cell antigen receptors more likely to be detected across multiple metastases or across time. This predictability was generalizable across other immune-mediated disorders. This work lays a foundation for prioritizing antibody sequences for therapeutic targeting in cancer.

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.