Search

Abstract Sex differences and age-related changes in the human heart at the tissue, cell, and molecular level have been well-documented and many may be relevant for cardiovascular disease. However, how molecular programs within individual cell types vary across individuals by age and sex remains poorly characterized. To better understand this variation, we performed single-nucleus combinatorial indexing (sci) ATAC- and RNA-Seq in human heart samples from nine donors. We identify hundreds of differentially expressed genes by age and sex and find epigenetic signatures of variation in ATAC-Seq data in this discovery cohort. We then scale up our single-cell RNA-Seq analysis by combining our data with five recently published single nucleus RNA-Seq datasets of healthy adult hearts. We find variation such as metabolic alterations by sex and immune changes by age in differential expression tests, as well as alterations in abundance of cardiomyocytes by sex and neurons with age. In addition, we compare our adult-derived ATAC-Seq profiles to analogous fetal cell types to identify putative developmental-stage-specific regulatory factors. Finally, we train predictive models of cell-type-specific RNA expression levels utilizing ATAC-Seq profiles to link distal regulatory sequences to promoters, quantifying the predictive value of a simple TF-to-expression regulatory grammar and identifying cell-type-specific TFs. Our analysis represents the largest single-cell analysis of cardiac variation by age and sex to date and provides a resource for further study of healthy cardiac variation and transcriptional regulation at single-cell resolution.

Retinal degeneration in mammals causes permanent loss of vision, due to an inability to regenerate naturally. Some non-mammalian vertebrates show robust regeneration, via Muller glia (MG). We have recently made significant progress in stimulating adult mouse MG to regenerate functional neurons by transgenic expression of the proneural transcription factor Ascl1. While these results showed that MG can serve as an endogenous source of neuronal replacement, the efficacy of this process is limited. With the goal of improving this in mammals, we designed a small molecule screen using sci-Plex, a method to multiplex up to thousands of single-nucleus RNA-seq conditions into a single experiment. We used this technology to screen a library of 92 compounds, identified, and validated two that promote neurogenesis in vivo. Our results demonstrate that high-throughput single-cell molecular profiling can substantially improve the discovery process for molecules and pathways that can stimulate neural regeneration and further demonstrate the potential for this approach to restore vision in patients with retinal disease.

Over the last three decades, human genetics has gone from dissecting high-penetrance Mendelian diseases to discovering the vast and complex genetic etiology of common human diseases. In tackling this complexity, scientists have discovered the importance of numerous genetic processes – most notably functional regulatory elements – in the development and progression of these diseases. Simultaneously, scientists have increasingly used multiplex assays of variant effect to systematically phenotype the cellular consequences of millions of genetic variants. In this article, we argue that the context of genetic variants – at all scales, from other genetic variants and gene regulation to cell biology to organismal environment – are critical components of how we can employ genomics to interpret these variants, and ultimately treat these diseases. We describe approaches to extend existing experimental assays and computational approaches to examine and quantify the importance of this context, including through causal analytic approaches. Having a unified understanding of the molecular, physiological, and environmental processes governing the interpretation of genetic variants is sorely needed for the field, and this perspective argues for feasible approaches by which the combined interpretation of cellular, animal, and epidemiological data can yield that knowledge.

Abstract Single-cell chromatin accessibility has emerged as a powerful means of understanding the epigenetic landscape of diverse tissues and cell types, but profiling cells from many independent specimens is challenging and costly. Here we describe a novel approach, sciPlex-ATAC-seq, which uses unmodified DNA oligos as sample-specific nuclear labels, enabling the concurrent profiling of chromatin accessibility within single nuclei from virtually unlimited specimens or experimental conditions. We first demonstrate our method with a chemical epigenomics screen, in which we identify drug-altered distal regulatory sites predictive of compound- and dose-dependent effects on transcription. We then analyze cell type-specific chromatin changes in PBMCs from multiple donors responding to synthetic and allogeneic immune stimulation. We quantify stimulation-altered immune cell compositions and isolate the unique effects of allogeneic stimulation on chromatin accessibility specific to T-lymphocytes. Finally, we observe that impaired global chromatin decondensation often coincides with chemical inhibition of allogeneic T-cell activation.

Abstract Recent advances in multiplexed single‐cell transcriptomics experiments facilitate the high‐throughput study of drug and genetic perturbations. However, an exhaustive exploration of the combinatorial perturbation space is experimentally unfeasible. Therefore, computational methods are needed to predict, interpret, and prioritize perturbations. Here, we present the compositional perturbation autoencoder (CPA), which combines the interpretability of linear models with the flexibility of deep‐learning approaches for single‐cell response modeling. CPA learns to in silico predict transcriptional perturbation response at the single‐cell level for unseen dosages, cell types, time points, and species. Using newly generated single‐cell drug combination data, we validate that CPA can predict unseen drug combinations while outperforming baseline models. Additionally, the architecture's modularity enables incorporating the chemical representation of the drugs, allowing the prediction of cellular response to completely unseen drugs. Furthermore, CPA is also applicable to genetic combinatorial screens. We demonstrate this by imputing in silico 5,329 missing combinations (97.6% of all possibilities) in a single‐cell Perturb‐seq experiment with diverse genetic interactions. We envision CPA will facilitate efficient experimental design and hypothesis generation by enabling in silico response prediction at the single‐cell level and thus accelerate therapeutic applications using single‐cell technologies.

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial–mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

Using spike-in controls with current single cell RNA-seq platforms remains a challenge. Here, the authors use a mixture of short, unmodified DNA oligos as a normalization standard for sci-RNAseq to improve the detection of global transcriptome changes.

