Imprinted SARS-CoV-2-specific memory lymphocytes define hybrid immunity.

Immune memory is tailored by cues that lymphocytes perceive during priming. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic created a situation in which nascent memory could be tracked through additional antigen exposures. Both SARS-CoV-2 infection and vaccination induce multifaceted, functional immune memory, but together, they engender improved protection from disease, termed hybrid immunity. We therefore investigated how vaccine-induced memory is shaped by previous infection. We found that following vaccination, previously infected individuals generated more SARS-CoV-2 RBD-specific memory B cells and variant-neutralizing antibodies and a distinct population of IFN-γ and IL-10-expressing memory SARS-CoV-2 spike-specific CD4+ T cells than previously naive individuals. Although additional vaccination could increase humoral memory in previously naive individuals, it did not recapitulate the distinct CD4+ T cell cytokine profile observed in previously infected subjects. Thus, imprinted features of SARS-CoV-2-specific memory lymphocytes define hybrid immunity.

CD8+ T cells specific for an immunodominant SARS-CoV-2 nucleocapsid epitope cross-react with selective seasonal coronaviruses.

Efforts are being made worldwide to understand the immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for the coronavirus disease 2019 (COVID-19) pandemic, including the impact of T cell immunity and cross-recognition with seasonal coronaviruses. Screening of SARS-CoV-2 peptide pools revealed that the nucleocapsid (N) protein induced an immunodominant response in HLA-B7+ COVID-19-recovered individuals that was also detectable in unexposed donors. A single N-encoded epitope that was highly conserved across circulating coronaviruses drove this immunodominant response. In vitro peptide stimulation and crystal structure analyses revealed T cell-mediated cross-reactivity toward circulating OC43 and HKU-1 betacoronaviruses but not 229E or NL63 alphacoronaviruses because of different peptide conformations. T cell receptor (TCR) sequencing indicated that cross-reactivity was driven by private TCR repertoires with a bias for TRBV27 and a long CDR3ß loop. Our findings demonstrate the basis of selective T cell cross-reactivity for an immunodominant SARS-CoV-2 epitope and its homologs from seasonal coronaviruses, suggesting long-lasting protective immunity.

Consistent patterns of common species across tropical tree communities.

Trees structure the Earth's most biodiverse ecosystem, tropical forests. The vast number of tree species presents a formidable challenge to understanding these forests, including their response to environmental change, as very little is known about most tropical tree species. A focus on the common species may circumvent this challenge. Here we investigate abundance patterns of common tree species using inventory data on 1,003,805 trees with trunk diameters of at least 10 cm across 1,568 locations1-6 in closed-canopy, structurally intact old-growth tropical forests in Africa, Amazonia and Southeast Asia. We estimate that 2.2%, 2.2% and 2.3% of species comprise 50% of the tropical trees in these regions, respectively. Extrapolating across all closed-canopy tropical forests, we estimate that just 1,053 species comprise half of Earth's 800 billion tropical trees with trunk diameters of at least 10 cm. Despite differing biogeographic, climatic and anthropogenic histories7, we find notably consistent patterns of common species and species abundance distributions across the continents. This suggests that fundamental mechanisms of tree community assembly may apply to all tropical forests. Resampling analyses show that the most common species are likely to belong to a manageable list of known species, enabling targeted efforts to understand their ecology. Although they do not detract from the importance of rare species, our results open new opportunities to understand the world's most diverse forests, including modelling their response to environmental change, by focusing on the common species that constitute the majority of their trees.

Ubiquitination regulates ER-phagy and remodelling of endoplasmic reticulum.

The endoplasmic reticulum (ER) undergoes continuous remodelling via a selective autophagy pathway, known as ER-phagy1. ER-phagy receptors have a central role in this process2, but the regulatory mechanism remains largely unknown. Here we report that ubiquitination of the ER-phagy receptor FAM134B within its reticulon homology domain (RHD) promotes receptor clustering and binding to lipidated LC3B, thereby stimulating ER-phagy. Molecular dynamics (MD) simulations showed how ubiquitination perturbs the RHD structure in model bilayers and enhances membrane curvature induction. Ubiquitin molecules on RHDs mediate interactions between neighbouring RHDs to form dense receptor clusters that facilitate the large-scale remodelling of lipid bilayers. Membrane remodelling was reconstituted in vitro with liposomes and ubiquitinated FAM134B. Using super-resolution microscopy, we discovered FAM134B nanoclusters and microclusters in cells. Quantitative image analysis revealed a ubiquitin-mediated increase in FAM134B oligomerization and cluster size. We found that the E3 ligase AMFR, within multimeric ER-phagy receptor clusters, catalyses FAM134B ubiquitination and regulates the dynamic flux of ER-phagy. Our results show that ubiquitination enhances RHD functions via receptor clustering, facilitates ER-phagy and controls ER remodelling in response to cellular demands.

