Large-scale cranial window for in vivo mouse brain imaging utilizing fluoropolymer nanosheet and light-curable resin

Two-photon microscopy enables in vivo imaging of neuronal activity in mammalian brains at high resolution. However, two-photon imaging tools for stable, long-term, and simultaneous study of multiple brain regions in same mice are lacking. Here, we propose a method to create large cranial windows covering such as the whole parietal cortex and cerebellum in mice using fluoropolymer nanosheets covered with light-curable resin (termed the 'Nanosheet Incorporated into light-curable REsin' or NIRE method). NIRE method can produce cranial windows conforming the curved cortical and cerebellar surfaces, without motion artifacts in awake mice, and maintain transparency for >5 months. In addition, we demonstrate that NIRE method can be used for in vivo two-photon imaging of neuronal ensembles, individual neurons and subcellular structures such as dendritic spines. The NIRE method can facilitate in vivo large-scale analysis of heretofore inaccessible neural processes, such as the neuroplastic changes associated with maturation, learning and neural pathogenesis.

Advances in cancer DNA methylation analysis with methPLIER: use of non-negative matrix factorization and knowledge-based constraints to enhance biological interpretability

DNA methylation is an epigenetic modification that results in dynamic changes during ontogenesis and cell differentiation. DNA methylation patterns regulate gene expression and have been widely researched. While tools for DNA methylation analysis have been developed, most of them have focused on intergroup comparative analysis within a dataset; therefore, it is difficult to conduct cross-dataset studies, such as rare disease studies or cross-institutional studies. This study describes a novel method for DNA methylation analysis, namely, methPLIER, which enables interdataset comparative analyses. methPLIER combines Pathway Level Information Extractor (PLIER), which is a non-negative matrix factorization (NMF) method, with regularization by a knowledge matrix and transfer learning. methPLIER can be used to perform intersample and interdataset comparative analysis based on latent feature matrices, which are obtained via matrix factorization of large-scale data, and factor-loading matrices, which are obtained through matrix factorization of the data to be analyzed. We used methPLIER to analyze a lung cancer dataset and confirmed that the data decomposition reflected sample characteristics for recurrence-free survival. Moreover, methPLIER can analyze data obtained via different preprocessing methods, thereby reducing distributional bias among datasets due to preprocessing. Furthermore, methPLIER can be employed for comparative analyses of methylation data obtained from different platforms, thereby reducing bias in data distribution due to platform differences. methPLIER is expected to facilitate cross-sectional DNA methylation data analysis and enhance DNA methylation data resources.

Superconductivity in a van der Waals layered quasicrystal

Van der Waals layered transition-metal chalcogenides are drawing significant attention owing to their intriguing physical properties. This group of materials consists of abundant members with various elements, having a variety of different structures. However, they are all crystalline materials, and the physical properties of van der Waals layered quasicrystals have never been studied to date. Here, we report on the discovery of superconductivity in a van der Waals layered quasicrystal of Ta1.6Te. The electrical resistivity, magnetic susceptibility, and specific heat of the quasicrystal unambiguously validate the occurrence of bulk superconductivity at a transition temperature of ~1 K. This discovery can promote new research on assessing the physical properties of novel van der Waals layered quasicrystals as well as two-dimensional quasicrystals; moreover, it paves the way toward new frontiers of superconductivity in thermodynamically stable quasicrystals.

Deep generative design of RNA family sequences

RNA engineering has immense potential to drive innovation in biotechnology and medicine. Despite its importance, a versatile platform for the automated design of functional RNA is still lacking. Here, we propose RNA family sequence generator (RfamGen), a deep generative model that designs RNA family sequences in a data-efficient manner by explicitly incorporating alignment and consensus secondary structure information. RfamGen can generate novel and functional RNA family sequences by sampling points from a semantically rich and continuous representation. We have experimentally demonstrated the versatility of RfamGen using diverse RNA families. Furthermore, we confirmed the high success rate of RfamGen in designing functional ribozymes through a quantitative massively parallel assay. Notably, RfamGen successfully generates artificial sequences with higher activity than natural sequences. Overall, RfamGen significantly improves our ability to design functional RNA and opens up new potential for generative RNA engineering in synthetic biology.

