Niche-DE: niche-differential gene expression analysis in spatial transcriptomics data identifies context-dependent cell-cell interactions.

Existing methods for analysis of spatial transcriptomic data focus on delineating the global gene expression variations of cell types across the tissue, rather than local gene expression changes driven by cell-cell interactions. We propose a new statistical procedure called niche-differential expression (niche-DE) analysis that identifies cell-type-specific niche-associated genes, which are differentially expressed within a specific cell type in the context of specific spatial niches. We further develop niche-LR, a method to reveal ligand-receptor signaling mechanisms that underlie niche-differential gene expression patterns. Niche-DE and niche-LR are applicable to low-resolution spot-based spatial transcriptomics data and data that is single-cell or subcellular in resolution.

Combing Genome-Wide Association Studies and Single-Cell Analysis to Elucidate the Mechanisms of Kidney Disease: Proceedings of the Henry Shavelle Professorship.

Background: Kidney diseases pose a significant global health burden; there is an urgent need to deepen our understanding of their underlying mechanisms. Summary: This review focuses on new innovative approaches that merge genome-wide association studies (GWAS) and single-cell omics (including transcriptomics) in kidney disease research. We begin by detailing how GWAS has identified numerous genetic risk factors, offering valuable insight into disease susceptibility. Then, we explore the application of scRNA-seq, highlighting its ability to unravel how genetic variants influence cellular phenotypes. Through a synthesis of recent studies, we illuminate the synergy between these two powerful methodologies, demonstrating their potential in elucidating the complex etiology of kidney diseases. Moreover, we discuss how this integrative approach could pave the way for precise diagnostics and personalized treatments. Key Message: This review underscores the transformative potential of combining GWAS and scRNA-seq in the journey toward a deeper understanding of kidney diseases.

Decomposing a deterministic path to mesenchymal niche formation by two intersecting morphogen gradients.

Organ formation requires integrating signals to coordinate proliferation, specify cell fates, and shape tissue. Tracing these events and signals remains a challenge, as intermediate states across many critical transitions are unresolvable over real time and space. Here, we designed a unique computational approach to decompose a non-linear differentiation process into key components to resolve the signals and cell behaviors that drive a rapid transition, using the hair follicle dermal condensate as a model. Combining scRNA sequencing with genetic perturbation, we reveal that proliferative Dkk1+ progenitors transiently amplify to become quiescent dermal condensate cells by the mere spatiotemporal patterning of Wnt/ß-catenin and SHH signaling gradients. Together, they deterministically coordinate a rapid transition from proliferation to quiescence, cell fate specification, and morphogenesis. Moreover, genetically repatterning these gradients reproduces these events autonomously in "slow motion" across more intermediates that resolve the process. This analysis unravels two morphogen gradients that intersect to coordinate events of organogenesis.

The Genomic and Phenotypic Landscape of Ichthyosis: An Analysis of 1000 Kindreds.