Heteromeric clusters of ubiquitinated ER-shaping proteins drive ER-phagy.

Membrane-shaping proteins characterized by reticulon homology domains play an important part in the dynamic remodelling of the endoplasmic reticulum (ER). An example of such a protein is FAM134B, which can bind LC3 proteins and mediate the degradation of ER sheets through selective autophagy (ER-phagy)1. Mutations in FAM134B result in a neurodegenerative disorder in humans that mainly affects sensory and autonomic neurons2. Here we report that ARL6IP1, another ER-shaping protein that contains a reticulon homology domain and is associated with sensory loss3, interacts with FAM134B and participates in the formation of heteromeric multi-protein clusters required for ER-phagy. Moreover, ubiquitination of ARL6IP1 promotes this process. Accordingly, disruption of Arl6ip1 in mice causes an expansion of ER sheets in sensory neurons that degenerate over time. Primary cells obtained from Arl6ip1-deficient mice or from patients display incomplete budding of ER membranes and severe impairment of ER-phagy flux. Therefore, we propose that the clustering of ubiquitinated ER-shaping proteins facilitates the dynamic remodelling of the ER during ER-phagy and is important for neuronal maintenance.

How the conception of control influences our understanding of actions.

Wilful movement requires neural control. Commonly, neural computations are thought to generate motor commands that bring the musculoskeletal system - that is, the plant - from its current physical state into a desired physical state. The current state can be estimated from past motor commands and from sensory information. Modelling movement on the basis of this concept of plant control strives to explain behaviour by identifying the computational principles for control signals that can reproduce the observed features of movements. From an alternative perspective, movements emerge in a dynamically coupled agent-environment system from the pursuit of subjective perceptual goals. Modelling movement on the basis of this concept of perceptual control aims to identify the controlled percepts and their coupling rules that can give rise to the observed characteristics of behaviour. In this Perspective, we discuss a broad spectrum of approaches to modelling human motor control and their notions of control signals, internal models, handling of sensory feedback delays and learning. We focus on the influence that the plant control and the perceptual control perspective may have on decisions when modelling empirical data, which may in turn shape our understanding of actions.

Cobalt-electrocatalytic HAT for functionalization of unsaturated C-C bonds.

The study and application of transition metal hydrides (TMHs) has been an active area of chemical research since the early 1960s1, for energy storage, through the reduction of protons to generate hydrogen2,3, and for organic synthesis, for the functionalization of unsaturated C-C, C-O and C-N bonds4,5. In the former instance, electrochemical means for driving such reactivity has been common place since the 1950s6 but the use of stoichiometric exogenous organic- and metal-based reductants to harness the power of TMHs in synthetic chemistry remains the norm. In particular, cobalt-based TMHs have found widespread use for the derivatization of olefins and alkynes in complex molecule construction, often by a net hydrogen atom transfer (HAT)7. Here we show how an electrocatalytic approach inspired by decades of energy storage research can be made use of in the context of modern organic synthesis. This strategy not only offers benefits in terms of sustainability and efficiency but also enables enhanced chemoselectivity and distinct, tunable reactivity. Ten different reaction manifolds across dozens of substrates are exemplified, along with detailed mechanistic insights into this scalable electrochemical entry into Co-H generation that takes place through a low-valent intermediate.

Inner ear biomechanics reveals a Late Triassic origin for mammalian endothermy.

Endothermy underpins the ecological dominance of mammals and birds in diverse environmental settings1,2. However, it is unclear when this crucial feature emerged during mammalian evolutionary history, as most of the fossil evidence is ambiguous3-17. Here we show that this key evolutionary transition can be investigated using the morphology of the endolymph-filled semicircular ducts of the inner ear, which monitor head rotations and are essential for motor coordination, navigation and spatial awareness18-22. Increased body temperatures during the ectotherm-endotherm transition of mammal ancestors would decrease endolymph viscosity, negatively affecting semicircular duct biomechanics23,24, while simultaneously increasing behavioural activity25,26 probably required improved performance27. Morphological changes to the membranous ducts and enclosing bony canals would have been necessary to maintain optimal functionality during this transition. To track these morphofunctional changes in 56 extinct synapsid species, we developed the thermo-motility index, a proxy based on bony canal morphology. The results suggest that endothermy evolved abruptly during the Late Triassic period in Mammaliamorpha, correlated with a sharp increase in body temperature (5-9 °C) and an expansion of aerobic and anaerobic capacities. Contrary to previous suggestions3-14, all stem mammaliamorphs were most probably ectotherms. Endothermy, as a crucial physiological characteristic, joins other distinctive mammalian features that arose during this period of climatic instability28.

