Heterozygous loss-of-function SMC3 variants are associated with variable growth and developmental features.

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping.

Zebrafish models of candidate human epilepsy-associated genes provide evidence of hyperexcitability.

Hundreds of novel candidate human epilepsy-associated genes have been identified thanks to advancements in next-generation sequencing and large genome-wide association studies, but establishing genetic etiology requires functional validation. We generated a list of >2200 candidate epilepsy-associated genes, of which 81 were determined suitable for the generation of loss-of-function zebrafish models via CRISPR/Cas9 gene editing. Of those 81 crispants, 48 were successfully established as stable mutant lines and assessed for seizure-like swim patterns in a primary F2 screen. Evidence of seizure-like behavior was present in 5 (arfgef1, kcnd2, kcnv1, ubr5, wnt8b) of the 48 mutant lines assessed. Further characterization of those 5 lines provided evidence for epileptiform activity via electrophysiology in kcnd2 and wnt8b mutants. Additionally, arfgef1 and wnt8b mutants showed a decrease in the number of inhibitory interneurons in the optic tectum of larval animals. Furthermore, RNAseq revealed convergent transcriptional abnormalities between mutant lines, consistent with their developmental defects and hyperexcitable phenotypes. These zebrafish models provide strongest experimental evidence supporting the role of ARFGEF1, KCND2, and WNT8B in human epilepsy and further demonstrate the utility of this model system for evaluating candidate human epilepsy genes.

Heterozygous loss-of-function SMC3 variants are associated with variable and incompletely penetrant growth and developmental features.

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 13 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated a milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, some instead having intriguing symptomatologies with rational biological links to SMC3 including bone marrow failure, acute myeloid leukemia, and Coats retinal vasculopathy. Analyses of transcriptomic and epigenetic data suggest that SMC3 pLoF variants reduce SMC3 expression but do not result in a blood DNA methylation signature clustering with that of CdLS, and that the global transcriptional signature of SMC3 loss is model-dependent. Our finding of substantial population-scale LoF intolerance in concert with variable penetrance in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multi-layered genomic data paired with careful phenotyping.

Utility of Exome Sequencing for Diagnosis in Unexplained Pediatric-Onset Epilepsy.

Importance: Genomic advances inform our understanding of epilepsy and can be translated to patients as precision diagnoses that influence clinical treatment, prognosis, and counseling. Objective: To delineate the genetic landscape of pediatric epilepsy and clinical utility of genetic diagnoses for patients with epilepsy. Design, Setting, and Participants: This cohort study used phenotypic data from medical records and treating clinicians at a pediatric hospital to identify patients with unexplained pediatric-onset epilepsy. Exome sequencing was performed for 522 patients and available biological parents, and sequencing data were analyzed for single nucleotide variants (SNVs) and copy number variants (CNVs). Variant pathogenicity was assessed, patients were provided with their diagnostic results, and clinical utility was evaluated. Patients were enrolled from August 2018 to October 2021, and data were analyzed through December 2022. Exposures: Phenotypic features associated with diagnostic genetic results. Main Outcomes and Measures: Main outcomes included diagnostic yield and clinical utility. Diagnostic findings included variants curated as pathogenic, likely pathogenic (PLP), or diagnostic variants of uncertain significance (VUS) with clinical features consistent with the involved gene's associated phenotype. The proportion of the cohort with diagnostic findings, the genes involved, and their clinical utility, defined as impact on clinical treatment, prognosis, or surveillance, are reported. Results: A total of 522 children (269 [51.5%] male; mean [SD] age at seizure onset, 1.2 [1.4] years) were enrolled, including 142 children (27%) with developmental epileptic encephalopathy and 263 children (50.4%) with intellectual disability. Of these, 100 participants (19.2%) had identifiable genetic explanations for their seizures: 89 participants had SNVs (87 germline, 2 somatic mosaic) involving 69 genes, and 11 participants had CNVs. The likelihood of identifying a genetic diagnosis was highest in patients with intellectual disability (adjusted odds ratio [aOR], 2.44; 95% CI, 1.40-4.26), early onset seizures (aOR, 0.93; 95% CI, 0.88-0.98), and motor impairment (aOR, 2.19; 95% CI 1.34-3.58). Among 43 patients with apparently de novo variants, 2 were subsequently determined to have asymptomatic parents harboring mosaic variants. Of 71 patients who received diagnostic results and were followed clinically, 29 (41%) had documented clinical utility resulting from their genetic diagnoses. Conclusions and Relevance: These findings suggest that pediatric-onset epilepsy is genetically heterogeneous and that some patients with previously unexplained pediatric-onset epilepsy had genetic diagnoses with direct clinical implications.