IMPORTANCE: Ichthyoses are clinically and genetically heterogeneous disorders characterized by scaly skin. Despite decades of investigation identifying pathogenic variants in more than 50 genes, clear genotype-phenotype associations have been difficult to establish. OBJECTIVE: To expand the genotypic and phenotypic spectra of ichthyosis and delineate genotype-phenotype associations. DESIGN, SETTING, AND PARTICIPANTS: This cohort study recruited an international group of individuals with ichthyosis and describes characteristic and distinguishing features of common genotypes, including genotype-phenotype associations, during a 10-year period from June 2011 to July 2021. Participants of all ages, races, and ethnicities were included and were enrolled worldwide from referral centers and patient advocacy groups. A questionnaire to assess clinical manifestations was completed by those with a genetic diagnosis. MAIN OUTCOMES AND MEASURES: Genetic analysis of saliva or blood DNA, a phenotyping questionnaire, and standardized clinical photographs. Descriptive statistics, such as frequency counts, were used to describe the cases in the cohort. Fisher exact tests identified significant genotype-phenotype associations. RESULTS: Results were reported for 1000 unrelated individuals enrolled from around the world (mean [SD] age, 50.0 [34.0] years; 524 [52.4%] were female, 427 [42.7%] were male, and 49 [4.9%] were not classified); 75% were from the US, 12% from Latin America, 4% from Canada, 3% from Europe, 3% from Asia, 2% from Africa, 1% from the Middle East, and 1% from Australia and New Zealand. A total of 266 novel disease-associated variants in 32 genes were identified among 869 kindreds. Of these, 241 (91%) pathogenic variants were found through multiplex amplicon sequencing and 25 (9%) through exome sequencing. Among the 869 participants with a genetic diagnosis, 304 participants (35%) completed the phenotyping questionnaire. Analysis of clinical manifestations in these 304 individuals revealed that pruritus, hypohydrosis, skin pain, eye problems, skin odor, and skin infections were the most prevalent self-reported features. Genotype-phenotype association analysis revealed that the presence of a collodion membrane at birth (odds ratio [OR], 6.7; 95% CI, 3.0-16.7; P < .001), skin odor (OR, 2.8; 95% CI, 1.1-6.8; P = .02), hearing problems (OR, 2.9; 95% CI, 1.6-5.5; P < .001), eye problems (OR, 3.0; 95% CI, 1.5-6.0; P < .001), and alopecia (OR, 4.6; 95% CI, 2.4-9.0; P < .001) were significantly associated with TGM1 variants compared with other ichthyosis genotypes studied. Skin pain (OR, 6.8; 95% CI, 1.6-61.2; P = .002), odor (OR, 5.7; 95% CI, 2.0-19.7; P < .001), and infections (OR, 3.1; 95% CI, 1.4-7.7; P = .03) were significantly associated with KRT10 pathogenic variants compared with disease-associated variants in other genes that cause ichthyosis. Pathogenic variants were identified in 869 (86.9%) participants. Most of the remaining individuals had unique phenotypes, enabling further genetic discovery. CONCLUSIONS AND RELEVANCE: This cohort study expands the genotypic and phenotypic spectrum of ichthyosis, establishing associations between clinical manifestations and genotypes. Collectively, the findings may help improve clinical assessment, assist with developing customized management plans, and improve clinical course prognostication.

Putting Cell Competition under the Microscope.

In two recent issues of Nature, Ellis et al. (2019) and Liu et al. (2019) examine mechanisms that influence clonal dynamics in the epidermis during development and homeostasis, respectively. These studies provide insight into the selection process that determines which cells expand and which cells are lost.

Second-Hit, Postzygotic PMVK and MVD Mutations in Linear Porokeratosis.

Importance: Linear porokeratosis features linear and whorled configurations of keratotic papules and plaques, with coronoid lamellae present on histologic examination. Because linear porokeratosis manifests in the lines of Blaschko representing the dorsoventral migration patterns of keratinocyte precursors, it has been suggested that postzygotic somatic mutation underlies the disease. However, no genetic evidence has supported this hypothesis to date. Objective: To identify genetic mutations associated with linear porokeratosis. Design, Setting, and Participants: Paired whole-exome sequencing of affected skin and blood/saliva samples from 3 participants from 3 academic medical centers with clinical and histologic diagnoses of linear porokeratosis. Interventions or Exposures: Whole-exome sequencing of paired blood/saliva and affected tissue samples isolated from linear porokeratosis lesions. Main Outcomes and Measures: Germline and somatic genomic characteristics of participants with linear porokeratosis. Results: Of the 3 participants, 2 were male. Participant ages ranged from 5 to 20 years old. We found a combination of a novel germline mutation and a novel somatic mutation within affected tissue in all cases. One participant had a germline heterozygous PMVK c.329G>A mutation and a somatic copy-neutral loss of heterozygosity confined to the lesional skin, while a second had a germline heterozygous PMVK c.79G>T mutation and an additional PMVK c.379C>T mutation in the lesional skin. In a third participant, there was a germline splice-site mutation in MVD (c.70 + 5G>A) and a somatic deletion in MVD causing frameshift and premature codon termination within the lesional skin (c.811_815del, p.F271Afs*33 frameshift). Conclusions and Relevance: Our findings suggest that linear porokeratosis is associated with the presence of second-hit postzygotic mutations in the genes that encode enzymes within the mevalonate biosynthesis pathway, and provide further evidence that the mevalonate pathway may be a potential target for therapeutic intervention in porokeratosis.