Abstract Common alphacoronaviruses and human rhinoviruses (HRV) induce type I and III interferon (IFN) responses important to limiting viral replication in the airway epithelium. In contrast, highly pathogenic betacoronaviruses including SARS-CoV-2 may evade or antagonize RNA-induced IFN I/III responses. In airway epithelial cells (AECs) from children and older adults we compared IFN I/III responses to SARS-CoV-2 and HRV-16, and assessed whether pre-infection with HRV-16, or pretreatment with recombinant IFN-β or IFN-λ, modified SARS-CoV-2 replication. Bronchial AECs from children (ages 6–18 years) and older adults (ages 60–75 years) were differentiated ex vivo to generate organotypic cultures. In a biosafety level 3 (BSL-3) facility, cultures were infected with SARS-CoV-2 or HRV-16, and RNA and protein was harvested from cell lysates 96 h. following infection and supernatant was collected 48 and 96 h. following infection. In additional experiments cultures were pre-infected with HRV-16, or pre-treated with recombinant IFN-β1 or IFN-λ2 before SARS-CoV-2 infection. In a subset of experiments a range of infectious concentrations of HRV-16, SARS-CoV-2 WA-01, SARS-CoV-2 Delta variant, and SARS-CoV-2 Omicron variant were studied. Despite significant between-donor heterogeneity SARS-CoV-2 replicated 100 times more efficiently than HRV-16. IFNB1, INFL2, and CXCL10 gene expression and protein production following HRV-16 infection was significantly greater than following SARS-CoV-2. IFN gene expression and protein production were inversely correlated with SARS-CoV-2 replication. Treatment of cultures with recombinant IFNβ1 or IFNλ2, or pre-infection of cultures with HRV-16, markedly reduced SARS-CoV-2 replication. In addition to marked between-donor heterogeneity in IFN responses and viral replication, SARS-CoV-2 (WA-01, Delta, and Omicron variants) elicits a less robust IFN response in primary AEC cultures than does rhinovirus, and heterologous rhinovirus infection, or treatment with recombinant IFN-β1 or IFN-λ2, reduces SARS-CoV-2 replication, although to a lesser degree for the Delta and Omicron variants.

Here, the authors use single cell RNA-sequencing to generate an atlas of signaling interactions regulating embryonic hematopoietic stem cell (HSC) development and apply this knowledge to engineer a niche sufficient to support HSC maturation in vitro.

A major cause of human deafness and vestibular dysfunction is permanent loss of the mechanosensory hair cells of the inner ear. In non-mammalian vertebrates such as zebrafish, regeneration of missing hair cells can occur throughout life. While a comparative approach has the potential to reveal the basis of such differential regenerative ability, the degree to which the inner ears of fish and mammals share common hair cells and supporting cell types remains unresolved. Here, we perform single-cell RNA sequencing of the zebrafish inner ear at embryonic through adult stages to catalog the diversity of hair cells and non-sensory supporting cells. We identify a putative progenitor population for hair cells and supporting cells, as well as distinct hair and supporting cell types in the maculae versus cristae. The hair cell and supporting cell types differ from those described for the lateral line system, a distributed mechanosensory organ in zebrafish in which most studies of hair cell regeneration have been conducted. In the maculae, we identify two subtypes of hair cells that share gene expression with mammalian striolar or extrastriolar hair cells. In situ hybridization reveals that these hair cell subtypes occupy distinct spatial domains within the three macular organs, the utricle, saccule, and lagena, consistent with the reported distinct electrophysiological properties of hair cells within these domains. These findings suggest that primitive specialization of spatially distinct striolar and extrastriolar hair cells likely arose in the last common ancestor of fish and mammals. The similarities of inner ear cell type composition between fish and mammals validate zebrafish as a relevant model for understanding inner ear-specific hair cell function and regeneration.

Abstract Background Mammalian development is associated with extensive changes in gene expression, chromatin accessibility, and nuclear structure. Here, we follow such changes associated with mouse embryonic stem cell differentiation and X inactivation by integrating, for the first time, allele-specific data from these three modalities obtained by high-throughput single-cell RNA-seq, ATAC-seq, and Hi-C. Results Allele-specific contact decay profiles obtained by single-cell Hi-C clearly show that the inactive X chromosome has a unique profile in differentiated cells that have undergone X inactivation. Loss of this inactive X-specific structure at mitosis is followed by its reappearance during the cell cycle, suggesting a “bookmark” mechanism. Differentiation of embryonic stem cells to follow the onset of X inactivation is associated with changes in contact decay profiles that occur in parallel on both the X chromosomes and autosomes. Single-cell RNA-seq and ATAC-seq show evidence of a delay in female versus male cells, due to the presence of two active X chromosomes at early stages of differentiation. The onset of the inactive X-specific structure in single cells occurs later than gene silencing, consistent with the idea that chromatin compaction is a late event of X inactivation. Single-cell Hi-C highlights evidence of discrete changes in nuclear structure characterized by the acquisition of very long-range contacts throughout the nucleus. Novel computational approaches allow for the effective alignment of single-cell gene expression, chromatin accessibility, and 3D chromosome structure. Conclusions Based on trajectory analyses, three distinct nuclear structure states are detected reflecting discrete and profound simultaneous changes not only to the structure of the X chromosomes, but also to that of autosomes during differentiation. Our study reveals that long-range structural changes to chromosomes appear as discrete events, unlike progressive changes in gene expression and chromatin accessibility.

