ABLs and TMKs are co-receptors for extracellular auxin.

Extracellular perception of auxin, an essential phytohormone in plants, has been debated for decades. Auxin-binding protein 1 (ABP1) physically interacts with quintessential transmembrane kinases (TMKs) and was proposed to act as an extracellular auxin receptor, but its role was disputed because abp1 knockout mutants lack obvious morphological phenotypes. Here, we identified two new auxin-binding proteins, ABL1 and ABL2, that are localized to the apoplast and directly interact with the extracellular domain of TMKs in an auxin-dependent manner. Furthermore, functionally redundant ABL1 and ABL2 genetically interact with TMKs and exhibit functions that overlap with those of ABP1 as well as being independent of ABP1. Importantly, the extracellular domain of TMK1 itself binds auxin and synergizes with either ABP1 or ABL1 in auxin binding. Thus, our findings discovered auxin receptors ABL1 and ABL2 having functions overlapping with but distinct from ABP1 and acting together with TMKs as co-receptors for extracellular auxin.

Coherent interfaces govern direct transformation from graphite to diamond.

Understanding the direct transformation from graphite to diamond has been a long-standing challenge with great scientific and practical importance. Previously proposed transformation mechanisms1-3, based on traditional experimental observations that lacked atomistic resolution, cannot account for the complex nanostructures occurring at graphite-diamond interfaces during the transformation4,5. Here we report the identification of coherent graphite-diamond interfaces, which consist of four basic structural motifs, in partially transformed graphite samples recovered from static compression, using high-angle annular dark-field scanning transmission electron microscopy. These observations provide insight into possible pathways of the transformation. Theoretical calculations confirm that transformation through these coherent interfaces is energetically favoured compared with those through other paths previously proposed1-3. The graphite-to-diamond transformation is governed by the formation of nanoscale coherent interfaces (diamond nucleation), which, under static compression, advance to consume the remaining graphite (diamond growth). These results may also shed light on transformation mechanisms of other carbon materials and boron nitride under different synthetic conditions.

Accurate and efficient integrative reference-informed spatial domain detection for spatial transcriptomics.

Spatially resolved transcriptomics (SRT) studies are becoming increasingly common and large, offering unprecedented opportunities in mapping complex tissue structures and functions. Here we present integrative and reference-informed tissue segmentation (IRIS), a computational method designed to characterize tissue spatial organization in SRT studies through accurately and efficiently detecting spatial domains. IRIS uniquely leverages single-cell RNA sequencing data for reference-informed detection of biologically interpretable spatial domains, integrating multiple SRT slices while explicitly considering correlations both within and across slices. We demonstrate the advantages of IRIS through in-depth analysis of six SRT datasets encompassing diverse technologies, tissues, species and resolutions. In these applications, IRIS achieves substantial accuracy gains (39-1,083%) and speed improvements (4.6-666.0) in moderate-sized datasets, while representing the only method applicable for large datasets including Stereo-seq and 10x Xenium. As a result, IRIS reveals intricate brain structures, uncovers tumor microenvironment heterogeneity and detects structural changes in diabetes-affected testis, all with exceptional speed and accuracy.

TMK-based cell-surface auxin signalling activates cell-wall acidification.

The phytohormone auxin controls many processes in plants, at least in part through its regulation of cell expansion1. The acid growth hypothesis has been proposed to explain auxin-stimulated cell expansion for five decades, but the mechanism that underlies auxin-induced cell-wall acidification is poorly characterized. Auxin induces the phosphorylation and activation of the plasma membrane H+-ATPase that pumps protons into the apoplast2, yet how auxin activates its phosphorylation remains unclear. Here we show that the transmembrane kinase (TMK) auxin-signalling proteins interact with plasma membrane H+-ATPases, inducing their phosphorylation, and thereby promoting cell-wall acidification and hypocotyl cell elongation in Arabidopsis. Auxin induced interactions between TMKs and H+-ATPases in the plasma membrane within seconds, as well as TMK-dependent phosphorylation of the penultimate threonine residue on the H+-ATPases. Our genetic, biochemical and molecular evidence demonstrates that TMKs directly phosphorylate plasma membrane H+-ATPase and are required for auxin-induced H+-ATPase activation, apoplastic acidification and cell expansion. Thus, our findings reveal a crucial connection between auxin and plasma membrane H+-ATPase activation in regulating apoplastic pH changes and cell expansion through TMK-based cell surface auxin signalling.