Intrinsic factors responsible for brittle versus ductile nature of refractory high-entropy alloys

Refractory high-entropy alloys (RHEAs) are of interest for ultrahigh-temperature applications. To overcome their drawbacks - low-temperature brittleness and poor creep strength at high temperatures - improved fundamental understanding is needed. Using experiments, theory, and modeling, we investigated prototypical body-centered cubic (BCC) RHEAs, TiZrHfNbTa and VNbMoTaW. The former is compressible to 77 K, whereas the latter is not below 298 K. Hexagonal close-packed (HCP) elements in TiZrHfNbTa lower its dislocation core energy, increase lattice distortion, and lower its shear modulus relative to VNbMoTaW whose elements are all BCC. Screw dislocations dominate TiZrHfNbTa plasticity, but equal numbers of edges and screws exist in VNbTaMoW. Dislocation cores are compact in VNbTaMoW and extended in TiZrHfNbTa, and different macroscopic slip planes are activated in the two RHEAs, which we attribute to the concentration of HCP elements. Our findings demonstrate how ductility and strength can be controlled through the ratio of HCP to BCC elements in RHEAs.

Transomics2cytoscape: an automated software for interpretable 2.5-dimensional visualization of trans-omic networks

Biochemical network visualization is one of the essential technologies for mechanistic interpretation of omics data. In particular, recent advances in multi-omics measurement and analysis require the development of visualization methods that encompass multiple omics data. Visualization in 2.5 dimension (2.5D visualization), which is an isometric view of stacked X-Y planes, is a convenient way to interpret multi-omics/trans-omics data in the context of the conventional layouts of biochemical networks drawn on each of the stacked omics layers. However, 2.5D visualization of trans-omics networks is a state-of-the-art method that primarily relies on time-consuming human efforts involving manual drawing. Here, we present an R Bioconductor package 'transomics2cytoscape' for automated visualization of 2.5D trans-omics networks. We confirmed that transomics2cytoscape could be used for rapid visualization of trans-omics networks presented in published papers within a few minutes. Transomics2cytoscape allows for frequent update/redrawing of trans-omics networks in line with the progress in multi-omics/trans-omics data analysis, thereby enabling network-based interpretation of multi-omics data at each research step. The transomics2cytoscape source code is available at https://github.com/ecell/transomics2cytoscape .

In vitro induction of patterned branchial arch-like aggregate from human pluripotent stem cells

Early patterning of neural crest cells (NCCs) in the craniofacial primordium is important for subsequent development of proper craniofacial structures. However, because of the complexity of the environment of developing tissues, surveying the early specification and patterning of NCCs is difficult. In this study, we develop a simplified in vitro 3D model using human pluripotent stem cells to analyze the early stages of facial development. In this model, cranial NCC-like cells spontaneously differentiate from neural plate border-like cells into maxillary arch-like mesenchyme after a long-term culture. Upon the addition of EDN1 and BMP4, these aggregates are converted into a mandibular arch-like state. Furthermore, temporary treatment with EDN1 and BMP4 induces the formation of spatially separated domains expressing mandibular and maxillary arch markers within a single aggregate. These results suggest that this in vitro model is useful for determining the mechanisms underlying cell fate specification and patterning during early facial development.

Ran-GTP assembles a specialized spindle structure for accurate chromosome segregation in medaka early embryos