Abstract Single‐cell RNA sequencing has emerged as a powerful tool for resolving cellular states associated with normal and maligned developmental processes. Here, we used scRNA‐seq to examine the cell cycle states of expanding human neural stem cells (hNSCs). From these data, we constructed a cell cycle classifier that identifies traditional cell cycle phases and a putative quiescent‐like state in neuroepithelial‐derived cell types during mammalian neurogenesis and in gliomas. The Neural G0 markers are enriched with quiescent NSC genes and other neurodevelopmental markers found in non‐dividing neural progenitors. Putative glioblastoma stem‐like cells were significantly enriched in the Neural G0 cell population. Neural G0 cell populations and gene expression are significantly associated with less aggressive tumors and extended patient survival for gliomas. Genetic screens to identify modulators of Neural G0 revealed that knockout of genes associated with the Hippo/Yap and p53 pathways diminished Neural G0 in vitro, resulting in faster G1 transit, down‐regulation of quiescence‐associated markers, and loss of Neural G0 gene expression. Thus, Neural G0 represents a dynamic quiescent‐like state found in neuroepithelial‐derived cells and gliomas.

Abstract Microscopy is a powerful tool for characterizing complex cellular phenotypes, but linking these phenotypes to genotype or RNA expression at scale remains challenging. Here, we present Visual Cell Sorting, a method that physically separates hundreds of thousands of live cells based on their visual phenotype. Automated imaging and phenotypic analysis directs selective illumination of Dendra2, a photoconvertible fluorescent protein expressed in live cells; these photoactivated cells are then isolated using fluorescence‐activated cell sorting. First, we use Visual Cell Sorting to assess hundreds of nuclear localization sequence variants in a pooled format, identifying variants that improve nuclear localization and enabling annotation of nuclear localization sequences in thousands of human proteins. Second, we recover cells that retain normal nuclear morphologies after paclitaxel treatment, and then derive their single‐cell transcriptomes to identify pathways associated with paclitaxel resistance in cancers. Unlike alternative methods, Visual Cell Sorting depends on inexpensive reagents and commercially available hardware. As such, it can be readily deployed to uncover the relationships between visual cellular phenotypes and internal states, including genotypes and gene expression programs.

Dimensionality reduction is often used to visualize expression profiling data in order to find relationships among cells. Here, the authors use Uniform Manifold Approximation and Projection (UMAP) on published expression data of gene deletions of S. cerevisiae to find novel protein interactions.

Summary: During embryogenesis, waves of hematopoietic progenitors develop from hemogenic endothelium (HE) prior to the emergence of self-renewing hematopoietic stem cells (HSCs). Although previous studies have shown that yolk-sac-derived erythromyeloid progenitors and HSCs emerge from distinct populations of HE, it remains unknown whether the earliest lymphoid-competent progenitors, multipotent progenitors, and HSCs originate from common HE. In this study, we demonstrate by clonal assays and single-cell transcriptomics that rare HE with functional HSC potential in the early murine embryo are distinct from more abundant HE with multilineage hematopoietic potential that fail to generate HSCs. Specifically, HSC-competent HE are characterized by expression of CXCR4 surface marker and by higher expression of genes tied to arterial programs regulating HSC dormancy and self-renewal. Taken together, these findings suggest a revised model of developmental hematopoiesis in which the initial populations of multipotent progenitors and HSCs arise independently from HE with distinct phenotypic and transcriptional properties.

Summary: To study the development of the human retina, we use single-cell RNA sequencing (RNA-seq) at key fetal stages and follow the development of the major cell types as well as populations of transitional cells. We also analyze stem cell (hPSC)-derived retinal organoids; although organoids have a very similar cellular composition at equivalent ages as the fetal retina, there are some differences in gene expression of particular cell types. Moreover, the inner retinal lamination is disrupted at more advanced stages of organoids compared with fetal retina. To determine whether the disorganization in the inner retina is due to the culture conditions, we analyze retinal development in fetal retina maintained under similar conditions. These retinospheres develop for at least 6 months, displaying better inner retinal lamination than retinal organoids. Our single-cell RNA sequencing (scRNA-seq) comparisons of fetal retina, retinal organoids, and retinospheres provide a resource for developing better in vitro models for retinal disease. : To determine how closely iPSC-derived retinal organoids model development, Sridhar et al. use scRNA-seq to compare organoids with fetal retina. Despite some defects in inner retinal lamination, organoids closely mimic fetal development in timing and cell composition. Our data provide a resource for improving organoid models for retinal disease. Keywords: organoids, stem cells, neurogenesis, lamination

Summary: Induced pluripotent stem cell (iPSC)-derived hematopoietic cells represent a highly attractive source for cell and gene therapy. Given the longevity, plasticity, and self-renewal potential of distinct macrophage subpopulations, iPSC-derived macrophages (iPSC-Mφ) appear of particular interest in this context. We here evaluated the airway residence, plasticity, and therapeutic efficacy of iPSC-Mφ in a murine model of hereditary pulmonary alveolar proteinosis (herPAP). We demonstrate that single pulmonary macrophage transplantation (PMT) of 2.5–4 × 106 iPSC-Mφ yields efficient airway residence with conversion of iPSC-Mφ to an alveolar macrophage (AMφ) phenotype characterized by a distinct surface marker and gene expression profile within 2 months. Moreover, PMT significantly improves alveolar protein deposition and other critical herPAP disease parameters. Thus, our data indicate iPSC-Mφ as a source of functional macrophages displaying substantial plasticity and therapeutic potential that upon pulmonary transplantation will integrate into the lung microenvironment, adopt an AMφ phenotype and gene expression pattern, and profoundly ameliorate pulmonary disease phenotypes. : Mucci and colleagues demonstrate marked plasticity of iPSC-derived macrophages and rapid adoption of an alveolar macrophage phenotype upon pulmonary transplantation, as well as profound therapeutic efficacy of iPSC-derived macrophages in the context of the severe lung disease pulmonary alveolar proteinosis. Key words: iPSC, hematopoiesis, macrophages, lung, cell therapy