STAREG: Statistical replicability analysis of high throughput experiments with applications to spatial transcriptomic studies.

Replicable signals from different yet conceptually related studies provide stronger scientific evidence and more powerful inference. We introduce STAREG, a statistical method for replicability analysis of high throughput experiments, and apply it to analyze spatial transcriptomic studies. STAREG uses summary statistics from multiple studies of high throughput experiments and models the the joint distribution of p-values accounting for the heterogeneity of different studies. It effectively controls the false discovery rate (FDR) and has higher power by information borrowing. Moreover, it provides different rankings of important genes. With the EM algorithm in combination with pool-adjacent-violator-algorithm (PAVA), STAREG is scalable to datasets with millions of genes without any tuning parameters. Analyzing two pairs of spatially resolved transcriptomic datasets, we are able to make biological discoveries that otherwise cannot be obtained by using existing methods.

mtPGS: Leverage multiple correlated traits for accurate polygenic score construction.

Accurate polygenic scores (PGSs) facilitate the genetic prediction of complex traits and aid in the development of personalized medicine. Here, we develop a statistical method called multi-trait assisted PGS (mtPGS), which can construct accurate PGSs for a target trait of interest by leveraging multiple traits relevant to the target trait. Specifically, mtPGS borrows SNP effect size similarity information between the target trait and its relevant traits to improve the effect size estimation on the target trait, thus achieving accurate PGSs. In the process, mtPGS flexibly models the shared genetic architecture between the target and the relevant traits to achieve robust performance, while explicitly accounting for the environmental covariance among them to accommodate different study designs with various sample overlap patterns. In addition, mtPGS uses only summary statistics as input and relies on a deterministic algorithm with several algebraic techniques for scalable computation. We evaluate the performance of mtPGS through comprehensive simulations and applications to 25 traits in the UK Biobank, where in the real data mtPGS achieves an average of 0.90%-52.91% accuracy gain compared to the state-of-the-art PGS methods. Overall, mtPGS represents an accurate, fast, and robust solution for PGS construction in biobank-scale datasets.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution.

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Spatial transcriptomics reveals profound subclonal heterogeneity and T-cell dysfunction in extramedullary myeloma.

Extramedullary disease (EMD) is a high-risk feature of multiple myeloma (MM) and remains a poor prognostic factor even in the era of novel immunotherapies. Here we applied spatial transcriptomics (tomo-seq [n=2] and 10X Visium [n=12]), and single-cell RNA sequencing (scRNAseq [n=3]) to a set of 14 EMD biopsies to dissect the three-dimensional architecture of tumor cells and their microenvironment. Overall, the infiltrating immune and stromal cells showed both intra- and inter-patient variation with no uniform distribution over the lesion. We observed substantial heterogeneity at the copy number level within plasma cells, including the emergence of new subclones in circumscribed areas of the tumor, consistent with genomic instability. We further identified spatial expression differences of GPRC5D and TNFRSF17, two important antigens for bispecific antibody therapy. EMD masses were infiltrated by various immune cells, including T-cells. Notably, exhausted TIM3+/PD-1+ T-cells diffusely co-localized with MM cells, whereas functional and activated CD8+ T-cells showed a focal infiltration pattern along with M1 macrophages in otherwise tumor-free regions. This segregation of fit and exhausted T-cells was resolved in the case of response to T-cell engaging bispecific antibodies. MM cells and microenvironment cells were embedded in a complex network that influenced immune activation and angiogenesis, and oxidative phosphorylation represented the major metabolic program within EMD lesions. In summary, spatial transcriptomics has revealed a multicellular ecosystem in EMD with checkpoint inhibition and dual targeting as potential new therapeutic avenues.

Inflammasome Activation Triggers Caspase-1-Mediated Cleavage of cGAS to Regulate Responses to DNA Virus Infection.