Despite drastic cellular changes during cleavage, a mitotic spindle assembles in each blastomere to accurately segregate duplicated chromosomes. Mechanisms of mitotic spindle assembly have been extensively studied using small somatic cells. However, mechanisms of spindle assembly in large vertebrate embryos remain little understood. Here, we establish functional assay systems in medaka (Oryzias latipes) embryos by combining CRISPR knock-in with auxin-inducible degron technology. Live imaging reveals several unexpected features of microtubule organization and centrosome positioning that achieve rapid, accurate cleavage. Importantly, Ran-GTP assembles a dense microtubule network at the metaphase spindle center that is essential for chromosome segregation in early embryos. This unique spindle structure is remodeled into a typical short, somatic-like spindle after blastula stages, when Ran-GTP becomes dispensable for chromosome segregation. We propose that despite the presence of centrosomes, the chromosome-derived Ran-GTP pathway has essential roles in functional spindle assembly in large, rapidly dividing vertebrate early embryos, similar to acentrosomal spindle assembly in oocytes.

Neuronal DSCAM regulates the peri-synaptic localization of GLAST in Bergmann glia for functional synapse formation

In the central nervous system, astrocytes enable appropriate synapse function through glutamate clearance from the synaptic cleft; however, it remains unclear how astrocytic glutamate transporters function at peri-synaptic contact. Here, we report that Down syndrome cell adhesion molecule (DSCAM) in Purkinje cells controls synapse formation and function in the developing cerebellum. Dscam-mutant mice show defects in CF synapse translocation as is observed in loss of function mutations in the astrocytic glutamate transporter GLAST expressed in Bergmann glia. These mice show impaired glutamate clearance and the delocalization of GLAST away from the cleft of parallel fibre (PF) synapse. GLAST complexes with the extracellular domain of DSCAM. Riluzole, as an activator of GLAST-mediated uptake, rescues the proximal impairment in CF synapse formation in Purkinje cell-selective Dscam-deficient mice. DSCAM is required for motor learning, but not gross motor coordination. In conclusion, the intercellular association of synaptic and astrocyte proteins is important for synapse formation and function in neural transmission.

Coherent electronic coupling in quantum dot solids induces cooperative enhancement of nonlinear optoelectronic responses

Synchronized dynamics of quantum dot (QD) ensembles are essential for generating ultrafast and giant optical responses beyond those of individual QDs. Increasing the strength of the direct electronic coupling between QDs is a key strategy for the realization of cooperative quantum phenomena. Here, we observe a quantum cooperative effect on nonlinear photocurrents caused by the coherent electronic coupling in semiconductor QD solids. We measure quantum interference signals cooperatively generated in QD solids. We control the inter-QD distance with atomic precision using bidentate ligands that strongly link the QDs. The harmonic quantum interference signals are strongly enhanced when shortening the molecular length of the ligand. Furthermore, we clarify that the coherence length of multiexcitons extends to neighbouring QDs. This finding is direct evidence that multiexciton coherent tunnelling assists the ultrafast exciton delocalization. Cooperative enhancement in QD solids may find application in advanced quantum optoelectronics.

Neuronal sensorimotor integration guiding salt concentration navigation in Caenorhabditis elegans

Animals navigate their environment by manipulating their movements and adjusting their trajectory which requires a sophisticated integration of sensory data with their current motor status. Here, we utilize the nematode Caenorhabditis elegans to explore the neural mechanisms of processing the sensory and motor information for navigation. We developed a microfluidic device which allows animals to freely move their heads while receiving temporal NaCl stimuli. We found that C. elegans regulates neck bending direction in response to temporal NaCl concentration changes in a way which is consistent with a C. elegans' navigational strategy which regulates traveling direction toward preferred NaCl concentrations. Our analysis also revealed that the activity of a neck motor neuron is significantly correlated with neck bending and activated by the decrease in NaCl concentration in a phase-dependent manner. By combining the analysis of behavioral and neural response to NaCl stimuli and optogenetic perturbation experiments, we revealed that NaCl decrease during ventral bending activates the neck motor neuron which counteracts ipsilateral bending. Simulations further suggest that this phase-dependent response of neck motor neurons can facilitate curving toward preferred salt concentrations.

Alkali metal bilayer intercalation in graphene