TWIST1 (TW) is a bHLH transcription factor (TF) and master regulator of the epithelial‐to‐mesenchymal transition (EMT). In vitro, TW promotes mesenchymal change, invasion, and self‐renewal in glioblastoma (GBM) cells. However, the potential therapeutic relevance of TW has not been established through loss‐of‐function studies in human GBM cell xenograft models. The effects of TW loss of function (gene editing and knockdown) on inhibition of tumorigenicity of U87MG and GBM4 glioma stem cells were tested in orthotopic xenograft models and conditional knockdown in established flank xenograft tumors. RNAseq and the analysis of tumors investigated putative TW‐associated mechanisms. Multiple bioinformatic tools revealed significant alteration of ECM, membrane receptors, signaling transduction kinases, and cytoskeleton dynamics leading to identification of PI3K/AKT signaling. We experimentally show alteration of AKT activity and periostin (POSTN) expression in vivo and/or in vitro. For the first time, we show that effect of TW knockout inhibits AKT activity in U87MG cells in vivo independent of PTEN mutation. The clinical relevance of TW and candidate mechanisms was established by analysis of the TCGA and ENCODE databases. TW expression was associated with decreased patient survival and LASSO regression analysis identified POSTN as one of top targets of TW in human GBM. While we previously demonstrated the role of TW in promoting EMT and invasion of glioma cells, these studies provide direct experimental evidence supporting protumorigenic role of TW independent of invasion in vivo and the therapeutic relevance of targeting TW in human GBM. Further, the role of TW driving POSTN expression and AKT signaling suggests actionable targets, which could be leveraged to mitigate the oncogenic effects of TW in GBM.

Summary: In vitro differentiation of pluripotent cells into β cells is a promising alternative to cadaveric-islet transplantation as a cure for type 1 diabetes (T1D). During the directed differentiation of human embryonic stem cells (hESCS) by exogenous factors, numerous genes that affect the differentiation process are turned on and off autonomously. Manipulating these reactions could increase the efficiency of differentiation and provide a more complete control over the final composition of cell populations. To uncover in vitro autonomous responses, we performed single-cell RNA sequencing on hESCs as they differentiate in spherical clusters. We observed that endocrine cells and their progenitors exist beside one another in separate compartments that activate distinct genetic pathways. WNT pathway inhibition in the endocrine domain of the differentiating clusters reveals a necessary role for the WNT inhibitor APC during islet formation in vivo. Accordingly, WNT inhibition in vitro causes an increase in the proportion of differentiated endocrine cells. : In vitro differentiation of pluripotent cells into β cells is a promising alternative to cadaveric islet transplantation as a cure for type 1 diabetes. Sharon et al. use scRNA-seq to identify the cell populations that form during the process and uncover a role for WNT pathway inhibition during endocrine differentiation. Keywords: Single-cell RNA sequencing, human embryonic stem cells, hESCs, in vitro differentiation, β cells, insulin, diabetes, pancreas development, WNT pathway, APC

Thyroid hormone (TH) regulates diverse developmental events and can drive disparate cellular outcomes. In zebrafish, TH has opposite effects on neural crest derived pigment cells of the adult stripe pattern, limiting melanophore population expansion, yet increasing yellow/orange xanthophore numbers. To learn how TH elicits seemingly opposite responses in cells having a common embryological origin, we analyzed individual transcriptomes from thousands of neural crest-derived cells, reconstructed developmental trajectories, identified pigment cell-lineage specific responses to TH, and assessed roles for TH receptors. We show that TH promotes maturation of both cell types but in distinct ways. In melanophores, TH drives terminal differentiation, limiting final cell numbers. In xanthophores, TH promotes accumulation of orange carotenoids, making the cells visible. TH receptors act primarily to repress these programs when TH is limiting. Our findings show how a single endocrine factor integrates very different cellular activities during the generation of adult form.

Viral infection can dramatically alter a cell’s transcriptome. However, these changes have mostly been studied by bulk measurements on many cells. Here we use single-cell mRNA sequencing to examine the transcriptional consequences of influenza virus infection. We find extremely wide cell-to-cell variation in the productivity of viral transcription – viral transcripts comprise less than a percent of total mRNA in many infected cells, but a few cells derive over half their mRNA from virus. Some infected cells fail to express at least one viral gene, but this gene absence only partially explains variation in viral transcriptional load. Despite variation in viral load, the relative abundances of viral mRNAs are fairly consistent across infected cells. Activation of innate immune pathways is rare, but some cellular genes co-vary in abundance with the amount of viral mRNA. Overall, our results highlight the complexity of viral infection at the level of single cells.

Single-cell RNA-sequencing methods are now robust and economically practical and are becoming a powerful tool for high-throughput, high-resolution transcriptomic analysis of cell states and dynamics. Single-cell approaches circumvent the averaging artifacts associated with traditional bulk population data, yielding new insights into the cellular diversity underlying superficially homogeneous populations. Thus far, single-cell RNA-sequencing has already shown great effectiveness in unraveling complex cell populations, reconstructing developmental trajectories, and modeling transcriptional dynamics. Ongoing technical improvements to single-cell RNA-sequencing throughput and sensitivity, the development of more sophisticated analytical frameworks for single-cell data, and an increasing array of complementary single-cell assays all promise to expand the usefulness and potential applications of single-cell transcriptomic profiling.