Viral infection triggers host innate immune responses that result in the production of various cytokines including type I interferons (IFN), activation of inflammasomes, and programmed cell death of the infected cells. Tight control of inflammatory cytokine production is crucial for the triggering of an effective immune response that can resolve the infection without causing host pathology. In examining the inflammatory response of Asc-/- and Casp1-/- macrophages, we found that deficiency in these molecules resulted in increased IFN production upon DNA virus infection, but not RNA virus challenge. Investigation of the underlying mechanism revealed that upon canonical and non-canonical inflammasome activation, caspase-1 interacted with cyclic GMP-AMP (cGAMP) synthase (cGAS), cleaving it and dampening cGAS-STING-mediated IFN production. Deficiency in inflammasome signaling enhanced host resistance to DNA virus in vitro and in vivo, and this regulatory role extended to other inflammatory caspases. Thus, inflammasome activation dampens cGAS-dependent signaling, suggesting cross-regulation between intracellular DNA-sensing pathways.

Activation of STAT6 by STING is critical for antiviral innate immunity.

STAT6 plays a prominent role in adaptive immunity by transducing signals from extracellular cytokines. We now show that STAT6 is required for innate immune signaling in response to virus infection. Viruses or cytoplasmic nucleic acids trigger STING (also named MITA/ERIS) to recruit STAT6 to the endoplasmic reticulum, leading to STAT6 phosphorylation on Ser(407) by TBK1 and Tyr(641), independent of JAKs. Phosphorylated STAT6 then dimerizes and translocates to the nucleus to induce specific target genes responsible for immune cell homing. Virus-induced STAT6 activation is detected in all cell-types tested, in contrast to the cell-type specific role of STAT6 in cytokine signaling, and Stat6(-/-) mice are susceptible to virus infection. Thus, STAT6 mediates immune signaling in response to both cytokines at the plasma membrane, and virus infection at the endoplasmic reticulum.

Hierarchically structured diamond composite with exceptional toughness.

The well known trade-off between hardness and toughness (resistance to fracture) makes simultaneous improvement of both properties challenging, especially in diamond. The hardness of diamond can be increased through nanostructuring strategies1,2, among which the formation of high-density nanoscale twins - crystalline regions related by symmetry - also toughens diamond2. In materials other than diamond, there are several other promising approaches to enhancing toughness in addition to nanotwinning3, such as bio-inspired laminated composite toughening4-7, transformation toughening8 and dual-phase toughening9, but there has been little research into such approaches in diamond. Here we report the structural characterization of a diamond composite hierarchically assembled with coherently interfaced diamond polytypes (different stacking sequences), interwoven nanotwins and interlocked nanograins. The architecture of the composite enhances toughness more than nanotwinning alone, without sacrificing hardness. Single-edge notched beam tests yield a toughness up to five times that of synthetic diamond10, even greater than that of magnesium alloys. When fracture occurs, a crack propagates through diamond nanotwins of the 3C (cubic) polytype along {111} planes, via a zigzag path. As the crack encounters regions of non-3C polytypes, its propagation is diffused into sinuous fractures, with local transformation into 3C diamond near the fracture surfaces. Both processes dissipate strain energy, thereby enhancing toughness. This work could prove useful in making superhard materials and engineering ceramics. By using structural architecture with synergetic effects of hardening and toughening, the trade-off between hardness and toughness may eventually be surmounted.

Identification of G-quadruplex-interacting proteins in living cells using an artificial G4-targeting biotin ligase.

G-quadruplexes (G4s) are noncanonical nucleic acid structures pivotal to cellular processes and disease pathways. Deciphering G4-interacting proteins is imperative for unraveling G4's biological significance. In this study, we developed a G4-targeting biotin ligase named G4PID, meticulously assessing its binding affinity and specificity both in vitro and in vivo. Capitalizing on G4PID, we devised a tailored approach termed G-quadruplex-interacting proteins specific biotin-ligation procedure (PLGPB) to precisely profile G4-interacting proteins. Implementing this innovative strategy in live cells, we unveiled a cohort of 149 potential G4-interacting proteins, which exhibiting multifaceted functionalities. We then substantiate the directly binding affinity of 7 candidate G4-interacting-proteins (SF3B4, FBL, PP1G, BCL7C, NDUV1, ILF3, GAR1) in vitro. Remarkably, we verified that splicing factor 3B subunit 4 (SF3B4) binds preferentially to the G4-rich 3' splice site and the corresponding splicing sites are modulated by the G4 stabilizer PDS, indicating the regulating role of G4s in mRNA splicing procedure. The PLGPB strategy could biotinylate multiple proteins simultaneously, which providing an opportunity to map G4-interacting proteins network in living cells.

A bisulfite-assisted and ligation-based qPCR amplification technology for locus-specific pseudouridine detection at base resolution.