Alkali metal (AM) intercalation between graphene layers holds promise for electronic manipulation and energy storage, yet the underlying mechanism remains challenging to fully comprehend despite extensive research. In this study, we employ low-voltage scanning transmission electron microscopy (LV-STEM) to visualize the atomic structure of intercalated AMs (potassium, rubidium, and cesium) in bilayer graphene (BLG). Our findings reveal that the intercalated AMs adopt bilayer structures with hcp stacking, and specifically a C6M2C6 composition. These structures closely resemble the bilayer form of fcc (111) structure observed in AMs under high-pressure conditions. A negative charge transferred from bilayer AMs to graphene layers of approximately 1~1.5×1014 e-/cm-2 was determined by electron energy loss spectroscopy (EELS), Raman, and electrical transport. The bilayer AM is stable in BLG and graphite superficial layers but absent in the graphite interior, primarily dominated by single-layer AM intercalation. This hints at enhancing AM intercalation capacity by thinning the graphite material.

The nucleolar shell provides anchoring sites for DNA untwisting

DNA underwinding (untwisting) is a crucial step in transcriptional activation. DNA underwinding occurs between the site where torque is generated by RNA polymerase (RNAP) and the site where the axial rotation of DNA is constrained. However, what constrains DNA axial rotation in the nucleus is yet unknown. Here, we show that the anchorage to the nuclear protein condensates constrains DNA axial rotation for DNA underwinding in the nucleolus. In situ super-resolution imaging of underwound DNA reveal that underwound DNA accumulates in the nucleolus, a nuclear condensate with a core-shell structure. Specifically, underwound DNA is distributed in the nucleolar core owing to RNA polymerase I (RNAPI) activities. Furthermore, underwound DNA in the core decreases when nucleolar shell components are prevented from binding to their recognition structure, G-quadruplex (G4). Taken together, these results suggest that the nucleolar shell provides anchoring sites that constrain DNA axial rotation for RNAPI-driven DNA underwinding in the core. Our findings will contribute to understanding how nuclear protein condensates make up constraints for the site-specific regulation of DNA underwinding and transcription.

Replicating shear-mediated self-assembly of spider silk through microfluidics

The development of artificial spider silk with properties similar to native silk has been a challenging task in materials science. In this study, we use a microfluidic device to create continuous fibers based on recombinant MaSp2 spidroin. The strategy incorporates ion-induced liquid-liquid phase separation, pH-driven fibrillation, and shear-dependent induction of β-sheet formation. We find that a threshold shear stress of approximately 72 Pa is required for fiber formation, and that β-sheet formation is dependent on the presence of polyalanine blocks in the repetitive sequence. The MaSp2 fiber formed has a β-sheet content (29.2%) comparable to that of native dragline with a shear stress requirement of 111 Pa. Interestingly, the polyalanine blocks have limited influence on the occurrence of liquid-liquid phase separation and hierarchical structure. These results offer insights into the shear-induced crystallization and sequence-structure relationship of spider silk and have significant implications for the rational design of artificially spun fibers.

Automatic feature engineering for catalyst design using small data without prior knowledge of target catalysis

The empirical aspect of descriptor design in catalyst informatics, particularly when confronted with limited data, necessitates adequate prior knowledge for delving into unknown territories, thus presenting a logical contradiction. This study introduces a technique for automatic feature engineering (AFE) that works on small catalyst datasets, without reliance on specific assumptions or pre-existing knowledge about the target catalysis when designing descriptors and building machine-learning models. This technique generates numerous features through mathematical operations on general physicochemical features of catalytic components and extracts relevant features for the desired catalysis, essentially screening numerous hypotheses on a machine. AFE yields reasonable regression results for three types of heterogeneous catalysis: oxidative coupling of methane (OCM), conversion of ethanol to butadiene, and three-way catalysis, where only the training set is swapped. Moreover, through the application of active learning that combines AFE and high-throughput experimentation for OCM, we successfully visualize the machine's process of acquiring precise recognition of the catalyst design. Thus, AFE is a versatile technique for data-driven catalysis research and a key step towards fully automated catalyst discoveries.