Ancient DNA and deep population structure in sub-Saharan African foragers.

Multiple lines of genetic and archaeological evidence suggest that there were major demographic changes in the terminal Late Pleistocene epoch and early Holocene epoch of sub-Saharan Africa1-4. Inferences about this period are challenging to make because demographic shifts in the past 5,000 years have obscured the structures of more ancient populations3,5. Here we present genome-wide ancient DNA data for six individuals from eastern and south-central Africa spanning the past approximately 18,000 years (doubling the time depth of sub-Saharan African ancient DNA), increase the data quality for 15 previously published ancient individuals and analyse these alongside data from 13 other published ancient individuals. The ancestry of the individuals in our study area can be modelled as a geographically structured mixture of three highly divergent source populations, probably reflecting Pleistocene interactions around 80-20 thousand years ago, including deeply diverged eastern and southern African lineages, plus a previously unappreciated ubiquitous distribution of ancestry that occurs in highest proportion today in central African rainforest hunter-gatherers. Once established, this structure remained highly stable, with limited long-range gene flow. These results provide a new line of genetic evidence in support of hypotheses that have emerged from archaeological analyses but remain contested, suggesting increasing regionalization at the end of the Pleistocene epoch.

Towards the understanding of molecular motors and its relationship with local unfolding.

Molecular motors are machines essential for life since they convert chemical energy into mechanical work. However, the precise mechanism by which nucleotide binding, catalysis, or release of products is coupled to the work performed by the molecular motor is still not entirely clear. This is due, in part, to a lack of understanding of the role of force in the mechanical-structural processes involved in enzyme catalysis. From a mechanical perspective, one promising hypothesis is the Haldane-Pauling hypothesis which considers the idea that part of the enzymatic catalysis is strain-induced. It suggests that enzymes cannot be efficient catalysts if they are fully complementary to the substrates. Instead, they must exert strain on the substrate upon binding, using enzyme-substrate energy interaction (binding energy) to accelerate the reaction rate. A novel idea suggests that during catalysis, significant strain energy is built up, which is then released by a local unfolding/refolding event known as 'cracking'. Recent evidence has also shown that in catalytic reactions involving conformational changes, part of the heat released results in a center-of-mass acceleration of the enzyme, raising the possibility that the heat released by the reaction itself could affect the enzyme's integrity. Thus, it has been suggested that this released heat could promote or be linked to the cracking seen in proteins such as adenylate kinase (AK). We propose that the energy released as a consequence of ligand binding/catalysis is associated with the local unfolding/refolding events (cracking), and that this energy is capable of driving the mechanical work.

A phase 3 randomized trial of mavorixafor, a CXCR4 antagonist, for WHIM syndrome.

ABSTRACT: We investigated efficacy and safety of mavorixafor, an oral CXCR4 antagonist, in participants with warts, hypogammaglobulinemia, infections, and myelokathexis (WHIM) syndrome, a rare immunodeficiency caused by CXCR4 gain-of-function variants. This randomized (1:1), double-blind, placebo-controlled, phase 3 trial enrolled participants aged ≥12 years with WHIM syndrome and absolute neutrophil count (ANC) ≤0.4 × 103/µL. Participants received once-daily mavorixafor or placebo for 52 weeks. The primary end point was time (hours) above ANC threshold ≥0.5 × 103/µL (TATANC; over 24 hours). Secondary end points included TAT absolute lymphocyte count ≥1.0 × 103/µL (TATALC; over 24 hours); absolute changes in white blood cell (WBC), ANC, and absolute lymphocyte count (ALC) from baseline; annualized infection rate; infection duration; and total infection score (combined infection number/severity). In 31 participants (mavorixafor, n = 14; placebo, n = 17), mavorixafor least squares (LS) mean TATANC was 15.0 hours and 2.8 hours for placebo (P < .001). Mavorixafor LS mean TATALC was 15.8 hours and 4.6 hours for placebo (P < .001). Annualized infection rates were 60% lower with mavorixafor vs placebo (LS mean 1.7 vs 4.2; nominal P = .007), and total infection scores were 40% lower (7.4 [95% confidence interval [CI], 1.6-13.2] vs 12.3 [95% CI, 7.2-17.3]). Treatment with mavorixafor reduced infection frequency, severity, duration, and antibiotic use. No discontinuations occurred due to treatment-emergent adverse events (TEAEs); no related serious TEAEs were observed. Overall, mavorixafor treatment demonstrated significant increases in LS mean TATANC and TATALC, reduced infection frequency, severity/duration, and was well tolerated. The trial was registered at www.clinicaltrials.gov as #NCT03995108.