Spatial transcriptomics reveals a cnidarian segment polarity program in Nematostella vectensis.

During early animal evolution, the emergence of axially polarized segments was central to the diversification of complex bilaterian body plans. Nevertheless, precisely how and when segment polarity pathways arose remains obscure. Here, we demonstrate the molecular basis for segment polarization in developing larvae of the sea anemone Nematostella vectensis. Utilizing spatial transcriptomics, we first constructed a 3D gene expression atlas of developing larval segments. Capitalizing on accurate in silico predictions, we identified Lbx and Uncx, conserved homeodomain-containing genes that occupy opposing subsegmental domains under the control of both bone morphogenetic protein (BMP) signaling and the Hox-Gbx cascade. Functionally, Lbx mutagenesis eliminated all molecular evidence of segment polarization at the larval stage and caused an aberrant mirror-symmetric pattern of retractor muscles (RMs) in primary polyps. These results demonstrate the molecular basis for segment polarity in a non-bilaterian animal, suggesting that polarized metameric structures were present in the Cnidaria-Bilateria common ancestor over 600 million years ago.

Spatial transcriptomics reveals a conserved segment polarity program that governs muscle patterning in Nematostella vectensis.

During early animal evolution, the emergence of axially-polarized segments was central to the diversification of complex bilaterian body plans. Nevertheless, precisely how and when segment polarity pathways arose remains obscure. Here we demonstrate the molecular basis for segment polarization in developing larvae of the pre-bilaterian sea anemone Nematostella vectensis . Utilizing spatial transcriptomics, we first constructed a 3-D gene expression atlas of developing larval segments. Capitalizing on accurate in silico predictions, we identified Lbx and Uncx, conserved homeodomain-containing genes that occupy opposing subsegmental domains under the control of both BMP signaling and the Hox-Gbx cascade. Functionally, Lbx mutagenesis eliminated all molecular evidence of segment polarization at larval stage and caused an aberrant mirror-symmetric pattern of retractor muscles in primary polyps. These results demonstrate the molecular basis for segment polarity in a pre-bilaterian animal, suggesting that polarized metameric structures were present in the Cnidaria-Bilateria common ancestor over 600 million years ago. Highlights: Nematostella endomesodermal tissue forms metameric segments and displays a transcriptomic profile similar to that observed in bilaterian mesoderm Construction of a comprehensive 3-D gene expression atlas enables systematic dissection of segmental identity in endomesoderm Lbx and Uncx , two conserved homeobox-containing genes, establish segment polarity in Nematostella The Cnidarian-Bilaterian common ancestor likely possessed the genetic toolkit to generate polarized metameric structures.

The SAGA core module is critical during Drosophila oogenesis and is broadly recruited to promoters.

