Anterior-posterior (A-P) elongation of the palate is a critical aspect of integrated midfacial morphogenesis. Reciprocal epithelial-mesenchymal interactions drive secondary palate elongation that is coupled to the periodic formation of signaling centers within the rugae growth zone (RGZ). However, the relationship between RGZ-driven morphogenetic processes, the differentiative dynamics of underlying palatal bone mesenchymal precursors, and the segmental organization of the upper jaw has remained enigmatic. A detailed ontogenetic study of these relationships is important because palatal segment growth is a critical aspect of normal midfacial growth, can produce dysmorphology when altered, and is a likely basis for evolutionary differences in upper jaw morphology. We completed a combined whole mount gene expression and morphometric analysis of normal murine palatal segment growth dynamics and resulting upper jaw morphology. Our results demonstrated that the first formed palatal ruga (ruga 1), found just posterior to the RGZ, maintained an association with important nasal, neurovascular and palatal structures throughout early midfacial development. This suggested that these features are positioned at a proximal source of embryonic midfacial directional growth. Our detailed characterization of midfacial morphogenesis revealed a one-to-one relationship between palatal segments and upper jaw bones during the earliest stages of palatal elongation. Growth of the maxillary anlage within the anterior secondary palate is uniquely coupled to RGZ-driven morphogenesis. This may help drive the unequaled proportional elongation of the anterior secondary palate segment prior to palatal shelf fusion. Our results also demonstrated that the future maxillary-palatine suture, approximated by the position of ruga 1 and consistently associated with the palatine anlage, formed predominantly via the posterior differentiation of the maxilla within the expanding anterior secondary palate. Our ontogenetic analysis provides a novel and detailed picture of the earliest spatiotemporal dynamics of intramembranous midfacial skeletal specification and differentiation within the context of the surrounding palatal segment A-P elongation and associated rugae formation.

Otolith organs in the inner ear and neuromasts in the fish lateral-line harbor two populations of hair cells oriented to detect stimuli in opposing directions. The underlying mechanism is highly conserved: the transcription factor EMX2 is regionally expressed in just one hair cell population and acts through the receptor GPR156 to reverse cell orientation relative to the other population. In mouse and zebrafish, loss of Emx2 results in sensory organs that harbor only one hair cell orientation and are not innervated properly. In zebrafish, Emx2 also confers hair cells with reduced mechanosensory properties. Here, we leverage mouse and zebrafish models lacking GPR156 to determine how detecting stimuli of opposing directions serves vestibular function, and whether GPR156 has other roles besides orienting hair cells. We find that otolith organs in Gpr156 mouse mutants have normal zonal organization and normal type I-II hair cell distribution and mechano-electrical transduction properties. In contrast, gpr156 zebrafish mutants lack the smaller mechanically evoked signals that characterize Emx2-positive hair cells. Loss of GPR156 does not affect orientation-selectivity of afferents in mouse utricle or zebrafish neuromasts. Consistent with normal otolith organ anatomy and afferent selectivity, Gpr156 mutant mice do not show overt vestibular dysfunction. Instead, performance on two tests that engage otolith organs is significantly altered – swimming and off-vertical-axis rotation. We conclude that GPR156 relays hair cell orientation and transduction information downstream of EMX2, but not selectivity for direction-specific afferents. These results clarify how molecular mechanisms that confer bi-directionality to sensory organs contribute to function, from single hair cell physiology to animal behavior.

Schlemm's canal (SC) is central in intraocular pressure regulation but requires much characterization. It has distinct inner and outer walls, each composed of Schlemm's canal endothelial cells (SECs) with different morphologies and functions. Recent transcriptomic studies of the anterior segment added important knowledge, but were limited in power by SEC numbers or did not focus on SC. To gain a more comprehensive understanding of SC biology, we performed bulk RNA sequencing on C57BL/6J SC, blood vessel, and lymphatic endothelial cells from limbal tissue (~4500 SECs). We also analyzed mouse limbal tissues by single-cell and single-nucleus RNA sequencing (C57BL/6J and 129/Sj strains), successfully sequencing 903 individual SECs. Together, these datasets confirm that SC has molecular characteristics of both blood and lymphatic endothelia with a lymphatic phenotype predominating. SECs are enriched in pathways that regulate cell-cell junction formation pointing to the importance of junctions in determining SC fluid permeability. Importantly, and for the first time, our analyses characterize 3 molecular classes of SECs, molecularly distinguishing inner wall from outer wall SECs and discovering two inner wall cell states that likely result from local environmental differences. Further, and based on ligand and receptor expression patterns, we document key interactions between SECs and cells of the adjacent trabecular meshwork (TM) drainage tissue. Also, we present cell type expression for a collection of human glaucoma genes. These data provide a new molecular foundation that will enable the functional dissection of key homeostatic processes mediated by SECs as well as the development of new glaucoma therapeutics.

Blood donor age and sex affect glycolysis in stored RBCs from 13,029 volunteers;Ancestry, genetic polymorphisms in PFKP, HK1, CD38/BST1 influence RBC glycolysis;Modeled PFKP effects relate to preventing loss of the total AXP pool in stored RBCs;ATP and hypoxanthine are biomarkers of hemolysis in vitro and in vivo.