HKDC1, a target of TFEB, is essential to maintain both mitochondrial and lysosomal homeostasis, preventing cellular senescence

Mitochondrial and lysosomal functions are intimately linked and are critical for cellular homeostasis, as evidenced by the fact that cellular senescence, aging, and multiple prominent diseases are associated with concomitant dysfunction of both organelles. However, it is not well understood how the two important organelles are regulated. Transcription factor EB (TFEB) is the master regulator of lysosomal function and is also implicated in regulating mitochondrial function; however, the mechanism underlying the maintenance of both organelles remains to be fully elucidated. Here, by comprehensive transcriptome analysis and subsequent chromatin immunoprecipitation-qPCR, we identified hexokinase domain containing 1 (HKDC1), which is known to function in the glycolysis pathway as a direct TFEB target. Moreover, HKDC1 was upregulated in both mitochondrial and lysosomal stress in a TFEB-dependent manner, and its function was critical for the maintenance of both organelles under stress conditions. Mechanistically, the TFEB-HKDC1 axis was essential for PINK1 (PTEN-induced kinase 1)/Parkin-dependent mitophagy via its initial step, PINK1 stabilization. In addition, the functions of HKDC1 and voltage-dependent anion channels, with which HKDC1 interacts, were essential for the clearance of damaged lysosomes and maintaining mitochondria-lysosome contact. Interestingly, HKDC1 regulated mitophagy and lysosomal repair independently of its prospective function in glycolysis. Furthermore, loss function of HKDC1 accelerated DNA damage-induced cellular senescence with the accumulation of hyperfused mitochondria and damaged lysosomes. Our results show that HKDC1, a factor downstream of TFEB, maintains both mitochondrial and lysosomal homeostasis, which is critical to prevent cellular senescence.

Lysophagy protects against propagation of α-synuclein aggregation through ruptured lysosomal vesicles

The neuron-to-neuron propagation of misfolded α-synuclein (αSyn) aggregates is thought to be key to the pathogenesis of synucleinopathies. Recent studies have shown that extracellular αSyn aggregates taken up by the endosomal-lysosomal system can rupture the lysosomal vesicular membrane; however, it remains unclear whether lysosomal rupture leads to the transmission of αSyn aggregation. Here, we applied cell-based αSyn propagation models to show that ruptured lysosomes are the pathway through which exogenous αSyn aggregates transmit aggregation, and furthermore, this process was prevented by lysophagy, i.e., selective autophagy of damaged lysosomes. αSyn aggregates accumulated predominantly in lysosomes, causing their rupture, and seeded the aggregation of endogenous αSyn, initially around damaged lysosomes. Exogenous αSyn aggregates induced the accumulation of LC3 on lysosomes. This LC3 accumulation was not observed in cells in which a key regulator of autophagy, RB1CC1/FIP200, was knocked out and was confirmed as lysophagy by transmission electron microscopy. Importantly, RB1CC1/FIP200-deficient cells treated with αSyn aggregates had increased numbers of ruptured lysosomes and enhanced propagation of αSyn aggregation. Furthermore, various types of lysosomal damage induced using lysosomotropic reagents, depletion of lysosomal enzymes, or more toxic species of αSyn fibrils also exacerbated the propagation of αSyn aggregation, and impaired lysophagy and lysosomal membrane damage synergistically enhanced propagation. These results indicate that lysophagy prevents exogenous αSyn aggregates from escaping the endosomal-lysosomal system and transmitting aggregation to endogenous cytosolic αSyn via ruptured lysosomal vesicles. Our findings suggest that the progression and severity of synucleinopathies are associated with damage to lysosomal membranes and impaired lysophagy.

Structural insights into the complex of oncogenic KRas4BG12V and Rgl2, a RalA/B activator

