Txnip Enhances Fitness of Dnmt3a-Mutant Hematopoietic Stem Cells via p21.

Clonal hematopoiesis (CH) refers to the age-related expansion of specific clones in the blood system, and manifests from somatic mutations acquired in hematopoietic stem cells (HSCs). Most CH variants occur in the gene DNMT3A, but while DNMT3A-mutant CH becomes almost ubiquitous in aging humans, a unifying molecular mechanism to illuminate how DNMT3A-mutant HSCs outcompete their counterparts is lacking. Here, we used interferon gamma (IFNγ) as a model to study the mechanisms by which Dnmt3a mutations increase HSC fitness under hematopoietic stress. We found Dnmt3a-mutant HSCs resist IFNγ-mediated depletion, and IFNγ-signaling is required for clonal expansion of Dnmt3a-mutant HSCs in vivo. Mechanistically, DNA hypomethylation-associated overexpression of Txnip in Dnmt3a-mutant HSCs leads to p53 stabilization and upregulation of p21. This preserves the functional potential of Dnmt3a-mutant HSCs through increased quiescence and resistance to IFNγ-induced apoptosis. These data identify a previously undescribed mechanism to explain increased fitness of DNMT3A-mutant clones under hematopoietic stress. SIGNIFICANCE: DNMT3A mutations are common variants in clonal hematopoiesis, and recurrent events in blood cancers. Yet the mechanisms by which these mutations provide hematopoietic stem cells a competitive advantage as a precursor to malignant transformation remain unclear. Here, we use inflammatory stress to uncover molecular mechanisms leading to this fitness advantage.See related commentary by De Dominici and DeGregori, p. 178. This article is highlighted in the In This Issue feature, p. 171.

OCT4 cooperates with distinct ATP-dependent chromatin remodelers in naïve and primed pluripotent states in human.

Understanding the molecular underpinnings of pluripotency is a prerequisite for optimal maintenance and application of embryonic stem cells (ESCs). While the protein-protein interactions of core pluripotency factors have been identified in mouse ESCs, their interactome in human ESCs (hESCs) has not to date been explored. Here we mapped the OCT4 interactomes in naïve and primed hESCs, revealing extensive connections to mammalian ATP-dependent nucleosome remodeling complexes. In naïve hESCs, OCT4 is associated with both BRG1 and BRM, the two paralog ATPases of the BAF complex. Genome-wide location analyses and genetic studies reveal that these two enzymes cooperate in a functionally redundant manner in the transcriptional regulation of blastocyst-specific genes. In contrast, in primed hESCs, OCT4 cooperates with BRG1 and SOX2 to promote chromatin accessibility at ectodermal genes. This work reveals how a common transcription factor utilizes differential BAF complexes to control distinct transcriptional programs in naïve and primed hESCs.

AIAP: A Quality Control and Integrative Analysis Package to Improve ATAC-seq Data Analysis.

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a technique widely used to investigate genome-wide chromatin accessibility. The recently published Omni-ATAC-seq protocol substantially improves the signal/noise ratio and reduces the input cell number. High-quality data are critical to ensure accurate analysis. Several tools have been developed for assessing sequencing quality and insertion size distribution for ATAC-seq data; however, key quality control (QC) metrics have not yet been established to accurately determine the quality of ATAC-seq data. Here, we optimized the analysis strategy for ATAC-seq and defined a series of QC metrics for ATAC-seq data, including reads under peak ratio (RUPr), background (BG), promoter enrichment (ProEn), subsampling enrichment (SubEn), and other measurements. We incorporated these QC tests into our recently developed ATAC-seq Integrative Analysis Package (AIAP) to provide a complete ATAC-seq analysis system, including quality assurance, improved peak calling, and downstream differential analysis. We demonstrated a significant improvement of sensitivity (20%-60%) in both peak calling and differential analysis by processing paired-end ATAC-seq datasets using AIAP. AIAP is compiled into Docker/Singularity, and it can be executed by one command line to generate a comprehensive QC report. We used ENCODE ATAC-seq data to benchmark and generate QC recommendations, and developed qATACViewer for the user-friendly interaction with the QC report. The software, source code, and documentation of AIAP are freely available at https://github.com/Zhang-lab/ATAC-seq_QC_analysis.