Genetically heterogeneous UM-HET3 mice born in 2020 were used to test possible lifespan effects of alpha-ketoglutarate (AKG), 2,4-dinitrophenol (DNP), hydralazine (HYD), nebivolol (NEBI), 16α-hydroxyestriol (OH_Est), and sodium thiosulfate (THIO), and to evaluate the effects of canagliflozin (Cana) when started at 16 months of age. OH_Est produced a 15% increase (p = 0.0001) in median lifespan in males but led to a significant (7%) decline in female lifespan. Cana, started at 16 months, also led to a significant increase (14%, p = 0.004) in males and a significant decline (6%, p = 0.03) in females. Cana given to mice at 6 months led, as in our previous study, to an increase in male lifespan without any change in female lifespan, suggesting that this agent may lead to female-specific late-life harm. We found that blood levels of Cana were approximately 20-fold higher in aged females than in young males, suggesting a possible mechanism for the sex-specific disparities in its effects. NEBI was also found to produce a female-specific decline (4%, p = 0.03) in lifespan. None of the other tested drugs provided a lifespan benefit in either sex. These data bring to 7 the list of ITP-tested drugs that induce at least a 10% lifespan increase in one or both sexes, add a fourth drug with demonstrated mid-life benefits on lifespan, and provide a testable hypothesis that might explain the sexual dimorphism in lifespan effects of the SGLT2 inhibitor Cana.

Systemic lupus erythematosus (SLE) is the most common form of lupus; however, the development of new drugs for SLE has been limited, with only two drugs approved in the past 60 years. Recently, the inhibition of glycolysis and disruption of N-linked glycosylation via 2- deoxyglucose (2DG), a competitive inhibitor of glucose-6-phosphate, has shown promising results in reducing autoimmune pathology. The precise mechanism of action, however, remains unclear. To elucidate how 2DG longitudinally modulates pathways in the context of SLE, BXSB.Yaa IL15-/-, CD8-/- mice were treated orally with 2DG for either 96-hours or 4-weeks.Subsequently, the transcriptional profiles of 9 tissues (heart, hippocampus, hypothalamus, kidney, liver, prefrontal cortex, skeletal muscle, small intestine, and spleen) were analyzed using unsupervised clustering and weighted gene correlation network analysis (WGCNA), ANOVA, and overrepresentation analysis. Principal component analysis showed distinct clustering based on tissue type. Through a three-step filtering procedure, we identified three response patterns to 2DG treatment: acute, long-term, or independent of dosage length. The three brain tissues had an increased immune response in both innate and adaptive immunity when treated with short-term 2DG. Five tissues, including kidney, responded to long-term 2DG treatment. Notably, 2DG predominantly affected pathways related to biosynthesis and metabolism of amino acids and genes related to glomerular nephropathy, suggesting its potential in preserving kidney function that is degraded in lupus pathogenesis. In addition, seven tissues (heart, hippocampus, hypothalamus, liver, prefrontal cortex, skeletal muscle, and spleen) were affected by 2DG regardless of dosage length. Although 2DG broadly effected pathways, each gene cluster within a tissue showed unique alterations based on time and/or treatment. Furthermore, each set of altered pathways could be classified as a downstream response to 2DG’s two major mechanisms of action: inhibition of glycolysis or disruption of N- linked glycosylation. These results show that 2DG has a systemic impact in altering immunity, amino acid metabolism, as well as other biological functions with tissue-specific variations. This study provides insights into the potential of 2DG as a therapeutic agent for addressing immune dysregulation, preserving kidney function, and treating lupus.

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic–epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic–epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

Inhibitory G alpha (GNAI or Gαi) proteins are critical for the polarized morphogenesis of sensory hair cells and for hearing. The extent and nature of their actual contributions remains unclear, however, as previous studies did not investigate all GNAI proteins and included non-physiological approaches. Pertussis toxin can downregulate functionally redundant GNAI1, GNAI2, GNAI3, and GNAO proteins, but may also induce unrelated defects. Here, we directly and systematically determine the role(s) of each individual GNAI protein in mouse auditory hair cells. GNAI2 and GNAI3 are similarly polarized at the hair cell apex with their binding partner G protein signaling modulator 2 (GPSM2), whereas GNAI1 and GNAO are not detected. In Gnai3 mutants, GNAI2 progressively fails to fully occupy the sub-cellular compartments where GNAI3 is missing. In contrast, GNAI3 can fully compensate for the loss of GNAI2 and is essential for hair bundle morphogenesis and auditory function. Simultaneous inactivation of Gnai2 and Gnai3 recapitulates for the first time two distinct types of defects only observed so far with pertussis toxin: (1) a delay or failure of the basal body to migrate off-center in prospective hair cells, and (2) a reversal in the orientation of some hair cell types. We conclude that GNAI proteins are critical for hair cells to break planar symmetry and to orient properly before GNAI2/3 regulate hair bundle morphogenesis with GPSM2.