Single-Cell Analysis Reveals a Hair Follicle Dermal Niche Molecular Differentiation Trajectory that Begins Prior to Morphogenesis.

Delineating molecular and cellular events that precede appendage morphogenesis has been challenging due to the inability to distinguish quantitative molecular differences between cells that lack histological distinction. The hair follicle (HF) dermal condensate (DC) is a cluster of cells critical for HF development and regeneration. Events that presage emergence of this distinctive population are poorly understood. Using unbiased single-cell RNA sequencing and in vivo methods, we infer a sequence of transcriptional states through which DC cells pass that begins prior to HF morphogenesis. Our data indicate that Wnt/ß-catenin signaling is required to progress into an intermediate stage that precedes quiescence and differentiation. Further, we provide evidence that quiescent DC cells are recent progeny of selectively proliferating cells present prior to morphogenesis and that are later identified in the peri-DC zone during DC expansion. Together, these findings provide an inferred path of molecular states that lead to DC cell differentiation.

Palmoplantar Keratoderma in Costello Syndrome Responsive to Acitretin.

Costello syndrome (CS) is a multisystem congenital disorder characterized by coarse facial features, cardiac defects, intellectual disability, and predisposition to malignancies. Dermatologic findings can include cutaneous papillomas, skin redundancy, acanthosis nigricans, and keratosis pilaris. Palmoplantar keratoderma (PPK) is present in approximately 76% of patients with CS, with disabling functional consequences in severe cases. We report a case of CS with severe PPK that improved dramatically with systemic administration of acitretin 0.3 mg/kg/day.

Somatic Mutations in NEK9 Cause Nevus Comedonicus.

Acne vulgaris (AV) affects most adolescents, and of those affected, moderate to severe disease occurs in 20%. Comedones, follicular plugs consisting of desquamated keratinocytes and sebum, are central to its pathogenesis. Despite high heritability in first-degree relatives, AV genetic determinants remain incompletely understood. We therefore employed whole-exome sequencing (WES) in nevus comedonicus (NC), a rare disorder that features comedones and inflammatory acne cysts in localized, linear configurations. WES identified somatic NEK9 mutations, each affecting highly conserved residues within its kinase or RCC1 domains, in affected tissue of three out of three NC-affected subjects. All mutations are gain of function, resulting in increased phosphorylation at Thr210, a hallmark of NEK9 kinase activation. We found that comedo formation in NC is marked by loss of follicular differentiation markers, expansion of keratin-15-positive cells from localization within the bulge to the entire sub-bulge follicle and cyst, and ectopic expression of keratin 10, a marker of interfollicular differentiation not present in normal follicles. These findings suggest that NEK9 mutations in NC disrupt normal follicular differentiation and identify NEK9 as a potential regulator of follicular homeostasis.

Mutational landscape of MCPyV-positive and MCPyV-negative Merkel cell carcinomas with implications for immunotherapy.