Over 150 types of chemical modifications have been identified in RNA to date, with pseudouridine (Ψ) being one of the most prevalent modifications in RNA. Ψ plays vital roles in various biological processes, and precise, base-resolution detection methods are fundamental for deep analysis of its distribution and function. In this study, we introduced a novel base-resolution Ψ detection method named pseU-TRACE. pseU-TRACE relied on the fact that RNA containing Ψ underwent a base deletion after treatment of bisulfite (BS) during reverse transcription, which enabled efficient ligation of two probes complementary to the cDNA sequence on either side of the Ψ site and successful amplification in subsequent real-time quantitative PCR (qPCR), thereby achieving selective and accurate Ψ detection. Our method accurately and sensitively detected several known Ψ sites in 28S, 18S, 5.8S, and even mRNA. Moreover, pseU-TRACE could be employed to measure the Ψ fraction in RNA and explore the Ψ metabolism of different pseudouridine synthases (PUSs), providing valuable insights into the function of Ψ. Overall, pseU-TRACE represents a reliable, time-efficient and sensitive Ψ detection method.

Catalytic-assembly of programmable atom equivalents.

Escape from metastable states in self-assembly of colloids is an intractable problem. Unlike the commonly adopted approach of thermal annealing, the recently developed enthalpy-mediated strategy provided a different option to address this dilemma in a dynamically controllable manner at room temperature. However, it required a complex catalytic-assembly DNA strand-displacement circuitry to mediate interaction between multiple components. In this work, we present a simple but effective way to achieve catalytic-assembly of DNA-functionalized colloidal nanoparticles, i.e., programmable atom equivalents, in a far-from-equilibrium system. A removable molecule named "catassembler" that acts as a catalyst was employed to rectify imperfect linkages and help the system escape from metastability without affecting the assembled framework. Notably, catalytic efficiency of the catassembler can be effectively improved by changing the seesaw catassembler in toehold length design or numbers of the repeat units. Leveraging this tractable catalytic-assembly approach, different ordered architectures were easily produced by directly mixing all reactants, as in chemical reactions. By switching bonding identities, solid-solid phase transformations between different colloidal crystals were achieved. This work opens up an avenue for programming colloid assembly in a far-from-equilibrium system.

Methods for mediation analysis with high-dimensional DNA methylation data: Possible choices and comparisons.

Epigenetic researchers often evaluate DNA methylation as a potential mediator of the effect of social/environmental exposures on a health outcome. Modern statistical methods for jointly evaluating many mediators have not been widely adopted. We compare seven methods for high-dimensional mediation analysis with continuous outcomes through both diverse simulations and analysis of DNAm data from a large multi-ethnic cohort in the United States, while providing an R package for their seamless implementation and adoption. Among the considered choices, the best-performing methods for detecting active mediators in simulations are the Bayesian sparse linear mixed model (BSLMM) and high-dimensional mediation analysis (HDMA); while the preferred methods for estimating the global mediation effect are high-dimensional linear mediation analysis (HILMA) and principal component mediation analysis (PCMA). We provide guidelines for epigenetic researchers on choosing the best method in practice and offer suggestions for future methodological development.

SPA, a Stigma-style-transmitting tract Physical microenvironment Assay for investigating mechano-signaling in pollen tubes.

Accurate sensing and responding to physical microenvironment are crucial for cell function and survival, but the underlying molecular mechanisms remain elusive. Pollen tube (PT) provides a perfect single-cell model for studying mechanobiology since it's naturally subjected to complex mechanical instructions from the pistil during invasive growth. Recent reports have revealed discrepant PT behaviors between in vivo and flat, two-dimensional in vitro cultures. Here, we established the Stigma-style-transmitting tract (TT) Physical microenvironment Assay (SPA) to recapitulate pressure changes in the pistil. This biomimetic assay has enabled us to swiftly identify highly redundant genes, GEF8/9/11/12/13, as new regulators for maintaining PTs integrity during style-to-TT emergence. In contrast to normal growth on solid medium, SPA successfully phenocopied gef8/9/11/12/13 PT in vivo growth-arrest deficiency. Our results suggest the existence of distinct signaling pathways regulating in vivo and in vitro PT integrity maintenance, underscoring the necessity of faithfully mimicking the physical microenvironment for studying plant cell biology.

ExPRSweb: An online repository with polygenic risk scores for common health-related exposures.