Changes in gene expression play an important role in evolution, yet the molecular mechanisms underlying regulatory evolution are poorly understood. Here we compare genome-wide binding of the six transcription factors that initiate segmentation along the anterior-posterior axis in embryos of two closely related species: Drosophila melanogaster and Drosophila yakuba. Where we observe binding by a factor in one species, we almost always observe binding by that factor to the orthologous sequence in the other species. Levels of binding, however, vary considerably. The magnitude and direction of the interspecies differences in binding levels of all six factors are strongly correlated, suggesting a role for chromatin or other factor-independent forces in mediating the divergence of transcription factor binding. Nonetheless, factor-specific quantitative variation in binding is common, and we show that it is driven to a large extent by the gain and loss of cognate recognition sequences for the given factor. We find only a weak correlation between binding variation and regulatory function. These data provide the first genome-wide picture of how modest levels of sequence divergence between highly morphologically similar species affect a system of coordinately acting transcription factors during animal development, and highlight the dominant role of quantitative variation in transcription factor binding over short evolutionary distances.

The thermodynamics of molecular recognition by proteins is a central determinant of complex biochemistry. For over a half-century, detailed cryogenic structures have provided deep insight into the energetic contributions to ligand binding by proteins. More recently, a dynamical proxy based on NMR-relaxation methods has revealed an unexpected richness in the contributions of conformational entropy to the thermodynamics of ligand binding. Here, we report the pressure dependence of fast internal motion within the ribonuclease barnase and its complex with the protein barstar. In what we believe is a first example, we find that protein dynamics are conserved along the pressure-binding thermodynamic cycle. The femtomolar affinity of the barnase-barstar complex exists despite a penalty by -TΔS conf of +11.7 kJ/mol at ambient pressure. At high pressure, however, the overall change in side-chain dynamics is zero, and binding occurs with no conformational entropy penalty, suggesting an important role of conformational dynamics in the adaptation of protein function to extreme environments. Distinctive clustering of the pressure sensitivity is observed in response to both pressure and binding, indicating the presence of conformational heterogeneity involving less efficiently packed alternative conformation(s). The structural segregation of dynamics observed in barnase is striking and shows how changes in both the magnitude and the sign of regional contributions of conformational entropy to the thermodynamics of protein function are possible., Competing Interests: A.J.W. is a founding member of Daedalus Innovations, LLC (Aston, PA, USA), a manufacturer of high-pressure NMR apparatus., (© 2022 The Author(s).)

This roundtable discussion is a creative contribution to this special issue on 'aging masculinities' based on questions posed by editor Stephen Katz to three leading senior scholars in the critical gerontological field in the United States. W. Andrew Achenbaum is a renowned historian who has devoted his career to writing about the relevance of past politics, cultures, and knowledges of aging to comprehending our current dilemmas. Thomas R. Cole is an acclaimed historical scholar and mentor to generations of Humanities researchers across the globe. His work also includes film, literature, ethics, and spirituality. Brian de Vries is a social gerontologist whose bold research on LGBTQ aging is a powerful voice in critiquing the multiple forms of discrimination, violence and hardships, and denied rights and life-chances imposed by hetero-patriarchal regimes in later life. That these men have shaped and been shaped by their work and advocacy is the key theme that inspires our conversation., (Copyright © 2022 Elsevier Inc. All rights reserved.)

To advance methods for bioaccumulation assessment of organic substances in air-breathing organisms, the present study developed an in vitro approach for screening neutral hydrophobic organic substances for their bioaccumulation potential in air-breathing organisms consisting of (1) depletion assays for chemicals in rat liver S9 subcellular fractions, (2) in vitro-in vivo extrapolation, and (3) whole-organism bioaccumulation modeling to assess the biomagnification potential of neutral organic substances in the rat. Testing of the in vitro method on 14 test chemicals of potentially biomagnifying substances showed that the bioassays could be conducted with a high level of reproducibility and that in vitro-derived elimination rate constants were in good agreement with in vivo-determined elimination rate constants in the rat. Exploring the potential of the in vitro approach for screening organic chemicals for bioaccumulation in air-breathing organisms indicated that chemical substances that exhibit a depletion rate constant in the S9 in vitro bioassay ≥0.3 h -1 are not expected to biomagnify in rats independent of their octanol-water partitioning coefficient (K OW ) or octanol-air partitioning coefficient (K OA ). The high level of reproducibility achieved in the test, combined with the good agreement between in vitro-derived and in vivo-determined depuration rates, suggests that the in vitro approach in combination with a K OA - and K OW -based screening approach has good potential for screening chemicals in commerce for their bioaccumulation potential in air-breathing organisms in a cost-effective and expedient manner, especially if the bioassay can be automated. Environ Toxicol Chem 2022;41:2565-2579. © 2022 SETAC., (© 2022 SETAC.)