Merkel cell carcinoma (MCC) is a rare but highly aggressive cutaneous neuroendocrine carcinoma, associated with the Merkel cell polyomavirus (MCPyV) in 80% of cases. To define the genetic basis of MCCs, we performed exome sequencing of 49 MCCs. We show that MCPyV-negative MCCs have a high mutation burden (median of 1121 somatic single nucleotide variants (SSNVs) per-exome with frequent mutations in RB1 and TP53 and additional damaging mutations in genes in the chromatin modification (ASXL1, MLL2, and MLL3), JNK (MAP3K1 and TRAF7), and DNA-damage pathways (ATM, MSH2, and BRCA1). In contrast, MCPyV-positive MCCs harbor few SSNVs (median of 12.5 SSNVs/tumor) with none in the genes listed above. In both subgroups, there are rare cancer-promoting mutations predicted to activate the PI3K pathway (HRAS, KRAS, PIK3CA, PTEN, and TSC1) and to inactivate the Notch pathway (Notch1 and Notch2). TP53 mutations appear to be clinically relevant in virus-negative MCCs as 37% of these tumors harbor potentially targetable gain-of-function mutations in TP53 at p.R248 and p.P278. Moreover, TP53 mutational status predicts death in early stage MCC (5-year survival in TP53 mutant vs wild-type stage I and II MCCs is 20% vs. 92%, respectively; P = 0.0036). Lastly, we identified the tumor neoantigens in MCPyV-negative and MCPyV-positive MCCs. We found that virus-negative MCCs harbor more tumor neoantigens than melanomas or non-small cell lung cancers (median of 173, 65, and 111 neoantigens/sample, respectively), two cancers for which immune checkpoint blockade can produce durable clinical responses. Collectively, these data support the use of immunotherapies for virus-negative MCCs.

A Somatic p.G45E GJB2 Mutation Causing Porokeratotic Eccrine Ostial and Dermal Duct Nevus.

IMPORTANCE: Recent data demonstrated somatic mutations in GJB2 that were present in affected porokeratotic eccrine ostial and dermal duct nevus (PEODDN) tissue but absent in unaffected skin. Recognizing that PEODDN lesions can also appear in individuals with keratitis-ichthyosis-deafness syndrome and finding somatic mutations in their cohort, the authors concluded that somatic GJB2 mutation may cause PEODDN. By using whole-exome sequencing, we show that somatic GJB2 mutation alone is sufficient to cause PEODDN. OBSERVATIONS: We performed whole-exome sequencing of paired blood and affected tissue samples isolated from a PEODDN lesion of a primary school-aged female patient with bands of hyperkeratotic-affected skin on the upper and lower extremities and trunk, and identified a single, protein-damaging p.Gly45Glu GJB2 mutation present in tissue samples but not in blood samples. CONCLUSION AND RELEVANCE: Our results prove that somatic GJB2 mutation is sufficient to cause PEODDN. Dominantly inherited GJB2 mutations, including the p.Gly45Glu found in our case, have been shown to cause the severe multisystem disorder keratitis-ichthyosis-deafness syndrome. GJB2 encodes connexin 26, a gap junction protein, which permits intercellular ion and macromolecule flux. Individuals with somatic mosaicism are at risk for transmitting systemic disease to their offspring, and all individuals with PEODDN lesions should be counseled regarding the risk of having a child with keratitis-ichthyosis-deafness syndrome.

Anthrax lethal factor cleaves mouse nlrp1b in both toxin-sensitive and toxin-resistant macrophages.

Anthrax lethal factor (LF) is the protease component of anthrax lethal toxin (LT). LT induces pyroptosis in macrophages of certain inbred mouse and rat strains, while macrophages from other inbred strains are resistant to the toxin. In rats, the sensitivity of macrophages to toxin-induced cell death is determined by the presence of an LF cleavage sequence in the inflammasome sensor Nlrp1. LF cleaves rat Nlrp1 of toxin-sensitive macrophages, activating caspase-1 and inducing cell death. Toxin-resistant macrophages, however, express Nlrp1 proteins which do not harbor the LF cleavage site. We report here that mouse Nlrp1b proteins are also cleaved by LF. In contrast to the situation in rats, sensitivity and resistance of Balb/cJ and NOD/LtJ macrophages does not correlate to the susceptibility of their Nlrp1b proteins to cleavage by LF, as both proteins are cleaved. Two LF cleavage sites, at residues 38 and 44, were identified in mouse Nlrp1b. Our results suggest that the resistance of NOD/LtJ macrophages to LT, and the inability of the Nlrp1b protein expressed in these cells to be activated by the toxin are likely due to polymorphisms other than those at the LF cleavage sites.

Anthrax lethal factor cleavage of Nlrp1 is required for activation of the inflammasome.

NOD-like receptor (NLR) proteins (Nlrps) are cytosolic sensors responsible for detection of pathogen and danger-associated molecular patterns through unknown mechanisms. Their activation in response to a wide range of intracellular danger signals leads to formation of the inflammasome, caspase-1 activation, rapid programmed cell death (pyroptosis) and maturation of IL-1ß and IL-18. Anthrax lethal toxin (LT) induces the caspase-1-dependent pyroptosis of mouse and rat macrophages isolated from certain inbred rodent strains through activation of the NOD-like receptor (NLR) Nlrp1 inflammasome. Here we show that LT cleaves rat Nlrp1 and this cleavage is required for toxin-induced inflammasome activation, IL-1 ß release, and macrophage pyroptosis. These results identify both a previously unrecognized mechanism of activation of an NLR and a new, physiologically relevant protein substrate of LT.

SHIP-1 increases early oxidative burst and regulates phagosome maturation in macrophages.

Although the inositol phosphatase SHIP-1 is generally thought to inhibit signaling for Fc receptor-mediated phagocytosis, the product of its activity, phosphatidylinositol 3,4 bisphosphate (PI(3,4)P(2)), has been implicated in activation of the NADPH oxidase. This suggests that SHIP-1 positively regulates the generation of reactive oxygen species after phagocytosis. To examine how SHIP-1 activity contributes to Fc receptor-mediated phagocytosis, we measured and compared phospholipid dynamics, membrane trafficking, and the oxidative burst in macrophages from SHIP-1-deficient and wild-type mice. SHIP-1-deficient macrophages showed significantly elevated ratios of PI(3,4,5)P(3) to PI(3,4)P(2) on phagosomal membranes. Imaging reactive oxygen intermediate activities in phagosomes revealed decreased early NADPH oxidase activity in SHIP-1-deficient macrophages. SHIP-1 deficiency also altered later stages of phagosome maturation, as indicated by the persistent elevation of PI(3)P and the early localization of Rab5a to phagosomes. These direct measurements of individual organelles indicate that phagosomal SHIP-1 enhances the early oxidative burst through localized alteration of the membrane 3'-phosphoinositide composition.

Differential association of phosphatidylinositol 3-kinase, SHIP-1, and PTEN with forming phagosomes.

In macrophages, enzymes that synthesize or hydrolyze phosphatidylinositol (3,4,5)-trisphosphate [PI(3,4,5)P(3)] regulate Fcgamma receptor-mediated phagocytosis. Inhibition of phosphatidylinositol 3-kinase (PI3K) or overexpression of the lipid phosphatases phosphatase and tensin homologue (PTEN) and Src homology 2 domain-containing inositol phosphatase (SHIP-1), which hydrolyze PI(3,4,5)P(3) to phosphatidylinositol 4,5-bisphosphate and phosphatidylinositol 3,4-bisphosphate [PI(3,4)P(2)], respectively, inhibit phagocytosis in macrophages. To examine how these enzymes regulate phagosome formation, the distributions of yellow fluorescent protein (YFP) chimeras of enzymes and pleckstrin homology (PH) domains specific for their substrates and products were analyzed quantitatively. PTEN-YFP did not localize to phagosomes, suggesting that PTEN regulates phagocytosis globally within the macrophage. SHIP1-YFP and p85-YFP were recruited to forming phagosomes. SHIP1-YFP sequestered to the leading edge and dissociated from phagocytic cups earlier than did p85-cyan fluorescent protein, indicating that SHIP-1 inhibitory activities are restricted to the early stages of phagocytosis. PH domain chimeras indicated that early during phagocytosis, PI(3,4,5)P(3) was slightly more abundant than PI(3,4)P(2) at the leading edge of the forming cup. These results support a model in which phagosomal PI3K generates PI(3,4,5)P(3) necessary for later stages of phagocytosis, PTEN determines whether those late stages can occur, and SHIP-1 regulates when and where they occur by transiently suppressing PI(3,4,5)P(3)-dependent activities necessary for completion of phagocytosis.