Lysine L-lactylation is the dominant lactylation isomer induced by glycolysis.

Lysine L-lactylation (Kl-la) is a novel protein posttranslational modification (PTM) driven by L-lactate. This PTM has three isomers: Kl-la, N-ε-(carboxyethyl)-lysine (Kce) and D-lactyl-lysine (Kd-la), which are often confused in the context of the Warburg effect and nuclear presence. Here we introduce two methods to differentiate these isomers: a chemical derivatization and high-performance liquid chromatography analysis for efficient separation, and isomer-specific antibodies for high-selectivity identification. We demonstrated that Kl-la is the primary lactylation isomer on histones and dynamically regulated by glycolysis, not Kd-la or Kce, which are observed when the glyoxalase system was incomplete. The study also reveals that lactyl-coenzyme A, a precursor in L-lactylation, correlates positively with Kl-la levels. This work not only provides a methodology for distinguishing other PTM isomers, but also highlights Kl-la as the primary responder to glycolysis and the Warburg effect.

Rapid mechanosensitive migration and dispersal of newly divided mesenchymal cells aid their recruitment into dermal condensates.

Embryonic mesenchymal cells are dispersed within an extracellular matrix but can coalesce to form condensates with key developmental roles. Cells within condensates undergo fate and morphological changes and induce cell fate changes in nearby epithelia to produce structures including hair follicles, feathers, or intestinal villi. Here, by imaging mouse and chicken embryonic skin, we find that mesenchymal cells undergo much of their dispersal in early interphase, in a stereotyped process of displacement driven by 3 hours of rapid and persistent migration followed by a long period of low motility. The cell division plane and the elevated migration speed and persistence of newly born mesenchymal cells are mechanosensitive, aligning with tissue tension, and are reliant on active WNT secretion. This behaviour disperses mesenchymal cells and allows daughters of recent divisions to travel long distances to enter dermal condensates, demonstrating an unanticipated effect of cell cycle subphase on core mesenchymal behaviour.

Splicing factor YBX1 mediates persistence of JAK2-mutated neoplasms.

Janus kinases (JAKs) mediate responses to cytokines, hormones and growth factors in haematopoietic cells1,2. The JAK gene JAK2 is frequently mutated in the ageing haematopoietic system3,4 and in haematopoietic cancers5. JAK2 mutations constitutively activate downstream signalling and are drivers of myeloproliferative neoplasm (MPN). In clinical use, JAK inhibitors have mixed effects on the overall disease burden of JAK2-mutated clones6,7, prompting us to investigate the mechanism underlying disease persistence. Here, by in-depth phosphoproteome profiling, we identify proteins involved in mRNA processing as targets of mutant JAK2. We found that inactivation of YBX1, a post-translationally modified target of JAK2, sensitizes cells that persist despite treatment with JAK inhibitors to apoptosis and results in RNA mis-splicing, enrichment for retained introns and disruption of the transcriptional control of extracellular signal-regulated kinase (ERK) signalling. In combination with pharmacological JAK inhibition, YBX1 inactivation induces apoptosis in JAK2-dependent mouse and primary human cells, causing regression of the malignant clones in vivo, and inducing molecular remission. This identifies and validates a cell-intrinsic mechanism whereby differential protein phosphorylation causes splicing-dependent alterations of JAK2-ERK signalling and the maintenance of JAK2V617F malignant clones. Therapeutic targeting of YBX1-dependent ERK signalling in combination with JAK2 inhibition could thus eradicate cells harbouring mutations in JAK2.

Integrating genomic features for non-invasive early lung cancer detection.

