Structural perturbation of chromatin domains with multiple developmental regulators can severely impact gene regulation and development.

Chromatin domain boundaries delimited by CTCF motifs can restrict the range of enhancer action. However, disruption of domain structure often results in mild gene dysregulation and thus predicting the impact of boundary rearrangements on animal development remains challenging. Here, we tested whether structural perturbation of a chromatin domain with multiple developmental regulators can result in more acute gene dysregulation and severe developmental phenotypes. We targeted clusters of CTCF motifs in a domain of the mouse genome containing three FGF ligand genes - Fgf3, Fgf4, and Fgf15 - that regulate several developmental processes. Deletion of the 23.9kb cluster that defines the centromeric boundary of this domain resulted in ectopic interactions of the FGF genes with enhancers located across the deleted boundary that are active in the developing brain. This caused strong induction of FGF expression and perinatal lethality with encephalocele and orofacial cleft phenotypes. Heterozygous boundary deletion was sufficient to cause these fully penetrant phenotypes, and strikingly, loss of a single CTCF motif within the cluster also recapitulated ectopic FGF expression and caused encephalocele. However, such phenotypic sensitivity to perturbation of domain structure did not extend to all CTCF clusters of this domain, nor to all developmental processes controlled by these three FGF genes - for example, the ability to undergo lineage specification in the blastocyst and pre-implantation development were not affected. By tracing the impact of different chromosomal rearrangements throughout mouse development, we start to uncover the determinants of phenotypic robustness and sensitivity to perturbation of chromatin boundaries. Our data show how small sequence variants at certain domain boundaries can have a surprisingly outsized effect and must be considered as potential sources of gene dysregulation during development and disease.

The role of LEDGF in transcription is exploited by HIV-1 to position integration.

HIV-1 integration occurs across actively transcribed genes due to the interaction of integrase with host chromatin factor LEDGF. Although LEDGF was originally isolated as a co-activator that stimulates promoter activity in purified systems, this role is inconsistent with LEDGF-mediated integration across gene bodies and with data indicating LEDGF is a histone chaperone that promotes transcriptional elongation. We found LEDGF is enriched in pronounced peaks that match the enrichments of H3K4me3 and RNA Pol II at transcription start sites (TSSs) of active promoters. Our genome-wide chromatin mapping revealed that MLL1 had a dominant role in recruiting LEDGF to promoters and the presence of LEDGF recruits RNA Pol II. Enrichment of LEDGF at TSSs correlates strongly with levels of integration across the transcribed sequences, indicating that LEDGF at TSSs contributed to integration across gene bodies. Although the N-terminal Pro-Trp-Trp-Pro (PWWP) domain of LEDGF interacts with nucleosomes containing H3K36me3, a modification thought to recruit LEDGF to chromatin, we found H3K36me3 does not contribute to gene specificity of integration. These data support a dual role model of LEDGF where it is tethered to promoters by MLL1 and recruits RNA Pol II. Subsequently, LEDGF travels across genes to effect HIV-1 integration. Our data also provides a mechanistic context for the contribution made by LEDGF to MLL1-based infant acute leukemia and acute myeloid leukemia in adults.

Messenger RNA transport on lysosomal vesicles maintains axonal mitochondrial homeostasis and prevents axonal degeneration.

In neurons, RNA granules are transported along the axon for local translation away from the soma. Recent studies indicate that some of this transport involves hitchhiking of RNA granules on lysosome-related vesicles. In the present study, we leveraged the ability to prevent transport of these vesicles into the axon by knockout of the lysosome-kinesin adaptor BLOC-one-related complex (BORC) to identify a subset of axonal mRNAs that depend on lysosome-related vesicles for transport. We found that BORC knockout causes depletion of a large group of axonal mRNAs mainly encoding ribosomal and mitochondrial/oxidative phosphorylation proteins. This depletion results in mitochondrial defects and eventually leads to axonal degeneration in human induced pluripotent stem cell (iPSC)-derived and mouse neurons. Pathway analyses of the depleted mRNAs revealed a mechanistic connection of BORC deficiency with common neurodegenerative disorders. These results demonstrate that mRNA transport on lysosome-related vesicles is critical for the maintenance of axonal homeostasis and that its failure causes axonal degeneration.

Loss of Ezh2 in the medial ganglionic eminence alters interneuron fate, cell morphology and gene expression profiles.