About a quarter of total human cancers carry mutations in Ras isoforms. Accumulating evidence suggests that small GTPases, RalA, and RalB, and their activators, Ral guanine nucleotide exchange factors (RalGEFs), play an essential role in oncogenic Ras-induced signalling. We studied the interaction between human KRas4B and the Ras association (RA) domain of Rgl2 (Rgl2RA), one of the RA-containing RalGEFs. We show that the G12V oncogenic KRas4B mutation changes the interaction kinetics with Rgl2RA The crystal structure of the KRas4BG12V: Rgl2RA complex shows a 2:2 heterotetramer where the switch I and switch II regions of each KRasG12V interact with both Rgl2RA molecules. This structural arrangement is highly similar to the HRasE31K:RALGDSRA crystal structure and is distinct from the well-characterised Ras:Raf complex. Interestingly, the G12V mutation was found at the dimer interface of KRas4BG12V with its partner. Our study reveals a potentially distinct mode of Ras:effector complex formation by RalGEFs and offers a possible mechanistic explanation for how the oncogenic KRas4BG12V hyperactivates the RalA/B pathway.

RENGE infers gene regulatory networks using time-series single-cell RNA-seq data with CRISPR perturbations

Single-cell RNA-seq analysis coupled with CRISPR-based perturbation has enabled the inference of gene regulatory networks with causal relationships. However, a snapshot of single-cell CRISPR data may not lead to an accurate inference, since a gene knockout can influence multi-layered downstream over time. Here, we developed RENGE, a computational method that infers gene regulatory networks using a time-series single-cell CRISPR dataset. RENGE models the propagation process of the effects elicited by a gene knockout on its regulatory network. It can distinguish between direct and indirect regulations, which allows for the inference of regulations by genes that are not knocked out. RENGE therefore outperforms current methods in the accuracy of inferring gene regulatory networks. When used on a dataset we derived from human-induced pluripotent stem cells, RENGE yielded a network consistent with multiple databases and literature. Accurate inference of gene regulatory networks by RENGE would enable the identification of key factors for various biological systems.

An archetype and scaling of developmental tissue dynamics across species

Morphometric studies have revealed the existence of simple geometric relationships among various animal shapes. However, we have little knowledge of the mathematical principles behind the morphogenetic dynamics that form the organ/body shapes of different species. Here, we address this issue by focusing on limb morphogenesis in Gallus gallus domesticus (chicken) and Xenopus laevis (African clawed frog). To compare the deformation dynamics between tissues with different sizes/shapes as well as their developmental rates, we introduce a species-specific rescaled spatial coordinate and a common clock necessary for cross-species synchronization of developmental times. We find that tissue dynamics are well conserved across species under this spacetime coordinate system, at least from the early stages of development through the phase when basic digit patterning is established. For this developmental period, we also reveal that the tissue dynamics of both species are mapped with each other through a time-variant linear transformation in real physical space, from which hypotheses on a species-independent archetype of tissue dynamics and morphogenetic scaling are proposed.

Structural and thermodynamic insights into antibody light chain tetramer formation through 3D domain swapping

Overexpression of antibody light chains in small plasma cell clones can lead to misfolding and aggregation. On the other hand, the formation of amyloid fibrils from antibody light chains is related to amyloidosis. Although aggregation of antibody light chain is an important issue, atomic-level structural examinations of antibody light chain aggregates are sparse. In this study, we present an antibody light chain that maintains an equilibrium between its monomeric and tetrameric states. According to data from X-ray crystallography, thermodynamic and kinetic measurements, as well as theoretical studies, this antibody light chain engages in 3D domain swapping within its variable region. Here, a pair of domain-swapped dimers creates a tetramer through hydrophobic interactions, facilitating the revelation of the domain-swapped structure. The negative cotton effect linked to the β-sheet structure, observed around 215 nm in the circular dichroism (CD) spectrum of the tetrameric variable region, is more pronounced than that of the monomer. This suggests that the monomer contains less β-sheet structures and exhibits greater flexibility than the tetramer in solution. These findings not only clarify the domain-swapped structure of the antibody light chain but also contribute to controlling antibody quality and advancing the development of future molecular recognition agents and drugs.