Altered neuronal physiology, development, and function associated with a common chromosome 15 duplication involving CHRNA7.

BACKGROUND: Copy number variants (CNVs) linked to genes involved in nervous system development or function are often associated with neuropsychiatric disease. While CNVs involving deletions generally cause severe and highly penetrant patient phenotypes, CNVs leading to duplications tend instead to exhibit widely variable and less penetrant phenotypic expressivity among affected individuals. CNVs located on chromosome 15q13.3 affecting the alpha-7 nicotinic acetylcholine receptor subunit (CHRNA7) gene contribute to multiple neuropsychiatric disorders with highly variable penetrance. However, the basis of such differential penetrance remains uncharacterized. Here, we generated induced pluripotent stem cell (iPSC) models from first-degree relatives with a 15q13.3 duplication and analyzed their cellular phenotypes to uncover a basis for the dissimilar phenotypic expressivity. RESULTS: The first-degree relatives studied included a boy with autism and emotional dysregulation (the affected proband-AP) and his clinically unaffected mother (UM), with comparison to unrelated control models lacking this duplication. Potential contributors to neuropsychiatric impairment were modeled in iPSC-derived cortical excitatory and inhibitory neurons. The AP-derived model uniquely exhibited disruptions of cellular physiology and neurodevelopment not observed in either the UM or unrelated controls. These included enhanced neural progenitor proliferation but impaired neuronal differentiation, maturation, and migration, and increased endoplasmic reticulum (ER) stress. Both the neuronal migration deficit and elevated ER stress could be selectively rescued by different pharmacologic agents. Neuronal gene expression was also dysregulated in the AP, including reduced expression of genes related to behavior, psychological disorders, neuritogenesis, neuronal migration, and Wnt, axonal guidance, and GABA receptor signaling. The UM model instead exhibited upregulated expression of genes in many of these same pathways, suggesting that molecular compensation could have contributed to the lack of neurodevelopmental phenotypes in this model. However, both AP- and UM-derived neurons exhibited shared alterations of neuronal function, including increased action potential firing and elevated cholinergic activity, consistent with increased homomeric CHRNA7 channel activity. CONCLUSIONS: These data define both diagnosis-associated cellular phenotypes and shared functional anomalies related to CHRNA7 duplication that may contribute to variable phenotypic penetrance in individuals with 15q13.3 duplication. The capacity for pharmacological agents to rescue some neurodevelopmental anomalies associated with diagnosis suggests avenues for intervention for carriers of this duplication and other CNVs that cause related disorders.

A genetic screen for regulators of muscle morphogenesis in Drosophila.

The mechanisms that determine the final topology of skeletal muscles remain largely unknown. We have been developing Drosophila body wall musculature as a model to identify and characterize the pathways that control muscle size, shape, and orientation during embryogenesis. Our working model argues muscle morphogenesis is regulated by (1) extracellular guidance cues that direct muscle cells toward muscle attachment sites, and (2) contact-dependent interactions between muscles and tendon cells. While we have identified several pathways that regulate muscle morphogenesis, our understanding is far from complete. Here, we report the results of a recent EMS-based forward genetic screen that identified a myriad of loci not previously associated with muscle morphogenesis. We recovered new alleles of known muscle morphogenesis genes, including back seat driver, kon-tiki, thisbe, and tumbleweed, arguing our screen had the depth and precision to uncover myogenic genes. We also identified new alleles of spalt-major, barren, and patched that presumably disrupt independent muscle morphogenesis pathways. Equally as important, our screen shows that at least 11 morphogenetic loci remain to be mapped and characterized. Our screen has developed exciting new tools to study muscle morphogenesis, which may provide future insights into the mechanisms that regulate skeletal muscle topology.

Deconstructing Stepwise Fate Conversion of Human Fibroblasts to Neurons by MicroRNAs.

Cell-fate conversion generally requires reprogramming effectors to both introduce fate programs of the target cell type and erase the identity of starting cell population. Here, we reveal insights into the activity of microRNAs miR-9/9∗ and miR-124 (miR-9/9∗-124) as reprogramming agents that orchestrate direct conversion of human fibroblasts into motor neurons by first eradicating fibroblast identity and promoting uniform transition to a neuronal state in sequence. We identify KLF-family transcription factors as direct target genes for miR-9/9∗-124 and show their repression is critical for erasing fibroblast fate. Subsequent gain of neuronal identity requires upregulation of a small nuclear RNA, RN7SK, which induces accessibilities of chromatin regions and neuronal gene activation to push cells to a neuronal state. Our study defines deterministic components in the microRNA-mediated reprogramming cascade.