Introduction: Enhancer of zeste homolog 2 (Ezh2) is responsible for trimethylation of histone 3 at lysine 27 (H3K27me3), resulting in repression of gene expression. Here, we explore the role of Ezh2 in forebrain GABAergic interneuron development. Methods: We removed Ezh2 in the MGE by generating Nkx2-1Cre;Ezh2 conditional knockout mice. We then characterized changes in MGE-derived interneuron fate and electrophysiological properties in juvenile mice, as well as alterations in gene expression, chromatin accessibility and histone modifications in the MGE. Results: Loss of Ezh2 increases somatostatin-expressing (SST+) and decreases parvalbumin-expressing (PV+) interneurons in the forebrain. We observe fewer MGE-derived interneurons in the first postnatal week, indicating reduced interneuron production. Intrinsic electrophysiological properties in SST+ and PV+ interneurons are normal, but PV+ interneurons display increased axonal complexity in Ezh2 mutant mice. Single nuclei multiome analysis revealed differential gene expression patterns in the embryonic MGE that are predictive of these cell fate changes. Lastly, CUT&Tag analysis revealed that some genomic loci are particularly resistant or susceptible to shifts in H3K27me3 levels in the absence of Ezh2, indicating differential selectivity to epigenetic perturbation. Discussion: Thus, loss of Ezh2 in the MGE alters interneuron fate, morphology, and gene expression and regulation. These findings have important implications for both normal development and potentially in disease etiologies.

Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen Borrelia burgdorferi.

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. Here, we use several RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. We identify complex gene arrangements and operons, untranslated regions and small RNAs. We predict intrinsic terminators and experimentally test examples of Rho-dependent transcription termination. Remarkably, 63% of RNA 3' ends map upstream of or internal to open reading frames (ORFs), including genes involved in the unique infectious cycle of B. burgdorferi. We suggest these RNAs result from premature termination, processing and regulatory events such as cis-acting regulation. Furthermore, the polyamine spermidine globally influences the generation of truncated mRNAs. Collectively, our findings provide insights into transcription termination and uncover an abundance of potential RNA regulators in B. burgdorferi.

Germline loss-of-function PAM variants are enriched in subjects with pituitary hypersecretion.

Introduction: Pituitary adenomas (PAs) are common, usually benign tumors of the anterior pituitary gland which, for the most part, have no known genetic cause. PAs are associated with major clinical effects due to hormonal dysregulation and tumoral impingement on vital brain structures. PAM encodes a multifunctional protein responsible for the essential C-terminal amidation of secreted peptides. Methods: Following the identification of a loss-of-function variant (p.Arg703Gln) in the peptidylglycine a-amidating monooxygenase (PAM) gene in a family with pituitary gigantism, we investigated 299 individuals with sporadic PAs and 17 familial isolated PA kindreds for PAM variants. Genetic screening was performed by germline and tumor sequencing and germline copy number variation (CNV) analysis. Results: In germline DNA, we detected seven heterozygous, likely pathogenic missense, truncating, and regulatory SNVs. These SNVs were found in sporadic subjects with growth hormone excess (p.Gly552Arg and p.Phe759Ser), pediatric Cushing disease (c.-133T>C and p.His778fs), or different types of PAs (c.-361G>A, p.Ser539Trp, and p.Asp563Gly). The SNVs were functionally tested in vitro for protein expression and trafficking by Western blotting, splicing by minigene assays, and amidation activity in cell lysates and serum samples. These analyses confirmed a deleterious effect on protein expression and/or function. By interrogating 200,000 exomes from the UK Biobank, we confirmed a significant association of the PAM gene and rare PAM SNVs with diagnoses linked to pituitary gland hyperfunction. Conclusion: The identification of PAM as a candidate gene associated with pituitary hypersecretion opens the possibility of developing novel therapeutics based on altering PAM function.

Integrated spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The differentiation of osteoblasts and the subsequent formation of bone marks an important terminal phase in palate formation that leads to the separation of the oral and nasal cavities. While the developmental events that precede palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms that lead to the bony union of fusing palatal shelves. Herein, the timeline of osteogenic transcriptional programming is unveiled in the embryonic palate by way of integrated bulk, single-cell, and spatially resolved RNA-seq analyses. We define spatially restricted expression patterns of key marker genes, both regulatory and structural, that are differentially expressed during palatal fusion, including the identification of several novel genes ( Deup1, Dynlrb2, Lrrc23 ) spatially restricted in expression to the palate, providing a relevant framework for future studies that identify new candidate genes for cleft palate anomalies in humans as well as the timing of mammalian embryonic palatal osteogenesis.

