Alterations to the remote control of Shh gene expression cause congenital abnormalities

A Bayesian Change Point Model for Dynamic Alternative Transcription Start Site Usage During Cellular Differentiation

ABSTRACT An alternative transcription start site (ATSS) is a major driving force for increasing the complexity of transcripts in human tissues. As a transcriptional regulatory mechanism, ATSS has biological significance. Many studies have confirmed that ATSS plays an important role in diseases and cell development and differentiation. However, exploration of its dynamic mechanisms remains insufficient. Identifying ATSS change points during cell differentiation is critical for elucidating potential dynamic mechanisms. For relative ATSS usage as percentage data, the existing methods lack sensitivity to detect the change point for ATSS longitudinal data. In addition, some methods have strict requirements for data distribution and cannot be applied to deal with this problem. In this study, the Bayesian change point detection model was first constructed using reparameterization techniques for two parameters of a beta distribution for the percentage data type, and the posterior distributions of parameters and change points were obtained using Markov Chain Monte Carlo (MCMC) sampling. With comprehensive simulation studies, the performance of the Bayesian change point detection model is found to be consistently powerful and robust across most scenarios with different sample sizes and beta distributions. Second, differential ATSS events in the real data, whose change points were identified using our method, were clustered according to their change points. Last, for each change point, pathway and transcription factor motif analyses were performed on its differential ATSS events. The results of our analyses demonstrated the effectiveness of the Bayesian change point detection model and provided biological insights into cell differentiation.

POSTRE: a tool to predict the pathological effects of human structural variants

Understanding the pathological impact of non-coding genetic variation is a major challenge in medical genetics. Accumulating evidences indicate that a significant fraction of genetic alterations, including structural variants (SVs), can cause human disease by altering the function of non-coding regulatory elements, such as enhancers. In the case of SVs, described pathomechanisms include changes in enhancer dosage and long-range enhancer-gene communication. However, there is still a clear gap between the need to predict and interpret the medical impact of non-coding variants, and the existence of tools to properly perform these tasks. To reduce this gap, we have developed POSTRE (Prediction Of STRuctural variant Effects), a computational tool to predict the pathogenicity of SVs implicated in a broad range of human congenital disorders. By considering disease-relevant cellular contexts, POSTRE identifies SVs with either coding or long-range pathological consequences with high specificity and sensitivity. Furthermore, POSTRE not only identifies pathogenic SVs, but also predicts the disease-causative genes and the underlying pathological mechanism (e.g, gene deletion, enhancer disconnection, enhancer adoption, etc.). POSTRE is available at https://github.com/vicsanga/Postre.

Function and Constraint in Enhancer Sequences with Multiple Evolutionary Origins

Thousands of human gene regulatory enhancers are composed of sequences with multiple evolutionary origins. These evolutionarily "complex" enhancers consist of older "core" sequences and younger "derived" sequences. However, the functional relationship between the sequences of different evolutionary origins within complex enhancers is poorly understood. We evaluated the function, selective pressures, and sequence variation across core and derived components of human complex enhancers. We find that both components are older than expected from the genomic background, and complex enhancers are enriched for core and derived sequences of similar evolutionary ages. Both components show strong evidence of biochemical activity in massively parallel report assays. However, core and derived sequences have distinct transcription factor (TF)-binding preferences that are largely similar across evolutionary origins. As expected, given these signatures of function, both core and derived sequences have substantial evidence of purifying selection. Nonetheless, derived sequences exhibit weaker purifying selection than adjacent cores. Derived sequences also tolerate more common genetic variation and are enriched compared with cores for expression quantitative trait loci associated with gene expression variability in human populations. In conclusion, both core and derived sequences have strong evidence of gene regulatory function, but derived sequences have distinct constraint profiles, TF-binding preferences, and tolerance to variation compared with cores. We propose that the step-wise integration of younger derived with older core sequences has generated regulatory substrates with robust activity and the potential for functional variation. Our analyses demonstrate that synthesizing study of enhancer evolution and function can aid interpretation of regulatory sequence activity and functional variation across human populations.

Enhancer-gene specificity in development and disease