Cysteine proteases comprise an important class of drug targets, especially for infectious diseases such as Chagas disease (cruzain) and COVID-19 (3CL protease, cathepsin L). Peptide aldehydes have proven to be potent inhibitors for all of these proteases. However, the intrinsic, high electrophilicity of the aldehyde group is associated with safety concerns and metabolic instability, limiting the use of aldehyde inhibitors as drugs. We have developed a novel class of self-masked aldehyde inhibitors (SMAIs) for cruzain, the major cysteine protease of the causative agent of Chagas disease- Trypanosoma cruzi . These SMAIs exerted potent, reversible inhibition of cruzain ( K i * = 18-350 nM) while apparently protecting the free aldehyde in cell-based assays. We synthesized prodrugs of the SMAIs that could potentially improve their pharmacokinetic properties. We also elucidated the kinetic and chemical mechanism of SMAIs and applied this strategy to the design of anti-SARS-CoV-2 inhibitors.

Conserved homology 1 (C1) domains are peripheral zinc finger domains that are responsible for recruiting their host signaling proteins, including Protein Kinase C (PKC) isoenzymes, to diacylglycerol-containing lipid membranes. In this work, we investigated the reactivity of the C1 structural zinc sites, using the cysteine-rich C1B regulatory region of the PKCα isoform as a paradigm. The choice of Cd 2+ as a probe was prompted by previous findings that xenobiotic metal ions modulate PKC activity. Using solution NMR and UV-vis spectroscopy, we found that Cd 2+ spontaneously replaced Zn 2+ in both structural sites of the C1B domain, with the formation of all-Cd and mixed Zn/Cd protein species. The Cd 2+ substitution for Zn 2+ preserved the C1B fold and function, as probed by its ability to interact with a potent tumor-promoting agent. Both Cys 3 His metal-ion sites of C1B have higher affinity to Cd 2+ than Zn 2+ , but are thermodynamically and kinetically inequivalent with respect to the metal ion replacement, despite the identical coordination spheres. We find that even in the presence of the oxygen-rich sites presented by the neighboring peripheral membrane-binding C2 domain, the thiol-rich sites can successfully compete for the available Cd 2+ . Our results indicate that Cd 2+ can target the entire membrane-binding regulatory region of PKCs, and that the competition between the thiol- and oxygen-rich sites will likely determine the activation pattern of PKCs., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Cole and Igumenova.)

Context: Kisspeptin-neurokinin B (NKB)-dynorphin neurons are critical regulators of the hypothalamic-pituitary-gonadal axis. NKB and dynorphin are hypothesized to influence the frequency of GnRH pulses, whereas kisspeptin is hypothesized to be a generator of the GnRH pulse. How these neuropeptides interact remains unclear., Objective: To probe the role of NKB in GnRH pulse generation and to determine the interactions between NKB, kisspeptin, and dynorphin in humans and mice with a complete absence of NKB., Design: Case/control., Setting: Academic medical center., Participants: Members of a consanguineous family bearing biallelic loss-of-function mutations in the gene encoding NKB and NKB-deficient mice., Interventions: Frequent blood sampling to characterize neuroendocrine profile and administration of kisspeptin, GnRH, and naloxone, a nonspecific opioid receptor antagonist used to block dynorphin., Main Outcome Measures: LH pulse characteristics., Results: Humans lacking NKB demonstrate slow LH pulse frequency, which can be increased by opioid antagonism. Mice lacking NKB also demonstrate impaired LH secretion, which can be augmented with an identical pharmacologic manipulation. Both mice and humans with NKB deficiency respond to exogenous kisspeptin., Conclusion: The preservation of LH pulses in the absence of NKB and dynorphin signaling suggests that both peptides are dispensable for GnRH pulse generation and kisspeptin responsiveness. However, NKB and dynorphin appear to have opposing roles in the modulation of GnRH pulse frequency., (Copyright © 2019 Endocrine Society.)

The field of maternal-fetal intervention is rapidly progressing and with it comes new and often complex ethical considerations that must be addressed. The purpose of this article is to review the ethical issues that arise in maternal-fetal intervention. We will provide two clinical scenarios and discuss the ethical issues related to each scenario and how they were addressed. We will also provide a list of recommended resources that any institutional offering maternal-fetal intervention should have in place to meet the ethical obligations of such work., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Cd 2+ and Pb 2+ are xenobiotic heavy metal ions that use ionic mimicry to interfere with the cellular function of biomacromolecules. Using a combination of SAXS, electron microscopy, FRET, and solution NMR spectroscopy, we demonstrate that treatment with Cd 2+ and Pb 2+ causes self-assembly of protein kinase C regulatory domains that peripherally associate with membranes. The self-assembly process successfully competes with ionic mimicry and is mediated by conserved protein regions that are distinct from the canonical Ca 2+ -binding motifs of protein kinase C. The ability of protein oligomers to interact with anionic membranes is enhanced compared to the monomeric species. Our findings suggest that metal-ion-dependent peripheral membrane domains can be utilized for generating protein-metal-ion nanoclusters and serve as biotemplates for the design of sequestration agents.

Congenital hypogonadotropic hypogonadism (CHH) is a rare genetic form of isolated gonadotropin-releasing hormone (GnRH) deficiency caused by mutations in > 30 genes. Fibroblast growth factor receptor 1 ( FGFR1 ) is the most frequently mutated gene in CHH and is implicated in GnRH neuron development and maintenance. We note that a CHH FGFR1 mutation (p.L342S) decreases signaling of the metabolic regulator FGF21 by impairing the association of FGFR1 with β-Klotho (KLB), the obligate co-receptor for FGF21. We thus hypothesized that the metabolic FGF21/KLB/FGFR1 pathway is involved in CHH Genetic screening of 334 CHH patients identified seven heterozygous loss-of-function KLB mutations in 13 patients (4%). Most patients with KLB mutations (9/13) exhibited metabolic defects. In mice, lack of Klb led to delayed puberty, altered estrous cyclicity, and subfertility due to a hypothalamic defect associated with inability of GnRH neurons to release GnRH in response to FGF21. Peripheral FGF21 administration could indeed reach GnRH neurons through circumventricular organs in the hypothalamus. We conclude that FGF21/KLB/FGFR1 signaling plays an essential role in GnRH biology, potentially linking metabolism with reproduction., (© 2017 The Authors. Published under the terms of the CC BY 4.0 license.)