Tissue-specific usage of transposable element-derived promoters in mouse development.

BACKGROUND: Transposable elements (TEs) are a significant component of eukaryotic genomes and play essential roles in genome evolution. Mounting evidence indicates that TEs are highly transcribed in early embryo development and contribute to distinct biological functions and tissue morphology. RESULTS: We examine the epigenetic dynamics of mouse TEs during the development of five tissues: intestine, liver, lung, stomach, and kidney. We found that TEs are associated with over 20% of open chromatin regions during development. Close to half of these accessible TEs are only activated in a single tissue and a specific developmental stage. Most accessible TEs are rodent-specific. Across these five tissues, 453 accessible TEs are found to create the transcription start sites of downstream genes in mouse, including 117 protein-coding genes and 144 lincRNA genes, 93.7% of which are mouse-specific. Species-specific TE-derived transcription start sites are found to drive the expression of tissue-specific genes and change their tissue-specific expression patterns during evolution. CONCLUSION: Our results suggest that TE insertions increase the regulatory potential of the genome, and some TEs have been domesticated to become a crucial component of gene and regulate tissue-specific expression during mouse tissue development.

Comparison of differential accessibility analysis strategies for ATAC-seq data.

ATAC-seq is widely used to measure chromatin accessibility and identify open chromatin regions (OCRs). OCRs usually indicate active regulatory elements in the genome and are directly associated with the gene regulatory network. The identification of differential accessibility regions (DARs) between different biological conditions is critical in determining the differential activity of regulatory elements. Differential analysis of ATAC-seq shares many similarities with differential expression analysis of RNA-seq data. However, the distribution of ATAC-seq signal intensity is different from that of RNA-seq data, and higher sensitivity is required for DARs identification. Many different tools can be used to perform differential analysis of ATAC-seq data, but a comprehensive comparison and benchmarking of these methods is still lacking. Here, we used simulated datasets to systematically measure the sensitivity and specificity of six different methods. We further discussed the statistical and signal density cut-offs in the differential analysis of ATAC-seq by applying them to real data. Batch effects are very common in high-throughput sequencing experiments. We illustrated that batch-effect correction can dramatically improve sensitivity in the differential analysis of ATAC-seq data. Finally, we developed a user-friendly package, BeCorrect, to perform batch effect correction and visualization of corrected ATAC-seq signals in a genome browser.

Transcriptomic analysis of bone and fibrous tissue morphogenesis during digit tip regeneration in the adult mouse.

Humans have limited regenerative potential of musculoskeletal tissues following limb or digit loss. The murine digit has been used to study mammalian regeneration, where stem/progenitor cells (the "blastema") completely regenerate the digit tip after distal, but not proximal, amputation. However, the molecular mechanisms responsible for this response remain to be determined. Here, we evaluated the spatiotemporal formation of bone and fibrous tissues after level-dependent amputation of the murine terminal phalanx and quantified the transcriptome of the repair tissue. Distal (regenerative) and proximal (non-regenerative) amputations showed significant differences in temporal gene expression and tissue regrowth over time. Genes that direct skeletal system development and limb morphogenesis are transiently upregulated during blastema formation and differentiation, including distal Hox genes. Overall, our results suggest that digit tip regeneration is controlled by a gene regulatory network that recapitulates aspects of limb development, and that failure to activate this developmental program results in fibrotic wound healing.

The Reissner Fiber Is Highly Dynamic In Vivo and Controls Morphogenesis of the Spine.