The Spt/Ada-Gcn5 Acetyltransferase (SAGA) coactivator complex has multiple modules with different enzymatic and non-enzymatic functions. How each module contributes to gene expression is not well understood. During Drosophila oogenesis, the enzymatic functions are not equally required, which may indicate that different genes require different enzymatic functions. An analogy for this phenomenon is the handyman principle: while a handyman has many tools, which tool he uses depends on what requires maintenance. Here we analyzed the role of the non-enzymatic core module during Drosophila oogenesis, which interacts with TBP. We show that depletion of SAGA-specific core subunits blocked egg chamber development at earlier stages than depletion of enzymatic subunits. These results, as well as additional genetic analyses, point to an interaction with TBP and suggest a differential role of SAGA modules at different promoter types. However, SAGA subunits co-occupied all promoter types of active genes in ChIP-seq and ChIP-nexus experiments, and the complex was not specifically associated with distinct promoter types in the ovary. The high-resolution genomic binding profiles were congruent with SAGA recruitment by activators upstream of the start site, and retention on chromatin by interactions with modified histones downstream of the start site. Our data illustrate that a distinct genetic requirement for specific components may conceal the fact that the entire complex is physically present and suggests that the biological context defines which module functions are critical.

A mouse-specific retrotransposon drives a conserved Cdk2ap1 isoform essential for development.

Retrotransposons mediate gene regulation in important developmental and pathological processes. Here, we characterized the transient retrotransposon induction during preimplantation development of eight mammals. Induced retrotransposons exhibit similar preimplantation profiles across species, conferring gene regulatory activities, particularly through long terminal repeat (LTR) retrotransposon promoters. A mouse-specific MT2B2 retrotransposon promoter generates an N-terminally truncated Cdk2ap1ΔN that peaks in preimplantation embryos and promotes proliferation. In contrast, the canonical Cdk2ap1 peaks in mid-gestation and represses cell proliferation. This MT2B2 promoter, whose deletion abolishes Cdk2ap1ΔN production, reduces cell proliferation and impairs embryo implantation, is developmentally essential. Intriguingly, Cdk2ap1ΔN is evolutionarily conserved in sequence and function yet is driven by different promoters across mammals. The distinct preimplantation Cdk2ap1ΔN expression in each mammalian species correlates with the duration of its preimplantation development. Hence, species-specific transposon promoters can yield evolutionarily conserved, alternative protein isoforms, bestowing them with new functions and species-specific expression to govern essential biological divergence.

Transcript assembly improves expression quantification of transposable elements in single-cell RNA-seq data.

Transposable elements (TEs) are an integral part of the host transcriptome. TE-containing noncoding RNAs (ncRNAs) show considerable tissue specificity and play important roles during development, including stem cell maintenance and cell differentiation. Recent advances in single-cell RNA-seq (scRNA-seq) revolutionized cell type-specific gene expression analysis. However, effective scRNA-seq quantification tools tailored for TEs are lacking, limiting our ability to dissect TE expression dynamics at single-cell resolution. To address this issue, we established a TE expression quantification pipeline that is compatible with scRNA-seq data generated across multiple technology platforms. We constructed TE-containing ncRNA references using bulk RNA-seq data and showed that quantifying TE expression at the transcript level effectively reduces noise. As proof of principle, we applied this strategy to mouse embryonic stem cells and successfully captured the expression profile of endogenous retroviruses in single cells. We further expanded our analysis to scRNA-seq data from early stages of mouse embryogenesis. Our results illustrated the dynamic TE expression at preimplantation stages and revealed 146 TE-containing ncRNA transcripts with substantial tissue specificity during gastrulation and early organogenesis.

A Role for FACT in RNA Polymerase II Promoter-Proximal Pausing.

FACT (facilitates chromatin transcription) is an evolutionarily conserved histone chaperone that was initially identified as an activity capable of promoting RNA polymerase II (Pol II) transcription through nucleosomes in vitro. In this report, we describe a global analysis of FACT function in Pol II transcription in Drosophila. We present evidence that loss of FACT has a dramatic impact on Pol II elongation-coupled processes including histone H3 lysine 4 (H3K4) and H3K36 methylation, consistent with a role for FACT in coordinating histone modification and chromatin architecture during Pol II transcription. Importantly, we identify a role for FACT in the maintenance of promoter-proximal Pol II pausing, a key step in transcription activation in higher eukaryotes. These findings bring to light a broader role for FACT in the regulation of Pol II transcription.