Enhancers control the establishment of spatiotemporal gene expression patterns throughout development. Over the past decade, the development of new technologies has improved our capacity to link enhancers with their target genes based on their colocalization within the same topological domains. However, the mechanisms that regulate how enhancers specifically activate some genes but not others within a given domain remain unclear. In this Review, we discuss recent insights into the factors controlling enhancer specificity, including the genetic composition of enhancers and promoters, the linear and 3D distance between enhancers and their target genes, and cell-type specific chromatin landscapes. We also discuss how elucidating the molecular principles of enhancer specificity might help us to better understand and predict the pathological consequences of human genetic, epigenetic and structural variants.

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

Gene-teratogen interactions influence the penetrance of birth defects by altering Hedgehog signaling strength

Birth defects result from interactions between genetic and environmental factors, but the mechanisms remain poorly understood. We find that mutations and teratogens interact in predictable ways to cause birth defects by changing target cell sensitivity to Hedgehog (Hh) ligands. These interactions converge on a membrane protein complex, the MMM complex, that promotes degradation of the Hh transducer Smoothened (SMO). Deficiency of the MMM component MOSMO results in elevated SMO and increased Hh signaling, causing multiple birth defects. In utero exposure to a teratogen that directly inhibits SMO reduces the penetrance and expressivity of birth defects in Mosmo-/- embryos. Additionally, tissues that develop normally in Mosmo-/- embryos are refractory to the teratogen. Thus, changes in the abundance of the protein target of a teratogen can change birth defect outcomes by quantitative shifts in Hh signaling. Consequently, small molecules that re-calibrate signaling strength could be harnessed to rescue structural birth defects.

Enhancer redundancy in development and disease

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

Transcriptional response of human articular chondrocytes treated with fibronectin fragments: an in vitro model of the osteoarthritis phenotype

OBJECTIVE

Fibronectin is a matrix protein that is fragmented during cartilage degradation in osteoarthritis (OA). Treatment of chondrocytes with fibronectin fragments (FN-f) has been used to model OA in vitro, but the system has not been fully characterized. This study sought to define the transcriptional response of chondrocytes to FN-f, and directly compare it to responses traditionally observed in OA.

DESIGN

Normal human femoral chondrocytes isolated from tissue donors were treated with either FN-f or PBS (control) for 3, 6, or 18 h. RNA-seq libraries were compared between time-matched FN-f and control samples in order to identify changes in gene expression over time. Differentially expressed genes were compared to a published OA gene set and used for pathway, transcription factor motif, and kinome analysis.

RESULTS

FN-f treatment resulted in 3,914 differentially expressed genes over the time course. Genes that are up- or downregulated in OA were significantly up- (P < 0.00001) or downregulated (P < 0.0004) in response to FN-f. Early response genes were involved in proinflammatory pathways, whereas many late response genes were involved in ferroptosis. The promoters of upregulated genes were enriched for NF-κB, AP-1, and IRF motifs. Highly upregulated kinases included CAMK1G, IRAK2, and the uncharacterized kinase DYRK3, while growth factor receptors TGFBR2 and FGFR2 were downregulated.

CONCLUSIONS

FN-f treatment of normal human articular chondrocytes recapitulated many key aspects of the OA chondrocyte phenotype. This in vitro model is promising for future OA studies, especially considering its compatibility with genomics and genome-editing techniques.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

A multidisciplinary review of triphalangeal thumb

Despite being a rare congenital limb anomaly, triphalangeal thumb is a subject of research in various scientific fields, providing new insights in clinical research and evolutionary biology. The findings of triphalangeal thumb can be predictive for other congenital anomalies as part of an underlying syndrome. Furthermore, triphalangeal thumb is still being used as a model in molecular genetics to study gene regulation by long-range regulatory elements. We present a review that summarizes a number of scientifically relevant topics that involve the triphalangeal thumb phenotype. Future initiatives involving multidisciplinary teams collaborating in the field of triphalangeal thumb research can lead to a better understanding of the pathogenesis and molecular mechanisms of this condition as well as other congenital upper limb anomalies.

Evidence that alternative transcriptional initiation is largely nonadaptive