This essay reads the protagonist of Cormac McCarthy's No Country for Old Men (2005), Sheriff Ed Tom Bell, as an exemplar of problems that contemporary aging men face when they look to ahead to the so-called Fourth Age. As the plot unfolds, Bell is an aging, increasingly ineffectual cowboy lawman who retires, renounces the violence that sustained his male dominance, and loses the moral certainty that ensured his identity. Like Bell, most old men struggle with four interrelated challenges as they move along the ever-lengthening journey of life: relevance, masculinity, love, and meaning.

Ca 2+ -dependent conserved-region 2 (C2) domains target their host signaling proteins to anionic membranes. The Ca 2+ -binding event is a prerequisite for membrane association. Here, we investigate multiscale metal-ion-dependent dynamics of the C2 domain of protein kinase Cα (C2α) using NMR spectroscopy. Interactions with metal ions attenuate microsecond-timescale motions of the loop regions, indicating that preorganization of the metal-binding loops occurs before membrane insertion. Binding of a full complement of Ca 2+ ions has a profound effect on the millisecond-timescale dynamics of the N- and C-terminal regions of C2α. We propose that Ca 2+ binding allosterically destabilizes the terminal regions of C2α and thereby facilitates the conformational rearrangement necessary for full membrane insertion and activation of protein kinase Cα., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

We report the protocol for heterologous expression and purification of the N-terminal regulatory region of two Protein Kinase C (PKC)(1) isozymes, one conventional and one novel. Previous studies of these domains relied almost exclusively on the fusion constructs with high-molecular-weight solubility fusion partners such as GST and MBP. We developed experimental procedures that enabled us to overcome challenges associated with the amphiphilic character of the regulatory domain and generate sufficient quantities of fusion partner-free proteins for biophysical work. The key features of the protocol are the identity of the cleavable fusion partner, expression conditions, growth medium additives, introduction of mutation/solubility tags, and incorporation of osmolytes. The protein yields are sufficient for cost-effective production of isotopically enriched proteins for NMR work and biophysical studies in general. Our work opens up an avenue for the structural studies of these challenging proteins with high amphiphilic character., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Conventional and novel isoenzymes of PKC are activated by the membrane-embedded second messenger diacylglycerol (DAG) through its interactions with the C1 regulatory domain. The affinity of C1 domains to DAG varies considerably among PKCs. To gain insight into the origin of differential DAG affinities, we conducted high-resolution NMR studies of C1B domain from PKCδ (C1Bδ) and its W252Y variant. The W252Y mutation was previously shown to render C1Bδ less responsive to DAG (Dries, D. R., Gallegos, L. L., and Newton, A. C. (2007) A single residue in the C1 domain sensitizes novel protein kinase C isoforms to cellular diacylglycerol production. J. Biol. Chem. 282, 826-830) and thereby emulate the behavior of C1B domains from conventional PKCs that have a conserved Tyr at the equivalent position. Our data revealed that W252Y mutation did not perturb the conformation of C1Bδ in solution but significantly reduced its propensity to partition into a membrane-mimicking environment in the absence of DAG. Using detergent micelles doped with a paramagnetic lipid, we determined that both the residue identity at position 252 and complexation with diacylglycerol influence the geometry of C1Bδ-micelle interactions. In addition, we identified the C-terminal helix α1 of C1Bδ as an interaction site with the head groups of phosphatidylserine, a known activator of PKCδ. Taken together, our studies (i) reveal the identities of C1Bδ residues involved in interactions with membrane-mimicking environment, DAG, and phosphatidylserine, as well as the affinities associated with each event and (ii) suggest that the initial ligand-independent membrane recruitment of C1B domains, which is greatly facilitated by the interfacial partitioning of Trp-252, is responsible, at least in part, for the differential DAG affinities., (© 2014 by The American Society for Biochemistry and Molecular Biology, Inc.)

Context: Pheochromocytomas and paragangliomas are notable for a high frequency of inherited cases, many of which present as apparently sporadic tumors., Objective: The objective of this study was to establish a comprehensive next generation sequencing (NGS)-based strategy for the diagnosis of patients with pheochromocytoma and paraganglioma by testing simultaneously for mutations in MAX, RET, SDHA, SDHB, SDHC, SDHD, SDHAF2, TMEM127, and VHL., Design: After the methodology for the assay was designed and established, it was validated on DNA samples with known genotype and then patients were studied prospectively., Setting: The study was performed in a diagnostic genetics laboratory., Patients: DNA samples from 205 individuals affected with adrenal or extraadrenal pheochromocytoma/head and neck paraganglioma (PPGL/HNPGL) were analyzed. A proof-of-principle study was performed using 85 samples known to contain a variant in 1 or more of the genes to be tested, followed by prospective analysis of an additional 120 samples., Main Outcome Measures: We assessed the ability to use an NGS-based method to perform comprehensive analysis of genes implicated in inherited PPGL/HNPGL., Results: The proof-of-principle study showed that the NGS assay and analysis gave a sensitivity of 98.7%. A pathogenic mutation was identified in 16.6% of the prospective analysis cohort of 120 patients., Conclusions: A comprehensive NGS-based strategy for the analysis of genes associated with predisposition to PPGL and HNPGL was established, validated, and introduced into diagnostic service. The new assay provides simultaneous analysis of 9 genes and allows more rapid and cost-effective mutation detection than the previously used conventional Sanger sequencing-based methodology.