Absence of TRIC-B from type XIV Osteogenesis Imperfecta osteoblasts alters cell adhesion and mitochondrial function - A multi-omics study.

Osteogenesis Imperfecta (OI) is a heritable collagen-related bone dysplasia characterized by bone fractures, growth deficiency and skeletal deformity. Type XIV OI is a recessive OI form caused by null mutations in TMEM38B, which encodes the ER membrane intracellular cation channel TRIC-B. Previously, we showed that absence of TMEM38B alters calcium flux in the ER of OI patient osteoblasts and fibroblasts, which further disrupts collagen synthesis and secretion. How the absence of TMEM38B affects osteoblast function is still poorly understood. Here we further investigated the role of TMEM38B in human osteoblast differentiation and mineralization. TMEM38B-null osteoblasts showed altered expression of osteoblast marker genes and decreased mineralization. RNA-Seq analysis revealed that cell-cell adhesion was one of the most downregulated pathways in TMEM38B-null osteoblasts, with further validation by real-time PCR and Western blot. Gap and tight junction proteins were also decreased by TRIC-B absence, both in patient osteoblasts and in calvarial osteoblasts of Tmem38b-null mice. Disrupted cell adhesion decreased mutant cell proliferation and cell cycle progression. An important novel finding was that TMEM38B-null osteoblasts had elongated mitochondria with altered fusion and fission markers, MFN2 and DRP1. In addition, TMEM38B-null osteoblasts exhibited a significant increase in superoxide production in mitochondria, further supporting mitochondrial dysfunction. Together these results emphasize the novel role of TMEM38B/TRIC-B in osteoblast differentiation, affecting cell-cell adhesion processes, gap and tight junction, proliferation, cell cycle, and mitochondrial function.

Cerebrospinal Fluid Protein Biomarker Discovery in CLN3.