Cerebrospinal fluid (CSF) physiology is important for the development and homeostasis of the central nervous system, and its disruption has been linked to scoliosis in zebrafish [1, 2]. Suspended in the CSF is an extracellular structure called the Reissner fiber, which extends from the brain through the central canal of the spinal cord. Zebrafish scospondin-null mutants are unable to assemble a Reissner fiber and fail to form a straight body axis during embryonic development [3]. Here, we describe hypomorphic missense mutations of scospondin, which allow Reissner fiber assembly and extension of a straight axis. However, during larval development, these mutants display progressive Reissner fiber disassembly, which is concomitant with the emergence of axial curvatures and scoliosis in adult animals. Using a scospondin-GFP knockin zebrafish line, we demonstrate several dynamic properties of the Reissner fiber in vivo, including embryonic fiber assembly, the continuous rostral to caudal movement of the fiber within the brain and central canal, and subcommissural organ (SCO)-spondin-GFP protein secretion from the floor plate to merge with the fiber. Finally, we show that disassembly of the Reissner fiber is also associated with the progression of axial curvatures in distinct scoliosis mutant zebrafish models. Together, these data demonstrate a critical role for the Reissner fiber for the maintenance of a straight body axis and spine morphogenesis in adult zebrafish. Our study establishes a framework for future investigations to address the cellular effectors responsible for Reissner-fiber-dependent regulation of axial morphology. VIDEO ABSTRACT.

Geminin is required for Hox gene regulation to pattern the developing limb.

Development of the complex structure of the vertebrate limb requires carefully orchestrated interactions between multiple regulatory pathways and proteins. Among these, precise regulation of 5' Hox transcription factor expression is essential for proper limb bud patterning and elaboration of distinct limb skeletal elements. Here, we identified Geminin (Gmnn) as a novel regulator of this process. A conditional model of Gmnn deficiency resulted in loss or severe reduction of forelimb skeletal elements, while both the forelimb autopod and hindlimb were unaffected. 5' Hox gene expression expanded into more proximal and anterior regions of the embryonic forelimb buds in this Gmnn-deficient model. A second conditional model of Gmnn deficiency instead caused a similar but less severe reduction of hindlimb skeletal elements and hindlimb polydactyly, while not affecting the forelimb. An ectopic posterior SHH signaling center was evident in the anterior hindlimb bud of Gmnn-deficient embryos in this model. This center ectopically expressed Hoxd13, the HOXD13 target Shh, and the SHH target Ptch1, while these mutant hindlimb buds also had reduced levels of the cleaved, repressor form of GLI3, a SHH pathway antagonist. Together, this work delineates a new role for Gmnn in modulating Hox expression to pattern the vertebrate limb.

Derivation of trophoblast stem cells from naïve human pluripotent stem cells.

Naïve human pluripotent stem cells (hPSCs) provide a unique experimental platform of cell fate decisions during pre-implantation development, but their lineage potential remains incompletely characterized. As naïve hPSCs share transcriptional and epigenomic signatures with trophoblast cells, it has been proposed that the naïve state may have enhanced predisposition for differentiation along this extraembryonic lineage. Here we examined the trophoblast potential of isogenic naïve and primed hPSCs. We found that naïve hPSCs can directly give rise to human trophoblast stem cells (hTSCs) and undergo further differentiation into both extravillous and syncytiotrophoblast. In contrast, primed hPSCs do not support hTSC derivation, but give rise to non-self-renewing cytotrophoblasts in response to BMP4. Global transcriptome and chromatin accessibility analyses indicate that hTSCs derived from naïve hPSCs are similar to blastocyst-derived hTSCs and acquire features of post-implantation trophectoderm. The derivation of hTSCs from naïve hPSCs will enable elucidation of early mechanisms that govern normal human trophoblast development and associated pathologies.