Alternative transcriptional initiation (ATI) refers to the frequent observation that one gene has multiple transcription start sites (TSSs). Although this phenomenon is thought to be adaptive, the specific advantage is rarely known. Here, we propose that each gene has one optimal TSS and that ATI arises primarily from imprecise transcriptional initiation that could be deleterious. This error hypothesis predicts that (i) the TSS diversity of a gene reduces with its expression level; (ii) the fractional use of the major TSS increases, but that of each minor TSS decreases, with the gene expression level; and (iii) cis-elements for major TSSs are selectively constrained, while those for minor TSSs are not. By contrast, the adaptive hypothesis does not make these predictions a priori. Our analysis of human and mouse transcriptomes confirms each of the three predictions. These and other findings strongly suggest that ATI predominantly results from molecular errors, requiring a major revision of our understanding of the precision and regulation of transcription.

The transcription of a gene may start from one of several transcription start sites, a phenomenon known as alternative transcriptional initiation. Contrary to common belief, this study shows that variation of the transcription start site of a given gene is nonadaptive and is largely attributable to transcriptional initiation error that is typically deleterious.

Multiple surveys of transcriptional initiation showed that mammalian genes typically have multiple transcription start sites such that transcription is initiated from any one of these sites. Many researchers believe that this phenomenon is adaptive because it allows production of multiple transcripts, from the same gene, that potentially vary in function or post-transcriptional regulation. Nevertheless, it is also possible that each gene has only one optimal transcription start site and that alternative transcriptional initiation arises primarily from molecular errors that are slightly deleterious. This error hypothesis makes a series of predictions about the amount of transcription start site diversity per gene, relative uses of the various start sites of a gene, among-tissue and across-species differences in start site usage, and the evolutionary conservation of cis-regulatory elements of various start sites, all of which are verified in our analyses of genome-wide transcription start site data from the human and mouse. These findings strongly suggest that alternative transcriptional initiation largely reflects molecular errors instead of molecular adaptations and require a rethink of the precision and regulation of transcription.

Classification of congenital upper limb anomalies: towards improved communication, diagnosis, and discovery

Recently the International Federation of Societies for Surgery of the Hand replaced the Swanson scheme for classifying congenital upper limb anomalies with the Oberg, Manske, Tonkin (OMT) classification. This review explores the reasons for this change after nearly 50 years of using the Swanson classification. In particular, it documents the state of our understanding regarding genetics and limb development at the time Swanson generated his classification. It also describes the continued progress in clinical genetics and developmental biology. Such progress drives the need to embrace and incorporate these changes within a new classification scheme; one that will improve communication, diagnosis, and support further discovery of the pathogenesis of congenital hand anomalies.

De novo mutations in regulatory elements in neurodevelopmental disorders

We previously estimated that 42% of patients with severe developmental disorders carry pathogenic de novo mutations in coding sequences. The role of de novo mutations in regulatory elements affecting genes associated with developmental disorders, or other genes, has been essentially unexplored. We identified de novo mutations in three classes of putative regulatory elements in almost 8,000 patients with developmental disorders. Here we show that de novo mutations in highly evolutionarily conserved fetal brain-active elements are significantly and specifically enriched in neurodevelopmental disorders. We identified a significant twofold enrichment of recurrently mutated elements. We estimate that, genome-wide, 1-3% of patients without a diagnostic coding variant carry pathogenic de novo mutations in fetal brain-active regulatory elements and that only 0.15% of all possible mutations within highly conserved fetal brain-active elements cause neurodevelopmental disorders with a dominant mechanism. Our findings represent a robust estimate of the contribution of de novo mutations in regulatory elements to this genetically heterogeneous set of disorders, and emphasize the importance of combining functional and evolutionary evidence to identify regulatory causes of genetic disorders.

Alternative intronic promoters in development and disease