Radiologic screening of high-risk adults reduces lung-cancer-related mortality1,2; however, a small minority of eligible individuals undergo such screening in the United States3,4. The availability of blood-based tests could increase screening uptake. Here we introduce improvements to cancer personalized profiling by deep sequencing (CAPP-Seq)5, a method for the analysis of circulating tumour DNA (ctDNA), to better facilitate screening applications. We show that, although levels are very low in early-stage lung cancers, ctDNA is present prior to treatment in most patients and its presence is strongly prognostic. We also find that the majority of somatic mutations in the cell-free DNA (cfDNA) of patients with lung cancer and of risk-matched controls reflect clonal haematopoiesis and are non-recurrent. Compared with tumour-derived mutations, clonal haematopoiesis mutations occur on longer cfDNA fragments and lack mutational signatures that are associated with tobacco smoking. Integrating these findings with other molecular features, we develop and prospectively validate a machine-learning method termed 'lung cancer likelihood in plasma' (Lung-CLiP), which can robustly discriminate early-stage lung cancer patients from risk-matched controls. This approach achieves performance similar to that of tumour-informed ctDNA detection and enables tuning of assay specificity in order to facilitate distinct clinical applications. Our findings establish the potential of cfDNA for lung cancer screening and highlight the importance of risk-matching cases and controls in cfDNA-based screening studies.

Asynchronous carbon sink saturation in African and Amazonian tropical forests.

Structurally intact tropical forests sequestered about half of the global terrestrial carbon uptake over the 1990s and early 2000s, removing about 15 per cent of anthropogenic carbon dioxide emissions1-3. Climate-driven vegetation models typically predict that this tropical forest 'carbon sink' will continue for decades4,5. Here we assess trends in the carbon sink using 244 structurally intact African tropical forests spanning 11 countries, compare them with 321 published plots from Amazonia and investigate the underlying drivers of the trends. The carbon sink in live aboveground biomass in intact African tropical forests has been stable for the three decades to 2015, at 0.66 tonnes of carbon per hectare per year (95 per cent confidence interval 0.53-0.79), in contrast to the long-term decline in Amazonian forests6. Therefore the carbon sink responses of Earth's two largest expanses of tropical forest have diverged. The difference is largely driven by carbon losses from tree mortality, with no detectable multi-decadal trend in Africa and a long-term increase in Amazonia. Both continents show increasing tree growth, consistent with the expected net effect of rising atmospheric carbon dioxide and air temperature7-9. Despite the past stability of the African carbon sink, our most intensively monitored plots suggest a post-2010 increase in carbon losses, delayed compared to Amazonia, indicating asynchronous carbon sink saturation on the two continents. A statistical model including carbon dioxide, temperature, drought and forest dynamics accounts for the observed trends and indicates a long-term future decline in the African sink, whereas the Amazonian sink continues to weaken rapidly. Overall, the uptake of carbon into Earth's intact tropical forests peaked in the 1990s. Given that the global terrestrial carbon sink is increasing in size, independent observations indicating greater recent carbon uptake into the Northern Hemisphere landmass10 reinforce our conclusion that the intact tropical forest carbon sink has already peaked. This saturation and ongoing decline of the tropical forest carbon sink has consequences for policies intended to stabilize Earth's climate.

Valuing improvements in the ecological integrity of local and regional waters using the biological condition gradient.

Scientific knowledge related to quantifying the monetized benefits for landscape-wide water quality improvements does not meet current regulatory and benefit-cost analysis needs in the United States. In this study we addressed this knowledge gap by incorporating the Biological Condition Gradient (BCG) as a water quality metric into a stated preference survey capable of estimating the total economic value (use and nonuse) for aquatic ecosystem improvements. The BCG is grounded in ecological principles and generalizable and transferable across space. Moreover, as the BCG translates available data on biological condition into a score on a 6-point scale, it provides a simple metric that can be readily communicated to the public. We applied our BCG-based survey instrument to households across the Upper Mississippi, Ohio, and Tennessee river basins and report values for a range of potential improvements that vary by location, spatial scale, and the scope of the water quality change. We found that people are willing to pay twice as much for an improvement policy that targets their home watershed (defined as a four-digit hydrologic unit) versus a more distant one. We also found that extending the spatial scale of a local policy beyond the home watershed does not generate additional benefits to the household. Finally, our results suggest that nonuse sources of value (e.g., bequest value, intrinsic aesthetic value) are an important component of overall benefits.

Divergent roles for STAT4 in shaping differentiation of cytotoxic ILC1 and NK cells during gut inflammation.