The placenta is one of the most important human organs, but it is perhaps the least understood. The first decision the earliest human cells have to make, shortly after the egg is fertilized by a sperm, is whether to become part of the embryo or part of the placenta. This choice happens before a pregnancy even implants into the uterus. The cells that commit to becoming the embryo transform into 'naïve pluripotent' cells, capable of becoming any cell in the body. Those that commit to becoming the placenta transform into 'trophectoderm' cells, capable of becoming the two types of cell in the placenta. Placental cells either invade into the uterus to anchor the placenta or produce hormones to support the pregnancy. Once a pregnancy implants into the uterus, the naïve pluripotent cells in the embryo become 'primed'. This prevents them from becoming cells of the placenta, and it poses a problem for placental research. In 2018, scientists in Japan reported conditions for growing trophectoderm cells in the laboratory, where they are known as "trophoblast stem cells". These cells were capable of transforming into specialized placental cells, but needed first to be isolated from the human embryo or placenta itself. Dong et al. now show how to reprogram other pluripotent cells grown in the laboratory to produce trophoblast stem cells. The first step was to reset primed pluripotent cells to put them back into a naïve state. Then, Dong et al. exposed the cells to the same concoction of nutrients and chemicals used in the 2018 study. This fluid triggered a transformation in the naïve pluripotent cells; they started to look like trophoblast stem cells, and they switched on genes normally active in trophectoderm cells. To test whether these cells had the same properties as trophoblast stem cells, Dong et al. gave them chemical signals to see if they could mature into placental cells. The stem cells were able to transform into both types of placental cell, either invading through a three-dimensional gel that mimics the wall of the uterus or making pregnancy hormones. There is a real need for a renewable supply of placental cells in pregnancy research. Animal placentas are not the same as human ones, so it is not possible to learn everything about human pregnancy from animal models. A renewable supply of trophoblast stem cells could aid in studying how the placenta forms and why this process sometimes goes wrong. This could help researchers to better understand miscarriage, pre-eclampsia and other conditions that affect the growth of an unborn baby. In the future, it may even be possible to make custom trophoblast stem cells to study the specific fertility issues of an individual.

Skeletal Muscle Regeneration in Advanced Diabetic Peripheral Neuropathy.

BACKGROUND: Decreased lean muscle mass in the lower extremity in diabetic peripheral neuropathy (DPN) is thought to contribute to altered joint loading, immobility, and disability. However, the mechanism behind this loss is unknown and could derive from a reduction in size of myofibers (atrophy), destruction of myofibers (degeneration), or both. Degenerative changes require participation of muscle stem (satellite) cells to regenerate lost myofibers and restore lean mass. Determining the degenerative state and residual regenerative capacity of DPN muscle will inform the utility of regeneration-targeted therapeutic strategies. METHODS: Biopsies were acquired from 2 muscles in 12 individuals with and without diabetic neuropathy undergoing below-knee amputation surgery. Biopsies were subdivided for histological analysis, transcriptional profiling, and satellite cell isolation and culture. RESULTS: Histological analysis revealed evidence of ongoing degeneration and regeneration in DPN muscles. Transcriptional profiling supports these findings, indicating significant upregulation of regeneration-related pathways. However, regeneration appeared to be limited in samples exhibiting the most severe structural pathology as only extremely small, immature regenerated myofibers were found. Immunostaining for satellite cells revealed a significant decrease in their relative frequency only in the subset with severe pathology. Similarly, a reduction in fusion in cultured satellite cells in this group suggests impairment in regenerative capacity in severe DPN pathology. CONCLUSION: DPN muscle exhibited features of degeneration with attempted regeneration. In the most severely pathological muscle samples, regeneration appeared to be stymied and our data suggest that this may be partly due to intrinsic dysfunction of the satellite cell pool in addition to extrinsic structural pathology (eg, nerve damage). CLINICAL RELEVANCE: Restoration of DPN muscle function for improved mobility and physical activity is a goal of surgical and rehabilitation clinicians. Identifying myofiber degeneration and compromised regeneration as contributors to dysfunction suggests that adjuvant cell-based therapies may improve clinical outcomes.

Transcriptional profiling of intramembranous and endochondral ossification after fracture in mice.

Bone fracture repair represents an important clinical challenge with nearly 1 million non-union fractures occurring annually in the U.S. Gene expression differs between non-union and healthy repair, suggesting there is a pattern of gene expression that is indicative of optimal repair. Despite this, the gene expression profile of fracture repair remains incompletely understood. In this work, we used RNA-seq of two well-established murine fracture models to describe gene expression of intramembranous and endochondral bone formation. We used top differentially expressed genes, enriched gene ontology terms and pathways, callus cellular phenotyping, and histology to describe and contrast these bone formation processes across time. Intramembranous repair, as modeled by ulnar stress fracture, and endochondral repair, as modeled by femur full fracture, exhibited vastly different transcriptional profiles throughout repair. Stress fracture healing had enriched differentially expressed genes associated with bone repair and osteoblasts, highlighting the strong osteogenic repair process of this model. Interestingly, the PI3K-Akt signaling pathway was one of only a few pathways uniquely enriched in stress fracture repair. Full fracture repair involved a higher level of inflammatory and immune cell related genes than did stress fracture repair. Full fracture repair also differed from stress fracture in a robust downregulation of ion channel genes following injury, the role of which in fracture repair is unclear. This study offers a broad description of gene expression in intramembranous and endochondral ossification across several time points throughout repair and suggests several potentially intriguing genes, pathways, and cells whose role in fracture repair requires further study.