Objectives: Research studies have reported that about a third of individuals with phaeochromocytoma/paraganglioma (PPGL) have an inherited predisposition, although the frequency of specific mutations can vary between populations. We evaluated VHL, SDHB and SDHD mutation testing in cohorts of patients with non-syndromic PPGL and head and neck paraganglioma (HNPGL)., Design: Prospective, observational evaluation of NHS practice., Patients: Individuals with PPGL/HNPGL referred to a supraregional genetics testing service over a 10-year period., Measurements: Clinical (age, tumour site, malignancy, etc.), mutation frequencies and characteristics., Results: A total of 501 probands with PPGL (n = 413) or HNPGL (n = 88) were studied. Thirty-one percent of patients with PPGL presented had a pathogenic mutation in SDHB, SDHD or VHL. Mutation detection rates were highest in those with a positive family history (62%), malignancy (53%), multiple tumours (33%) or PGL (44%). Twenty-eight percent of individuals with a single sporadic phaeochromocytoma had a mutation. Overall, 63% of patients with HNPGL had a mutation (92% of those with a family history, 89% of those with multicentric tumours and 34% of those with a single sporadic HNPGL). Penetrance was calculated in 121 SDHB mutation-positive probands and 187 of their mutation-positive relatives. Most relatives were asymptomatic and lifetime penetrance in non-proband SDHB mutation carriers was <50%., Conclusions: Practice-based evaluations of genetic testing in PPGL reveal high mutation detection rates. Although clinical criteria can be used to prioritize mutation testing, mutations were detected in 'low risk groups' indicating a need for comprehensive and inexpensive genetic testing strategies for PPGL and HNPGL., (© 2012 John Wiley & Sons Ltd.)

Context: P450 oxidoreductase deficiency (PORD) is a unique congenital adrenal hyperplasia variant that manifests with glucocorticoid deficiency, disordered sex development (DSD), and skeletal malformations. No comprehensive data on genotype-phenotype correlations in Caucasian patients are available., Objective: The objective of the study was to establish genotype-phenotype correlations in a large PORD cohort., Design: The design of the study was the clinical, biochemical, and genetic assessment including multiplex ligation-dependent probe amplification (MLPA) in 30 PORD patients from 11 countries., Results: We identified 23 P450 oxidoreductase (POR) mutations (14 novel) including an exonic deletion and a partial duplication detected by MLPA. Only 22% of unrelated patients carried homozygous POR mutations. p.A287P was the most common mutation (43% of unrelated alleles); no other hot spot was identified. Urinary steroid profiling showed characteristic PORD metabolomes with variable impairment of 17α-hydroxylase and 21-hydroxylase. Short cosyntropin testing revealed adrenal insufficiency in 89%. DSD was present in 15 of 18 46,XX and seven of 12 46,XY individuals. Homozygosity for p.A287P was invariably associated with 46,XX DSD but normal genitalia in 46,XY individuals. The majority of patients with mild to moderate skeletal malformations, assessed by a novel scoring system, were compound heterozygous for missense mutations, whereas nearly all patients with severe malformations carried a major loss-of-function defect on one of the affected alleles., Conclusions: We report clinical, biochemical, and genetic findings in a large PORD cohort and show that MLPA is a useful addition to POR mutation analysis. Homozygosity for the most frequent mutation in Caucasians, p.A287P, allows for prediction of genital phenotype and moderate malformations. Adrenal insufficiency is frequent, easily overlooked, but readily detected by cosyntropin testing.

Hypertension affects one billion people and is a principal reversible risk factor for cardiovascular disease. Pseudohypoaldosteronism type II (PHAII), a rare Mendelian syndrome featuring hypertension, hyperkalaemia and metabolic acidosis, has revealed previously unrecognized physiology orchestrating the balance between renal salt reabsorption and K(+) and H(+) excretion. Here we used exome sequencing to identify mutations in kelch-like 3 (KLHL3) or cullin 3 (CUL3) in PHAII patients from 41 unrelated families. KLHL3 mutations are either recessive or dominant, whereas CUL3 mutations are dominant and predominantly de novo. CUL3 and BTB-domain-containing kelch proteins such as KLHL3 are components of cullin-RING E3 ligase complexes that ubiquitinate substrates bound to kelch propeller domains. Dominant KLHL3 mutations are clustered in short segments within the kelch propeller and BTB domains implicated in substrate and cullin binding, respectively. Diverse CUL3 mutations all result in skipping of exon 9, producing an in-frame deletion. Because dominant KLHL3 and CUL3 mutations both phenocopy recessive loss-of-function KLHL3 mutations, they may abrogate ubiquitination of KLHL3 substrates. Disease features are reversed by thiazide diuretics, which inhibit the Na-Cl cotransporter in the distal nephron of the kidney; KLHL3 and CUL3 are expressed in this location, suggesting a mechanistic link between KLHL3 and CUL3 mutations, increased Na-Cl reabsorption, and disease pathogenesis. These findings demonstrate the utility of exome sequencing in disease gene identification despite the combined complexities of locus heterogeneity, mixed models of transmission and frequent de novo mutation, and establish a fundamental role for KLHL3 and CUL3 in blood pressure, K(+) and pH homeostasis.