Natural killer (NK) cells and type 1 innate lymphoid cells (ILC1) require signal transducer and activator of transcription 4 (STAT4) to elicit rapid effector responses and protect against pathogens. By combining genetic and transcriptomic approaches, we uncovered divergent roles for STAT4 in regulating effector differentiation of these functionally related cell types. Stat4 deletion in Ncr1-expressing cells led to impaired NK cell terminal differentiation as well as to an unexpected increased generation of cytotoxic ILC1 during intestinal inflammation. Mechanistically, Stat4-deficient ILC1 exhibited upregulation of gene modules regulated by STAT5 in vivo and an aberrant effector differentiation upon in vitro stimulation with IL-2, used as a prototypical STAT5 activator. Moreover, STAT4 expression in NCR+ innate lymphocytes restrained gut inflammation in the dextran sulfate sodium-induced colitis model limiting pathogenic production of IL-13 from adaptive CD4+ T cells in the large intestine. Collectively, our data shed light on shared and distinctive mechanisms of STAT4-regulated transcriptional control in NK cells and ILC1 required for intestinal inflammatory responses.

Higher HIV-1 evolutionary rate is associated with cytotoxic T lymphocyte escape mutations in infants.

Escape from cytotoxic T lymphocyte (CTL) responses toward HIV-1 Gag and Nef has been associated with reduced control of HIV-1 replication in adults. However, less is known about CTL-driven immune selection in infants as longitudinal studies of infants are limited. Here, 1,210 gag and 1,264 nef sequences longitudinally collected within 15 months after birth from 14 HIV-1 perinatally infected infants and their mothers were analyzed. The number of transmitted founder (T/F) viruses and associations between virus evolution, selection, CTL escape, and disease progression were determined. The analyses indicated that a paraphyletic-monophyletic relationship between the mother-infant sequences was common (80%), and that the HIV-1 infection was established by a single T/F virus in 10 of the 12 analyzed infants (83%). Furthermore, most HIV-1 CTL escape mutations among infants were transmitted from the mothers and did not revert during the first year of infection. Still, immune-driven selection was observed at approximately 3 months after HIV-1 infection in infants. Moreover, virus populations with CTL escape mutations in gag evolved faster than those without, independently of disease progression rate. These findings expand the current knowledge of HIV-1 transmission, evolution, and CTL escape in infant HIV-1 infection and are relevant for the development of immune-directed interventions in infants.IMPORTANCEDespite increased coverage in antiretroviral therapy for the prevention of perinatal transmission, paediatric HIV-1 infection remains a significant public health concern, especially in areas of high HIV-1 prevalence. Understanding HIV-1 transmission and the subsequent virus adaptation from the mother to the infant's host environment, as well as the viral factors that affect disease outcome, is important for the development of early immune-directed interventions for infants. This study advances our understanding of vertical HIV-1 transmission, and how infant immune selection pressure is shaping the intra-host evolutionary dynamics of HIV-1.

Pycallingcards: an integrated environment for visualizing, analyzing, and interpreting Calling Cards data.

MOTIVATION: Unraveling the transcriptional programs that control how cells divide, differentiate, and respond to their environments requires a precise understanding of transcription factors' (TFs) DNA-binding activities. Calling cards (CC) technology uses transposons to capture transient TF binding events at one instant in time and then read them out at a later time. This methodology can also be used to simultaneously measure TF binding and mRNA expression from single-cell CC and to record and integrate TF binding events across time in any cell type of interest without the need for purification. Despite these advantages, there has been a lack of dedicated bioinformatics tools for the detailed analysis of CC data. RESULTS: We introduce Pycallingcards, a comprehensive Python module specifically designed for the analysis of single-cell and bulk CC data across multiple species. Pycallingcards introduces two innovative peak callers, CCcaller and MACCs, enhancing the accuracy and speed of pinpointing TF binding sites from CC data. Pycallingcards offers a fully integrated environment for data visualization, motif finding, and comparative analysis with RNA-seq and ChIP-seq datasets. To illustrate its practical application, we have reanalyzed previously published mouse cortex and glioblastoma datasets. This analysis revealed novel cell-type-specific binding sites and potential sex-linked TF regulators, furthering our understanding of TF binding and gene expression relationships. Thus, Pycallingcards, with its user-friendly design and seamless interface with the Python data science ecosystem, stands as a critical tool for advancing the analysis of TF functions via CC data. AVAILABILITY AND IMPLEMENTATION: Pycallingcards can be accessed on the GitHub repository: https://github.com/The-Mitra-Lab/pycallingcards.