Complex traits are influenced by genetic risk factors, lifestyle, and environmental variables, so-called exposures. Some exposures, e.g., smoking or lipid levels, have common genetic modifiers identified in genome-wide association studies. Because measurements are often unfeasible, exposure polygenic risk scores (ExPRSs) offer an alternative to study the influence of exposures on various phenotypes. Here, we collected publicly available summary statistics for 28 exposures and applied four common PRS methods to generate ExPRSs in two large biobanks: the Michigan Genomics Initiative and the UK Biobank. We established ExPRSs for 27 exposures and demonstrated their applicability in phenome-wide association studies and as predictors for common chronic conditions. Especially the addition of multiple ExPRSs showed, for several chronic conditions, an improvement compared to prediction models that only included traditional, disease-focused PRSs. To facilitate follow-up studies, we share all ExPRS constructs and generated results via an online repository called ExPRSweb.

METRO: Multi-ancestry transcriptome-wide association studies for powerful gene-trait association detection.

Integrative analysis of genome-wide association studies (GWASs) and gene expression studies in the form of a transcriptome-wide association study (TWAS) has the potential to better elucidate the molecular mechanisms underlying disease etiology. Here we present a method, METRO, that can leverage gene expression data collected from multiple genetic ancestries to enhance TWASs. METRO incorporates expression prediction models constructed in different genetic ancestries through a likelihood-based inference framework, producing calibrated p values with substantially improved TWAS power. We illustrate the benefits of METRO in both simulations and applications to seven complex traits and diseases obtained from four GWASs. These GWASs include two of primarily European ancestry (n = 188,577 and 339,226) and two of primarily African ancestry (n = 42,752 and 23,827). In the real data applications, we leverage gene expression data measured on 1,032 African Americans and 801 European Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study to identify a substantially larger number of gene-trait associations as compared to existing TWAS approaches. The benefits of METRO are most prominent in applications to GWASs of African ancestry where the sample size is much smaller than GWASs of European ancestry and where a more powerful TWAS method is crucial. Among the identified associations are high-density lipoprotein-associated genes including PLTP and PPARG that are critical for maintaining lipid homeostasis and the type II diabetes-associated gene MAPT that supports microtubule-associated protein tau as a key component underlying impaired insulin secretion.

Mutant p53 Gains Its Function via c-Myc Activation upon CDK4 Phosphorylation at Serine 249 and Consequent PIN1 Binding.

TP53 missense mutations significantly influence the development and progression of various human cancers via their gain of new functions (GOF) through different mechanisms. Here we report a unique mechanism underlying the GOF of p53-R249S (p53-RS), a p53 mutant frequently detected in human hepatocellular carcinoma (HCC) that is highly related to hepatitis B infection and aflatoxin B1. A CDK inhibitor blocks p53-RS's nuclear translocation in HCC, whereas CDK4 interacts with p53-RS in the G1/S phase of the cells, phosphorylates it, and enhances its nuclear localization. This is coupled with binding of a peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (PIN1) to p53-RS, but not the p53 form with mutations of four serines/threonines previously shown to be crucial for PIN1 binding. As a result, p53-RS interacts with c-Myc and enhances c-Myc-dependent rDNA transcription key for ribosomal biogenesis. These results unveil a CDK4-PIN1-p53-RS-c-Myc pathway as a novel mechanism for the GOF of p53-RS in HCC.

Transcriptome-wide profiling of A-to-I RNA editing by Slic-seq.

Adenosine-to-inosine (A-to-I) RNA editing is a post-transcriptional processing event involved in diversifying the transcriptome and is responsible for various biological processes. In this context, we developed a new method based on the highly selective cleavage activity of Endonuclease V against Inosine and the universal activity of sodium periodate against all RNAs to enrich the inosine-containing RNA and accurately identify the editing sites. We validated the reliability of our method in human brain in both Alu and non-Alu elements. The conserved sites of A-to-I editing in human cells (HEK293T, HeLa, HepG2, K562 and MCF-7) primarily occurs in the 3'UTR of the RNA, which are highly correlated with RNA binding and protein binding. Analysis of the editing sites between the human brain and mouse brain revealed that the editing of exons is more conserved than that in other regions. This method was applied to three neurological diseases (Alzheimer's, epilepsy and ageing) of mouse brain, reflecting that A-to-I editing sites significantly decreased in neuronal activity genes.