Reporter-ChIP-nexus reveals strong contribution of the Drosophila initiator sequence to RNA polymerase pausing.

RNA polymerase II (Pol II) pausing is a general regulatory step in transcription, yet the stability of paused Pol II varies widely between genes. Although paused Pol II stability correlates with core promoter elements, the contribution of individual sequences remains unclear, in part because no rapid assay is available for measuring the changes in Pol II pausing as a result of altered promoter sequences. Here, we overcome this hurdle by showing that ChIP-nexus captures the endogenous Pol II pausing on transfected plasmids. Using this reporter-ChIP-nexus assay in Drosophila cells, we show that the pausing stability is influenced by downstream promoter sequences, but that the strongest contribution to Pol II pausing comes from the initiator sequence, in which a single nucleotide, a G at the +2 position, is critical for stable Pol II pausing. These results establish reporter-ChIP-nexus as a valuable tool to analyze Pol II pausing.

Author Correction: Suppression of m6A reader Ythdf2 promotes hematopoietic stem cell expansion.

In the initial published version of this article, there was an inadvertent omission from the Acknowledgements that this work was supported by Stowers Institute for Medical Research (SIMR-1004) and NIH National Cancer Institute grant to University of Kansas Cancer Center (P30 CA168524). This omission does not affect the description of the results or the conclusions of this work.

Suppression of m6A reader Ythdf2 promotes hematopoietic stem cell expansion.

Transplantation of hematopoietic stem cells (HSCs) from human umbilical cord blood (hUCB) holds great promise for treating a broad spectrum of hematological disorders including cancer. However, the limited number of HSCs in a single hUCB unit restricts its widespread use. Although extensive efforts have led to multiple methods for ex vivo expansion of human HSCs by targeting single molecules or pathways, it remains unknown whether it is possible to simultaneously manipulate the large number of targets essential for stem cell self-renewal. Recent studies indicate that N6-methyladenosine (m6A) modulates the expression of a group of mRNAs critical for stem cell-fate determination by influencing their stability. Among several m6A readers, YTHDF2 is recognized as promoting targeted mRNA decay. However, the physiological functions of YTHDF2 in adult stem cells are unknown. Here we show that following the conditional knockout (KO) of mouse Ythdf2 the numbers of functional HSC were increased without skewing lineage differentiation or leading to hematopoietic malignancies. Furthermore, knockdown (KD) of human YTHDF2 led to more than a 10-fold increase in the ex vivo expansion of hUCB HSCs, a fivefold increase in colony-forming units (CFUs), and more than an eightfold increase in functional hUCB HSCs in the secondary serial of a limiting dilution transplantation assay. Mapping of m6A in RNAs from mouse hematopoietic stem and progenitor cells (HSPCs) as well as from hUCB HSCs revealed its enrichment in mRNAs encoding transcription factors critical for stem cell self-renewal. These m6A-marked mRNAs were recognized by Ythdf2 and underwent decay. In Ythdf2 KO HSPCs and YTHDF2 KD hUCB HSCs, these mRNAs were stabilized, facilitating HSC expansion. Knocking down one of YTHDF2's key targets, Tal1 mRNA, partially rescued the phenotype. Our study provides the first demonstration of the function of YTHDF2 in adult stem cell maintenance and identifies its important role in regulating HSC ex vivo expansion by regulating the stability of multiple mRNAs critical for HSC self-renewal, thus identifying potential for future clinical applications.

Capicua controls Toll/IL-1 signaling targets independently of RTK regulation.