Approximately 20,000 mammalian genes are estimated to encode between 250 thousand and 1 million different proteins. This enormous diversity of the mammalian proteome is caused by the ability of a single-gene locus to encode multiple protein isoforms. Protein isoforms encoded by one gene locus can be functionally distinct, and they can even have antagonistic functions. One of the mechanisms involved in creating this proteome complexity is alternative promoter usage. Alternative intronic promoters are located downstream from their canonical counterparts and drive the expression of alternative RNA isoforms that lack upstream exons. These upstream exons can encode some important functional domains, and proteins encoded by alternative mRNA isoforms can be thus functionally distinct from the full-length protein encoded by canonical mRNA isoforms. Since any misbalance of functionally distinct protein isoforms is likely to have detrimental consequences for the cell and the whole organism, their expression must be precisely regulated. Misregulation of alternative intronic promoters is frequently associated with various developmental defects and diseases including cancer, and it is becoming increasingly clear that this phenomenon deserves more attention.

Sonic Hedgehog Signaling in Limb Development

The gene encoding the secreted protein Sonic hedgehog (Shh) is expressed in the polarizing region (or zone of polarizing activity), a small group of mesenchyme cells at the posterior margin of the vertebrate limb bud. Detailed analyses have revealed that Shh has the properties of the long sought after polarizing region morphogen that specifies positional values across the antero-posterior axis (e.g., thumb to little finger axis) of the limb. Shh has also been shown to control the width of the limb bud by stimulating mesenchyme cell proliferation and by regulating the antero-posterior length of the apical ectodermal ridge, the signaling region required for limb bud outgrowth and the laying down of structures along the proximo-distal axis (e.g., shoulder to digits axis) of the limb. It has been shown that Shh signaling can specify antero-posterior positional values in limb buds in both a concentration- (paracrine) and time-dependent (autocrine) fashion. Currently there are several models for how Shh specifies positional values over time in the limb buds of chick and mouse embryos and how this is integrated with growth. Extensive work has elucidated downstream transcriptional targets of Shh signaling. Nevertheless, it remains unclear how antero-posterior positional values are encoded and then interpreted to give the particular structure appropriate to that position, for example, the type of digit. A distant cis-regulatory enhancer controls limb-bud-specific expression of Shh and the discovery of increasing numbers of interacting transcription factors indicate complex spatiotemporal regulation. Altered Shh signaling is implicated in clinical conditions with congenital limb defects and in the evolution of the morphological diversity of vertebrate limbs.

Progressive Loss of Function in a Limb Enhancer during Snake Evolution

The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP.

Regulation of Hedgehog Signalling Inside and Outside the Cell

The hedgehog (Hh) signalling pathway is conserved throughout metazoans and plays an important regulatory role in both embryonic development and adult homeostasis. Many levels of regulation exist that control the release, reception, and interpretation of the hedgehog signal. The fatty nature of the Shh ligand means that it tends to associate tightly with the cell membrane, and yet it is known to act as a morphogen that diffuses to elicit pattern formation. Heparan sulfate proteoglycans (HSPGs) play a major role in the regulation of Hh distribution outside the cell. Inside the cell, the primary cilium provides an important hub for processing the Hh signal in vertebrates. This review will summarise the current understanding of how the Hh pathway is regulated from ligand production, release, and diffusion, through to signal reception and intracellular transduction.

Use of a Conditional Ubr5 Mutant Allele to Investigate the Role of an N-End Rule Ubiquitin-Protein Ligase in Hedgehog Signalling and Embryonic Limb Development

Hedgehog (Hh) signalling is a potent regulator of cell fate and function. While much is known about the events within a Hh-stimulated cell, far less is known about the regulation of Hh-ligand production. Drosophila Hyperplastic Discs (Hyd), a ubiquitin-protein ligase, represents one of the few non-transcription factors that independently regulates both hh mRNA expression and pathway activity. Using a murine embryonic stem cell system, we revealed that shRNAi of the mammalian homologue of hyd, Ubr5, effectively prevented retinoic-acid-induced Sonic hedgehog (Shh) expression. We next investigated the UBR5:Hh signalling relationship in vivo by generating and validating a mouse bearing a conditional Ubr5 loss-of-function allele. Conditionally deleting Ubr5 in the early embryonic limb-bud mesenchyme resulted in a transient decrease in Indian hedgehog ligand expression and decreased Hh pathway activity, around E13.5. Although Ubr5-deficient limbs and digits were, on average, shorter than control limbs, the effects were not statistically significant. Hence, while loss of UBR5 perturbed Hedgehog signalling in the developing limb, there were no obvious morphological defects. In summary, we report the first conditional Ubr5 mutant mouse and provide evidence for a role for UBR5 in influencing Hh signalling, but are uncertain to whether the effects on Hedgehog signaling were direct (cell autonomous) or indirect (non-cell-autonomous). Elaboration of the cellular/molecular mechanism(s) involved may help our understanding on diseases and developmental disorders associated with aberrant Hh signalling.