Microautophagy regulated by STK38 and GABARAPs is essential to repair lysosomes and prevent aging

Lysosomes are degradative organelles and signaling hubs that maintain cell and tissue homeostasis, and lysosomal dysfunction is implicated in aging and reduced longevity. Lysosomes are frequently damaged, but their repair mechanisms remain unclear. Here, we demonstrate that damaged lysosomal membranes are repaired by microautophagy (a process termed "microlysophagy") and identify key regulators of the first and last steps. We reveal the AGC kinase STK38 as a novel microlysophagy regulator. Through phosphorylation of the scaffold protein DOK1, STK38 is specifically required for the lysosomal recruitment of the AAA+ ATPase VPS4, which terminates microlysophagy by promoting the disassembly of ESCRT components. By contrast, microlysophagy initiation involves non-canonical lipidation of ATG8s, especially the GABARAP subfamily, which is required for ESCRT assembly through interaction with ALIX. Depletion of STK38 and GABARAPs accelerates DNA damage-induced cellular senescence in human cells and curtails lifespan in C. elegans, respectively. Thus, microlysophagy is regulated by STK38 and GABARAPs and could be essential for maintaining lysosomal integrity and preventing aging.

Secretory GFP reconstitution labeling of neighboring cells interrogates cell-cell interactions in metastatic niches

Cancer cells inevitably interact with neighboring host tissue-resident cells during the process of metastatic colonization, establishing a metastatic niche to fuel their survival, growth, and invasion. However, the underlying mechanisms in the metastatic niche are yet to be fully elucidated owing to the lack of methodologies for comprehensively studying the mechanisms of cell-cell interactions in the niche. Here, we improve a split green fluorescent protein (GFP)-based genetically encoded system to develop secretory glycosylphosphatidylinositol-anchored reconstitution-activated proteins to highlight intercellular connections (sGRAPHIC) for efficient fluorescent labeling of tissue-resident cells that neighbor on and putatively interact with cancer cells in deep tissues. The sGRAPHIC system enables the isolation of metastatic niche-associated tissue-resident cells for their characterization using a single-cell RNA sequencing platform. We use this sGRAPHIC-leveraged transcriptomic platform to uncover gene expression patterns in metastatic niche-associated hepatocytes in a murine model of liver metastasis. Among the marker genes of metastatic niche-associated hepatocytes, we identify Lgals3, encoding galectin-3, as a potential pro-metastatic factor that accelerates metastatic growth and invasion.

Periportal hepatocyte proliferation at midgestation governs maternal glucose homeostasis in mice

The maternal liver is challenged by metabolic demands throughout pregnancy. However, hepatocyte dynamics and their physiological significance in pregnancy remain unclear. Here, we show in mice that hepatocyte proliferation is spatiotemporally regulated in each liver lobular zone during pregnancy, with transient proliferation of periportal and pericentral hepatocytes during mid and late gestation, respectively. Using adeno-associated virus (AAV)-8-mediated expression of the cell cycle inhibitor p21 in hepatocytes, we show that inhibition of hepatocyte proliferation during mid, but not late, gestation impairs liver growth. Transcriptionally, genes involved in glucose/glycogen metabolism are downregulated in late pregnancy when midgestational hepatocyte proliferation is attenuated. In addition, hepatic glycogen storage is abolished, with concomitant elevated blood glucose concentrations, glucose intolerance, placental glycogen deposition, and fetal overgrowth. Laser capture microdissection and RNA-seq analysis of each liver lobular zone show zone-specific changes in the transcriptome during pregnancy and identify genes that are periportally expressed at midgestation, including the hyaluronan-mediated motility receptor (Hmmr). Knockdown of Hmmr in hepatocytes by AAV8-shHmmr suppresses periportal hepatocyte proliferation at midgestation and induces impaired hepatic glycogen storage, glucose intolerance, placental glycogen deposition and fetal overgrowth. Our results suggest that periportal hepatocyte proliferation during midgestation is critical for maternal glycogen metabolism and fetal size.