The HMG-box protein Capicua (Cic) is a conserved transcriptional repressor that functions downstream of receptor tyrosine kinase (RTK) signaling pathways in a relatively simple switch: In the absence of signaling, Cic represses RTK-responsive genes by binding to nearly invariant sites in DNA, whereas activation of RTK signaling down-regulates Cic activity, leading to derepression of its targets. This mechanism controls gene expression in both Drosophila and mammals, but whether Cic can also function via other regulatory mechanisms remains unknown. Here, we characterize an RTK-independent role of Cic in regulating spatially restricted expression of Toll/IL-1 signaling targets in Drosophila embryogenesis. We show that Cic represses those targets by binding to suboptimal DNA sites of lower affinity than its known consensus sites. This binding depends on Dorsal/NF-κB, which translocates into the nucleus upon Toll activation and binds next to the Cic sites. As a result, Cic binds to and represses Toll targets only in regions with nuclear Dorsal. These results reveal a mode of Cic regulation unrelated to the well-established RTK/Cic depression axis and implicate cooperative binding in conjunction with low-affinity binding sites as an important mechanism of enhancer regulation. Given that Cic plays a role in many developmental and pathological processes in mammals, our results raise the possibility that some of these Cic functions are independent of RTK regulation and may depend on cofactor-assisted DNA binding.

Paused RNA polymerase II inhibits new transcriptional initiation.

RNA polymerase II (Pol II) pauses downstream of the transcription initiation site before beginning productive elongation. This pause is a key component of metazoan gene expression regulation. Some promoters have a strong disposition for Pol II pausing and often mediate faster, more synchronous changes in expression. This requires multiple rounds of transcription and thus cannot rely solely on pause release. However, it is unclear how pausing affects the initiation of new transcripts during consecutive rounds of transcription. Using our recently developed ChIP-nexus method, we find that Pol II pausing inhibits new initiation. We propose that paused Pol II helps prevent new initiation between transcription bursts, which may reduce noise.

Droplet topology control of Janus emulsion prepared in one-step high energy mixing.

Droplet topology of a Janus emulsion in a vegetable oil (VO)/silicone oil (SO)/Tween 80 aqueous solution (Aq) system prepared in a one-step high energy mixing was investigated, mainly by image observation. Quantitative information of the topology was analyzed referring to the curvature of VO/SO interface, the location of contact plane, and the volume ratio of VO/SO within individual droplets. The results show that the "stable Janus emulsion" region in the phase map enlarges with surfactant concentration. The average volume ratio of two oil lobes within an individual Janus droplet agrees with the emulsion composition in the "stable Janus emulsion" region, which means that the droplet topology can be controlled by the emulsion preparation process within realistic limits. The volume ratio of VO/SO within individual droplet ranges from about 0.54 to 0.17 in the VO/SO/3 wt% Tween 80(Aq) system, beyond which separate VO and SO droplets are observed. The topology of a Janus droplet is found to be determined by both the contact angle of three liquids in the contact line and the location of the contact plane. The contact angle of the oil cap is determined by the interfacial tension referring to the local equilibrium. The location of the contact plane is the dominant factor determining the volume ratio of two oil lobes. Composition change in the emulsion results in the corresponding tune of the location of contact plane and subsequently, the volume ratio of two oils within Janus droplets.

A global change in RNA polymerase II pausing during the Drosophila midblastula transition.

Massive zygotic transcription begins in many organisms during the midblastula transition when the cell cycle of the dividing egg slows down. A few genes are transcribed before this stage but how this differential activation is accomplished is still an open question. We have performed ChIP-seq experiments on tightly staged Drosophila embryos and show that massive recruitment of RNA polymerase II (Pol II) with widespread pausing occurs de novo during the midblastula transition. However, ∼100 genes are strongly occupied by Pol II before this timepoint and most of them do not show Pol II pausing, consistent with a requirement for rapid transcription during the fast nuclear cycles. This global change in Pol II pausing correlates with distinct core promoter elements and associates a TATA-enriched promoter with the rapid early transcription. This suggests that promoters are differentially used during the zygotic genome activation, presumably because they have distinct dynamic properties. DOI:http://dx.doi.org/10.7554/eLife.00861.001.