Syndromic CLN3-Batten is a fatal, pediatric, neurodegenerative disease caused by variants in CLN3, which encodes the endolysosomal transmembrane CLN3 protein. No approved treatment for CLN3 is currently available. The protracted and asynchronous disease presentation complicates the evaluation of potential therapies using clinical disease progression parameters. Biomarkers as surrogates to measure the progression and effect of potential therapeutics are needed. We performed proteomic discovery studies using cerebrospinal fluid (CSF) samples from 28 CLN3-affected and 32 age-similar non-CLN3 individuals. Proximal extension assay (PEA) of 1467 proteins and untargeted data-dependent mass spectrometry [MS; MassIVE FTP server (ftp://MSV000090147@massive.ucsd.edu)] were used to generate orthogonal lists of protein marker candidates. At an adjusted p-value of <0.1 and threshold CLN3/non-CLN3 fold-change ratio of 1.5, PEA identified 54 and MS identified 233 candidate biomarkers. Some of these (NEFL, CHIT1) have been previously linked with other neurologic conditions. Others (CLPS, FAM217B, QRICH2, KRT16, ZNF333) appear to be novel. Both methods identified 25 candidate biomarkers, including CHIT1, NELL1, and ISLR2 which had absolute fold-change ratios >2. NELL1 and ISLR2 regulate axonal development in neurons and are intriguing new candidates for further investigation in CLN3. In addition to identifying candidate proteins for CLN3 research, this study provides a comparison of two large-scale proteomic discovery methods in CSF.

Characterization of serotonin-5-HTR1E signaling pathways and its role in cell survival.

5-Hydroxytryptamine receptor 1E (5-HTR1E) is reported to activate cyclic AMP (cAMP) and extracellular-signal related kinases (ERK) pathways via its ligands and binding partners, but the detailed mechanism underlying the serotonin-induced 5-HTR1E signaling is still not known. In the present study, we determined the cellular regulators of ERK and cAMP signaling pathways in response to serotonin-induced 5-HTR1E activation in 5-HTR1E overexpressing HEK293 cells. We found that Pertussis Toxin (PTX) treatment completely reversed the effect of serotonin-5-HTR1E mediated signaling on cAMP and ERK pathways, confirming the involvement of a Gαi-linked cascade. We also observed that Gßγ and Gq were not associated with 5-HTR1E activation, while blocking protein kinase A (PKA) inhibited ERK signaling only, and had no effect on cAMP. Additionally, serotonin-stimulated ERK1/2 phosphorylation was similar in 5-HTR1E overexpressing, ß-arrestin-deficient HEK293 cells and is solely dependent on G protein signaling. siRNA mediated gene knockdown studies in SH-SY5Y cells revealed that the inhibition of 5-HTR1E reduced the expression of cMyc, Cyclin D1, Cyclin E and BCL2 genes which are related to cell cycle regulation and survival. MTT assays showed that 5-HTR1E knockdown in SHSY-5Y and U118 cells inhibited cell survival significantly. In addition to the signaling mechanism, we also performed RNA-seq analysis in 5-HTR1E overexpressing HEK293 cells and found that 5-HTR1E can regulate the expression of Receptor activity modifying protein 1 (RAMP1), Nuclear receptor 1 (NR4A1) and other Cyclin genes. These findings indicate that serotonin interaction with 5-HTR1E receptor simultaneously activates cAMP and ERK pathway in HEK293 cells and its expression is important for cell survival.

miRNA Expression Profiles of Mouse Round Spermatids in GRTH/DDX25-Mediated Spermiogenesis: mRNA-miRNA Network Analysis.

GRTH/DDX25 is a testis-specific DEAD-box family of RNA helicase, which plays an essential role in spermatogenesis and male fertility. There are two forms of GRTH, a 56 kDa non-phosphorylated form and a 61 kDa phosphorylated form (pGRTH). GRTH-KO and GRTH Knock-In (KI) mice with R242H mutation (lack pGRTH) are sterile with a spermatogenic arrest at step 8 of spermiogenesis due to failure of round spermatids (RS) to elongate. We performed mRNA-seq and miRNA-seq analysis on RS of WT, KI, and KO to identify crucial microRNAs (miRNAs) and mRNAs during RS development by establishing a miRNA-mRNA network. We identified increased levels of miRNAs such as miR146, miR122a, miR26a, miR27a, miR150, miR196a, and miR328 that are relevant to spermatogenesis. mRNA-miRNA target analysis on these DE-miRNAs and DE-mRNAs revealed miRNA target genes involved in ubiquitination process (Ube2k, Rnf138, Spata3), RS differentiation, and chromatin remodeling/compaction (Tnp1/2, Prm1/2/3, Tssk3/6), reversible protein phosphorylation (Pim1, Hipk1, Csnk1g2, Prkcq, Ppp2r5a), and acrosome stability (Pdzd8). Post-transcriptional and translational regulation of some of these germ-cell-specific mRNAs by miRNA-regulated translation arrest and/or decay may lead to spermatogenic arrest in KO and KI mice. Our studies demonstrate the importance of pGRTH in the chromatin compaction and remodeling process, which mediates the differentiation of RS into elongated spermatids through miRNA-mRNA interactions.

Identification of cerebral spinal fluid protein biomarkers in Niemann-Pick disease, type C1.

BACKGROUND: Niemann-Pick disease, type C1 (NPC1) is an ultrarare, recessive, lethal, lysosomal disease characterized by progressive cerebellar ataxia and cognitive impairment. Although the NPC1 phenotype is heterogeneous with variable age of onset, classical NPC1 is a pediatric disorder. Currently there are no therapies approved by the FDA and therapeutics trials for NPC1 are complicated by disease rarity, heterogeneity, and the relatively slow rate of neurological decline. Thus, identification of disease relevant biomarkers is necessary to provide tools that can support drug development efforts for this devastating neurological disease. METHODS: Proximal extension assays (O-link® Explore 1536) were used to compare cerebrospinal fluid (CSF) samples from individuals with NPC1 enrolled in a natural history study and non-NPC1 comparison samples. Relative expression levels of 1467 proteins were determined, and candidate protein biomarkers were identified by evaluating fold-change and adjusted Kruskal-Wallis test p-values. Selected proteins were orthogonally confirmed using ELISA. To gain insight into disease progression and severity we evaluated the altered protein expression with respect to clinically relevant phenotypic aspects: NPC Neurological Severity Score (NPC1 NSS), Annual Severity Increment Score (ASIS) and age of neurological onset. RESULTS: This study identified multiple proteins with altered levels in CSF from individuals with NPC1 compared to non-NPC1 samples. These included proteins previously shown to be elevated in NPC1 (NEFL, MAPT, CHIT1, CALB1) and additional proteins confirmed by orthogonal assays (PARK7, CALB2/calretinin, CHI3L1/YKL-40, MIF, CCL18 and ENO2). Correlations with clinically relevant phenotypic parameters demonstrated moderate negative (p = 0.0210, r = -0.41) and possible moderate positive (p = 0.0631, r = 0.33) correlation of CSF CALB2 levels with age of neurological onset and ASIS, respectively. CSF CHI3L1 levels showed a moderate positive (p = 0.0183, r = 0.40) correlation with the concurrent NPC1 NSS. A strong negative correlation (p = 0.0016, r = -0.648) was observed between CSF CCL18 and age of neurological onset for childhood/adolescent cases. CSF CCL18 levels also showed a strong positive correlation (p = 0.0017, r = 0.61) with ASIS. CONCLUSION: Our study identified and validated multiple proteins in CSF from individuals with NPC1 that are candidates for further investigation in a larger cohort. These analytes may prove to be useful as supportive data in therapeutic trials. TRIAL REGISTRATIONS: NCT00344331, NCT00001721, NCT02931682.

Characterization of serotonin-5-HTR1E signaling pathways and its role in cell survival.

5-Hydroxy tryptamine receptor 1E (5-HTR1E) is reported to activate cAMP and ERK pathways via its ligands and binding partners, but the detailed mechanism underlying the serotonin induced 5-HTR1E signaling is still not known. In the present study, we determined the cellular regulators of ERK and cAMP signaling pathways in response to serotonin induced 5-HTR1E activation in 5-HTR1E overexpressing HEK293 cells. We found that Pertussis Toxin (PTX) treatment completely reversed the effect of serotonin-5-HTR1E mediated signaling on cAMP and ERK pathways, confirming the involvement of a Gαi-linked cascade. We also observed that Gßγ and Gq were not associated with 5-HTR1E activation, while blocking PKA inhibited ERK signaling only, and had no effect on cAMP. Additionally, serotonin-stimulated ERK1/2 phosphorylation was similar in 5-HTR1E overexpressing, ß-arrestin-deficient HEK293 cells and is solely dependent on G protein signaling. siRNA mediated gene knockout studies in SH-SY5Y cells revealed that the inhibition of 5-HTR1E reduced the expression of cMyc, Cyclin D1, Cyclin E and BCL2 genes which are related to cell cycle regulation and survival. MTT assays showed that 5-HTR1E knockdown in SHSY-5Y and U118 cells inhibited cell survival significantly. In addition to the signaling mechanism, we also performed RNA-seq analysis in 5-HTR1E overexpressing HEK293 cells and found that 5-HTR1E can regulate the expression of Receptor activity modifying protein 1 ( RAMP1 ), Nuclear receptor 1 ( NR4A1 ) and other Cyclin genes. These findings indicate that serotonin interaction with 5-HTR1E receptor simultaneously activates cAMP and ERK pathway in HEK293 cells and its expression is important for cell survival.

Enhancer-promoter interactions can bypass CTCF-mediated boundaries and contribute to phenotypic robustness.

How enhancers activate their distal target promoters remains incompletely understood. Here we dissect how CTCF-mediated loops facilitate and restrict such regulatory interactions. Using an allelic series of mouse mutants, we show that CTCF is neither required for the interaction of the Sox2 gene with distal enhancers, nor for its expression. Insertion of various combinations of CTCF motifs, between Sox2 and its distal enhancers, generated boundaries with varying degrees of insulation that directly correlated with reduced transcriptional output. However, in both epiblast and neural tissues, enhancer contacts and transcriptional induction could not be fully abolished, and insertions failed to disrupt implantation and neurogenesis. In contrast, Sox2 expression was undetectable in the anterior foregut of mutants carrying the strongest boundaries, and these animals fully phenocopied loss of SOX2 in this tissue. We propose that enhancer clusters with a high density of regulatory activity can better overcome physical barriers to maintain faithful gene expression and phenotypic robustness.

Germline loss-of-function PAM variants are enriched in subjects with pituitary hypersecretion.

Pituitary adenomas (PAs) are common, usually benign tumors of the anterior pituitary gland which, for the most part, have no known genetic cause. PAs are associated with major clinical effects due to hormonal dysregulation and tumoral impingement on vital brain structures. Following the identification of a loss-of-function variant (p.Arg703Gln) in the PAM gene in a family with pituitary gigantism, we investigated 299 individuals with sporadic PAs and 17 familial isolated pituitary adenomas kindreds for PAM variants. PAM encodes a multifunctional protein responsible for the essential C-terminal amidation of secreted peptides. Genetic screening was performed by germline and tumor sequencing and germline copy number variation (CNV) analysis. No germline CNVs or somatic single nucleotide variants (SNVs) were identified. We detected seven likely pathogenic heterozygous missense, truncating, and regulatory SNVs. These SNVs were found in sporadic subjects with GH excess (p.Gly552Arg and p.Phe759Ser), pediatric Cushing disease (c.-133T>C and p.His778fs), or with different types of PAs (c.-361G>A, p.Ser539Trp, and p.Asp563Gly). The SNVs were functionally tested in vitro for protein expression and trafficking by Western blotting, for splicing by minigene assays, and for amidation activity in cell lysates and serum samples. These analyses confirmed a deleterious effect on protein expression and/or function. By interrogating 200,000 exomes from the UK Biobank, we confirmed a significant association of the PAM gene and rare PAM SNVs to diagnoses linked to pituitary gland hyperfunction. Identification of PAM as a candidate gene associated with pituitary hypersecretion opens the possibility of developing novel therapeutics based on altering PAM function.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen B. burgdorferi.

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. We employed complementary RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. By systematically mapping B. burgdorferi RNA ends at single nucleotide resolution, we delineated complex gene arrangements and operons and mapped untranslated regions (UTRs) and small RNAs (sRNAs). We experimentally tested modes of B. burgdorferi transcription termination and compared our findings to observations in E. coli , P. aeruginosa , and B. subtilis . We discovered 63% of B. burgdorferi RNA 3' ends map upstream or internal to open reading frames (ORFs), suggesting novel mechanisms of regulation. Northern analysis confirmed the presence of stable 5' derived RNAs from mRNAs encoding gene products involved in the unique infectious cycle of B. burgdorferi . We suggest these RNAs resulted from premature termination and regulatory events, including forms of cis- acting regulation. For example, we documented that the polyamine spermidine globally influences the generation of truncated mRNAs. In one case, we showed that high spermidine concentrations increased levels of RNA fragments derived from an mRNA encoding a spermidine import system, with a concomitant decrease in levels of the full- length mRNA. Collectively, our findings revealed new insight into transcription termination and uncovered an abundance of potential RNA regulators.

Profiles of Wnt pathway gene expression during tooth morphogenesis.

Mouse and human genetic studies indicate key roles of the Wnt10a ligand in odontogenesis. Previous studies have identified effectors and regulators of the Wnt signaling pathway actively expressed during key stages of tooth morphogenesis. However, limitations in multiplexing and spatial resolution hindered a more comprehensive analysis of these signaling molecules. Here, profiling of transcriptomes using fluorescent multiplex in situ hybridization and single-cell RNA-sequencing (scRNA-seq) provide robust insight into the synchronized expression patterns of Wnt10a, Dkk1, and Sost simultaneously during tooth development. First, we identified Wnt10a transcripts restricted to the epithelium at the stage of tooth bud morphogenesis, contrasting that of Sost and Dkk1 localization to the dental mesenchyme. By embryonic day 15.5 (E15.5), a marked shift of Wnt10a expression from dental epithelium to mesenchyme was noted, while Sost and Dkk1 expression remained enriched in the mesenchyme. By postnatal day 0 (P0), co-localization patterns of Wnt10a, Dkk1, and Sost were observed in both terminally differentiating and secreting odontoblasts of molars and incisors. Interestingly, Wnt10a exhibited robust expression in fully differentiated ameloblasts at the developing cusp tip of both molars and incisors, an observation not previously noted in prior studies. At P7 and 14, after the mineralization of dentin and enamel, Wnt10a expression was limited to odontoblasts. Meanwhile, Wnt modulators showed reduced or absent signals in molars. In contrast, strong signals persisted in ameloblasts (for Wnt10a) and odontoblasts (for Wnt10a, Sost, and Dkk1) towards the proximal end of incisors, near the cervical loop. Our scRNA-seq analysis used CellChat to further contextualize Wnt pathway-mediated communication between cells by examining ligand-receptor interactions among different clusters. The co-localization pattern of Wnt10a, Dkk1, and Sost in both terminally differentiating and secreting odontoblasts of molars and incisors potentially signifies the crucial ligand-modulator interaction along the gradient of cytodifferentiation starting from each cusp tip towards the apical region. These data provide cell type-specific insight into the role of Wnt ligands and mediators during epithelial-mesenchymal interactions in odontogenesis.

SMAD3 mutation in LDS3 causes bone fragility by impairing the TGF-ß pathway and enhancing osteoclastogenesis.

Loss-of-function mutations in SMAD3 cause Loeys-Dietz syndrome type 3 (LDS3), a rare autosomal-dominant connective tissue disorder characterized by vascular pathology and skeletal abnormalities. Dysregulation of TGF-ß/SMAD signaling is associated with abnormal skeletal features and bone fragility. To date, histomorphometric and ultrastructural characteristics of bone with SMAD3 mutations have not been reported in humans and the exact mechanism by which SMAD3 mutations cause the LDS3 phenotype is poorly understood. Here, we investigated bone histomorphometry and matrix mineralization in human bone with a SMAD3 mutation and explored the associated cellular defect in the TGF-ß/SMAD pathway in vitro. The index patient had recurrent fractures, mild facial dysmorphism, arachnodactyly, pectus excavatum, chest asymmetry and kyphoscoliosis. Bone histomorphometry revealed markedly reduced cortical thickness (-68 %), trabecular thickness (-32 %), bone formation rate (-50 %) and delayed mineralization. Quantitative backscattered electron imaging demonstrated undermineralized bone matrix with increased heterogeneity in mineralization. The patient's SMAD3 mutation (c.200 T > G; p.I67S), when expressed from plasmid vectors in HEK293 cells, showed reduced phosphorylation and transcription factor activity compared to normal control and SMAD3 (p.S264Y), a gain-of-function mutation, somatic mosaicism of which causes melorheostosis. Transfection study of the patients' SMAD3 (p.I67S) mutation displayed lower luciferase reporter activity than normal SMAD3 and reduced expression of TGF-ß signaling target genes. Patient fibroblasts also demonstrated impaired SMAD3 protein stability. Osteoclastogenic differentiation significantly increased and osteoclast-associated genes, including ACP5 (encoding TRAP), ATP6V0D2, and DCSTAMP, were up-regulated in CD14 (+) peripheral blood mononuclear cells (PBMCs) with the SMAD3 (p.I67S) mutation. Upregulation of osteoclastogenic genes was associated with decreased expression of TGF-ß signaling target genes. We conclude that bone with the SMAD3 (p.I67S) mutation features reduced bone formation, and our functional studies revealed decreased SMAD3 activation and protein stability as well as increased osteoclastogenesis. These findings enhance our understanding of the pathophysiology of LDS3 caused by SMAD3 mutations. Emerging therapies targeting in the TGF-ß/SMAD pathway also raise hope for treatment of LDS3.

An epigenome atlas of neural progenitors within the embryonic mouse forebrain.

A comprehensive characterization of epigenomic organization in the embryonic mouse forebrain will enhance our understanding of neurodevelopment and provide insight into mechanisms of neurological disease. Here we collected single-cell chromatin accessibility profiles from four distinct neurogenic regions of the embryonic mouse forebrain using single nuclei ATAC-Seq (snATAC-Seq). We identified thousands of differentially accessible peaks, many restricted to distinct progenitor cell types or brain regions. We integrated snATAC-Seq and single cell transcriptome data to characterize changes of chromatin accessibility at enhancers and promoters with associated transcript abundance. Multi-modal integration of histone modifications (CUT&Tag and CUT&RUN), promoter-enhancer interactions (Capture-C) and high-order chromatin structure (Hi-C) extended these initial observations. This dataset reveals a diverse chromatin landscape with region-specific regulatory mechanisms and genomic interactions in distinct neurogenic regions of the embryonic mouse brain and represents an extensive public resource of a 'ground truth' epigenomic landscape at this critical stage of neurogenesis.