FeatSNP: An Interactive Database for Brain-Specific Epigenetic Annotation of Human SNPs.

FeatSNP is an online tool and a curated database for exploring 81 million common SNPs' potential functional impact on the human brain. FeatSNP uses the brain transcriptomes of the human population to improve functional annotation of human SNPs by integrating transcription factor binding prediction, public eQTL information, and brain specific epigenetic landscape, as well as information of Topologically Associating Domains (TADs). FeatSNP supports both single and batched SNP searching, and its interactive user interface enables users to explore the functional annotations and generate publication-quality visualization results. FeatSNP is freely available on the internet at FeatSNP.org with all major web browsers supported.

Kdm6b regulates context-dependent hematopoietic stem cell self-renewal and leukemogenesis.

The histone demethylase KDM6B (JMJD3) is upregulated in blood disorders, suggesting that it may have important pathogenic functions. Here we examined the function of Kdm6b in hematopoietic stem cells (HSC) to evaluate its potential as a therapeutic target. Loss of Kdm6b lead to depletion of phenotypic and functional HSCs in adult mice, and Kdm6b is necessary for HSC self-renewal in response to inflammatory and proliferative stress. Loss of Kdm6b leads to a pro-differentiation poised state in HSCs due to the increased expression of the AP-1 transcription factor complex (Fos and Jun) and immediate early response (IER) genes. These gene expression changes occurred independently of chromatin modifications. Targeting AP-1 restored function of Kdm6b-deficient HSCs, suggesting that Kdm6b regulates this complex during HSC stress response. We also show Kdm6b supports developmental context-dependent leukemogenesis for T-cell acute lymphoblastic leukemia (T-ALL) and M5 acute myeloid leukemia (AML). Kdm6b is required for effective fetal-derived T-ALL and adult-derived AML, but not vice versa. These studies identify a crucial role for Kdm6b in regulating HSC self-renewal in different contexts, and highlight the potential of KDM6B as a therapeutic target in different hematopoietic malignancies.

Cellular and molecular characterization of multiplex autism in human induced pluripotent stem cell-derived neurons.

Background: Autism spectrum disorder (ASD) is a neurodevelopmental disorder with pronounced heritability in the general population. This is largely attributable to the effects of polygenic susceptibility, with inherited liability exhibiting distinct sex differences in phenotypic expression. Attempts to model ASD in human cellular systems have principally involved rare de novo mutations associated with ASD phenocopies. However, by definition, these models are not representative of polygenic liability, which accounts for the vast share of population-attributable risk. Methods: Here, we performed what is, to our knowledge, the first attempt to model multiplex autism using patient-derived induced pluripotent stem cells (iPSCs) in a family manifesting incremental degrees of phenotypic expression of inherited liability (absent, intermediate, severe). The family members share an inherited variant of uncertain significance (VUS) in GPD2, a gene that was previously associated with developmental disability but here is insufficient by itself to cause ASD. iPSCs from three first-degree relatives and an unrelated control were differentiated into both cortical excitatory (cExN) and cortical inhibitory (cIN) neurons, and cellular phenotyping and transcriptomic analysis were conducted. Results: cExN neurospheres from the two affected individuals were reduced in size, compared to those derived from unaffected related and unrelated individuals. This reduction was, at least in part, due to increased apoptosis of cells from affected individuals upon initiation of cExN neural induction. Likewise, cIN neural progenitor cells from affected individuals exhibited increased apoptosis, compared to both unaffected individuals. Transcriptomic analysis of both cExN and cIN neural progenitor cells revealed distinct molecular signatures associated with affectation, including the misregulation of suites of genes associated with neural development, neuronal function, and behavior, as well as altered expression of ASD risk-associated genes. Conclusions: We have provided evidence of morphological, physiological, and transcriptomic signatures of polygenic liability to ASD from an analysis of cellular models derived from a multiplex autism family. ASD is commonly inherited on the basis of additive genetic liability. Therefore, identifying convergent cellular and molecular phenotypes resulting from polygenic and monogenic susceptibility may provide a critical bridge for determining which of the disparate effects of rare highly deleterious mutations might also apply to common autistic syndromes.