Towards understanding pre-mRNA splicing mechanisms and the role of SR proteins

Alternative pre-mRNA splicing provides a source of vast protein diversity by removing non-coding sequences (introns) and accurately linking different exonic regions in the correct reading frame. The regulation of alternative splicing is essential for various cellular functions in both pathological and physiological conditions. In eukaryotic cells, this process is commonly used to increase proteomic diversity and to control gene expression either co- or post-transcriptionally. Alternative splicing occurs within a megadalton-sized, multi-component machine consisting of RNA and proteins; during the splicing process, this complex undergoes dynamic changes via RNA-RNA, protein-protein and RNA-protein interactions. Co-transcriptional splicing functionally integrates the transcriptional machinery, thereby enabling the two processes to influence one another, whereas post-transcriptional splicing facilitates the coupling of RNA splicing with post-splicing events. This review addresses the structural aspects of spliceosomes and the mechanistic implications of their stepwise assembly on the regulation of pre-mRNA splicing. Moreover, the role of phosphorylation-based, signal-induced changes in the regulation of the splicing process is demonstrated.

Advances in the Molecular Genetics of Non-syndromic Syndactyly

Syndactyly, webbing of adjacent digits with or without bony fusion, is one of the most common hereditary limb malformations. It occurs either as an isolated abnormality or as a component of more than 300 syndromic anomalies. There are currently nine types of phenotypically diverse nonsyndromic syndactyly. Non-syndromic syndactyly is usually inherited as an autosomal dominant trait, although the more severe presenting types and subtypes may show autosomal recessive or X-linked pattern of inheritance. The phenotype appears to be not only caused by a main gene, but also dependant on genetic background and subsequent signaling pathways involved in limb formation. So far, the principal genes identified to be involved in congenital syndactyly are mainly involved in the zone of polarizing activity and sonic hedgehog pathway. This review summarizes the recent progress made in the molecular genetics, including known genes and loci responsible for non-syndromic syndactyly, and the signaling pathways those genetic factors involved in, as well as clinical features and animal models. We hope our review will contribute to the understanding of underlying pathogenesis of this complicated disorder and have implication on genetic counseling.

Alternative mRNA transcription, processing, and translation: insights from RNA sequencing

The human transcriptome comprises >80,000 protein-coding transcripts and the estimated number of proteins synthesized from these transcripts is in the range of 250,000 to 1 million. These transcripts and proteins are encoded by less than 20,000 genes, suggesting extensive regulation at the transcriptional, post-transcriptional, and translational level. Here we review how RNA sequencing (RNA-seq) technologies have increased our understanding of the mechanisms that give rise to alternative transcripts and their alternative translation. We highlight four different regulatory processes: alternative transcription initiation, alternative splicing, alternative polyadenylation, and alternative translation initiation. We discuss their transcriptome-wide distribution, their impact on protein expression, their biological relevance, and the possible molecular mechanisms leading to their alternative regulation. We conclude with a discussion of the coordination and the interdependence of these four regulatory layers.

A novel mutation (g.106737G>T) in zone of polarizing activity regulatory sequence (ZRS) causes variable limb phenotypes in Werner mesomelia

Werner mesomelia is characterized by a sequence variation in the specific region (position 404) of the enhancer ZRS of SHH. The phenotype comprises variable mesomelia, abnormalities of the thumb and great toe and supernumerary digits. We describe extensive variation in limb phenotype in a large family and report on a novel sequence variation NG_009240.1: g.106737G>T (traditional nomenclature: ZRS404G>T) in the ZRS within the LMBR1 gene. The newly recognized clinical features in this family include small thenar eminence, sandal gap, broad first metatarsals, mesoaxial polydactyly, and postaxial polydactyly. We provide information on 12 affected family members. We review the literature on how a sequence variation in ZRS may cause such diverse phenotypes.