Gon4l regulates notochord boundary formation and cell polarity underlying axis extension by repressing adhesion genes.

Anteroposterior (AP) axis extension during gastrulation requires embryonic patterning and morphogenesis to be spatiotemporally coordinated, but the underlying genetic mechanisms remain poorly understood. Here we define a role for the conserved chromatin factor Gon4l, encoded by ugly duckling (udu), in coordinating tissue patterning and axis extension during zebrafish gastrulation through direct positive and negative regulation of gene expression. Although identified as a recessive enhancer of impaired axis extension in planar cell polarity (PCP) mutants, udu functions in a genetically independent, partially overlapping fashion with PCP signaling to regulate mediolateral cell polarity underlying axis extension in part by promoting notochord boundary formation. Gon4l limits expression of the cell-cell and cell-matrix adhesion molecules EpCAM and Integrinα3b, excesses of which perturb the notochord boundary via tension-dependent and -independent mechanisms, respectively. By promoting formation of this AP-aligned boundary and associated cell polarity, Gon4l cooperates with PCP signaling to coordinate morphogenesis along the AP embryonic axis.

Regulatory networks specifying cortical interneurons from human embryonic stem cells reveal roles for CHD2 in interneuron development.

Cortical interneurons (cINs) modulate excitatory neuronal activity by providing local inhibition. During fetal development, several cIN subtypes derive from the medial ganglionic eminence (MGE), a transient ventral telencephalic structure. While altered cIN development contributes to neurodevelopmental disorders, the inaccessibility of human fetal brain tissue during development has hampered efforts to define molecular networks controlling this process. Here, we modified protocols for directed differentiation of human embryonic stem cells, obtaining efficient, accelerated production of MGE-like progenitors and MGE-derived cIN subtypes with the expected electrophysiological properties. We defined transcriptome changes accompanying this process and integrated these data with direct transcriptional targets of NKX2-1, a transcription factor controlling MGE specification. This analysis defined NKX2-1-associated genes with enriched expression during MGE specification and cIN differentiation, including known and previously unreported transcription factor targets with likely roles in MGE specification, and other target classes regulating cIN migration and function. NKX2-1-associated peaks were enriched for consensus binding motifs for NKX2-1, LHX, and SOX transcription factors, suggesting roles in coregulating MGE gene expression. Among the NKX2-1 direct target genes with cIN-enriched expression was CHD2, which encodes a chromatin remodeling protein mutated to cause human epilepsies. Accordingly, CHD2 deficiency impaired cIN specification and altered later electrophysiological function, while CHD2 coassociated with NKX2-1 at cis-regulatory elements and was required for their transactivation by NKX2-1 in MGE-like progenitors. This analysis identified several aspects of gene-regulatory networks underlying human MGE specification and suggested mechanisms by which NKX2-1 acts with chromatin remodeling activities to regulate gene expression programs underlying cIN development.

SRmapper: a fast and sensitive genome-hashing alignment tool.

UNLABELLED: Modern sequencing instruments have the capability to produce millions of short reads every day. The large number of reads produced in conjunction with variations between reads and reference genomic sequences caused both by legitimate differences, such as single-nucleotide polymorphisms and insertions/deletions (indels), and by sequencer errors make alignment a difficult and computationally expensive task, and many reads cannot be aligned. Here, we introduce a new alignment tool, SRmapper, which in tests using real data can align 10s of billions of base pairs from short reads to the human genome per computer processor day. SRmapper tolerates a higher number of mismatches than current programs based on Burrows-Wheeler transform and finds about the same number of alignments in 2-8× less time depending on read length (with higher performance gain for longer read length). The current version of SRmapper aligns both single and pair-end reads in base space fastq format and outputs alignments in Sequence Alignment/Map format. SRmapper uses a probabilistic approach to set a default number of mismatches allowed and determines alignment quality. SRmapper's memory footprint (∼2.5 GB) is small enough that it can be run on a computer with 4 GB of random access memory for a genome the size of a human. Finally, SRmapper is designed so that its function can be extended to finding small indels as well as long deletions and chromosomal translocations in future versions. AVAILABILITY: http://www.umsl.edu/∼wongch/software.html.